plot.HTSCluster: Visualize results from clustering using a Poisson mixture model

Description

A function to visualize the clustering results obtained from a Poisson mixture model.

Usage

## S3 method for class 'HTSCluster':
plot(x, file.name = FALSE, 
    graphs = c("map", "map.bycluster", "lambda"), data=NA, ...)
## S3 method for class 'HTSClusterWrapper':
plot(x, file.name = FALSE,
    graphs = c("capushe", "ICL", "BIC"), capushe.validation=NA, ...)

Arguments

An object of class "HTSCluster" or "HTSClusterWrapper"

file.name

Optional file name if plots are to be saved in a PDF file.

graphs

Type of graph to be included in plots. May be equal to "map", "may.bycluster", "weighted.histograms", and/or "lambda" for objects of class "HTSCluster" and c("ICL", "BIC") f

capushe.validation

Optional number of clusters to use for capushe validation (should be less than the maximum number of clusters specificed in the "HTSClusterWrapper" object).

data

(n x q) matrix of observed counts for n observations and q variables (only required for the plotting of weighted histograms)

...

Additional arguments (mainly useful for plotting)

Details

For objects of class "HTSCluster", the plotting function provides the possibility for the following visualizations: 1) A histogram of maximum conditional probabilities across all clusters. 2) Per-cluster boxplots of maximum conditional probabilities. 3) Weighted histograms of observation profiles (with weights equal to the corresponding conditional probability for each observation in each cluster), plotted independently for each variable. Fitted densities after fitting the Poisson mixture model are overlaid in red. 4) A global view of $\ensuremath\boldsymbol{\lambda}$ and $\ensuremath\boldsymbol{\pi}$ values for the selected model. When the number of conditions <= 2,="" bar="" heights="" represent="" the="" value="" of="" $\ensuremath\boldsymbol{\lambda}_k$="" for="" each="" cluster,="" and="" width="" corresponds="" to="" $\ensuremath\boldsymbol{\pi_k}$.="" <="" p="">

For objects of class "HTSClusterWrapper", the plotting function provides the possibility for one or all of the following visualizations:

1) ICL plot for all fitted models.

2) BIC plot for all fitted models. 5) Capushe diagnostic plots.

References

Andrea Rau, Gilles Celeux, Marie-Laure Martin-Magniette, and Cathy Maugis-Rabusseau (2011). Clustering high-throughput sequencing data with Poisson mixture models. Technical report RR-7786, Inria Saclay -- Ile-de-France.

Examples

Run this code

set.seed(12345)

## Simulate data as shown in Rau et al. (2011)
## Library size setting "A", high cluster separation
## n = 2000 observations
simulate <- PoisMixSim(n = 200, libsize = "A", separation = "high")
y <- simulate$y
conds <- simulate$conditions

## Run the PMM-II model for g = 3
## "TC" library size estimate, EM algorithm
run <- PoisMixClus(y, g = 3, lib.size = TRUE,
    lib.type = "TC", conds = conds, init.type = "small-em")

## Visualization of results (not run):
## plot(run)