Hi, I am using edger to analyze RNA-seq data. I have generated count table for each genes under each condition. However, my goal is not to look at differential expression of specific genes. Instead, I want to look at differential expression of some gene families. For example, I have gene A1, A2,...,An under gene family A . I want to see whether gene family A is differentially expressed. My question is whether the expression of gene family A is simply the addition of expression of genes A1 to An regarding to normalized read counts. By the way, the gene lengths are slightly different for genes within same family. What is in my mind is that I am quite sure FPKM is addible but not read counts. Thus, if read count is not addible, what should I do to analyze DE of gene families. Even if read count is addible, can I simply generate a new table for gene families and follow the edger guide as genes. Thank you for your time.
For an edgeR analysis, fry() is quicker and better than roast(). It's actually equivalent to using roast() with nrot=Inf.
Thanks for the correction. I've done most of my gene set testing with limma, so I'm less familiar with the options for edgeR.