为什么我不能使用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,],...):整数溢出 - 使用sum(as.numeric(。))
你没有使用boot()
作为文件记录(诚然,这是非常复杂的)。来自?boot
:
在所有其他情况下,“统计”必须至少采用两个参数。传递的第一个参数将始终是原始数据。第二个是定义引导样本的索引,频率或权重的向量。
我想你想要:
bsum <- function(x,i) sum(x[i])
bs <- boot(rnorm(150000), bsum, R = 1000)
我没有花时间弄清楚boot()
在你的情况下实际上做了什么 - 几乎肯定不是你想要的。