Biostar Beta. Not for public use.
Considering that the basal expression level is not 0 for the calculation of the differential expression
0
Entering edit mode
18 months ago
Pin.Bioinf • 240
Malaga

Hello, I am doing a RNA-Seq DE analysis with DESeq2 and I was told to not consider the gene basal expression level as 0, so, how would I do that step? Are they referring to the relevel function in DESeq2? or setting the control samples as the targets of the comparison? Or is there more to it that i am not aware of?

Also, would results differ a lot if I do the analysis as indicated in the DESeq2 tutorial and not do the "not considering basal expression level as 0" step ?

Thank you

ADD COMMENTlink
0
Entering edit mode

I think this has to do with calculating fold changes and log changes. With basal levels of 0, you would end up dividing my zero and/or taking logs of 0s, both of which will give you infinite values.

ADD REPLYlink
0
Entering edit mode

and I was told to not consider the gene basal expression level as 0

Who told you?

Are you referring to zero raw counts or 0 normalised counts?

ADD REPLYlink
0
Entering edit mode

My biologist fellow. Normalized counts

ADD REPLYlink
0
Entering edit mode

You should obtain the logged counts and then covert these to the Z scale. Then work on expressed / not expressed based on the Z scores.

ADD REPLYlink
0
Entering edit mode

Logged counts? You mean the log2 of my normalized data and then transform to Z scale?

ADD REPLYlink
0
Entering edit mode

Yep - that's it.

ADD REPLYlink

Login before adding your answer.

Similar Posts
Loading Similar Posts
Powered by the version 2.1