[Bioc-devel] Splitting design matrix contrast pivoting logic into a separate function?
Ryan C. Thompson
rct at thompsonclan.org
Thu Aug 29 23:16:15 CEST 2013
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?
More information about the Bioc-devel