[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


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



More information about the Bioc-devel mailing list