Interaction of fixed effects in DESeq2
1
0
Entering edit mode
Silva • 0
@silva-7695
Last seen 9.5 years ago
Finland

Hi,

probably a silly question but is there a difference between these two models in DESeq2:

counts ~ Condition : Generation

counts ~ Condition * Generation

They both model the interaction between condition and generation, am I correct?

However, when I try to run the upper model, I get an error message:

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.

The model with the asterisk works just fine.

Thanks!

 

 

 

deseq2 interactions • 1.7k views
ADD COMMENT
1
Entering edit mode
@mikelove
Last seen 3 days ago
United States

From ?formula:

"In addition to + and :, a number of other operators are useful in model formulae. The * operator denotes factor crossing: a*b interpreted as a+b+a:b."

So it doesn't really make sense to have only an interaction term (~ a:b). Interaction terms are usually meaningful in designs with main effect terms.

ADD COMMENT

Login before adding your answer.

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