Извлечь стандартные ошибки из glm

Question

Я сделал glm, и я просто хочу извлечь стандартные ошибки каждого коэффициента. Я видел в Интернете функцию se.coef(), но он не работает, он возвращает "Error: could not find function "se.coef"".

Answer 1

Информация, которую вы ищете, хранится в объекте coefficients, который был возвращен summary(). Вы можете извлечь его таким образом: 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 

Взгляните на names(summary(glm.D93)) для краткого обзора всего того, что возвращается. Более подробную информацию можно найти, проверив summary.glm, если вы хотите увидеть конкретные расчеты, которые происходят, хотя этот уровень детализации, вероятно, не нужен каждый раз, если вы не указали 0 статистики.

Answer 2

Другой способ:

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

Answer 3

se.coef() на самом деле работает. Но это не в базовом пакете: он находится в пакете {arm}: http://www.inside-r.org/packages/cran/arm/docs/se.ranef