nineplot <- function(y, ...)
{
# NINEPLOT
# Normal probability plot surrounded by reference plots
# GKS  Jul 96, 29 Dec 96
#
	oldpar <- par(mfrow = c(3, 3), mar = rep(1, 4), col = 1)
	on.exit(par(oldpar))
	doqqplot <- function(x, p, ...)
	{
		q <- qnorm(p)[order(order(x))]
		.S(plot(range(q), range(x), type = "n", axes = F, ...), "plot")
		box(...)
		points(q, x, ...)
		qqline(x, ...)
	}
	y <- y[!is.na(y)]
	n <- length(y)
	p <- ppoints(n)
	for(i in 1:4) {
		z <- rnorm(n)
		doqqplot(z, p, col = 1)
	}
	doqqplot(y, p, col = 2)
	for(i in 1:4) {
		z <- rnorm(n)
		doqqplot(z, p, col = 1)
	}
	invisible()
}