2012-11-10 48 views
0

我的目標是可視化固定效果模型的交互項。我試過包'visreg'和'ggplot2',但似乎不可能。如何可視化r中面板模型的交互項?

我該怎麼做?

在此先感謝。

示例數據集和麪板型號:

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")) 
+0

你究竟是什麼意思與「可視化的條件效應」?您是否想要將模型擬合成數據?一年固定效應?公司的固定效應? – RoyalTS

+0

@RoyalTS - 抱歉不清楚,我想將交互項(資本:lag(output,1))可視化。 – tfr

+0

如果有人感興趣:我使用線性模型(lm)併爲id變量添加了虛擬變量。使用'visreg'的可視化效果很好。 – tfr

回答

1

在當前的CRAN版本sjPlot中,sjp.int函數不支持plm對象。

但是,我在當前的GitHub-version中實現了它。

考慮從你上面的例子,下面的代碼:下面的圖形

sjp.int(m1)

生產:

enter image description here

我不知道是否配備了所有類型的設計/模型plmsjp.int一起工作,如果您發現任何錯誤,請隨時提出問題(at GitHub)。

相關問題