删除R中具有因子值的交互项

我有以下型号：

> summary(mymodel.gml6) 

Call:
glm(formula = anumber ~ poly(coun, 2, raw = TRUE) + pharm+ 
    pharm:patus, family = poisson, data = mydata)

Deviance Residuals: 
    Min       1Q   Median       3Q      Max  
-2.4805  -1.7070  -0.7171   0.4482  12.6264  

Coefficients:
                                 Estimate Std. Error z value Pr(>|z|)    
(Intercept)                 -0.2786022  0.0570305  -4.885 1.03e-06 ***
poly(coun, 2, raw = TRUE)1   0.2217527  0.0162746  13.626  < 2e-16 ***
poly(coun, 2, raw = TRUE)2  -0.0156164  0.0009538 -16.372  < 2e-16 ***
pharmyes                     0.3945798  0.0343844  11.476  < 2e-16 ***
pharmno:patusyes             0.0178374  0.0352272   0.506    0.613    
pharmyes:patusyes            0.2206909  0.0311678   7.081 1.43e-12 ***

---

问题和答案略有不同，因为它们侧重于删除一个完整的交互术语。

R不允许您这样做，因为不应该这样做。交互作用由两个变量表示；你需要两者都保留，或者两者都不保留。与其中一个变量相关联的行上报告的p值是否显著并不重要，因为您将通过指定相同模型的不同方式，获得相同数据的不同显著性模式。要测试交互，需要同时测试这两个变量。在我看来，您的模型中似乎有一些奇怪的地方，但是没有足够的信息显示出来，尽管如此，您应该能够执行

drop1（mymodel.glm6，test=“LRT”）

以在R中进行测试。

R不允许您这样做，因为它不应该这样做。交互作用由两个变量表示；你需要两者都保留，或者两者都不保留。与其中一个变量相关联的行上报告的p值是否显著并不重要，因为您将通过指定相同模型的不同方式，获得相同数据的不同显著性模式。要测试交互，需要同时测试这两个变量。在我看来，您的模型中似乎有一些奇怪的地方，但是没有提供足够的信息来查看，尽管如此，您应该能够执行

drop1（mymodel.glm6，test=“LRT”）

以在R中进行测试。

您可以要求管理员将其移动到SO（stackoverflow）。那些人在那里做更多的编程工作。你可以要求管理员将此移到SO（stackoverflow）。那些人在那里做的比我多。

mymodel<-update(mymodel, .~.-pharmno:patusyes)