Add panel=points as an extra parameter.
Replace
if (partial.resid)
points(xx, pres[, i], cex = cex, pch = pch, col = col.res, ...)
by
if (partial.resid)
panel(xx, pres[, i], cex = cex, pch = pch, col = col.res, ...)
Also, I think, add col.smooth, span and iter as named parameters; see below.
This allows, e.g., the following
par(mfcol=c(2, 2))
data(mtcars)
cars.lm <- lm(mpg~disp+hp, data=mtcars)
termplot(cars.lm, partial.resid=T, panel=panel.smooth, col.smooth="blue")
termplot(cars.lm, partial.resid=T, panel=panel.smooth, col.terms=0,
col.smooth="blue")
(col.terms=0 hides the fitted response line or curve.) These partial regression plots can be highly insightful in a regression context.
Inclusion of col.smooth, span and iter as named parameters is not strictly necessary, but it does avert complaints from plot() when they are passed as unnamed parameters.
(I foreshadowed this when I submitted problem #6327, and offered a fix
that, in default of anything better, I assume will be made or has been made.)
John Maindonald email: [EMAIL PROTECTED] phone : +61 2 (6125)3473 fax : +61 2(6125)5549 Centre for Bioinformation Science, Room 1194, John Dedman Mathematical Sciences Building (Building 27) Australian National University, Canberra ACT 0200.
______________________________________________ [EMAIL PROTECTED] mailing list https://www.stat.math.ethz.ch/mailman/listinfo/r-devel
