Density heatmap — irGSEA.densityheatmap • irGSEA

Easy to show the data distribution by density heatmap

irGSEA.densityheatmap(
  object = NULL,
  method = NULL,
  show.geneset = NULL,
  group.by = NULL,
  heatmap_width = 12,
  heatmap_height = 12,
  cluster.levels = NULL,
  ...
)

Arguments

object: A Seurat after perform irGSEA.score
method: A character. It should be one of the followling : AUCell, UCell, singscore, ssgsea.
show.geneset: A character. It should be one of the rownames of enrichment score matrix.
group.by: Default ident when it is set to NULL. You can specify other column of metadata.
heatmap_width: Width of the whole heatmap (including heatmap components), default 12.
heatmap_height: Heigh of the whole heatmap (including heatmap components), default 12.
cluster.levels: A vector equal to the number of clusters.
...: More parameters pass to densityHeatmap

Value

density heatmap

Examples


if (FALSE) {
# load PBMC dataset by R package SeuratData
library(Seurat)
library(SeuratData)
# download 3k PBMCs from 10X Genomics
InstallData("pbmc3k")
data("pbmc3k.final")
pbmc3k.final <- SeuratObject::UpdateSeuratObject(pbmc3k.final)

# Seurat object
pbmc3k.final <- irGSEA.score(object = pbmc3k.final, assay = "RNA",
slot = "data", msigdb = T, species = "Homo sapiens",
category = "H", geneid = "symbol",
method = c("AUCell", "UCell", "singscore", "ssgsea"), kcdf = 'Gaussian')

irGSEA.densityheatmap.plot1 <- irGSEA.densityheatmap(object = pbmc3k.final,
method = "UCell", show.geneset = "HALLMARK-INFLAMMATORY-RESPONSE")
irGSEA.densityheatmap.plot2 <- irGSEA.densityheatmap(object = pbmc3k.final,
method = "ssgsea", show.geneset = "HALLMARK-IL6-JAK-STAT3-SIGNALING")

}