r - 从固定效应模型 (plm) 绘制交互项

Question

我正在尝试从固定效应模型中绘制交互效应。出于某种原因，我总是收到相同的错误消息：Error in crossprod(beta, t(X)) : non-conformable arguments

首先，我认为这只是我的数据（可能是多重共线性），但是在使用 plm 包中包含的数据时会发生同样的问题。这是一个例子：

library(plm)
library(sjPlot)

data("Grunfeld", package="plm")
grun.fe <- plm(inv~value+capital, data = Grunfeld, model = "within")

plot_model(grun.fe, type = "pred", terms = c("value","capital")) 

事实上，如果我想使用任何估计的固定效应模型（例如使用predict）进行任何类型的预测，我会收到相同的错误消息。随机效应和池化面板模型工作正常，只有固定效应。

我想知道这是否与 plm 存储固定效果结果矩阵的方式有关，但在 plm 手册中找不到任何帮助。相同的任务在 Stata 中使用 xtreg 运行良好，例如使用边距和边距图。任何帮助将非常感激！

标签： rinteractionpanel-dataplm