Bookmark and Share

Notice: On March 31, it was announced that Statalist is moving from an email list to a forum. The old list will shut down at the end of May, and its replacement, statalist.org is already up and running.


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

RE: st: Summarizing properties of other group members meeting specified conditions


From   Erik Aadland <erikaadland@hotmail.com>
To   <statalist@hsphsun2.harvard.edu>
Subject   RE: st: Summarizing properties of other group members meeting specified conditions
Date   Tue, 28 Jun 2011 11:27:45 +0000

I modified the suggested code a little bit, moving the cond_2 == 1 to the replace line, and now it almost works as intended.
 
The modified code is:
 
sort project_id org_id 
 
gen maxwins = . 
su org_id, meanonly 
quietly forvalues i = 1/`r(max)' { 
    gen include = 1 if org_id != `i' & cond_1 == 1 
    egen x = max(var_2 * include), by(project_id) 
    replace maxwins = x if org_id == `i' & cond_2 == 1 
    drop include x 
} 

However, I do get a strange result when I have more than one org_id with cond_1 == 1 in the same project_id.
 
Here is the resulting output:
 
project_id    org_id    cond_1    cond_2    cond_3    var_1    var_2    maxwins
1             1         0         1         0         3        0        3
1             2         1         0         0         0        3        .
1             3         0         0         1         0        1        . 
2             4         0         1         0         2        0        3
2             4         1         0         0         0        4        . 
2             5         1         0         0         0        3        . 
2             6         0         1         0         1        0        4
3             8         0         1         0         1        0        . 

 
Why does org_id 4 in project_id 2 get a maxwins value of 3 when the max var_2 value for org_id with cond_1 == 1 is 4?
And why does org_id 6 in project_id 2 get the correct maxwins value of 4 when org_id 4 does not?
 
Kind regards,
 
Erik Aadland.
 
 
 


----------------------------------------
> Date: Tue, 28 Jun 2011 11:32:53 +0100
> Subject: Re: st: Summarizing properties of other group members meeting specified conditions
> From: njcoxstata@gmail.com
> To: statalist@hsphsun2.harvard.edu
>
> I guess that your original typo remains in
>
> gen include 1 = if org_id != `i' & cond_1 == 1
>
> which should be more like
>
> gen include = 1 if org_id != `i' & cond_1 == 1
>
> Warning: I am not sure that I have grasped your verbal summary of what
> you want. My emphasis is on showing you Stata techniques.
>
> Nick
>
> On Tue, Jun 28, 2011 at 11:16 AM, Erik Aadland <erikaadland@hotmail.com> wrote:
> > Thank you so much for your help and input, Nick. I sincerely appreciate it.
> >
> > When testing the suggested code, after having sorted on project_id and org_id, I get an error message for the following line of code:
> >
> > quietly forvalues i = 1/`r(max)' {
> >
> > The error message is:
> >
> > 1 invalid name
> > r(198);
> >
> > How to get around this problem?
> >
> >
> > Kind regards,
> >
> > Erik Aadland.
> >
> >
> >
> >
> >> Date: Tue, 28 Jun 2011 08:20:57 +0100
> >> Subject: Re: st: Summarizing properties of other group members meeting specified conditions
> >> From: njcoxstata@gmail.com
> >> To: statalist@hsphsun2.harvard.edu
> >>
> >> There was another typo in the line I corrected (you can't have two
> >> -if-s) and it is not clear where your r(max) comes from. This may be
> >> closer, but I've not tested it.
> >>
> >> gen maxwins = .
> >> su org_id, meanonly
> >> quietly forvalues i = 1/`r(max)' {
> >> gen include 1 = if org_id != `i' & cond_1 == 1
> >> egen x = max(var_2 * include) if cond_2 == 1, by(project_id)
> >> replace maxwins = x if org_id == `i'
> >> drop include x
> >> }
> >>
> >>
> >>
> >> On Mon, Jun 27, 2011 at 3:26 PM, Erik Aadland <erikaadland@hotmail.com> wrote:
> >> > Dear statalist.
> >> >
> >> > After having read the following FAQ http://www.stata.com/support/faqs/data/members.html and experimented with solutions, I am still struggling.
> >> >
> >> > For each project_id, I am trying to generate a variable that for each org_id that meets a specified condition generates the highest observed value on a given variable for another org_id that meets a different specified condition. I think an example is in order:
> >> >
> >> > Example structure:
> >> >
> >> > project_id org_id cond_1 cond_2 cond_3 var_1 var_2
> >> > 1 1 0 1 0 3 0
> >> > 1 2 1 0 0 0 3
> >> > 1 3 0 0 1 0 1
> >> > 2 4 0 1 0 2 0
> >> > 2 4 1 0 0 0 4
> >> > 2 5 1 0 0 0 3
> >> > 2 6 0 1 0 1 0
> >> > 3 8 0 1 0 1 0
> >> >
> >> >
> >> > For instance, for each project_id I try generate a variable that for each org_id meeting cond_2 == 1 produces the highest value of var_2 for org_id that meets cond_1 == 1. In the case of project_id 1 above org_id 1 should get a score of 3, while org_id 2 and 3 should get a missing value. Some project_ids may have only one observation (consist of a single org_id) as in project_id 3 above, or may not have other org_ids represented that meet the specified condition. Then, it is of course not possible to generate a value on the variable.
> >> >
> >> > I have tried with the following code, but I cannot figure out how to specify cond_2 == 1, and it comes out all wrong (Stata does not even accept my code):
> >> >
> >> > sort project_id org_id
> >> > gen maxwins = .
> >> > quietly forvalues i = 1/`r(max)' {
> >> > gen include 1 = if org_id != `i' if cond_1 == 1
> >> > egen x = max(var_2 * include), by(project_id)
> >> > replace maxwins = x if org_id == `i'
> >> > drop include x
> >> > }
> >> >
> >> > Is it possible to do this without loops also?
> >>
> >> *
> >> * 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/
> >
>
> *
> * 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–2014 StataCorp LP   |   Terms of use   |   Privacy   |   Contact us   |   Site index