包"prais"
包含函数prais.winsten
,用于使用 Prais Winsten 估计器运行回归模型。但是,stargazer
和broom
包似乎都不适用于prais.winsten
函数的结果。
拟合表单模型后
pw<- prais.winsten(speed ~ dist, cars)
当我尝试使用观星者报告结果时,出现以下错误:
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]
有什么方法可以在stargazer
和broom
中使用prais.wisten
结果?
或者,有没有办法在glm
运行 Prais Winsten 估算器?
最简单的方法是编辑 Prais.Winsten 函数以返回一个 lm 对象。
getAnywhere(prais.winsten())
将返回该函数。然后,您可以将代码复制到脚本中,只需将返回函数从"结果"更改为"lm"。