[英]Confidence Interval of Sample Means using R
我的數據框包含500個樣本,每個樣本的大小為100。 下面是快照。 我需要計算平均值為90/95/99的置信區間。
head(Means_df)
Means
1 14997
2 11655
3 12471
4 12527
5 13810
6 13099
我正在使用以下代碼,但僅獲得一行的置信區間。 誰能幫我提供代碼?
tint <- matrix(NA, nrow = dim(Means_df)[2], ncol = 2)
for (i in 1:dim(Means_df)[2]) {
temp <- t.test(Means_df[, i], conf.level = 0.9)
tint[i, ] <- temp$conf.int
}
colnames(tint) <- c("lcl", "ucl")
Means_df
是具有500行和1列的數據幀。 因此
dim(Means_df)[2]
將給出值1
。
這就是為什么您只能獲得一個價值的原因。
通過使用dim(Means_df)[1]
或更好的nrow(Means_df)
代替dim(Means_df)[2]
解決此問題。
對於任何單一平均值,例如14997,您都無法在不知道數據方差或標准差的情況下計算95%-CI。 如果可以訪問每個樣本的標准偏差,則可以計算平均值的標准誤,並由此輕松計算95%-CI。 顯然,您缺少任務所需的信息。
聲明:本站的技術帖子網頁,遵循CC BY-SA 4.0協議,如果您需要轉載,請注明本站網址或者原文地址。任何問題請咨詢:yoyou2525@163.com.