I'm trying to replace the panel.smooth
for argument panel
in pairs()
with a regression line drawing function instead of lowess line, with no success.
I tried to create function reg
and place it for argument panel
but that does not work? Any fix?
reg <- function(x, y) abline(lm(y~x)) # made to draw regression line instead of lowess line
panel.cor <- function(x, y, digits = 2, prefix = "", cex.cor, ...)
{
usr <- par("usr"); on.exit(par(usr))
par(usr = c(0, 1, 0, 1))
r <- abs(cor(x, y))
txt <- format(c(r, 0.123456789), digits = digits)[1]
txt <- paste0(prefix, txt)
text(0.5, 0.5, txt, cex = 1.1, font = 4)
}
#EXAMPLE OF USE:
pairs(USJudgeRatings[1:5], panel = panel.smooth, # replace HERE for panel.smooth #
cex = 1.5, pch = 19, col = adjustcolor(4, .4), cex.labels = 2, font.labels = 2, lower.panel = panel.cor)
You never used your function reg
. I slightly modified reg
to take a color argument. Instead of using panel.smooth
(which gives the loess curve) I modified it to use the linear model and your reg function.
reg <- function(x, y, col) abline(lm(y~x), col=col)
panel.lm = function (x, y, col = par("col"), bg = NA, pch = par("pch"),
cex = 1, col.smooth = "red", span = 2/3, iter = 3, ...) {
points(x, y, pch = pch, col = col, bg = bg, cex = cex)
ok <- is.finite(x) & is.finite(y)
if (any(ok)) reg(x[ok], y[ok], col.smooth)
}
pairs(USJudgeRatings[1:5], panel = panel.lm,
cex = 1.5, pch = 19, col = adjustcolor(4, .4), cex.labels = 2,
font.labels = 2, lower.panel = panel.cor)