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?