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
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!