r - 推定からいくつかの因子間相互作用項を削除する

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ます。ただし、他のすべての相互作用を保持します。

要因を 1 つの列に格納された複数のダミー変数として理解しているので、論理的な問題はないと思いますか?

（どうやって）これを達成できますか？

[編集]
ワンステップソリューションは、私がそれを使用したいので素晴らしいでしょうstepAIC()

Answer 1

関数を使用しますupdate。

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