R/bdm_boxp.R

In bigMap: Big Data Mapping

# +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++
# The bigMap Package for R.

# Copyright (c) 2018, Joan Garriga <jgarriga@ceab.csic.es>, Frederic Bartumeus <fbartu@ceab.csic.es> (Blanes Centre for Advanced Studies, CEAB-CSIC).

# bigMap 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 3 of the License, or (at your option) any later version.

# bigMap 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 for more details.

# You should have received a copy of the GNU General Public License along with this program. If not, see http://www.gnu.org/licenses.
# +++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++++

#' Clustering statistics box-plot.
#'
#' @param bdm  A \var{bdm} instance as generated by \code{bdm.init()}.
#'
#' @param byVars A logical value. By default (\code{byVars = FALSE}) box-plots are grouped by cluster. With \code{byVars = TRUE} box-plots are grouped by input feature.
#'
#' @param layer The number of a layer (1 by default).
#'
#' @return None.
#'
#' @details If the number of clusters is large, only the first 25 clusters will be plotted. Note that the WTT algorithm numbers the clusters based on density value at the peak cell of the cluster. Thus, the numbering of the clusters is highly correlated with their relevance in terms of partial density. Therefore, in case of more than 25 clusters, the most relevant should always be included in the plot.
#'
#' @examples
#'
#' bdm.example()
#' bdm.boxp(exMap)
#' bdm.boxp(exMap, byVars = TRUE)

bdm.boxp <- function(bdm, byVars = F , layer = 1)
{
    if (is.null(bdm$wtt)) {
        return(message('+++ Error: up-stream step WTT not found'))
    }

    if (bdm$wtt[[layer]]$s > 25) {
        message('+++ Showing only 25 first clusters')
    }
    L <- bdm.labels(bdm, layer = layer)
    K <- sort(unique(bdm$wtt[[layer]]$C))

    parbdm.set(mar = c(1.25, 1.25, 1.25, 0.75))

    if (byVars)
    {
        layout(layout.get(ncol(bdm$data)))
        nulL <- lapply(seq(ncol(bdm$data)), function(j) {
                X <- sapply(K, function(k) bdm$data[which(L == k), j])
                boxplot(X, main = colnames(bdm$data)[j], names = paste('C.', K))
            })
        nullPlts <- max(layout.get(ncol(bdm$data))) -ncol(bdm$data) -1
    }
    else
    {
        layout(layout.get(length(K)))
        nulL <- lapply(K, function(k) {
				l <- which(L == k)
				if (length(l)) {
                	boxplot(bdm$data[l, ], main = paste('cluster', k))
				} else {
					cat('+++  cluster ', k, 'empty \n')
				}
            })
        nullPlts <- max(layout.get(length(K))) -length(K) -1
    }

    if (nullPlts > 1) nulL <- sapply(seq(nullPlts-1), function(p) plot.null())

    parbdm.def()

}