data(Haibara)
## put histograms on the diagonal
panel.hist <- function(x, ...)
{
usr <- par("usr")
par(usr = c(usr[1:2], 0, 1.3))
nB <- 15; nB1 <- nB + 1
xmin <- min(x, na.rm = TRUE)
xmax <- max(x, na.rm = TRUE)
w <- (xmax - xmin) / nB
breaks <- xmin
b <- xmin
for (i in 1:nB) {
b <- b + w
breaks <- c(breaks, b)
}
h <- hist(x, breaks = breaks, plot = FALSE)
y <- h$counts
y <- y / max(y)
rect(breaks[1:nB], 0, breaks[2:nB1], y, ...)
}
par(xaxs = "i", yaxs = "i", xaxt = "n", yaxt = "n")
pairs(Haibara, diag.panel = panel.hist, pch = 20, cex.labels = 1.5,
label.pos = 0.9, lower.panel = NULL)
Run the code above in your browser using DataLab