Thanks for the up front comment but I defined count before the forval command, otherwise the -forv- will not work.
Do you have any suggestions regarding my question (which I summarize below for convenience)? Can I use inrange for my problem, i.e with begin and end date? How to write it such that it works? Sorry, if this questions sound too easy for you but I try hard since days.
clear
input Id begin end state s2 res, auto
1 12 14 1 . 0
1 15 46 2 1 0
1 47 52 3 . 1
1 59 66 2 1 1
1 70 80 1 . 1
1 81 88 2 1 1
1 89 105 1 . 2
2 14 21 2 1 0
2 55 59 3 . 1
2 61 109 1 . 0
2 110 115 2 1 0
2 116 200 1 . 1
end
gen count = .
local N = _N
qui forval i = 1/`N' {
egen tag = tag(s2) ///
if Id == Id[`i'] & ///
inrange(begin, begin[`i'] - 30, end)
count if tag == 1
replace count = r(N) in `i'
drop tag
}
l, sepby(Id)
I want to know how often a person has been in a certain state (e.g. state
2) during the past 30 days before the current state. Important is that the
current state should not count.
I have data in which the time intervals are not of equal length. And to
keep things for the moment simple I defined a variable s2 which takes on
the value 1 if the state of interest occurs. The variable res shows my
desired results.
The next step would be to calculate the duration the person has been in s2
during the past 30 days. If there is a related solution I would be
extremely happy if you let me know this as well.
-----Ursprüngliche Nachricht-----
Von: "Martin Weiss" <[email protected]>
Gesendet: 31.05.09 15:50:53
An: <[email protected]>
Betreff: st: Re: Count occurrence and duration of states in panel data
<>
Two observations up front:
You do not have to use -:mylabel, auto- here. I used it a while back
to -input- a string that should have labels attached automatically, but your
vars are numeric.
In your -forv- you replace the variable "count" which does not exist. You
also want to make "Id" lowercase for your code to stand a chance...
HTH
Martin
_______________________
----- Original Message -----
From: <[email protected]>
To: <[email protected]>
Sent: Sunday, May 31, 2009 3:07 PM
Subject: st: Count occurrence and duration of states in panel data
>I want to know how often a person has been in a certain state (e.g. state
>2) during the past 30 days before the current state. Important is that the
>current state should not count.
>
> I have data in which the time intervals are not of equal length. And to
> keep things for the moment simple I defined a variable s2 which takes on
> the value 1 if the state of interest occurs. The variable res shows my
> desired results.
>
>
> clear
> input id begin end state s2 res :mylabel , auto
> 1 12 14 1 . 0
> 1 15 46 2 1 0
> 1 47 52 3 . 1
> 1 59 66 2 1 1
> 1 70 80 1 . 1
> 1 81 88 2 1 1
> 1 89 105 1 . 2
> 2 14 21 2 1 0
> 2 55 59 3 . 1
> 2 61 109 1 . 0
> 2 110 115 2 1 0
> 2 116 200 1 . 1
> end
>
> I followed some earlier questions at statalist and stata tips 39, 51 on
> inrange and events in intervals:
>
> gen count = .
>
> local N = _N
> qui forval i = 1/`N' {
> egen tag = tag(s2) ///
> if Id == Id[`i'] & ///
> inrange(begin, begin[`i'] - 30, end)
> count if tag == 1
> replace count = r(N) in `i'
> drop tag
> }
> l, sepby(Id)
>
> I surmise that the inrange command is not appropriate. As s2 should also
> count if even the begin date has been more than 30 days away but the
> person is in the state within the 30 day interval (with respect to the end
> date),
>
> Can inrange be adjusted such that it takes the begin date and the end date
> into account? Do I have to think about an entirely different approach?
>
> The next step would be to calculate the duration the person has been in s2
> during the past 30 days. If there is a related solution I would be
> extremely happy if you let me know this as well.
>
> Any ideas will be appreciated!
> Elke
> ____________________________________________________________
> Text: GRATIS für alle WEB.DE-Nutzer: Die maxdome Movie-FLAT!
> Jetzt freischalten unter http://movieflat.web.de
>
>
> *
> * 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/
____________________________________________________________
Text: GRATIS für alle WEB.DE-Nutzer: Die maxdome Movie-FLAT!
Jetzt freischalten unter http://movieflat.web.de
*
* 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/
