Hi, I am trying to use edgeR for analysis of gene expression across 4 RNASeq samples(different stations) (3 replicates each. The information I would like to have, is the expression (fold change) for each gene across all 4 stations. If I use the GLM method by selecting one intercept, I will only get for each gene, the fold change in 3 of the 4 stations. How can I retrive the expression (fold change) for the 4th station? Should I use another approach?
Thanks for the help, any suggestion is really appreciated.