考虑以下引导程序: 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)