Использование Stargazer с Zelig
Я пытаюсь использовать stargazer (версия 5.2) для вывода стандартных сводных () результатов регрессии, оцененной с помощью Zelig (версия 5.0.13). Тем не менее, я получаю ошибку
Error in envRefInferField(x, what, getClass(class(x)), selfEnv) :
‘result’ is not a valid field or method name for reference class “Zelig-ls”
при попытке передать оценочную модель в Stargazer.
Минимальный рабочий пример воспроизведения ошибки будет
library(Zelig)
library(stargazer)
data("swiss")
z <- zls$new()
z$zelig(Fertility~Education, data=swiss)
stargazer(z)
Я предполагаю, что мне нужно каким-то образом получить доступ к месту, где оценочная модель хранится в объекте zelig, но у меня возникают проблемы с его поиском.
1 ответ
Решение
Начиная с Zelig 5.0-13 и stargazer 5.2, stargazer, похоже, не в состоянии справиться с внутренней моделью Zelig. В качестве обходного пути вы можете просто привести его к "нормальному" lm, чтобы он работал:
stargazer(lm(z$zelig.out$z.out[[1]]), type='text')
Результат:
===============================================
Dependent variable:
---------------------------
Fertility
-----------------------------------------------
Education -0.862***
(0.145)
Constant 79.610***
(2.104)
-----------------------------------------------
Observations 47
R2 0.441
Adjusted R2 0.428
Residual Std. Error 9.446 (df = 45)
F Statistic 35.446*** (df = 1; 45)
===============================================
Note: *p<0.1; **p<0.05; ***p<0.01