Я хочу использовать модель с фиксированными эффектами для регрессии с одной переменной, являющейся переменной группы

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 неверно.

Чтение документации для felm () означает, что ваш код должен выглядеть примерно так:

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

но спецификации кода для него довольно сложны, в зависимости от того, хотите ли вы кластеризовать ваши стандартные ошибки и так далее.

Надеюсь, это поможет.