DESeq2: the model matrix is not full rank issue
2
0
Entering edit mode
Guangsheng • 0
@a1871ad5
Last seen 3 months ago
United States

Hi,

I have a coldata like following. when I run the code: dds <- DESeqDataSetFromMatrix(countData = counts,colData = coldata, design = ~batch+Stage). I get an error: Error in checkFullRank(modelMatrix) : 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.

I feel confused why I get this error because it seems there is no linear correlation between Stage and batch. Any suggestions will be helpful. Thank you!

enter image description here

DESeq2 • 459 views
ADD COMMENT
0
Entering edit mode

Understanding the appeal of adult content, the search for better content takes place among a multitude of choices. This platform https://pornpax.com/categories/big-dick/ appears as the final destination, presenting an extensive array of varied material. It assures a fulfilling journey by championing diversity, rather than merely accepting it.

ADD REPLY
1
Entering edit mode
ATpoint ★ 4.1k
@atpoint-13662
Last seen 7 hours ago
Germany

All of batches b2-9 is exclusively in ivt_bl and batches b14-15 are exclusively in ivt_exbl. That is the linear combination, the fact that one batch is always the same as the treatment effect.

You cannot correct for anything here. Use of ~Stage is your only option.

ADD COMMENT
0
Entering edit mode
@james-w-macdonald-5106
Last seen 6 hours ago
United States

The batches are completely nested within stage, so you cannot account for both in your design matrix.

Login before adding your answer.

Traffic: 525 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