как работать с prais.результаты winsten в stargazer and broom (r)


Пакет "prais" содержит функцию prais.winsten для запуска регрессионных моделей с оценкой Prais Winsten. Однако ни stargazer, ни broom пакеты, по-видимому, не работают с результатами функции prais.winsten.

После подгонки модели вида

pw<- prais.winsten(speed ~ dist, cars)

Когда я пытаюсь сообщить результаты с помощью stargazer, я получаю следующую ошибку:

stargazer(pw, out = "pw.html") 
Error: $ operator is invalid for atomic vectors

И аналогично очищая модель с помощью пакета broom я получаю

tidy(pw)
Error in names(object) <- nm : 
  'names' attribute [1] must be the same length as the vector [0]

Есть ли способ работать с результатами prais.wisten в stargazer и broom?

Или, в качестве альтернативы, есть ли способ запустить оценку Prais Winsten в glm?

1 2

1 ответ:

Самый простой способ сделать это-отредактировать Prais.Функция Winsten для возврата объекта lm.

getAnywhere(prais.winsten())вернет функцию. Затем вы можете скопировать код в свой скрипт и просто изменить функцию возврата с "results" на "lm".