我有一个在我的摘要中产生不同级别的glm,我想使用此摘要中的特定级别作为我的因素之一。
我的代码是
Model3<-glm(status ~ month + age + centre + treatment, family="binomial")
summary(Model3)
这是
Coefficients:
Estimate Std. Error z value Pr(>|z|)
(Intercept) -0.487717 0.315126 -1.548 0.12170
month1 0.606370 0.286407 2.117 0.03425 *
month2 0.402645 0.284686 1.414 0.15726
month3 0.606370 0.286407 2.117 0.03425 *
month4 0.362311 0.284486 1.274 0.20282
age -0.020200 0.006931 -2.915 0.00356 **
centre2 1.090830 0.189501 5.756 8.60e-09 ***
treatment1 0.822554 0.181614 4.529 5.92e-06 ***
我希望以某种方式能够从我的下一个glm中移除month2
,但是如果不将它们全部移除,我看不出这是怎么回事。