[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]

From |
"Justin B Echouffo Tcheugui" <jbe24@medschl.cam.ac.uk> |

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

Subject |
st: Multilevel modelling of survival data |

Date |
Mon, 16 Mar 2009 09:27:22 -0000 |

Dear Statalisters, I am a novice in multilevel modelling; I am trying to assess the effect of a preventive intervention on mortality data in a cluster randomised trial. The randomisation units were clinics in which my patients are nested. I am trying to assess the efcct of the intervention on mortality. I originally did a Cox regression using the following command: xi: stcox i.randomgroup, cluster(clinic). After reading the section about survival data in the book Multilevel and Longitudinal Modelling Using Stata by Sophia Rabe-Hesketh. I am now confused as to which command or function to use: Should I stick to my first analytical approach or use one of the following commands with the cloglog functions? cloglog event randomgroup yearsfollowup, eform vce(cluster clinic) xtcloglog event randomgroup yearsfollowup, eform i(clinic) gllamm event randomgroup yearsfollowup randomnum, eform i(clinic) link(cll) family(binom) adapt These commands gave me results that are different from Cox analysis. The results of xtclolog and gllamm with the cll link are similar. I understand that these commands are equivalent and allow the introduction of a random intercept for clinic. I also understand that the complementary log log link function follows if a proportional hazard model holds in continuous time and the survival times are interval censored. I am tempted to believe that my data is continuous-time data, but it is not interval censored as I know precisely the time to failure. Can I still use the command with the complementary log log function? If yes how to I include the time to failure in the model? I hope I have not completely mess up my model with years of follow-up to replace censoring interval. I hope my questions make sense. I look forward to hearing from you. Thanks Justin B. * * 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: Multilevel modelling of survival data***From:*Maarten buis <maartenbuis@yahoo.co.uk>

**References**:**st: RE: Return r(111) this time***From:*"Nick Cox" <n.j.cox@durham.ac.uk>

**Re: st: RE: Return r(111) this time***From:*Eva Poen <eva.poen@gmail.com>

**RE: st: RE: Return r(111) this time***From:*"Nick Cox" <n.j.cox@durham.ac.uk>

**Re: st: RE: Return r(111) this time***From:*Eva Poen <eva.poen@gmail.com>

**Re: st: RE: Return r(111) this time***From:*jjc.li@utoronto.ca

**Re: st: RE: Return r(111) this time***From:*Eva Poen <eva.poen@gmail.com>

**Re: st: RE: Return r(111) this time***From:*jjc.li@utoronto.ca

**RE: st: RE: Return r(111) this time***From:*"Nick Cox" <n.j.cox@durham.ac.uk>

**RE: st: RE: Return r(111) this time***From:*jjc.li@utoronto.ca

**Re: st: RE: Return r(111) this time***From:*Eva Poen <eva.poen@gmail.com>

**Re: st: RE: Return r(111) this time***From:*jjc.li@utoronto.ca

- Prev by Date:
**Re: st: Maximum length of a string expression in a program** - Next by Date:
**Re: st: do-file on dataset not currently used by Stata** - Previous by thread:
**Re: st: the command to specify one observation** - Next by thread:
**Re: st: Multilevel modelling of survival data** - Index(es):

© Copyright 1996–2016 StataCorp LP | Terms of use | Privacy | Contact us | What's new | Site index |