¿Por qué el error "ajuste estimado 'a' es NA" generado a partir del paquete de arranque R al calcular intervalos de confianza utilizando el método bca?

Question

Tengo un vector de números que he subido aquí (alrig.com/code/MyData.Rdata) usando dput.

Me gustaría obtener el bca ci así que escribí este código:

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

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

Pero cuando ejecuto lo siguiente, obtengo esto:

 > 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
 

¿Alguien puede ayudarme a descubrir qué está causando este error? ¡Gracias por tu ayuda!

Answer 1

Como se puede ver en su mensaje de error, boot.ci llamadas bca.ci. Debido a que el boot.out objeto de no suministrar L, el empírica de la influencia de los valores de la estadística que trata de calcular en los datos, bca.ci intenta calcular utilizando la empinf de la función y, a continuación, (como dice Michael) se utiliza para calcular la aceleración constante:

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

Pero con un pequeño número de repeticiones, empinf a veces falla y devuelve un vector de NA valores. El resultado es que usted no tiene los valores de L, a no puede ser calculado, y usted consigue su error. Como ocram dice, aumentando el número de boostrap repeticiones va a solucionar esto. Incluso duplicando R a la de 2000, probablemente debería hacerlo.