Valori delle statistiche di bootstrap

Question

Desidero ottenere i valori delle statistiche di bootstrap (originale, bias e errore) in un elenco separato, ma non riesco a capire come farlo.

Ecco un esempio:

> 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 

Ora, al posto del testo che voglio i valori effettivi. Apparentemente data$t0 è "originale" ma non vedo come ottenere i valori di bias e di errore.

---

Inoltre, dal momento che a digitare un nome di funzione vi dà il suo codice, ho digitato boost in R e copiato un frammento dal codice sorgente, e provato a cercare sul mio installazione R locale. Ma non ho trovato nulla. Perché, non dovresti prendere questo codice sorgente da una memoria locale?

Answer 1

Std.error e bias non vengono memorizzati come parte dell'oggetto di avvio. Si è calcolato al volo (vedi: https://stat.ethz.ch/pipermail/r-help/2011-July/284660.html)

Nel tuo caso, prova:

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