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

Tengo un vector de números que he subido aquí (... / 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 me sale 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!

Como puede ver en su mensaje de error, boot.cillamadas bca.ci. Debido a que el boot.outobjeto no se suministra L, los valores de influencia empírica para la estadística que está calculando en los datos, bca.ciintenta calcularlos usando la empinffunción y luego (como dice Michael) los usa para calcular la constante de aceleración:

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, a empinfveces falla y devuelve un vector de NAvalores. El resultado es que no tiene valores L, ano puede calcularse y obtiene su error. Como dice ocram, aumentar el número de repeticiones boostrap solucionará esto. Incluso duplicar Ra 2000 probablemente debería hacerlo.