我一直在研究一个项目，该项目研究印度两个身份（种姓）群体之间的工资差距。我正在使用Oaxaca软件包，它是R工具。在运行Oaxaca代码本身的同时，我无法提取聚集在“区域”级别上的“健壮”标准错误。
在正常的线性回归中，这是使用coeftest完成的：
coeftest(linear.model, vcov. = vcovHC(linear.model,type = "HC1")).
当前，我的瓦哈卡代码的格式为
R <- oaxaca(y ~ x1 + x2+|z, data = data)
该问题的解决方案是什么？