Hi Gordon,
I am currently working on some code that would benefit from having
direct access to the code in edgeR::glmLRT that handles reparametrizing
the design matrix so that the N linearly independent components of the
specified contrasts are represented in the first N columns of the design
matrix (i.e. the part of glmLRT that does all the QR decomposition
stuff). Would you consider splitting this logic into a separate function
that takes a design matrix and contrast matrix and returns the
equivalent reparametrized design matrix and vector of coefficients?
Thanks,
-Ryan Thompson
_______________________________________________
Bioc-devel@r-project.org mailing list
https://stat.ethz.ch/mailman/listinfo/bioc-devel