Hello sir, I have a question regarding interpretation of DESeq2 output. I have 30 samples with treatment1/treatment2/control information and metadata of age and sex. I want to see the treatment effect adjusting for age and sex, so I binned age for 4 categories and make DESeq object as below.
dds <- DESeqDataSetFromMatrix(countData=raw_count, colData=mapping, design = ~ AGE_CAT + SEX + treatment)
I run the DESeq(), and got the result specifying the contrast by `c("treatment", "treatment1", "treatment2")`.
I got the result table of some of genes are differentially expressed across group, this time treatment1 and treatment2)
My explanation to this result when I do contrast treatment1 vs treatment2 is that "About one specific gene, comparing between the full model which consider the age, sex and treatment and reduced model which consider the age and sex, Walt statistics seem to be significant, so there seems to be some effect of treatment considering age and sex".
Here is my question: Is this safe to interpret this differentialy expressed result as "This gene is differentially expressed between treatment1 and treatment2 adjusted for age and sex ( i.e. considered that if age and sex are same)" ?
Thank you very much in advance for your answer.