使用 'boot' 计算 R 中的自举平均值

Question

我想使用 R 命令“boot”来计算自举平均值。 但是，这不起作用，并且总是返回错误消息：

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

这是我正在使用的代码示例：

df <- data.frame(
  v1 <- c(0,159028,236220,9127,0,0),
  v2 <- c(13, 42, 56, 9, 77, 34)
)

boot(df$v1, statistic = mean, R = 100)

返回上述错误。

为什么这个命令不起作用？ 我该怎么做才能解决这个问题？

Answer 1

您的statistic应该是一个 function，它需要两个 arguments：第一个应该是要采样的数据，第二个应该是boot将用来对数据进行子集化的索引。 由于您直接将mean传递给statistic参数，因此索引向量将传递给mean.default的第二个参数，即trim 。 此参数需要一个数字，并在传递向量时抛出错误。

您需要创建一个小包装器boot来显示您希望它如何对数据进行子集化以进行采样：

boot(df$v1, statistic = function(x, inds) mean(x[inds]), R = 100)
#>
#> ORDINARY NONPARAMETRIC BOOTSTRAP
#>
#>
#> Call:
#> boot(data = df$v1, statistic = function(x, inds) mean(x[inds]), 
#>     R = 100)
#>
#>
#> Bootstrap Statistics :
#>     original   bias    std. error
#> t1* 67395.83 839.4583    39734.55