function (iseed = 27) 
{
    set.seed(iseed)
    mycolors = rep(c("blue", "black", "red", "green", "orange", 
        "blueviolet", "azure", "cyan"), 4)
    par(mfrow = c(5, 4), mar = c(1, 1, 1, 1))
    for (i in 1:20) {
        x = rnorm(9)
        y = rnorm(9)
        m = min(x, y)
        M = max(x, y)
        qqplot(x, y, xlim = c(m, M), ylim = c(m, M), pch = 16, 
            col = mycolors[i], axes = F, xlab = "", ylab = "")
        abline(0, 1)
        locx = m + 0.01 * (M - m)
        locy = M - 0.05 * (M - m)
        text(locx, locy, substitute(bar(y) - bar(x) == Dbar, 
            list(Dbar = format(signif(mean(y) - mean(x), 3)))), 
            adj = 0, cex = 1.4)
    }
}