有没有人对如何使用面板校正标准误差在 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 创建边际效应图?