从估算中删除一些因子交互项

Question

我想从估算中删除某些因素互动。以下是来自想象中的劳动力市场的生成数据的示例（我在此处上传：http://pastebin.com/raw.php?i=EcMEVqUC）

s <- source("http://pastebin.com/raw.php?i=EcMEVqUC")$value

lm(income ~ age + cit * prof, data=s)

在这个示例经济中，外国人不允许在公共部门工作，因此citforeign:profofficial是NA。因此，我想排除citforeign:profofficial的互动术语。但保留所有其他互动。

我理解因为存储在一列中的多个虚拟变量的因素我不认为存在逻辑问题吗？

（如何）我可以实现这个目标吗？

[编辑]
一步到位的解决方案会非常棒，因为我希望将其与stepAIC()

一起使用

Answer 1

使用函数update。

model1 <- lm(income ~ age + cit * prof, data=s)
model2 <- update(model1, . ~ . - citforeign:profofficial)