I would like to test for group-specific effects in my dataset. This contains a group, time and individual (ID) factors as is in the DESEQ2 Vignette under "Group-specific condition effects, individuals nested within groups". After adding a 'id.n' column to my dds, the vignette states:
"Now we can reassign our DESeqDataSet a design of ~ grp + grp:ind.n + grp:cnd, before we call DESeq. This new design will result in the following model matrix:"
It provides no information on reassigning the formula. I tried the following but the formula is linked to ibjects and this replaces it with language.
dds@design[] <- 'group + group:ind.n+group:time'
Thank you for you help.