Trích xuất các lỗi chuẩn từ glm

Question

Tôi đã thực hiện glm và tôi chỉ muốn trích xuất các lỗi chuẩn của từng hệ số. Tôi thấy trên internet chức năng se.coef() nhưng nó không hoạt động, nó trả về "Error: could not find function "se.coef"".

Answer 1

Thông tin bạn sau được lưu trữ trong đối tượng coefficients được trả lại bởi summary(). Bạn có thể giải nén nó thusly: summary(glm.D93)$coefficients[, 2]

#Example from ?glm 
counts <- c(18,17,15,20,10,20,25,13,12) 
outcome <- gl(3,1,9) 
treatment <- gl(3,3) 
print(d.AD <- data.frame(treatment, outcome, counts)) 
glm.D93 <- glm(counts ~ outcome + treatment, family=poisson()) 

#coefficients has the data of interest 
> summary(glm.D93)$coefficients 
       Estimate Std. Error  z value  Pr(>|z|) 
(Intercept) 3.044522e+00 0.1708987 1.781478e+01 5.426767e-71 
outcome2 -4.542553e-01 0.2021708 -2.246889e+00 2.464711e-02 
outcome3 -2.929871e-01 0.1927423 -1.520097e+00 1.284865e-01 
treatment2 1.337909e-15 0.2000000 6.689547e-15 1.000000e+00 
treatment3 1.421085e-15 0.2000000 7.105427e-15 1.000000e+00 

#So extract the second column 
> summary(glm.D93)$coefficients[, 2] 
(Intercept) outcome2 outcome3 treatment2 treatment3 
    0.1708987 0.2021708 0.1927423 0.2000000 0.2000000 

Hãy xem names(summary(glm.D93)) cho một bài đánh giá nhanh chóng của tất cả những gì được trả về. Bạn có thể tìm thêm chi tiết bằng cách xem summary.glm nếu bạn muốn xem các tính toán cụ thể đang diễn ra, mặc dù mức chi tiết có thể không cần thiết mỗi lần trừ khi bạn có < 3 thống kê.

Answer 2

Một cách khác:

sqrt(diag(vcov(glm.D93))) 

Answer 3

se.coef() thực hiện công việc. Nhưng nó không nằm trong gói cơ sở: nó nằm trong gói {arm}: http://www.inside-r.org/packages/cran/arm/docs/se.ranef