Предположим, я использую lfe пакет felm для учета очень большого количества фиксированных эффектов.
Y = b0 + b1*X1 + fixed effects
Y = b0 + b1*X1 + b2*X2 + b3*X3 + fixed effects
model1 = felm(Y~X1| fixedeffects, data = data)
model1 = felm(Y~X1 + X2 + X3| fixedeffects, data = data)
Стандартные lm приемы с объектами, такие как anova(), не работают.Я мог бы запрограммировать это сам, но мне просто любопытно.