我想在R中复制Stata do.file（面板模型），但不幸的是我最终得到了错误的标准误差估计。数据是专有的，所以我不能在这里发布。使用的Stata代码如下所示：
xtreg Y X, vce(cluster countrycodeid) fe nonest dfadj
使用fe表示固定效果，nonest表示面板未嵌套在群集中，而dfadj表示某种DF调整发生 - 不可能找出现在的那种。
我的R-Code看起来像这样，让我最终得到了正确的系数值：
model <- plm(Y~X+as.factor(year),data=panel,model="within",index=c("codeid","year"))
现在出现了困难的部分，即使在尝试了多种标准误差稳健估计方法之后，我还没有找到解决方案，例如广泛使用lmtest和不同程度自由变换方法。标准错误应该遵循国家 - 年对模式（由Stata代码中的变量countrycodeid捕获，其采用形式为codeid-year，因为对于某些不可用的变量，似乎缺少数据每月一次。
在使用不平衡面板和plm()包时，有没有人知道是否有特殊的技巧可以使用，可以使用哪种DF调整，以及是否有可能在coeftest()以国家/年为基础运作？
答案 0 :(得分：1)
这不是一个完整的答案。
Stata使用this post中描述的有限样本校正。我认为这可能会使您的标准错误更加接近。
此外，您可以通过发出help whatsnew9了解有关nonest / dfadj的更多信息。当指定cluster()选项时，Stata用于调整内部转换的VCE。除非指定dfadj，否则不再调整集群强大的VCE。您可能需要使用版本控制来复制旧的估计值。