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
Addition of plot=F argument to termplot
2 messages · Terry Therneau, Brian Ripley
On 19/10/2012 18:15, Terry Therneau wrote:
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?
bugs.r-project.org is the best place to lodge enhancement suggestions: it accepts attachments.
Brian D. Ripley, ripley at stats.ox.ac.uk Professor of Applied Statistics, http://www.stats.ox.ac.uk/~ripley/ University of Oxford, Tel: +44 1865 272861 (self) 1 South Parks Road, +44 1865 272866 (PA) Oxford OX1 3TG, UK Fax: +44 1865 272595