celdaHeatmap for celda Gene clustering model

# S3 method for celda_G
celdaHeatmap(counts, celda.mod, nfeatures = 25, ...)

Arguments

counts

Integer matrix. Rows represent features and columns represent cells. This matrix should be the same as the one used to generate `celda.mod`.

celda.mod

Celda object of class "celda_G".

nfeatures

Integer. Maximum number of features to select for each module. Default 25.

...

Additional parameters.

Examples

celda.sim = simulateCells("celda_G") celda.mod = celda_G(celda.sim$counts, L=celda.sim$L, nchains=1, max.iter=1)
#> --------------------------------------------------------------------
#> Starting Celda_G: Clustering genes.
#> --------------------------------------------------------------------
#> Thu Sep 06 12:55:31 2018 .. Initializing chain 1 with 'random' (seed=12345)
#> Thu Sep 06 12:55:31 2018 .... Determining if any gene clusters should be split.
#> Thu Sep 06 12:55:32 2018 .... Cluster 8 was reassigned and cluster 2 was split in two.
#> Thu Sep 06 12:55:32 2018 .... Completed iteration: 1 | logLik: -289530.128307233
#> Thu Sep 06 12:55:32 2018 .. Finished chain 1 with seed 12345
#> --------------------------------------------------------------------
#> Completed Celda_G. Total time: 0.6326349 secs
#> --------------------------------------------------------------------
celdaHeatmap(celda.sim$counts, celda.mod)