Я хотел бы реализовать разложение Оахака в R. Он используется, например, в. экономика труда, чтобы отличить объясненную дисперсию от необъяснимой дисперсии.
Это довольно легко сделать с невзвешенными данными, используя пакет Oaxaca (см.
предыдущее объяснение здесь для общего использования оахака) . Однако пакет Oaxaca в настоящее время не поддерживает взвешенные данные опросов, такие как текущий опрос населения.
"Survey" является наиболее популярным пакетом для работы с данными опросов в R, но в нем отсутствует возможность прямого выполнения разложения Оахака.
Ниже приведен пример, который отмечает очевидные ограничения двух пакетов:
# Note the lack of support for a "Oaxaca decomposition command":
library(survey)
data(api)
# The line below weights the data
dclus2<-svydesign(id=~dnum+snum, weights=~pw, data=apiclus2)
model0<-svyglm(I(sch.wide=="Yes")~ell+meals+mobility, design=dclus2,
family=quasibinomial())
# Note the lack of support for survey weights:
library(oaxaca)
data("chicago")
# The line below will not work if my data is a survey.design object (i.e. weighted data)
oaxaca.results <- oaxaca(ln.real.wage ~ age + female + LTHS + some.college
+ college + advanced.degree | foreign.born, data = chicago, R = 50)
Если кто-либо может сказать мне, как сделать пакет oaxaca совместимым с весами опроса или как реализовать команду разложения oaxaca на объекте survey.design, то это будет очень полезно. Есть указатели?