R中引导标准差的置信区间

Question

我是初学者。我想找到从引导程序生成的SD的CI，并且想知道我是否已经正确地完成了它。谢谢！

Repnumber <- 1000                    
group1.sd.values <- numeric(Repnumber)    

for (i in 1:Repnumber){
group1 = sample(Data, size=length(Data), replace=T)
group1.sd.values[i] = sd(group1)
}

group1.sd.upperCI <- quantile(group1.sd.values, probs = .975) 
group1.sd.lowerCI <- quantile(group1.sd.values, probs = .025)

Answer 1

以下是我在单个向量上进行自举的方法：

CI <- quantile(replicate(Repnumber, sd(sample(Data, rep = TRUE))), c(.025, .975))

这将为您提供长度为2的向量中间隔的上限值和下限值。

请注意，使用replicate通常应该比循环快得多。