引导统计值答案

【问题标题】：Values of bootstrap statistics引导统计值
【发布时间】：2012-11-21 11:24:46
【问题描述】：

我想将引导统计的值（原始、偏差和错误）放入一个单独的列表中 - 但我不知道该怎么做。

这是一个例子：

> 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 并从源代码中复制了一个 sn-p，并尝试在我的本地 R 安装中搜索它。但是什么也找不到。为什么，R 不应该从本地存储中获取源代码？

【问题讨论】：

标签： r statistics-bootstrap

【解决方案1】：

std.error 和bias 不作为引导对象的一部分存储。它是即时计算的（请参阅：https://stat.ethz.ch/pipermail/r-help/2011-July/284660.html）

在你的情况下，尝试：

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

【讨论】：

有道理，谢谢。你能告诉我为什么我在我的机器上找不到boot 的 R 代码吗？我在 R 中输入“boot”，从输出中选取一行（例如 lapply(seq_len(RR), fn)）并搜索它 - 在我的本地 R 安装中没有找到。