Notice: On April 23, 2014, Statalist moved from an email list to a forum, based at statalist.org.
From | <isabel.lechner@vetsuisse.unibe.ch> |
To | <statalist@hsphsun2.harvard.edu> |
Subject | RE: st: Group identifiers within group |
Date | Fri, 7 Jun 2013 14:08:11 +0000 |
Thanks Nick, that's exactly what I was looking for! And I guess lakes and glaciers still wouldn't do it, so maybe adding swiss cheese holes would help ;) Isabel -----Original Message----- From: owner-statalist@hsphsun2.harvard.edu [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Nick Cox Sent: Freitag, 7. Juni 2013 15:52 To: statalist@hsphsun2.harvard.edu Subject: Re: st: Group identifiers within group Quite so. I didn't read your mail carefully enough. bysort GroupID AnimalID Parity : gen level2 = _n == 1 by GroupID : replace level2 = sum(level2) Optimistic numbering scheme allows for just under one trillion animals in Switzerland. You would need to use the lakes and glaciers too.... Nick njcoxstata@gmail.com On 7 June 2013 14:06, <isabel.lechner@vetsuisse.unibe.ch> wrote: > Thanks for the suggestions! But the input "egen level2 = group(GroupID > AnimalID Parity)" just grouped the "GroupID", but not by Animal ID > anymore (see output values of level2 below) > > AnimalID Parity GroupID level2 > CH003381451070 3 1 1 > CH003381451070 4 1 2 > CH003381451070 4 1 2 > CH003381462076 2 2 3 > CH003381462076 3 2 4 > CH003381462076 3 2 4 > CH110003099108 10 3 5 > CH110003099108 10 3 5 > CH110003099108 11 3 6 > CH110003099108 11 3 6 > CH110003099108 11 3 6 > CH110003099108 12 3 7 > CH110003099108 12 3 7 > CH110003099108 12 3 7 > > Also the suggestion by Alexis "by AnimalID Parity: gen var1= _n" didn't give me the right values. > > AnimalID Parity GroupID var1 > CH003381451070 3 1 1 > CH003381451070 4 1 1 > CH003381451070 4 1 2 > CH003381462076 2 2 1 > CH003381462076 3 2 1 > CH003381462076 3 2 2 > CH110003099108 10 3 1 > CH110003099108 10 3 2 > CH110003099108 11 3 1 > CH110003099108 11 3 2 > CH110003099108 11 3 3 > CH110003099108 12 3 1 > CH110003099108 12 3 2 > CH110003099108 12 3 3 > > And here again what I'm trying to get, do you have any other suggestions? > > AnimalID Parity GroupID 2nd level? > CH003381451070 3 1 1 > CH003381451070 4 1 2 > CH003381451070 4 1 2 > CH003381462076 2 2 1 > CH003381462076 3 2 2 > CH003381462076 3 2 2 > CH110003099108 10 3 1 > CH110003099108 10 3 1 > CH110003099108 11 3 2 > CH110003099108 11 3 2 > CH110003099108 11 3 2 > CH110003099108 12 3 3 > CH110003099108 12 3 3 > CH110003099108 12 3 3 > > > > > > -----Original Message----- > From: owner-statalist@hsphsun2.harvard.edu > [mailto:owner-statalist@hsphsun2.harvard.edu] On Behalf Of Nick Cox > Sent: Freitag, 7. Juni 2013 14:17 > To: statalist@hsphsun2.harvard.edu > Subject: Re: st: Group identifiers within group > > -group()- here is a very old and very undocumented function. I would suggest now to StataCorp, and I think I have suggested before, that calling -group()- should be illegal except under version control. > > You are confusing it with the -egen- function of the same name. > > I think you seek > > egen level2 = group(GroupID AnimalID Parity) > > Best not to use informative variable names such as "var1". > > Nick > njcoxstata@gmail.com > > > On 7 June 2013 13:02, <isabel.lechner@vetsuisse.unibe.ch> wrote: >> Dear all, >> >> I'm trying to create a new variable with a group identifier, but on a 2nd level. For the variable "GroupID" I did the following: >> >> by AnimalID: gen GroupID = _n >> >> Is it possible to do the same but on the 2nd level (based on AnimalID AND Parity), so that I will get the values as shown in the last variable "2nd level"? >> I've tried: by TierID: gen var1 = group(par), but this doesn't seem to do the job. Data are sorted by AnimalID and then Parity. >> >> >> AnimalID Parity GroupID 2nd level? >> CH003381451070 3 1 1 >> CH003381451070 4 1 2 >> CH003381451070 4 1 2 >> CH003381462076 2 2 1 >> CH003381462076 3 2 2 >> CH003381462076 3 2 2 >> CH110003099108 10 3 1 >> CH110003099108 10 3 1 >> CH110003099108 11 3 2 >> CH110003099108 11 3 2 >> CH110003099108 11 3 2 >> CH110003099108 12 3 3 >> CH110003099108 12 3 3 >> CH110003099108 12 3 3 >> >> Thanks for your help! >> Isabel >> >> * >> * For searches and help try: >> * http://www.stata.com/help.cgi?search >> * http://www.stata.com/support/faqs/resources/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/faqs/resources/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/faqs/resources/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/faqs/resources/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/faqs/resources/statalist-faq/ * http://www.ats.ucla.edu/stat/stata/