考虑以下 bootstrap :
library(MASS)
library(boot)
# c)
set.seed(1)
boot.fn= function(data, index) mean(data[index])
output=boot(Boston$medv, boot.fn, 1000)
如果我们运行
print(output)
,我们得到
Call:
boot(data = Boston$medv, statistic = boot.fn, R = 1000)
Bootstrap Statistics :
original bias std. error
t1* 22.53281 0.008517589 0.4119374
但是,当我检查
output
时对象,我找不到表示引导统计数据的值。在哪里
original
,
bias
和
std. error
在实际
output
引导返回的对象?
请您参考如下方法:
它们由 print.boot
计算而不是存储在 boot
中目的。看 getAnywhere(print.boot)
详情。
您可以自己计算这些值,也可以使用 capture.output
.
对于您的示例:
#original:
output$t0
#bias:
mean(output$t)-output$t0
#se:
sd(output$t)