Bootstrap统计量的含义

Question

我希望能够将Bootstrap统计值（原始值、偏差和误差）获取到一个单独的列表中，但我无法弄清楚如何实现。

以下是一个示例：

> library(boot)
> set.seed(123)
> mean.fun <- function(data, idx) { mean(data[idx]) }
> data <- boot(data=rnorm(100), statistic=mean.fun, R=999)
> names(data)
 [1] "t0"        "t"         "R"         "data"     
 [5] "seed"      "statistic" "sim"       "call"     
 [9] "stype"     "strata"    "weights"  
> data
ORDINARY NONPARAMETRIC BOOTSTRAP


Call:
boot(data = rnorm(100), statistic = mean.fun, R = 999)


Bootstrap Statistics :
    original   bias    std. error
t1* 0.09040591 0.004751773  0.08823615

现在，我希望看到实际的数值，而不是文字。显然，data$t0 是“原始数据”，但我不知道如何获取偏差和误差的值。

---

另外，由于输入函数名称可以显示其代码，我在R中输入了boost并从源代码中复制了一小段代码，试图在本地的R安装中搜索它。但是找不到任何东西。为什么？难道R不应该从本地存储中获取源代码吗？

- andreister

Answer 1

标准误差和偏差没有作为引导对象的一部分存储。它是即时计算的（参见：https://stat.ethz.ch/pipermail/r-help/2011-July/284660.html）

在您的情况下，请尝试：

mean(data$t) - data$t0
sd(data$t)

- jokel

---

那很有道理，谢谢。您能否告诉我为什么我在我的计算机上找不到boot的R代码？我在R中键入了"boot"，从输出中选取一行（例如lapply(seq_len(RR),fn)）并进行搜索 - 在我的本地R安装中没有找到任何内容。 - andreister