Декомпозиция Оахаки в R [закрыто]

Закрыт . Этот вопрос должен быть более . В настоящее время он не принимает ответы.

Хотите улучшить этот вопрос? Обновите вопрос, чтобы он фокусировался только на одной проблеме, отредактировав это сообщение .

Закрыт 8 месяцев назад .

Я хотел бы провести декомпозицию Оахаки в R. Я полагаю, что он используется, например, в экономике труда, чтобы отличить объясненную дисперсию от необъяснимой. Мне не удалось найти подходящее решение в R, и я довольно неохотно создаю его сам (я бы, вероятно, все испортил).

В любом случае, процедура кратко объясняется здесь:

http: // en.wikipedia.org/wiki/Ronald_Oaxaca

Stata наделена довольно хорошим пакетом для этого, но Stata мне нелегко получить.

www.stata.com/meeting/5german/SINNING_stata_presentation.pdf

Обратите внимание: я также отправил сообщение в R-help, но оно не получило ответа. Я надеюсь, что можно и разместить в этом списке.

Заранее спасибо, stata.com/meeting/5german/SINNING_stata_presentation.pdf

Обратите внимание: я также разместил сообщение в R-help, но ответа на него не было. Я надеюсь, что можно и разместить в этом списке.

Заранее спасибо, stata.com/meeting/5german/SINNING_stata_presentation.pdf

Заранее спасибо, Расмус

Редактировать: Я сделал следующую функцию, которая, кажется, дает неправильные ответы (ура). Я попытался перейти по ссылке 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)
 }

r stata economics

задан Metrics 25 August 2013 в 16:54

0 ответов

Другие вопросы по тегам:

r stata economics

Декомпозиция Оахаки в R [закрыто]

0 ответов

Похожие вопросы: