Valeurs des statistiques bootstrap

Question

Je veux obtenir les valeurs des statistiques bootstrap (original, bias et error) dans une liste séparée - mais je n'arrive pas à comprendre comment faire cela.

Voici un exemple:

> 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 

Maintenant, au lieu de texte que je veux que les valeurs réelles. Apparemment data$t0 est le "original" mais je ne vois pas comment obtenir les valeurs de biais et d'erreur.

---

En outre, depuis la saisie d'un nom de fonction vous donne son code, je dactylographiées boost en R et copier un extrait du code source, et a essayé de chercher sur mon installation R local. Mais n'a rien trouvé. Pourquoi, ne devrait pas prendre ce code source à partir d'un stockage local?

Answer 1

Les erreurs std.error et bias ne sont pas stockées en tant que partie de l'objet de démarrage. Il est calculé à la volée (voir: https://stat.ethz.ch/pipermail/r-help/2011-July/284660.html)

Dans votre cas, essayez:

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