使用引导包引导大 n 时出错（错误：整数溢出）

Question

为什么我不能使用boot 包引导具有大 n 的统计信息？虽然，150,000 obs 并不大，所以我不知道为什么这不起作用。

示例

library(boot)

bs <- boot(rnorm(150000), sum, R = 1000)
bs

ORDINARY NONPARAMETRIC BOOTSTRAP


Call:
boot(data = rnorm(150000), statistic = sum, R = 1000)


Bootstrap Statistics :
WARNING: All values of t1* are NA

错误信息

在 statistic(data, i[r, ], ...) 中：整数溢出 - 使用 sum(as.numeric(.))

Answer 1

您没有使用文档中的boot()（诚然，这非常复杂）。来自?boot：

在所有其他情况下，“统计” 必须至少有两个参数。第一个参数通过 将始终是原始数据。第二个将是 定义 引导样本。

我想你想要：

bsum <- function(x,i) sum(x[i])
bs <- boot(rnorm(150000), bsum, R = 1000)

我还没有花时间弄清楚 boot() 在你的情况下实际上在做什么 - 虽然几乎可以肯定不是你想要的。