R中加权调查数据的瓦哈卡分解

Question

我想在R中实现Oaxaca分解。劳动经济学来区分解释性方差和无法解释的方差。

使用Oaxaca package (See previous explanation here for general oaxaca usage)处理未加权的数据非常容易。但是，瓦哈卡州一揽子计划目前不支持加权调查数据，例如当前人口调查。

"Survey"是用于处理R中调查数据的最受欢迎的软件包，但是它缺乏直接执行Oaxaca分解的功能。

下面是一个示例，指出了两个软件包的明显局限性：

# 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软件包与调查权重兼容，或者如何对Survey.design对象实施oaxaca分解命令，那么将不胜感激。有指针吗？

Answer 1

“ decr”库中有一些脚本可以对加权数据进行加权和分解。