web-dev-qa-db-ja.com

glmから標準エラーを抽出する

glmを実行し、各係数の標準誤差を抽出したいだけです。インターネットで関数se.coef()を見ましたが、機能せず、"Error: could not find function "se.coef""。

rextractglmstandard-error
14
user1096592

必要な情報は、summary()によって返されるcoefficientsオブジェクトに格納されます。このように抽出できます：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_をチェックして詳細を確認できますが、統計が3未満でない限り、そのレベルの詳細は毎回必要になるとは限りません。

33
Chase

別の方法：

sqrt(diag(vcov(glm.D93)))
22
Wojciech Sobala

se.coef（）は実際に機能します。ただし、基本パッケージには含まれていません。{arm}パッケージに含まれています： http://www.inside-r.org/packages/cran/arm/docs/se.ranef

4
Joel Chan