r - 如何在 r 中可视化面板模型的交互项？

Question

我的目标是可视化固定效应模型的交互项。我尝试了“visreg”和“ggplot2”包，但似乎不可能。

我怎样才能做到这一点？

提前致谢。

示例数据集和面板模型：

library(plm)
data(EmplUK)

m1 <- plm(emp ~ wage + capital + lag(output,1) + capital:lag(output,1), 
          data = EmplUK, model="within", effect="twoway", index=c("firm","year"))

Answer 1

在sjPlot的当前 CRAN 版本中，该sjp.int函数不支持plm对象。

但是，我已经在当前的GitHub-version中实现了它。

以上面的示例为例，以下代码：

sjp.int(m1)

产生以下情节：

我不确定是否适合所有类型的设计/模型plm，sjp.int如果您发现任何错误，请随时提出问题（在 GitHub 上）。

Answer 2

也许这里的包裹会派上用场？http://www.strengejacke.de/sjPlot/sjp.int/

CRAN的probemod软件包也可以提供帮助。