Bca yöntemi kullanılarak güven aralıkları hesaplanırken neden R önyükleme paketinden “tahmini düzenleme 'a' NA” hatası üretiliyor?

Ben dput kullanarak burada yüklediğim sayıların bir vektör (... / code / MyData.Rdata) var.

Ben bca ci almak istiyorum bu yüzden bu kodu yazdım:

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

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

Ama aşağıdakileri çalıştırdığımda bunu alıyorum:

> 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

Birisi bu hataya neyin neden olduğunu anlamama yardımcı olabilir mi? Yardımın için teşekkürler!

Hata mesajınızdan da görebileceğiniz gibi, boot.ciarayın bca.ci. Çünkü boot.outnesne arz etmez L, sen verilere hesaplanması ediyoruz istatistik için ampirik etkisi değerleri bca.ciçalışır kullanarak bunları hesaplamak için empinf(Michael diyor gibi) ivme sabitini hesaplamak için bunları kullanır sonra işlevini ve:

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

Ancak az sayıda çoğaltma ile empinfbazen başarısız olur ve bir değer vektörü döndürür NA. Sonuç olarak, için hiçbir değeriniz yoktur L, ahesaplanamaz ve hatanızı alırsınız. Ocram'ın dediği gibi, boostrap replikasyonlarının sayısını arttırmak bunu düzeltir. R2000'e kadar iki katına çıkması bile muhtemelen bunu yapmalıdır.