引导统计信息的值

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

导读：我想将bootstrap统计值(原始,偏差和错误)的值放到一个单独的列表中 – 但我无法弄清楚如何做到这一点. 这是一个例子： library(boot) set.seed(123) mean.fun - function(data,idx) { mean(data[idx]) } data - boot(data=rnorm(100),statistic=mean.fun,R=9

我想将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

现在,而不是文本,我想要实际值.显然数据$t0是“原始”但我不知道如何获得偏差和错误的值.

此外,由于键入一个函数名称为您提供了代码,我在R中键入了boost并从源代码中复制了一个片段,并尝试在我的本地R安装中搜索它.但找不到任何东西.为什么,R不应该从本地存储中获取源代码？

解决方法

std.error和bias不会存储为引导对象的一部分.它是即时计算的(见： https://stat.ethz.ch/pipermail/r-help/2011-July/284660.html)

在你的情况下,尝试：

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

（编辑：李大同）

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