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: Group identifiers within group
From
<[email protected]>
To
<[email protected]>
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: [email protected] [mailto:[email protected]] On Behalf Of Nick Cox
Sent: Freitag, 7. Juni 2013 15:52
To: [email protected]
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
[email protected]
On 7 June 2013 14:06, <[email protected]> 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: [email protected]
> [mailto:[email protected]] On Behalf Of Nick Cox
> Sent: Freitag, 7. Juni 2013 14:17
> To: [email protected]
> 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
> [email protected]
>
>
> On 7 June 2013 13:02, <[email protected]> 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/