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

【问题标题】：How to visualize an interaction term of a panel model in r?如何在 r 中可视化面板模型的交互项？
【发布时间】：2012-10-30 18:01:53
【问题描述】：

我的目标是可视化固定效应模型的交互项。我尝试了“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"))

【问题讨论】：

“可视化条件效果”到底是什么意思？您想可视化模型对数据的拟合吗？年度固定效应？公司固定效应？
@RoyalTS - 抱歉不清楚，我想可视化交互项 (capital:lag(output,1))。
如果有人感兴趣：我使用了线性模型 (lm) 并为 id 变量添加了假人。使用“visreg”的可视化效果很好。

标签： r visualization regression panel-data plm

【解决方案1】：

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

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

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

sjp.int(m1)

产生以下情节：

我不确定所有带有plm 的设计/模型是否适用于sjp.int，如果您发现任何错误，请随时提出问题 (at GitHub)。

【讨论】：

【解决方案2】：

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

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

【讨论】：