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,method = "spearman"))
txt <- format(c(r, 0.123456789), digits = digits)[1]
txt <- paste(prefix, txt, sep = "")
if(missing(cex.cor)) cex.cor <- 0.8/strwidth(txt)
text(0.5, 0.5, txt, cex = cex.cor * r)
}
pairs(lily, lower.panel = panel.smooth, upper.panel = panel.cor,col = "grey38",
pch = 19,lwd = 2,
main = "Glacier lily seedlings relationship to \nnumber of flowers and rockiness")
Run the code above in your browser using DataLab