如何打印glm系数,包括参考水平的所有因子水平?
summary(glm_obj) 仅打印偏离参考值的值。
我知道那些是0,但我需要这个来进行集成,即告诉其他系统可能发生的因子水平。
如果太简单,请原谅,我找不到任何地方。
谢谢
为了说明我面临的问题:
上面的模型描述只提供了鸢尾花属种和维吉尼亚种的系数,从模型本身的角度来看这是完全可以接受的,正如Dason所指出的那样。
但是,我需要与另一个应用程序共享该模型,该应用程序必须知道何种“物种”水平可预期(例如,在出现未研究过的新水平时发出警告)。
Summary()给出了相同的结果:
我知道那些是0,但我需要这个来进行集成,即告诉其他系统可能发生的因子水平。
如果太简单,请原谅,我找不到任何地方。
谢谢
为了说明我面临的问题:
> glm(Petal.Width~Species,data=iris)
Call: glm(formula = Petal.Width ~ Species, data = iris)
Coefficients:
(Intercept) Speciesversicolor Speciesvirginica
0.246 1.080 1.780
Degrees of Freedom: 149 Total (i.e. Null); 147 Residual
Null Deviance: 86.57
Residual Deviance: 6.157 AIC: -45.29`
上面的模型描述只提供了鸢尾花属种和维吉尼亚种的系数,从模型本身的角度来看这是完全可以接受的,正如Dason所指出的那样。
但是,我需要与另一个应用程序共享该模型,该应用程序必须知道何种“物种”水平可预期(例如,在出现未研究过的新水平时发出警告)。
Summary()给出了相同的结果:
> summary(glm(Petal.Width~Species,data=iris))
Call:
glm(formula = Petal.Width ~ Species, data = iris)
Deviance Residuals:
Min 1Q Median 3Q Max
-0.626 -0.126 -0.026 0.154 0.474
Coefficients:
Estimate Std. Error t value Pr(>|t|)
(Intercept) 0.24600 0.02894 8.50 1.96e-14 ***
Speciesversicolor 1.08000 0.04093 26.39 < 2e-16 ***
Speciesvirginica 1.78000 0.04093 43.49 < 2e-16 ***
---
Signif. codes: 0 ‘***’ 0.001 ‘**’ 0.01 ‘*’ 0.05 ‘.’ 0.1 ‘ ’ 1
(Dispersion parameter for gaussian family taken to be 0.04188163)
Null deviance: 86.5699 on 149 degrees of freedom
Residual deviance: 6.1566 on 147 degrees of freedom
AIC: -45.285
Number of Fisher Scoring iterations: 2