Как не статистик, я достиг своего лимита здесь:
Я пытаюсь подогнать модель Пуассона для данных панели (используя pglm) и хочу вычислить устойчивые стандартные ошибки (используя lmtest).
Мой код в настоящее время выглядит так:
#poisson model (panel with year fixed effects):
poisson_model <- pglm(y ~ a + b + c + factor(year), data = regression_data,
model = "pooling", family = poisson, index = c("ID", "year"))
#robust standard errors:
robust_SE_model <- coeftest(poisson_model, vcov. = vcovHC(poisson_model, type = "HC1"))
Этот код прекрасно работает для одной из моих других спецификаций модели, когда я подгоняю обычную модель панели с plm, но когда я пробую модель Пуассона с pglm, я получаю следующее сообщение об ошибке:
Error in terms.default(object) : no terms component nor attribute
Это из-за ограничения пакета lmtest или я здесь ошибаюсь? Я действительно надеюсь, что смогу решить эту проблему с помощью пакетов (не обязательно pglm и lmtest) и не буду погружаться в ручной расчет надежных ошибок.
Любая помощь высоко ценится!