This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Seurat: How to make heatmap of Log transformed expression levels

I'm new to single-cell sequencing analysis so I am sorry if this is a super simple question. I am working with a single-cell sequencing dataset and I'm trying to generate a heatmap of the average expression levels of select genes in select clusters. This is what I did:

avgexp <- AverageExpression(seurat.object, features = c("genes", "of", "interest"), return.seurat = TRUE) subset.1 <- subset(avgexp, idents = "clusters", "of", "interest")

The way I understand it , performing DoHeatmap(subset.1) plots scaled expression data, useful for telling apart "high expressing" and "low expressing" clusters for each gene. But if I want to make between-gene expression level comparisons I am thinking I would need the heatmap to plot the non-scaled log-transformed expression levels instead. I thought I could use DoHeatmap(subset.1, slot = "data") to do this, but according to the heatmap legend, the data is still scaled. Is there something I am missing here or not understanding?

Thank you!

seurat rna-seq r

0 answers

No answers yet.

Log in to answer this question.