Hello,
I'm a master student in a biological-related field and I'm using DESeq2 for differential gene expression analysis; when I evaluate the result of the output 3 genes results as very down-regulated, however if I look at the read counts there are 4 reads in just one of the 8 cells of untreated samples and no reads in the 26 cells of the treated samples. A similar thing happens for the other two genes. Reading the guide that you provide I found that no pre-filtering of the read counts is needed rather is better to give the reads as they are obtained to DESeq2; however it seems that the software has a bias due to the abovementioned results. I have not further investigated the problem, I just report it.
Data that I used are freely available on CommonMind Consortium: https://www.synapse.org/#!Synapse:syn11617751, related metadata are https://www.synapse.org/#!Synapse:syn11638462; as model I used Dx (treatment received by monkeys) and as covariates Sex and DLPFC_RNA_isolation_Batch as reported in https://www.ncbi.nlm.nih.gov/pubmed/27668389 (online methods related to monkeys).
Regards,
Francesco Biagi
Hello MIcheal,
sorry for the delay, yes they are all significant. Yes, I have already fixed the problem, this post was ment in order to make aware mantainers of that problem, because as they said in their manual this pre-filtering step should not be needed at all.. as you said: "..filtering is not a requirement for the method to work."
However thank you a lot for you help.
Best Regards,
Francesco
Thanks for the post. Could you also say what version of DESeq2 you use? Are you using a Wald test (the default)? Just simple two group comparison? I wouldn’t think that these would have very small pvalues at all.
Is this similar to the setup you describe? I'm trying to figure out what range of fitted parameters would give a small p-value for such a gene.