Created attachment 1376[details]
Updated function
We have a local mod that is used whenever none of the termplot options
is quite right. It is invoked 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")
----
It is also useful to overlay the plots from two different fits.
The changed to termplot code simply adds a new block of statements.
I've attached updated versions of the .R and .Rd files; all lines except my addition agree with the current devel version.
The best reason I can give for why you should do this: you will be able to defer requests for adding other new arguments!