root/MGET/Branches/Jason/PythonPackage/src/GeoEco/Statistics/ClevelandPlotForDataframe.r @ 973

# ClevelandPlotForDataframe.r
#
# Copyright (C) 2011 Jason J. Roberts
#
# This program is free software; you can redistribute it and/or
# modify it under the terms of the GNU General Public License
# as published by the Free Software Foundation; either version 2
# of the License, or (at your option) any later version.
#
# This program is distributed in the hope that it will be useful,
# but WITHOUT ANY WARRANTY; without even the implied warranty of
# MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
# GNU General Public License (available in the file LICENSE.TXT)
# for more details.
#
# You should have received a copy of the GNU General Public License
# along with this program; if not, write to the Free Software
# Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA  02110-1301, USA.

ClevelandPlot <- function(df, transforms=NULL)
{
    library(stats)
    
    # Apply the transforms (if any).
    
    labs <- names(df)

    if (!is.null(transforms))
        for (i in 1:length(names(df)))
        {
            name <- names(df)[[i]]
            if (name %in% names(transforms))
            {
                d <- df[[name]]
                df[[name]] <- eval(parse(text=as.character(transforms[[name]])))
                labs[i] <- deparse(eval(parse(text=paste("substitute(", transforms[[name]], ", list(d = as.name(\"", name, "\")))", sep=""))))
            }
        }

    # Open a graphics device for a suitable number of plots.

    numPlots <- length(names(df))
    par(mfrow=c(ceiling(numPlots / ceiling(sqrt(numPlots))), ceiling(sqrt(numPlots))), mar=c(2, 0, 2, 0))

    # Create the plots.

    for (i in 1:length(names(df)))
    {
        dotchart(df[[names(df)[[i]]]], lcolor=NULL, bg="black", main=labs[i])
    }
}

ClevelandPlotToFile <- function(filename, df, transforms, res=1000.0, width=3000.0, height=3000.0, pointSize=10.0, bg="white")
{
    if (tolower(substr(filename, nchar(filename) - 3, nchar(filename))) == ".emf")
        win.metafile(filename, width=width, height=height, pointsize=pointSize)
    else
        png(filename, res=res, width=width, height=height, pointsize=pointSize, bg=bg)

    tryCatch(ClevelandPlot(df, transforms), finally=dev.off())
}