我想在 R 中进行瓦哈卡分解。它用于例如我相信,劳动经济学可以区分可解释的方差和不可解释的方差。我一直无法在 R 中找到合适的解决方案,而且我不太愿意自己创建一个解决方案(我可能会把它搞砸)。
无论如何,这里简要解释一下该过程:
http://en.wikipedia.org/wiki/Ronald_Oaxaca http://en.wikipedia.org/wiki/Ronald_Oaxaca
Stata 对此有一个相当好的软件包,但 Stata 对我来说并不容易获得。
www.stata.com/meeting/5german/SINNING_stata_presentation.pdf http://www.stata.com/meeting/5german/SINNING_stata_presentation.pdf
请注意:我也在R-help上发了消息,但没有得到回复。我希望也可以在这个列表上发布。
提前致谢,
拉斯穆斯
编辑:我做了以下函数,这似乎会产生错误的答案(呃)。我尝试按照上面的 Stata 链接进行操作,但没有达到我的预期:)
oaxaca <- function (fsex,frace1,frace2) {
## First we make regresions
data1 <- subset(l2,sex==fsex & race==frace1)
data2 <- subset(l2,sex==fsex & race==frace2)
mindata1 <- subset(cbind(grade,exp,I(exp^2)),sex==fsex & race==frace1)
mindata2 <- subset(cbind(grade,exp,I(exp^2)),sex==fsex & race==frace2)
reg1 <- lm(log(wage)~grade+exp+I(exp^2), data=data1)
reg2 <- lm(log(wage)~grade+exp+I(exp^2), data=data2)
## DECOMPOSITION
################
## Variables
gap <- mean(log(wage[race==frace1 & sex==fsex]))-mean(log(wage[race==frace2 & sex==fsex]))
mean1 <- colMeans(mindata1)
mean2 <- colMeans(mindata2)
beta1 <- summary(reg1)$coefficients[,1]
beta2 <- summary(reg2)$coefficients[,1]
beta1incep <- summary(reg1)$coefficients[1,1]
beta2incep <- summary(reg2)$coefficients[1,1]
beta1coef <- summary(reg1)$coefficients[c(2,3,4),1]
beta2coef <- summary(reg2)$coefficients[c(2,3,4),1]
betastar <- .5*(beta1coef+beta2coef)
betastar2 <- (beta1+beta2)/2
expl <- sum((mean1-mean2)*beta1coef)
uexpl <- sum(mean2*(beta2coef-beta1coef))
pct=expl/gap
pct2=uexpl/gap
## output
out <- data.frame(Gap=gap,
Explained=expl,
Unexplained=uexpl,
Pct=pct*100)
return(out)
}