Extraiga errores estándar de glm

Question

Hice un glm y solo quiero extraer los errores estándar de cada coeficiente. Vi en Internet la función se.coef() pero no funciona, devuelve "Error: could not find function "se.coef"".

Answer 1

La información que busca se almacena en el objeto coefficients devuelto por summary(). Puede extraer de esta manera: 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 

Tome un vistazo a names(summary(glm.D93)) para una revisión rápida de todo lo que se devuelve. Se pueden encontrar más detalles visitando summary.glm si desea ver los cálculos específicos que están en marcha, aunque ese nivel de detalle probablemente no sea necesario todas las veces, a menos que tenga < 3 estadísticas.

Answer 2

Otra forma:

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

Answer 3

se.coef() realmente hace el trabajo. Pero no está en el paquete base: está en el paquete {arm}: http://www.inside-r.org/packages/cran/arm/docs/se.ranef