Warum wird der Fehler "geschätzte Anpassung 'a' ist NA" aus dem R-Boot-Paket generiert, wenn Konfidenzintervalle mit der bca-Methode berechnet werden?

Ich habe einen Vektor von Zahlen, den ich hier mit dput hochgeladen habe (... / code / MyData.Rdata).

Ich möchte das bca ci bekommen, also habe ich diesen Code geschrieben:

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

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

Aber wenn ich Folgendes ausführe, bekomme ich Folgendes:

> 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

Kann mir jemand helfen, herauszufinden, was diesen Fehler verursacht? Danke für deine Hilfe!

Wie Sie Ihrer Fehlermeldung entnehmen können, werden boot.ciAnrufe getätigt bca.ci. Da das boot.outObjekt nicht liefert L, bca.civersuchen die empirischen Einflusswerte für die Statistik, die Sie für die Daten berechnen, diese mithilfe der empinfFunktion zu berechnen und verwenden sie dann (wie Michael sagt) zur Berechnung der Beschleunigungskonstante:

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

Bei einer kleinen Anzahl von Replikationen empinfschlägt dies jedoch manchmal fehl und gibt einen Wertevektor zurück NA. Das Ergebnis ist, dass Sie keine Werte für haben L, anicht berechnet werden können und Ihren Fehler erhalten. Wie ocram sagt, wird dies durch Erhöhen der Anzahl der Boostrap-Replikationen behoben. Selbst eine Verdoppelung Rauf 2000 sollte es wahrscheinlich tun.