I am trying to compare between 10 groups each with two replicates (RNA-seq data). I wish to determine the genes which distinguish each individual apart from one another. So far looking through the edgeR manual I have not found a way to compare all of the groups at once in order to obtain values such as logFC, logCPM, PValue, and FDR. Do you have any suggestions of how to compare all of these groups to one another at once in edgeR instead of comparing in a pairwise manner?
In other words, is it possible to compare more than two groups at once using a negative binomial model?
Otherwise, do you have any suggestions of how to go about this problem?