R/plot.landscape.R

In hubbell: Functions for Hubbell's Grand Unified Theory

"plot.landscape" <-
function(x, main, ...)
{
    land <- x
    op <- par(no.readonly = TRUE)
    layout(matrix(c(1, 1, 2, 1, 1, 2), 2, 3, byrow = TRUE))
    d <- dim(land)
    nx <- d[1]
    ny <- d[2]
    J <- d[3]
    plot(1:nx, 1:ny, type="n", axes=FALSE, xlab="", ylab="")
    if (!missing(main))
        mtext(main)
    box()
    ijit <- jitter(rep(0,J))
    jjit <- jitter(rep(0,J))*100
    for (i in 1:nx)
        for (j in 1:ny) {
            plt <- land[i,j,]
            plt <- plt[plt != "#FFFFFF"]
            nplt <- length(plt)
            if (nplt) {
                ijit <- 10*jitter(rep(0,nplt))
                jjit <- 10*jitter(rep(0,nplt))
                points(rep(i,nplt) + ijit, rep(j,nplt) + jjit, 
                       col=plt, ...)
            }
        }
    rad <- rev(sort(table(as.vector(land))))
    plot(1:length(rad), rad, log="y", type="l", xlab="Rank", ylab="Abundance")
    points(1:length(rad), rad, col=names(rad), pch=16) 
    par(ylog=op$ylog)
    par(op)
    invisible()
}