[R] Survival statistics--displaying multiple plots
Gregory Pierce
pierce.gregory at gmail.com
Thu May 3 04:21:16 CEST 2007
Hello all!
I am once again analyzing patient survival data with chronic liver disease.
The severity of the liver disease is given by a number which is continuously
variable. I have referred to this number as "meld"--model for end stage
liver disease--which is the result of a mathematical calculation on
underlying laboratory values. So, for example, I can generate a Kaplan-Meier
plot of patients undergoing a TIPS procedure with the following:
>plot(survfit(Surv(days,status==1),subset(tips,meld<10))
where "tips" is my data set, "days" is the number of days alive, and meld is
the meld score.
What I would like to do is display the survival graphs of patients with
meld<10, 10<meld<20, and meld>20. I am unsure about how to go about this.
Any suggestions would be appreciated.
Greg
More information about the R-help
mailing list