Skip to content

creating log-log survival plots that are not inverted

1 message · Bob Green

#
I am hoping for some advice regarding how to obtain a log-log 
survival plot that is not in the inverse. On page 128 of Modelling 
survival data by Therneau & Grambsch there is the an example of the 
type of desired plot, with a log of the survival curve  by 
years.  Marc Schwartz has provided me with some reproducible code.

The code below produces inverted plots.

library(surv2sample)
data(gastric)
fit <- survfit ( Surv(time, event) ~ treatment, data = gastric)
#Default plot:
plot(fit)
plot(fit, fun = "cloglog")

library(Design)
survplot(fit, loglog = TRUE, conf = "none")

Any assistance is much appreciated,

regards

Bob