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

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

ОТВЕТЫ

Ответ 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, если вы хотите увидеть конкретные вычисления, которые происходят, хотя этот уровень детализации, вероятно, не нужен каждый раз, если вы не станете статистикой.

Ответ 2

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

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