search for: gamterms

I'm trying to plot pspline'd explanatory variables from coxph() models as per Therneau and Grambsch (Modeling Survival Data: Extending the Cox Model). They have s-plus functions for this at: http://www.mayo.edu/hsr/people/therneau/book/sfunction/gamterms.s http://www.mayo.edu/hsr/people/therneau/book/sfunction/plotterm.s I'd like to make these plots in R, but they make use of non-R functions, and I'm not sure how to adapt gamterms and plotterm to R. Have others done this, or does anyone have advice? Thanks, Jacob Etches University of To...

Brian, I used termplot(..., plot=FALSE) recently in R-devel: works like a charm. Thanks much for the update. Our in-house "gamterms" function, which this obviates, would also return the "constant" attribute from the underlying predict(..., type="terms") call. I have occasionally found this useful, and so it would be a worthwhile addition to termplot. Currently fit <- coxph(Surv(time, statu...