Pourquoi l'erreur «l'ajustement estimé« a »est NA» est-elle générée à partir du package de démarrage R lors du calcul des intervalles de confiance à l'aide de la méthode bca?

J'ai un vecteur de nombres que j'ai téléchargé ici (... / code / MyData.Rdata) en utilisant dput.

Je voudrais obtenir le bca ci donc j'ai écrit ce code:

my.mean <- function(dat, idx){
 return (mean(dat[idx], na.rm = TRUE))
} 

boot.out<-boot(data=my.data, statistic = my.mean, R=1000)

Mais quand je lance ce qui suit, j'obtiens ceci:

> boot.ci(boot.out)
Error in bca.ci(boot.out, conf, index[1L], L = L, t = t.o, t0 = t0.o,  : 
estimated adjustment 'a' is NA
In addition: Warning message:
In boot.ci(boot.out) : bootstrap variances needed for studentized intervals

Quelqu'un peut-il m'aider à déterminer la cause de cette erreur? Merci de votre aide!

Comme vous pouvez le voir dans votre message d'erreur, les boot.ciappels bca.ci. Parce que l' boot.outobjet ne fournit pas L, les valeurs d'influence empiriques pour la statistique que vous calculez sur les données, bca.ciessaient de les calculer en utilisant la empinffonction, puis (comme Michael le dit), il les utilise pour calculer la constante d'accélération:

L <- empinf(boot.out, index = index, t = t.o, ...)
a <- sum(L^3)/(6 * sum(L^2)^1.5)

Mais avec un petit nombre de réplications, empinféchoue parfois et renvoie un vecteur de NAvaleurs. Le résultat est que vous n'avez aucune valeur pour L, ane peut pas être calculé et vous obtenez votre erreur. Comme le dit ocram, l'augmentation du nombre de réplications boostrap corrigera ce problème. Même doubler Rà 2000 devrait probablement le faire.