r - Bootstrap 标准错误位于引导类中的什么位置？

Question

考虑以下 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引导返回的对象？

Answer 1

它们由 print.boot 计算而不是存储在 boot 中目的。看 getAnywhere(print.boot)详情。

您可以自己计算这些值，也可以使用 capture.output .

对于您的示例:

#original:
output$t0
#bias:
mean(output$t)-output$t0
#se: 
sd(output$t)