Survival Plot in R 2.10.0
I would like to produce a complimentary log-log survival plot with only the points appearing on the graph. I am using the code below, taken from the plot.survfit page of help for the the survival package (version 2.35-7). I am running in R 2.10.0 on Windows XP, and the list of packages following the error is loaded. Is there some specific 'type= ' syntax, or an additional parameter that I must specify?
leukemia.surv <- survfit(Surv(time, status) ~ x, data = aml) plot(leukemia.surv, lty = 2:3) ## works fine plot(leukemia.surv, lty = 2:3, fun="cloglog") ## works fine plot(leukemia.surv, type="p") ## generates error below:
Error in plot.default(tempx, tempy * yscale, type = "n", log = logax, : formal argument "type" matched by multiple actual arguments
search()
[1] ".GlobalEnv" "package:survival" "package:splines" [4] "package:stats" "package:graphics" "package:grDevices" [7] "package:utils" "package:datasets" "package:methods" [10] "Autoloads" "package:base" I can get a plot of points using plot(log(leukemia.surv$time),log(-log (leukemia.surv$surv)), type="p"), but then I lose the functionality of grid(). Ulimately what I want to do is overlay these plotted points with lines generated by various parametric survival models (dist=weibull, lognormal, etc). All help is appreciated. Thanks in advance!