如果我想估计具有(区域)固定效应的线性概率模型,这与仅运行固定效应回归是否一样?也许我被这门语言绊倒了。我的目标是估计婴儿奖金的效果。我的因变量是NEWBORN的二进制指标,而我感兴趣的主要自变量是接收婴儿奖金的指标。我控制年龄,年龄平方,教育程度,婚姻状况和家庭收入。
我应该使用:
## 1.) Linear Probability
LPM <- lm(newborn ~ treatment + age + age_sq + highest_education + marital_stat +
hh_income_log, data=fertility_15_45)
或
## 2.) FE Model
FE_model <- plm(newborn ~ treatment + age + age_sq + highest_education + marital_stat +
hh_income_log, data = fertility_15_45, index="region", model="within")
答案 0 :(得分:1)
您可能希望在LPM中添加一个区域虚拟对象以获得区域固定效果。 示例:
library(plm)
data(Cigar)
summary(plm(I(sales > 121.2) ~ price + pop, data=Cigar, model="within", index="state"))$coe
# Estimate Std. Error t-value Pr(>|t|)
# price -2.880255e-03 2.626505e-04 -10.9661107 7.519348e-27
# pop -6.922327e-06 1.311006e-05 -0.5280165 5.975758e-01
summary(lm(I(sales > 121.2) ~ 0 + price + pop + factor(state), data=Cigar))$coe[1:2, ]
# Estimate Std. Error t-value Pr(>|t|)
# price -2.880255e-03 2.626505e-04 -10.9661107 7.519348e-27
# pop -6.922327e-06 1.311006e-05 -0.5280165 5.975758e-01