I'm using the RMS package to fit and interpret a Cox-PH model:
############
require(rms)
f <- cph(S ~ rcs(x1, 3) + x2 + x3, method="efron", x=TRUE, y=TRUE,
surv=TRUE, time.inc=12)
f1 <- summary(f, x1=c(110, 130))
plot(f1, log=TRUE, main='')
############
I would like to display the HR for x1 for two different changes in the
predictor, say from 25th to the 50th, and 50th to 75th percentiles of
x1. Is it possible to do this in the same plot and _only_ for x1?
Thanks,
Thomas Speidel
______________________________________________
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.
