Bookmark and Share

Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.


[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

Re: st: loop


From   Chiara Mussida <[email protected]>
To   [email protected]
Subject   Re: st: loop
Date   Fri, 30 Mar 2012 10:08:41 +0200

Thanks Kit,

margins, by(married occupation)
gives me the product between the Betas of the estimated wage reg and
the mean value of characteristics for each occupation, for married and
unmarried. Right?

Chiara

On 29/03/2012, Christopher Baum <[email protected]> wrote:
> <>
> Chiara said
> I wrote a loop for wage regressions to obtain the predicted wages if
> men are paid as men at each occupation:
>
> reg lwage $x if fem==0
> predict pip if fem==0, xb
> forvalues k = 1 2 to 7 {
> 	predict pip`k' if fem==0 & occ`k'==1, xb
> }
> summarize pip, meanonly
> scalar xMbM=r(mean) /*Predicted wages if men are paid as men*/
> forvalues k = 1 2 to 7 {
> 	summarize pip`k', meanonly
> 	scalar xMbM`k'=r(mean) /*Predicted wages if men are paid as men at
> each occupation*/
> }
>
> I now want to get the xMbM (and also xFbF for females) by occupation.
> For each occupation I want the product between ythe mean individual
> charachteristics (xM) and the coefficients (bM). I tried with this
> loop, but i'm not sure it's the correct one:
>
> reg lwage $x if fem==0
> predict pip if occupation==1, xb
> forvalues k = 1  {
> 	predict pip`k' if fem==0 & occ`k'==1, xb
> }
> summarize pip, meanonly
> scalar xMbM=r(mean) /*Predicted wages if men are paid as men*/
> forvalues k = 1 {
> 	summarize pip`k', meanonly
> 	scalar xMbM`k'=r(mean) /*Predicted wages if men are paid as men at
> each occupation*/
> }
>
> More precisely, i'm not sure this wille give me each product between
> mean values of characteristics (x, for M or F) and estomated wage
> equation coefficients (b, for M or F) at each 1...7 occupation.
>
>
>
>
> No need for all this manual labor...
>
> --------------------
> webuse nlsw88,clear
> // lets treat marital status as equivalent of gender, since all people here
> are women
> // run regression over all cases so can generate pred wage by mar.stat. and
> occup.
> reg wage age collgrad south i.occupation
> margins, by(married occupation)
> marginsplot,graph(married) xlab(,angle(90))
> ---------------------
>
> Kit
>
>
> Kit Baum   |   Boston College Economics & DIW Berlin   |
> http://ideas.repec.org/e/pba1.html
>                              An Introduction to Stata Programming  |
> http://www.stata-press.com/books/isp.html
>   An Introduction to Modern Econometrics Using Stata  |
> http://www.stata-press.com/books/imeus.html
>
>
> *
> *   For searches and help try:
> *   http://www.stata.com/help.cgi?search
> *   http://www.stata.com/support/statalist/faq
> *   http://www.ats.ucla.edu/stat/stata/
>


-- 
Chiara Mussida
PhD candidate
Doctoral school of Economic Policy
Catholic University, Piacenza (Italy)
*
*   For searches and help try:
*   http://www.stata.com/help.cgi?search
*   http://www.stata.com/support/statalist/faq
*   http://www.ats.ucla.edu/stat/stata/


© Copyright 1996–2018 StataCorp LLC   |   Terms of use   |   Privacy   |   Contact us   |   Site index