From ripley at stats.ox.ac.uk  Fri Oct 19 19:20:57 2012
From: ripley at stats.ox.ac.uk (Prof Brian Ripley)
Date: Fri, 19 Oct 2012 18:20:57 +0100
Subject: [Rd] Addition of plot=F argument to termplot
In-Reply-To: <50818AAE.6030008@mayo.edu>
References: <50818AAE.6030008@mayo.edu>
Message-ID: <50818BF9.1010304@stats.ox.ac.uk>

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



[Rd] Addition of plot=F argument to termplot
Prof Brian Ripley