Deseq2 - how to save zero counts on normalization step
1
0
Entering edit mode
12 days ago
Sandhiya • 0

Hi all,

I have raw read counts [of rna-seq]. I need to normalize this data and have done it using DESeq2. But, I found that DESeq2 eliminates the read counts of genes when it's equal to 0. I don't want to lose this data as it might affect my result.

Is there any way to do the normalization without eliminating zero read counts?

Thanks in advance.

deseq2 normalization • 243 views
ADD COMMENT
0
Entering edit mode

The author suggest removing empty genes are you sure you know better than he? I think empty genes are not removed by default, so can't you just not remove them?

Genes with read count 0 will be ignored when normalizing. They have to be.

ADD REPLY
0
Entering edit mode

@swbarnes2 I don't compete OR get into argument with what authors say.. I'm writing my requirement and how I can find a solution for it.

ADD REPLY
1
Entering edit mode
12 days ago

DESeq2 does not remove any genes by default (0 counts or not). Your results will contain all genes, though those that don't pass independent filtering due to low expression will have adjusted p-values set to NA.

Similarly, 0's are not removed from the normalized counts, e.g. min(counts(dds, normalized = TRUE)) will still be 0.

ADD COMMENT
0
Entering edit mode

Thank you @jared.andrews..

I have solved this ..

ADD REPLY

Login before adding your answer.

Traffic: 2890 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6