对于plm对象，如何在R中获得多路聚类标准错误，其中聚类不在面板的时间/组ID级别？
包plm支持使用函数plm::vcovHC计算集群健壮的标准错误。不幸的是，此功能仅支持组的聚类或面板的时间ID。在某些情况下，人们会希望将标准误差聚集在与面板观察单位不同的水平上。一个例子是具有单个固定效应的回归，其中较高聚合水平的变量被用作独立变量。
2014年提出类似的question，建议使用自助功能。这个问题不同，因为（1）我想集中两个变量，而不仅仅是一个和（2）我宁愿不使用引导程序。
包multiwaycov提供了与want非常接近的内容，但遗憾的是它不支持plm对象。