Hello!

So currently i'm doing a research on gene expression particularly on miRNAs. I'm looking for the most differential expressed genes and the most similarly expressed from treated vs untreated samples comparison. I have conducted differential expression gene (DGE) analysis using NOISeq with the following script to complete and get the differential expressed genes:

de <- degenes(myresults, q = 0.95, M = NULL)

de_up <- degenes(myresults, q = 0.95, M = "up")

de_down <- degenes(myresults, q = 0.95, M = "down")

Now, i'm in need of help to find the most similarly expressed genes.

Based on my research DeSeq2 have a section/scripts that allows users to extract result table from DGE with alternatives hypothesis using 'results' function. And one of them is abs(log2FoldChange) < x. As shown in the link below

https://www.bioconductor.org/packages/devel/bioc/vignettes/DESeq2/inst/doc/DESeq2.html#tests-of-log2-fold-change-above-or-below-a-threshold

So my question is, can this script be applied to NOIseq as well? Or can the above script that I'm using for DGE be changed/adjusted to get similarly expressed?

I have been trying to search for other alternative but couldn't find any. I hope someone can help me out.

Thank you!