bootstrap标准错误在引导类中存在于何处？

发布时间：2020-12-17 20:39:17 所属栏目：安全来源：网络整理

导读：考虑以下引导程序： library(MASS)library(boot)# c)set.seed(1)boot.fn= function(data,index) mean(data[index])output=boot(Boston$medv,boot.fn,1000) 如果我们运行print(输出),我们得到 Call:boot(data = Boston$medv,statistic = boot.fn,R = 1000)Boo

考虑以下引导程序：

library(MASS)
library(boot)

# c)
set.seed(1)
boot.fn= function(data,index) mean(data[index])
output=boot(Boston$medv,boot.fn,1000)

如果我们运行print(输出),我们得到

Call:
boot(data = Boston$medv,statistic = boot.fn,R = 1000)


Bootstrap Statistics :
    original      bias    std. error
t1* 22.53281 0.008517589   0.4119374

但是,当我检查输出对象时,我找不到表示引导程序统计信息的值.哪里是原创,偏见和标准. boot返回的实际输出对象中的错误？

解决方法

它们由print.boot计算,不存储在引导对象中.查看getAnywhere(print.boot)了解详细信息.

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

对于你的例子：

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

（编辑：李大同）

【声明】本站内容均来自网络，其相关言论仅代表作者个人观点，不代表本站立场。若无意侵犯到您的权利，请及时与联系站长删除相关内容!