有没有人对如何使用面板校正的标准错误在R中制作边际效应图有任何建议？
为了估计R中面板校正的标准误差，我使用plm和lmtest包。
首先我用plm估计回归模型：
plm <- plm(Y ~ X*Z, data = a, model='pooling', index=c('cicode', 'year'))
然后我使用lmtest包通过coeftest传递这个plm对象：
plm2 <- coeftest(plm, vcov=function(x) vcovBK(x,type="HC1", cluster="time"))
plm2包含系数列表和面板校正的标准误差。我可以使用什么包和代码来创建使用plm2的边际效果图？