如何在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，如果发现任何错误，请随时提出问题（{{3} }）。

Answer 2

也许这里的套餐派上用场了？ http://www.strengejacke.de/sjPlot/sjp.int/

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