Hi
Have a look on how panel.smooth is defined:
panel.smooth
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))
lines(stats::lowess(x[ok], y[ok], f = span, iter = iter),
col = col.smooth, ...)
}
<bytecode: 0x000000000bfb2ce0>
<environment: namespace:graphics>
And change that to something like this:
panel.regression <- function (x, y, col = par("col"), bg = NA, pch =
par("pch"),
cex = 1, col.regres = "red", ...)
{
points(x, y, pch = pch, col = col, bg = bg, cex = cex)
ok <- is.finite(x) & is.finite(y)
if (any(ok))
abline(stats::lm(y[ok] ~ x[ok]), col = col.regres, ...)
}
Substitute panel.smooth with panel.regression in your pairs.panel() function
and you have it.
Yours sincerely / Med venlig hilsen
Frede Aakmann Tøgersen
Specialist, M.Sc., Ph.D.
Plant Performance & Modeling
Technology & Service Solutions
T +45 9730 5135
M +45 2547 6050
fr...@vestas.com
http://www.vestas.com
Company reg. name: Vestas Wind Systems A/S
This e-mail is subject to our e-mail disclaimer statement.
Please refer to www.vestas.com/legal/notice
If you have received this e-mail in error please contact the sender.
> -----Original Message-----
> From: r-help-boun...@r-project.org [mailto:r-help-boun...@r-project.org]
> On Behalf Of Shane Carey
> Sent: 25. april 2014 12:26
> To: r-help@r-project.org
> Subject: [R] Linear line on pairs plot
>
> Hi,
>
> Im trying to plot a linear line on the scatter plot using the pairs()
> function. At the moment the line is non linear. However, I want a linear
> line and the associated R value.
>
> Here is my current code:
>
> panel.cor.scale <- function(x, y, digits=2, prefix="", cex.cor)
> {
> usr <- par("usr"); on.exit(par(usr))
> par(usr = c(0, 1, 0, 1))
> r = (cor(x, y,use="pairwise"))
> txt <- format(c(r, 0.123456789), digits=digits)[1]
> txt <- paste(prefix, txt, sep="")
> if(missing(cex.cor)) cex <- 0.8/strwidth(txt)
> text(0.5, 0.5, txt, cex = cex * abs(r))
> }
>
>
> 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 = (cor(x, y,use="pairwise"))
> txt <- format(c(r, 0.123456789), digits=digits)[1]
> txt <- paste(prefix, txt, sep="")
> if(missing(cex.cor)) cex <- 0.8/strwidth(txt)
> text(0.5, 0.5, txt, cex = cex )
> }
>
>
> panel.hist <- function(x, ...)
> {
> usr <- par("usr"); on.exit(par(usr))
> par(usr = c(usr[1:2], 0, 1.5) )
> h <- hist(x, plot = FALSE)
> breaks <- h$breaks; nB <- length(breaks)
> y <- h$counts; y <- y/max(y)
> rect(breaks[-nB], 0, breaks[-1], y, col="cyan", ...)
> }
>
>
> pairs.panels <- function (x,y,smooth=TRUE,scale=FALSE)
> {if (smooth ){
> if (scale) {
>
> pairs(x,diag.panel=panel.hist,upper.panel=panel.cor.scale,lower.panel=pane
> l.smooth)
> }
> else
> {pairs(x,diag.panel=panel.hist,upper.panel=panel.cor,lower.panel=panel.sm
> ooth)
> } #else
> {pairs(x,diag.panel=panel.hist,upper.panel=panel.cor,lower.panel=panel.sm
> ooth)
> }
> else #smooth is not true
> { if (scale) {pairs(x,diag.panel=panel.hist,upper.panel=panel.cor.scale)
> } else {pairs(x,diag.panel=panel.hist,upper.panel=panel.cor) }
> } #end of else (smooth)
> } #end of function
>
> pairs.panels(iris[1:4])
>
> Thanks
>
> --
> Shane
>
> [[alternative HTML version deleted]]
>
> ______________________________________________
> R-help@r-project.org mailing list
> https://stat.ethz.ch/mailman/listinfo/r-help
> PLEASE do read the posting guide http://www.R-project.org/posting-
> guide.html
> and provide commented, minimal, self-contained, reproducible code.
______________________________________________
R-help@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/r-help
PLEASE do read the posting guide http://www.R-project.org/posting-guide.html
and provide commented, minimal, self-contained, reproducible code.
