Я работал над проектом, который рассматривает разрыв в заработной плате между двумя группами идентичностей (каст) в Индии. Я использую пакет Oaxaca, который является инструментом R для того же самого. Несмотря на то, что я запускаю сам код Oaxaca, я не могу извлечь устойчивые стандартные ошибки, сгруппированные на уровне округа.
В обычной линейной регрессии это делается с использованием coeftest:
coeftest(linear.model, vcov. = vcovHC(linear.model,type = "HC1")).
В настоящее время мой код оахака имеет вид
R <- oaxaca(y ~ x1 + x2+|z, data = data)
Каким может быть решение этой проблемы?