使用R的样本均值的置信区间

Question

我的数据框包含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")

Answer 1

Means_df是具有500行和1列的数据帧。 因此

dim(Means_df)[2]

将给出值1 。

这就是为什么您只能获得一个价值的原因。

通过使用dim(Means_df)[1]或更好的nrow(Means_df)代替dim(Means_df)[2]解决此问题。

Answer 2

对于任何单一平均值，例如14997，您都无法在不知道数据方差或标准差的情况下计算95％-CI。 如果可以访问每个样本的标准偏差，则可以计算平均值的标准误，并由此轻松计算95％-CI。 显然，您缺少任务所需的信息。