R help - Feb 2009 - log-minus-log plot

Dear experts,

I would like to know how to plot the log-minus-log plot for survival
analysis (to check the proportional assumption) in R.

Using the AML example.

fit <- survfit(Surv(time, status) ~ x, data=aml)
length(fit$surv) #20

as the length of fit$surv is shorter than aml$x and aml$time. I don't
know how to plot.

Thank you.

Regards,

CH


-- 
CH Chan
Research Assistant - KWH
http://www.macgrass.com

C.H. wrote:
> Dear experts,
> 
> I would like to know how to plot the log-minus-log plot for survival
> analysis (to check the proportional assumption) in R.
> 
> Using the AML example.
> 
> fit <- survfit(Surv(time, status) ~ x, data=aml)
> length(fit$surv) #20
> 
> as the length of fit$surv is shorter than aml$x and aml$time. I don't
> know how to plot.
> 
> Thank you.
> 
> Regards,
> 
> CH
> 
> 
library(Design)  # uses library(survival)
fit <- survfit(. . .)
survplot(fit, loglog=TRUE)
?survplot for details

-- 
Frank E Harrell Jr   Professor and Chair           School of Medicine
                      Department of Biostatistics   Vanderbilt University