Wyodrębnij standardowe błędy z glm

Question

Zrobiłem glm i chcę tylko wyodrębnić standardowe błędy każdego współczynnika. Widziałem w Internecie funkcję se.coef(), ale to nie działa, zwraca "Error: could not find function "se.coef"".

Answer 1

Informacje, których szukasz, są przechowywane w obiekcie coefficients zwróconym przez summary(). Można wyodrębnić go 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 

Spójrz na names(summary(glm.D93)) na szybki przegląd wszystkiego, co jest zwracany. Więcej szczegółów można znaleźć, sprawdzając summary.glm, jeśli chcesz zobaczyć konkretne obliczenia, które są w toku, choć ten poziom szczegółowości prawdopodobnie nie jest potrzebny za każdym razem, chyba że statystyki te są <.

Answer 2

Innym sposobem:

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

Answer 3

se.coef() faktycznie działa. Ale nie ma go w pakiecie podstawowym: jest w pakiecie {arm}: http://www.inside-r.org/packages/cran/arm/docs/se.ranef