DESeq and checkFullRank(modelMatrix)
1
0
Entering edit mode
Merlin ▴ 10
@merlin-15723
Last seen 4.9 years ago
Vancouver

Hello guys,

I am running DESeq2 to compare 4 different conditions and I have a problem with a commmand line.

dds = DESeqDataSetFromMatrix(countData = countMatrix, colData = coldata, design = ~ group + treatment +group:treatment)

although I found online many others people solved the same problem, I couldn't solve by myself. The error that I get is this

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.

Do you have suggestions on how to solve this?

Thank you

RNA DESeq2 • 1.0k views
ADD COMMENT
0
Entering edit mode
@mikelove
Last seen 13 hours ago
United States

The vignette has suggestions on this. There’s a section devoted to this design error.

ADD COMMENT
0
Entering edit mode

Also I notice that many of our past interactions involve me suggesting you to read certain parts of the vignette. At this point I think it would be reasonable for you to first check the entire vignette and the function help pages before posting a support site question.

ADD REPLY

Login before adding your answer.

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