I need to identify differential genes across 8 groups. I know this can be done using DESeq2 or EdgeR as these are better suited and equipped.
However, I was thinking if something like this piece of code could be tweaked. Here it works for two groups each with 3 replicates and generates a heatmap. The code is in R
pvals <- apply(dm, 1, function(x) wilcox.test(x~c(rep(0,3), rep(1,9)))$p.value) heatmap(dm[pvals<0.05,])
Thank you for your inputs and suggestions.