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]

From |
"Ciril Bosch-Rosa" <cirilbosch@gmail.com> |

To |
<statalist@hsphsun2.harvard.edu> |

Subject |
st: Interval Censored Data Formatting and Intcens |

Date |
Fri, 1 Jun 2012 14:12:40 -0700 |

Hello, I am trying to apply survival analysis in an environment where I can gather data in a very limited way. I have been looking at literature on how to manage the data but am not being successful at transmitting the structure of my data to Stata. It turns out that my data is very limited because I can only observe once in time each subject. I am taking this as an interval censored problem, where all subjects start being at risk at t_0 and then are censored until I look at some of them at (say) t_13, where they are either 1 (dead) or 0 (right censored). So, either the subject died at some before (and including) t_13 or he survived in which case I will never see him again (hence the right censoring). I looked at Cleves et al. but they are not clear on how I can tell Stata that an even happened while "I was not looking". They tell me how to tell stata I was not looking, but not how to tell I was not looking & something happened: id t_0 t_1 event 1 1 3 0 1 13 15 1 the example above I tell stata there is a censored interval from t=3 until t=13 and then I observe id=1 until he gets a failure at t=15, but I am not able to tell stata that they might have died in the t3-t13 interval: id t_0 t_1 event 1 1 3 0 1 13 13 1 2 1 3 0 2 13 13 0 In other words I observe you from t=1 to t=3 and then you get censored until I get one peek at t=13, where I can see if you died between 3 and 13 or not. I believe this can be done with ""intcens depvar1 depvar2 " where depvar are my t=3 and t=13 respectively, my problem though is that there are different interval censoring for different subjects id t_0 t_1 event 1 1 3 0 1 13 13 1 2 1 3 0 2 13 13 0 3 1 3 0 3 21 21 1 4 1 3 0 4 21 21 1 how can I tell intcens, "intcens depvar1 depvar2 if id=1|2 but depvar1 depvar3 if id=3|4" ?? Any help/advice is welcome Ciril PS If you are wondering in what environment I get this kind of data; this is a lab experiment where I apply a dose to groups of subjects and get instantaneous results, hence my observing groups of them at the same "time". My X-axis is not time, but rather "dosage" and I am trying to plot the survivor curve given the very limited observable outcomes for each dosage. * * 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/

**Follow-Ups**:**Re: st: Interval Censored Data Formatting and Intcens***From:*Steve Samuels <sjsamuels@gmail.com>

- Prev by Date:
**RE: st: Variable running totals** - Next by Date:
**st: Re: from variable names to values** - Previous by thread:
**RE: st: Variable running totals** - Next by thread:
**Re: st: Interval Censored Data Formatting and Intcens** - Index(es):