Thank you very much for your answer. I would like to construct for presentation purposes the HR(t), not the beta(t). How can I perform this?
I obtained the x-y values of the cox.zph plot time= as.numeric(as.character(rownames(cox.zph.object$y))) HR=exp(cox.zph.object$y[,1]) However when I plot the plot(time, HR) and fit a lowess line (to time, HR), the HR seems to jump over proportionally high (higher that what the cox.zph plot suggests) Am I making any error in my thinking? Many thanks JT -- View this message in context: http://r.789695.n4.nabble.com/How-to-understand-the-plotting-of-the-cox-zph-function-tp3788886p3793651.html Sent from the R help mailing list archive at Nabble.com. ______________________________________________ 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.
