gói Boot trong hỗ trợ đơn giản R

Question

Nếu tôi muốn sử dụng các boot() chức năng từ gói R của boot tính tầm quan trọng của hệ số Pearson mối tương quan giữa hai vectơ, tôi nên làm điều đó như thế này:

boot(re1, cor, R = 1000) 

nơi re1 là ma trận hai cột cho hai vectơ quan sát này? Tôi dường như không thể có được quyền này bởi vì cor của các vectơ này là 0.8, nhưng hàm trên trả về -0.2 là t0.

Answer 1

Chỉ cần nhấn mạnh ý tưởng chung về khởi động trong R, mặc dù @caracal đã trả lời câu hỏi của bạn thông qua nhận xét của anh ấy. Khi sử dụng boot, bạn cần phải có cấu trúc dữ liệu (thường là ma trận) có thể được lấy mẫu theo hàng. Việc tính toán thống kê của bạn thường được thực hiện trong một hàm nhận ma trận dữ liệu này và trả về số liệu thống kê lãi suất được tính sau khi lấy lại mẫu. Sau đó, bạn gọi số boot() xử lý chức năng này để R sao chép và thu thập kết quả theo định dạng có cấu trúc. Những kết quả này có thể được đánh giá bằng cách sử dụng boot.ci() lần lượt.

Dưới đây là hai ví dụ hoạt động với nghiên cứu low birth baby trong gói MASS.

require(MASS) 
data(birthwt) 
# compute CIs for correlation between mother's weight and birth weight 
cor.boot <- function(data, k) cor(data[k,])[1,2] 
cor.res <- boot(data=with(birthwt, cbind(lwt, bwt)), 
       statistic=cor.boot, R=500) 
cor.res 
boot.ci(cor.res, type="bca") 
# compute CI for a particular regression coefficient, e.g. bwt ~ smoke + ht 
fm <- bwt ~ smoke + ht 
reg.boot <- function(formula, data, k) coef(lm(formula, data[k,])) 
reg.res <- boot(data=birthwt, statistic=reg.boot, 
       R=500, formula=fm) 
boot.ci(reg.res, type="bca", index=2) # smoke