Giá trị của số liệu thống kê bootstrap

Question

Tôi muốn lấy các giá trị của số liệu thống kê bootstrap (bản gốc, sai lệch và lỗi) vào một danh sách riêng biệt - nhưng tôi không thể tìm ra cách để làm điều đó.

Dưới đây là một ví dụ:

> library(boot) 
> set.seed(123) 
> mean.fun <- function(data, idx) { mean(data[idx]) } 
> data <- boot(data=rnorm(100), statistic=mean.fun, R=999) 
> names(data) 
[1] "t0"  "t"   "R"   "data"  
[5] "seed"  "statistic" "sim"  "call"  
[9] "stype"  "strata" "weights" 
> data 
ORDINARY NONPARAMETRIC BOOTSTRAP 


Call: 
boot(data = rnorm(100), statistic = mean.fun, R = 999) 


Bootstrap Statistics : 
    original bias std. error 
t1* 0.09040591 0.004751773 0.08823615 

Bây giờ, thay vì văn bản tôi muốn các giá trị thực tế. Rõ ràng data$t0 là "bản gốc" nhưng tôi không thấy cách nhận giá trị sai lệch và sai lệch.

---

Ngoài ra, kể từ khi gõ một tên hàm cho phép bạn mã của nó, tôi gõ boost trong R và sao chép một đoạn từ mã nguồn, và cố gắng tìm kiếm nó trên cài đặt R địa phương của tôi. Nhưng không thể tìm thấy gì cả. Tại sao, không nên R lấy mã nguồn đó từ bộ nhớ cục bộ?

Answer 1

Std.error và bias không được lưu trữ như một phần của đối tượng khởi động. Nó được tính toán một cách nhanh chóng (xem: https://stat.ethz.ch/pipermail/r-help/2011-July/284660.html)

Trong trường hợp của bạn, hãy thử:

mean(data$t) - data$t0 
sd(data$t)