R中的Bootstrap置信区间

Question

我是新用户,使用启动包时遇到问题.我想做的就是使用自举来产生数字向量的均值周围的置信区间,例如:

x <- rnorm(100, 1, .5)

有小费吗？

Answer 1

以下不够吗？

library(boot)
x <- rnorm(100, 1, .5)
b <- boot(x, function(u,i) mean(u[i]), R = 999)
boot.ci(b, type = c("norm", "basic", "perc"))

@ cryptic0:在`？boot`中解释.`u`将是原始数据(在本例中为'x`),'i`是对应于boostrap样本的索引向量(对于不同的样本,该函数将被称为`R`次).如果数据是向量,则引导样本是"u [i]",如果是data.frame,则为"u [i,]".(当然,因为`u`和`i`是函数的形式参数,所以它们在此函数之外没有任何值.) (2认同)