交互术语中变量的绘图效果

Question

我想使用面板数据和有限元模型在交互术语中绘制变量的效果。

我的等式中有各种交互效果，例如这个：

FIXED1 <- plm(GDPPCgrowth ~ FDI * PRIVCR, data = dfp)

我只能找到lm的解决方案，但不能找到plm的解决方案。

因此，在x轴上应该有PRIVCR，在y轴上有外国直接投资对增长的影响。

感谢您的帮助！

莉莎

Answer 1

我不知道直接支持plm个对象的包。当您要求FE模型时，您可以采用LSDV方法进行FE并按lm进行估算，以获得与effects包一起使用的lm对象。以下是Grunfeld数据的示例：

library(plm)
library(effects)
data("Grunfeld", package = "plm")

mod_fe <- plm(inv ~ value + capital + value:capital, data = Grunfeld, model = "within")
Grunfeld[ , "firm"] <- factor(Grunfeld[ , "firm"]) # needs to be factor in the data NOT in the formula [required by package effects]
mod_lsdv <- lm(inv ~ value + capital + value:capital + firm, data = Grunfeld)
coefficients(mod_fe)    # estimates are the same
coefficients(mod_lsdv)  # estimates are the same

eff_obj <- effects::Effect(c("value", "capital"), mod_lsdv)
plot(eff_obj)