考虑一个由双边贸易数据和政策虚拟组成的大型小组（0级和1级，参考级别0）。
head(data)
home partner year gdphome gdppartner exports policy
Austria Albania 1988 133018182771 2126000000 8833653 0
Austria Albania 1989 132785154342 2335124988 10116497 0
Austria Albania 1990 166062376740 2101624963 12125119 0
Austria Albania 1991 173375508073 1139166646 8548989 0
Austria Albania 1992 194608183696 709452584 4949408 1
Austria Albania 1993 189921096652 1228071038 6107622 1
我正在尝试估算政策假人对出口的影响，在时间段内使用假人并放弃常数：
> reg <-lm(log(exports)~0+log(gdphome)+log(gdppartner)+policy
+factor(year), data=fulldata)
R第一次下降虚拟（1988），但包括政策虚拟的两个级别的估计：
> coefs <- coef(reg)[1:5]
> coefs
log(gdphome) log(gdppartner) policy0
1.216703971 1.097967926 -31.447622433
policy1 factor(year)1989
-30.251306284 -0.009447245
但我想要的只是政策假人的一个估计和时间假人的所有估计。当你放弃常数时，似乎R包括模型的第一个虚拟的所有级别。如何告诉R使用第二个假人？
我可以更改模型中的顺序，但我不想跳过oder输出中的时间虚拟对象上的所有系数来得到我对策略变量的估计。 另一种方法是包括一个常数，但是时间假人的估计值的含义会改变。