plot_ClussCluster: Plots the results of 'ClussCluster'

In ClussCluster: Simultaneous Detection of Clusters and Cluster-Specific Genes in High-Throughput Transcriptome Data

Description

Plots the number of signature genes against the tuning parameters if multiple tuning parameters are evaluated in the object. If only one is included, then plot_ClussCluster returns a venn diagram and a heatmap at this particular tuning parameter.

Usage

plot_ClussCluster(object, m = 10, snames = NULL, gnames = NULL, ...)

top.m.hm(object, m, snames = NULL, gnames = NULL, ...)

Arguments

object	An object that is obtained by applying the ClussCluster function to the data set.
m	The number of top signature genes selected to produce the heatmap.
snames	The names of the cells.
gnames	The names of the genes
...	Addtional parameters, sent to the method

Details

Takes the normalized and log transformed number of reads mapped to genes (e.g., log(RPKM+1) or log(TPM+1) where RPKM stands for Reads Per Kilobase of transcript per Million mapped reads and TPM stands for transcripts per million) but NOT centered.

If multiple tuning parameters are evaluated in the object, the number of signature genes is computed for each cluster and is plotted against the tuning parameters. Each color and line type corresponds to a cell type.

If only one tuning parameter is evaluated, two plots will be produced. One is the venn diagram of the cell-type-specific genes, the other is the heatmap of the data with the cells and top m signature genes. See more details in the paper.

Value

a ggplot2 object of the heatmap with top signature genes selected by ClussCluster

Examples

data(Hou_sim)
run.cc <- ClussCluster(Hou_sim$x, nclust = 3, ws = c(2.4, 5, 8.8))
plot_ClussCluster(run.cc, m = 5, snames=Hou$snames, gnames=Hou$gnames)

ClussCluster documentation built on July 2, 2019, 5:03 p.m.