[Rd] Addition of plot=F argument to termplot

From: Terry Therneau <therneau_at_mayo.edu>
Date: Fri, 19 Oct 2012 12:15:26 -0500

I have a suggested addition to termplot. We have a local mod that is used whenever none of the termplot options is quite right. It is used here almost daily for Cox models in order to put the y axis on a risk scale:

fit <- coxph(Surv(time, status) ~ ph.ecog + pspline(age), data=lung) zz <- termplot(fit, se=TRUE, plot=FALSE)

yy <- zz$age$y + outer(zz$age$se, c(0, -2, 2), '*') matplot(zz$age$x, exp(yy), log='y', type='l', lty=c(1,2,2),

     xlab="Age", ylab="Relative risk")

The return value is a list, each element of which is a dataframe. Another use is to overlay the fits from two different models. The changes to termplot.R and .Rd are small. Since we're not supposed to add attachments
to these emails, where to I send the updated files?

Terry T

R-devel_at_r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel Received on Fri 19 Oct 2012 - 17:15:05 GMT

This quarter's messages: by month, or sorted: [ by date ] [ by thread ] [ by subject ] [ by author ]

All messages

Archive maintained by Robert King, hosted by the discipline of statistics at the University of Newcastle, Australia.
Archive generated by hypermail 2.2.0, at Fri 19 Oct 2012 - 17:30:47 GMT.

Mailing list information is available at https://stat.ethz.ch/mailman/listinfo/r-devel. Please read the posting guide before posting to the list.

list of date sections of archive