DESeq2 - Test for linear combinations
1
0
Entering edit mode
xavi • 0
@xavi-15151
Last seen 5.9 years ago

Hi,

How could we test for linear combinations between variables added to the model before running "DESeq()"?

It would be nice if we could detect vars with lin. comb. and remove them beforehand.

We basically would like to protect against this error before calling "DESeq()":

"Error in DESeqDataSet(se, design = design, ignoreRank) :
  the model matrix is not full rank, so the model cannot be fit as specified.
  one or more variables or interaction terms in the design formula
  are linear combinations of the others and must be removed
Calls: DESeqDataSetFromMatrix -> DESeqDataSet"
Execution halted

Many thanks

 

deseq2 linear combination • 749 views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 8 hours ago
United States

You can look at the rank of the design matrix using model.matrix(design(dds), colData(dds)) and qr(x)$rank.

ADD COMMENT

Login before adding your answer.

Traffic: 506 users visited in the last hour
Help About
FAQ
Access RSS
API
Stats

Use of this site constitutes acceptance of our User Agreement and Privacy Policy.

Powered by the version 2.3.6