st: generating grouped relative survival estimates from stpm2

[Tim Evans <Tim.Evans@wmciu.nhs.uk>]

Hi all

I was hoping someone had had experience of estimating relative survival from stpm2. I am trying to take account of case mix by hospital and I can do this sort of ok, but what I need to do is extract the % relative survival as a result of correcting for certain factors.

My current unadjusted model is this:

Definition:
hosp1, hosp2, hosp3 hosp4 include binary variables that state whether a case belongs to that hospital or not and I have 4 hospitals

Model:
stpm2 hosp2 hosp3 hosp4, scale(hazard) df(5) tvc(hosp2 hosp3 hosp4) dftvc(3)
range temptime 0 5 200
predict rs1, survival timevar(temptime) zeros ci

predict rs2, survival timevar(temptime) at(hosp2 1) zeros ci

predict rs3, survival timevar(temptime) at(hosp3 1) zeros ci

predict rs4, survival timevar(temptime) at(hosp4 1) zeros ci

What I would to do is collapse the data to obtain the grouped 1 year survival for each hospital. After I get this working, I can start to add other variables into the model.

This is my first time on here, so apologies for any lack of information.

Dr T Evans
Cancer Registration Information Manager

West Midlands Cancer Intelligence Unit
Public Health Building
University of Birmingham
Birmingham
B15 2TT
 
Tel: 0121 414 4274

_DISCLAIMER:
This email and any attachments hereto contains proprietary information, some or all of which may be confidential or legally privileged. It is for the exclusive use of the intended recipient(s) only. If an addressing or transmission error has misdirected this e-mail and you are not the intended recipient(s), please notify the author by replying to this e-mail. If you are not the intended recipient you must not use, disclose, distribute, copy, print, or rely on this e-mail or any attachments, as this may be unlawful.


*
*   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/