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