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/

また、probemodCRAN のパッケージも役立ちます。