R使用bootstrap计算标准误差

Question

我有这个值数组:

> df
[1] 2 0 0 2 2 0 0 1 0 1 2 1 0 1 3 0 0 1 1 0 0 0 2 1 2 1 3 1 0 0 0 1 1 2 0 1 3
[38] 1 0 2 1 1 2 2 1 2 2 2 1 1 1 2 1 0 0 0 0 0 0 0 0 0 0 1 0 1 1 0 1 0 0 0 0 0
[75] 0 0 0 0 0 1 1 0 1 1 1 1 3 1 3 0 1 2 2 1 2 3 1 0 0 1

我想使用包boot来计算数据的标准误差。 http://www.ats.ucla.edu/stat/r/faq/boot.htm

所以，我用这个命令来追求:

library(boot)
boot(df, mean, R=10)

我得到了这个错误:

Error in mean.default(data, original, ...) : 
'trim' must be numeric of length one

有人可以帮我解决这个问题吗？谢谢

Answer 1

如果您正在引导均值，您可以执行以下操作:

set.seed(1)
library(boot)
x<-rnorm(100)
meanFunc <- function(x,i){mean(x[i])}
bootMean <- boot(x,meanFunc,100)
>bootMean

ORDINARY NONPARAMETRIC BOOTSTRAP


Call:
boot(data = x, statistic = meanFunc, R = 100)


Bootstrap Statistics :
     original      bias    std. error
t1* 0.1088874 0.002614105  0.07902184

如果您只是输入 mean 作为参数，您将得到如下错误:

bootMean <- boot(x,mean,100)
Error in mean.default(data, original, ...) : 
  'trim' must be numeric of length one