First, my apologies if this has been covered already. I thought for sure it would have been, but I can't find the relevant info with my searches.
It's pretty obvious from the DESeq2 vignette how to test whether a gene is differentially expressed (DE) and how to do so at various LFC thresholds, but what if I want to determine which genes are not DE? Is there a way to tease apart the genes for which the null hypothesis is not rejected because the gene's expression has actually not changed vs the genes that just have too low read counts to make a determination?
Said another way, if I want to create 3 groups of genes, those (1) highly likely to be DE, (2) highly unlikely to be DE, and (3) with insufficient power to make a determination, how do I do that? Is there something more sophisticated I can do with DESeq2 than just applying a mean count cutoff to define the third class? Maybe use the LFC thresholds in some way?