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: Re: calculating proportions


From   Ekaterina Hertog <[email protected]>
To   "[email protected]" <[email protected]>
Subject   Re: st: Re: calculating proportions
Date   Tue, 21 Feb 2012 18:20:07 +0000

Dear Nick,
thanks a lot for the pointer.
I must have had N_female defined earlier. Once I dropped all macros, it worked,
Ekaterina
On 21/02/2012 17:33, Nick Cox wrote:
There are various possible reasons for a problem here with what
Ekaterina did, including

1. No local macro N_female was defined by her earlier.

2. One such was defined by her but is not visible to Stata at that point.

The use of local macros looks fine in the code posted by Pathmes and Oliver.

Nicl

On Tue, Feb 21, 2012 at 5:22 PM, Ekaterina Hertog
<[email protected]>  wrote:
Dear Pathmes and Oliver,
thanks a lot, great code and it solves my problems! And I learned more than
3 things from it
Just a small hitch:
converting numbers to proportions syntax:

gen fap=fa/`N_female'

gives an error message: 'invalid syntax'.
I went around it by creating a variable which is a constant and equals a
total number of women in the population and using it instead of the local
macro. I clumsy way to do it.
katya

On 21/02/2012 16:45, Pathmeswaran wrote:
Thanks Oliver.

With your help I was able to learn 3 things on my first day in statalist!

   forvalues
   `r(N)'
   indentation in foreach/ forvalue loops

*************** example code begins **************************

version 8
clear all
set obs 7

gen byte id = _n

gen byte sex = 1
replace sex = 0 in 4/7
label def lbl_sex 0 "male" 1 "female"
label values sex lbl_sex

gen int income = 7
replace income = 10 in 3
replace income = 0 in 4
replace income = 7 in 5
replace income = 9 in 6
replace income = 11 in 7

list

* I have just copied Oliver's code

sort sex
gen idm=_n if sex==0
* idm is a unique number for each male

gen fa=.
label variable fa "Number of females with higher income"

gen fb=.
label variable fb "Number of females with lower income"

gen fe=.
label variable fe "Number of females with equal income"

* The variables fa, fb&    fe are generated individually for each male
in the following loop
quietly{
        count if sex == 0
        local N_male = `r(N)'

        count if sex == 1
        local N_female = `r(N)'

        forvalues i = 1/`N_male' {
                egen int feabove = rank(income) if sex==1 | (sex==0&
  idm==`i'), field
                replace fa= feabove - 1 if idm==`i'
                drop feabove

                egen int febelow = rank(income) if sex==1 | (sex==0&
  idm==`i'), track
                replace fb=  febelow - 1 if idm==`i'
                drop febelow

                replace fe=`N_female' - ( fa+ fb) if idm==`i'
        }
}

* Converting the numbers to proportions

gen fap=fa/`N_female'
label variable fap "Proportion of females with higher income"

gen fbp=fb/`N_female'
label variable fbp "Proportion of females with lower income"

gen fep=fe/`N_female'
label variable fep "Proportion of females with equal income"

**************** example code ends ***************************


*
*   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/

*
*   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