Плоские эффекты переменных в условиях взаимодействия

Question

Я хотел бы построить эффекты переменных в условиях взаимодействия, используя данные панели и модель FE.

У меня есть различные эффекты взаимодействия в моем уравнении, например, это один здесь:

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

Я могу только найти решения для фильма, но не для ППСТЗА.

Таким образом, по оси x должен быть PRIVCR, а по оси Y - влияние ПИИ на рост.

Благодарим за помощь!

Лиза

Answer 1

Я не знаю о пакете, который поддерживает plm объекты непосредственно. Поскольку вы запрашиваете FE-модели, вы можете просто взять подход LSDV для FE и сделать оценку на lm, чтобы получить объект lm, который работает с пакетом effects. Вот пример данных 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)