Entering edit mode

4.1 years ago

seta
★
1.5k

Hi all,

Please kindly advise me if using **TMM method for read count normalization and then DEGseq** for identification of differentially expressed genes is the right approach?
I never this tool, but as I searched this package shouldn't be used for the analysis because its statistical assumptions are wrong. Could you please provide me some references for it?

Thank you

Check out DESeq2

I asked about DEGseq, it differs from DEseq.

Ok. According to the vignette the

`DEGexp`

function expects raw counts, as such TMM normalised counts breaks that assumption. The reason I suggest DESeq2 is that it's a very well established package for analysing RNA Seq counts.Where did you find this information?

Thank you for your response. As far as I know, TMM is between-sample normalized count (not break the assumption that you mentioned, please correct me if I'm wrong), and, therefore, more suitable for DE analysis, do you agree with me? However, DEseq or edgeR has internally calculated TMM, but unfortunately, it is not clearly specified about DEGseq, have you any opinion about normalization way in DEGseq?

As I search, DEGseq assumes a Poisson distribution, whereas edgeR and DEseq2 assume a negative binomial (NB) distribution that is more appropriate than the Poisson distribution for an RNA seq analysis. Do you agree with that the Poisson distribution assumption is not suitable for RNA-seq analysis?