如何在R中绘制plm中固定效应模型的效应？

Question

在plm中，我无法得到固定效果模型的效果图.我试过使用effect()、predict()和各种包(如sjPlot等)。有没有一种方法来绘制它，特别是使用交互？

我总是收到这样的错误消息：

Error in mod.matrix %*% scoef : non-conformable arguments

Answer 1

试试fixef？例如，见下文：

plm_2 <- plm(wealth ~  Volatility, data = ds_panel,index=c("rho"), model = "within")
y1 <- fixef(plm_2)
x1 <- as.numeric(names(y1))
plot(y1~x1, pch = 20, ylab = "FE", xlab = expression(rho))