如何在不同级别而不是id或时间上对plm的标准错误进行聚类？

Question

我通过plm进行面板回归。 plm只能将标准错误归类为“组”或“时间”级别。但是我想集中在国家一级，而不是小组或时间上。 （根据我的数据，组是公司级别。）

我使用以下代码进行集群： summary(r, vcov = function(x) vcovHC(x, cluster = "group", method = "arellano")) 但是cluster的选项仅是组和时间。我想聚集在国家一级。