I have the following problem: I have my model (controlling for Batch, Age, Sex) and I am interested in the condition.
dds <- DESeqDataSetFromMatrix(countData = cts, colData = colData, design = ~ BATCH + AGE + SEX + CONDITION
Now I need to do some grouping the condition, following vignette with
dds$group <- factor(paste0(dds$genotype, dds$condition)) design(dds) <- ~ group
But doing it this way I remove the control for Batch, Age and Sex.
How can I control for all these factors and do the grouping as well?
Thank you very much!