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:49:04 +0000

I'll answer my own question.
Org_id 4 gets a maxvalue of 3 because it meets both cond_1 and cond_2.
 
The code therefore appears to be correct now, until I stand corrected again.
 
Thanks again.
 
 
All the best,
 
Erik.
 


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