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

I want to use a fixed effects model on a regression with one variable being the group variable

我正在使用 felm() 并且代码在所有模型上都是 运行ning……但我只需要它在状态上 运行……问题是 "Estimate the model using fixed effects (FE) at the state level"。使用 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 中是正确的。

看起来 felm() 是当你有多个分组变量时,但听起来你只使用一个分组变量来获得固定效果? (即状态)。

你应该得到与

相同的正确结果
mod3 <- lm(totfatrte ~ D81 + state, data = DRIVING)

此外,如果 stata 和 R 之间的系数或标准误差不一致,那并不一定意味着 R 是错误的。

阅读 documentation for felm() 表明您的代码应该更像这样:

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

但是它的代码规范非常复杂,具体取决于您是否要对标准错误进行聚类等等。

希望这对您有所帮助。