引導程序統計值

我想將引導程序統計信息（原始值，偏差和錯誤）的值放入單獨的列表中 - 但我無法弄清楚如何做到這一點。引導程序統計值

下面是一個例子：

> 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不應該從本地存儲中獲取源代碼？

來源

2012-12-03 andreister

std.error和bias不作爲引導對象的一部分存儲。它是實時計算（見：https://stat.ethz.ch/pipermail/r-help/2011-July/284660.html）

在你的情況下，嘗試：

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

來源

2012-12-03 16:31:08 jokel

這是有道理的，謝謝。你能否告訴我爲什麼我的機器上找不到'boot'的R代碼？我在R中輸入「boot」，從輸出中找到一行（例如'lapply（seq_len（RR），fn）'）並搜索它 - 在我的本地R安裝中找不到任何內容。 – andreister

引導程序統計值

回答

相關問題