考虑以下 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 , biasstd. error在实际 output引导返回的对象?

请您参考如下方法:

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

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

对于您的示例:

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


评论关闭
IT干货网

微信公众号号:IT虾米 (左侧二维码扫一扫)欢迎添加!