This is a test version of Biostars. For the public version, visit https://www.biostars.org.
Normalize Bimodal Distribution

I am analyzing a single cell dataset from SMART-SEQ having 300 cells. I am following Seurat's guided tutorial for clustering. I am confused about how to normalize the data, and I am using NormalizeData() from Seurat.

However, I do not see much difference between raw and normalized counts. The distribution still appears to be bi-modal; how should I proceed with the normalization? I there a workaround through which I can evaluate the normalized counts? It is a follow-up question from the previous post

Here's how the data looks, Various Normalization options

scrna-seq seurat normalization transcriptomics kevin-blighe

It looks like you plotted the same thing 4 times with different labels. In particular, your raw counts here looks very different than in the previous question (i would expect to see some counts near 0 and no counts near 10^6)

enter image description here

1 answer

Likely because count normalization has nothing to do with the plotted metric - it's just the number of RNA molecules detected per cell. Count normalization isn't going to change that value.

See this answer for more info.

Hi Jared, thanks for your response. I wonder if there is a way to evaluate whether the normalization has been done correctly or not. I can always plot counts vs. density, but is there any other measure or plot?

Thanks, LChart, appreciate your quick response.

Log in to answer this question.