R/ecospat.cor.panel.R

In ecospat: Spatial Ecology Miscellaneous Methods

# Adjusted by L. Mathys, 2006, modified by N.E. Zimmermann
 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 <- paste(prefix, txt, sep = "")
  if (missing(cex.cor)) 
    cex <- 0.8/strwidth(txt)
  text(0.5, 0.5, txt, cex = cex * r)
}


# Adjusted by L. Mathys, 2006
panel.hist <- function(x) {
  # usr <- par("usr")
  # on.exit(par(usr))
  par(usr = c(0, 1, 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 = "skyblue3")
}


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)
}



ecospat.cor.plot <- function(data) {
  pairs(data, lower.panel = panel.smooth, diag.panel = panel.hist, upper.panel = panel.cor)
  
}