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

# 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/
```