Entering edit mode
2.3 years ago
mohsamir2016
▴
30
Dear All, May be a naïve question, but buzzled me. I have a dataset where I created a DGElist object with TMM normalized count, and then estimated the tagwise dispersion from a design matrix.. code looks like :
dgeObj <- calcNormFactors(dgeObj)
design <- model.matrix(~Animal+Status)
y <- estimateDisp(dgeObj, design)
Tagwisedie<- y$tagwise.dispersion
From this code I obtained the 30 genes with the least Tagwise dispersion between two condition
I have also run DE analyses on the same dataset, and obtained the 30 genes that are the least DE between two condition
The question is : Should the 30 genes with least tagwise dispersion is the same as the 30 genes with least DE ?
In my hand it is not .. any comments ? or explanation
Thanks