我想在回归上使用固定效应模型，其中一个变量是组变量

Question

我正在使用felm（）并且代码正在所有模型上运行……但是我只需要它在状态下运行……问题问“在状态级别使用固定效果（FE）估计模型”。使用felm（）不能给我正确的结果，因为我不知道是否需要将状态包括为因变量（没有给我正确的答案）或如何指定一个变量需要为组变量（我假设这是获得准确结果的方法。

我尝试使用

plm(ind~depvar+state,data=data, model='within')

我尝试使用

felm(ind~depvar+state,data=data)


FELinMod3<-felm(DRIVING$totfatrte~DRIVING$D81+DRIVING$state, data=DRIVING)

FELinMod3<-plm(DRIVING$totfatrte~DRIVING$D81+DRIVING$state, data=DRIVING, model='within')

输出给我的系数与我在STATA中正确的系数不正确。

Answer 1

看起来felm()适用于具有多个分组变量的情况，但是听起来好像您只对固定效果使用了一个分组变量？ （即状态）。

您应该获得

相同的正确结果

mod3 <- lm(totfatrte ~ D81 + state, data = DRIVING)

此外，如果Stata和R之间的系数或标准误差不一致，则不一定表示R是错误的。

读取documentation for felm()表示您的代码应更像这样：

model3<-felm(totfatrte ~ D81 | state, data = DRIVING)

，但是根据是否要对标准错误进行聚类等，其代码规范非常复杂。

希望这会有所帮助。