为什么我不能使用boot
包来引导大n的统计信息?虽然,150,000只不大,所以我不知道为什么这不起作用。
示例
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
错误消息
在统计中(data,i [r,],...):整数溢出 - 使用 总和(as.numeric(。))
答案 0 :(得分:1)
您未按照文件记录使用boot()
(无可否认,这非常复杂)。来自?boot
:
在所有其他情况下'统计' 必须至少采取两个论点。第一个论点通过了 永远是原始数据。第二个是 指数,频率或权重的矢量 bootstrap样本。
我想你想要:
bsum <- function(x,i) sum(x[i])
bs <- boot(rnorm(150000), bsum, R = 1000)
我还没有花时间弄清楚boot()
实际在你的情况下做了什么 - 几乎肯定不是你想要的。