我想在 R 中复制一个 Stata do.file（面板模型），但不幸的是我最终得到了错误的标准误差估计。这些数据是专有的，所以我不能在这里发布。使用的 Stata 代码如下所示：
xtreg Y X, vce(cluster countrycodeid) fe nonest dfadj
对于fe固定效果，nonest表明面板没有嵌套在集群中，并且dfadj由于发生了某种 DF 调整这一事实 - 目前尚无法确定哪种类型。
我的 R 代码看起来像这样，让我最终得到正确的系数值：
model <- plm(Y~X+as.factor(year),data=panel,model="within",index=c("codeid","year"))
现在是困难的部分，到目前为止我还没有找到解决方案，即使在尝试了多种标准误差稳健估计方法之后，例如广泛使用lmtest各种自由度变换方法。标准错误应该遵循国家-年份对模式（由countrycodeidStata代码中的变量捕获，其形式为codeid-year，因为某些变量似乎缺少数据，而这些数据不能按月获得。
有谁知道在使用不平衡面板和包装时是否需要牢记特殊技巧，plm()可以使用哪种 DF 调整，以及是否有可能按coeftest()国家/地区对函数中的数据进行分组？