Hi, in Design package, a plot of survival probability vs. a covariate can be
generated by survplot() on a cph object using the folliowing code:
n <- 1000
set.seed(731)
age <- 50 + 12*rnorm(n)
label(age) <- "Age"
sex <- factor(sample(c('male','female'), n, TRUE))
cens <- 15*runif(n)
h <- .02*exp(.04*(age-50)+.8*(sex=='Female'))
dt <- -log(runif(n))/h
label(dt) <- 'Follow-up Time'
e <- ifelse(dt <= cens,1,0)
dt <- pmin(dt, cens)
units(dt) <- "Year"
dd <- datadist(age, sex)
options(datadist='dd')
S <- Surv(dt,e)
library(Design)
f <- cph(S ~ age, surv=TRUE,x=T,y=T)
plot(f,age=NA,time=5)
But the same code won't work if I used rms package:
detach(package:Design)
library(rms)
f <- cph(S ~ age, surv=TRUE,x=T,y=T)
plot(f,age=NA,time=5)
Error in xy.coords(x, y, xlabel, ylabel, log) :
'x' and 'y' lengths differ
Is there a way to plot the same graph using rms package. I like to use Frank
Harrell's new package rms and try to avoid using old Design package.
Thanks
John
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
