Por que o erro "ajuste estimado 'a' é NA" gerado no pacote de inicialização R ao calcular intervalos de confiança usando o método bca?

Eu tenho um vetor de números que eu enviei aqui (... / code / MyData.Rdata) usando dput.

Gostaria de obter o bca ci, então eu escrevi 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)

Mas quando executo o seguinte, recebo o seguinte:

> 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

Alguém pode me ajudar a descobrir o que está causando esse erro? Obrigado pela ajuda!

Como você pode ver na sua mensagem de erro, boot.cichama bca.ci. Como o boot.outobjeto não fornece L, os valores empíricos de influência para a estatística que você está calculando nos dados, bca.citentam calculá-los usando a empinffunção e, em seguida (como Michael diz), eles os usam para calcular a constante de aceleração:

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

Mas, com um pequeno número de repetições, empinfàs vezes falha e retorna um vetor de NAvalores. O resultado é que você não tem valores para L, anão pode ser calculado e obtém seu erro. Como o ocram diz, aumentar o número de replicações de boostrap corrigirá isso. Mesmo dobrar Rpara 2000 provavelmente deve fazê-lo.