R/plot.tessresid.R

In stppResid: Perform residual analysis on space-time point process models.

plot.tessresid <- function(x, ..., col.key = rev(heat.colors(100)), cutoffs = NULL)
{
  X <- x
	residuals <- X$residuals
	if(is.null(cutoffs))	{
		cutoffs <- seq(min(residuals) - 1e-09, max(residuals) + 1e-09, length.out = length(col.key)+1)
		cutoff.m <- data.frame(cbind(cutoffs[-length(cutoffs)], cutoffs[-1]))
		CL.res <- rep(0, length(residuals))
		for(j in 1:length(residuals)) {
			CL.res[j] <- col.key[which((cutoff.m[,2] >= residuals[j]) & (cutoff.m[,1] < residuals[j]))]	
		}
	}
	else {
		if((length(cutoffs) - 1) != length(col.key))
			stop("length of cutoffs should be 1 more than length of col.key")
		if(!is.character(col.key))
			stop("col.key must be character vector of colors in hexadecimal")
		cutoff.m <- data.frame(cbind(cutoffs[-length(cutoffs)], cutoffs[-1]))
		CL.res <- rep(0, length(residuals))
		for(j in 1:length(residuals)) {
			CL.res[j] <- col.key[which((cutoff.m[,2] >= residuals[j]) & (cutoff.m[,1] < residuals[j]))]	
		}
	}
	layout(matrix(c(1,2), ncol=1), heights=c(2,.5))
	par(mar=c(4,4,4,4), bty="n")
	plot(X[[1]]$x, X[[1]]$y, type = "n", xlab = "x", ylab = "y", xlim = X[[1]]$xcoord, ylim = X[[1]]$ycoord)
	for(i in 1:length(X[[2]])) {
		temp <- X[[2]][i]
		x <- temp[[1]]$x
		y <- temp[[1]]$y
		p <- polygon(x, y, col = CL.res[i])	
	}
	points(X[[1]]$x, X[[1]]$y, ...)
	par(mar=c(1, 1, 1.5, 1))
	key <- (0:length(CL.res))/length(CL.res)
	plot(NULL, ylim=c(-3,0), xlim=c(-0.2, 1.2), type="n", axes=F, xlab="", ylab="", main="")
	image(key, -2:0, matrix(rep(key,2), nrow=length(CL.res)+1, byrow=F), add=T, col=col.key)
	for(i in seq(0, 1, length=10)) lines(c(i,i), c(-2,0), col=gray(0.3), lty=2)
	text(0, -2.5, round(min(cutoffs), 3), cex=1)
	text(0.25, -2.5, round(quantile(cutoffs, 0.25), 3), cex=1)
	text(0.5, -2.5, round(median(cutoffs), 3), cex=1)
	text(0.75, -2.5, round(quantile(cutoffs, 0.75), 3), cex=1)
	text(1, -2.5, round(max(cutoffs), 3), cex=1)
	mtext("Voronoi tessellation residuals", 3)
}