当前位置 : 主页 > 网页制作 > Bootstarp >

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

来源:互联网 收集:自由互联 发布时间:2021-06-12
考虑以下引导程序: 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, statisti
考虑以下引导程序:

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)
网友评论