Tại sao lỗi điều chỉnh ước tính 'a' là NA NA được tạo từ gói khởi động R khi tính khoảng tin cậy bằng phương pháp bca?

Tôi có một vectơ số mà tôi đã tải lên ở đây (... / code / MyData.Rdata) bằng cách sử dụng dput.

Tôi muốn lấy bca ci vì vậy tôi đã viết mã này:

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

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

Nhưng khi tôi chạy như sau tôi nhận được điều này:

> 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

Ai đó có thể giúp tôi tìm ra những gì gây ra lỗi này? Cảm ơn sự giúp đỡ của bạn!

Như bạn có thể thấy từ thông báo lỗi của bạn, boot.cicác cuộc gọi bca.ci. Bởi vì boot.outđối tượng không cung cấp L, các giá trị ảnh hưởng theo kinh nghiệm cho thống kê mà bạn đang tính toán trên dữ liệu, bca.cicố gắng tính toán chúng bằng empinfhàm và sau đó (như Michael nói), nó sử dụng chúng để tính hằng số gia tốc:

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

Nhưng với một số lượng nhỏ các bản sao, empinfđôi khi thất bại và trả về một vectơ NAgiá trị. Kết quả là bạn không có giá trị nào L, akhông thể tính được và bạn gặp lỗi. Như ocram nói, việc tăng số lượng bản sao boostrap sẽ khắc phục điều này. Thậm chí tăng gấp đôi Rlên 2000 có lẽ nên làm điều đó.