Внедрение 2SLS с фиксированными эффектами и эндогенным термином взаимодействия И весами в R

Question

Я хочу смоделировать следующее в R:

outcome = beta Var1 + beta Var2+ beta Var1:Var2+ controls + county FE + year FE

У меня есть два инструмента для Var1. Я также должен весить население графства. Что мне нужно, это запустить регрессию 2SLS с двумя инструментами для Var1, с фиксированными эффектами округа и года, и все они будут взвешены по населению округа.

Пакет felm, похоже, не позволяет мне использовать термин взаимодействия. plm не реализовал веса для анализа 2SLS, поэтому я не могу использовать этот пакет и вес по населению округа.

Мои вопросы:

1. Может кто-нибудь порекомендовать? пакет, который позволяет мне использовать термин взаимодействия, включая мои два фиксированных эффекта И вес по населению округа?

2. Если нет, то есть ли простой способ исправить мои стандартные ошибки 2SLS если я использую felm, но использую предсказанные значения из регрессии первого этапа вместо собственных вычислений felm 2SLS?

Любая помощь приветствуется! Спасибо!

ETA: Я тоже только что попробовал ivreg из пакета AER. Я использовал эту команду:

test <- ivreg(data = mydata, outcome~ Var1 + Var2 + 
                Var1:Var2 + Control1+ Control2 + Control3 + Control4 +
                Control5 + Control6 + Control7+ 
                Control8 + as.factor(FIPS) + as.factor(Year) | Var2+ 
                Control1 + Control2 + Control3 + Control4 +
                Control5 + Control6 + Control7 + 
                Control8 + Instrument1+ Instrument2, 
                weights = mydata$population)

Это привело меня к ошибке:

In ivreg.fit(X, Y, Z, weights, offset, ...) :
  more regressors than instruments