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""
.Wyodrębnij standardowe błędy z glm
Odpowiedz
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ą <.
Czy w obiekcie 'glm.D93' przechowywane są błędy standardowe? Nie mogłem tego odczytać za pomocą 'str()'. Czy też 'summary()' jawnie obliczyć błędy? –
@ mindless.panda - AFAIK są obliczane bezpośrednio przez 'summary.glm'. Jeśli wpiszesz funkcję do konsoli sans '()', a następnie przewiń w dół o około 25 linii, zobaczysz, gdzie jest ona obliczana. – Chase
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
Może pomóc w przygotowaniu danych i przykładowego kodu. – screechOwl