I would like to use edgeR's GLM functionality to determine the influence of several factors on a gene's expression levels. For performance reasons I would like to estimate the dispersion using the entire count table, but then determine the fit just for the gene I am interested in.
Is there any way to do this, i.e. call glmFit() on just one row of the table? I already tried to create a new CDS for just the gene and copy over the dispersion from the full CDS, but that does not seem to work.
Is there any way to do this, i.e. call glmFit() on just one row of the table? I already tried to create a new CDS for just the gene and copy over the dispersion from the full CDS, but that does not seem to work.
Comment