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, status) ~ pspline(age) + sex + ns(wt.loss),
data=lung)
zed <- termplot(fit, se=TRUE, plot=FALSE)
returns a list with components zed$age, zed$sex, zed$wt.loss. The constant could be added
as another element of the list or as an attribute, I don't have an opinion either way so
have not suggested a patch. You may have a reason for preferring one or the other.
Clearly this is not critical for version 3.0 release.
I sent this to you since you impliemented the plot=FALSE change, cc to the list in case
someone else is appropriate.
For those on the list, the recent change has three nice features:
a. Use of predict(.., type='terms') is a nuisance because the result is in data set
order rather than x order, a lines() call becomes a scribble. This has reduced each term
to the set of unique x values, in order, just what you need for a plot.
b. In the coxph example above I use plot(zed$age$x, exp(zed$age$y), log='y') to get
a better y-axis label. For all the developers, this is a nice way to deflect requests for
yet-another-plotting-option addition to termplot.
c. Easy to overlay results from two separate fits.
Terry T.
______________________________________________
R-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-devel
