Good morning.
I have something that I am confused about when trying to use DESeq2 for my analysis. I read the tutorial from the developer of DESeq2, but I just wanted to confirm if I am understanding how LRT works. In this scenario, I have a dataset with conditions (control, cases). I want to test the differential gene expression while controlling for age and gender. This is the code that I have:
dds <- DESeqDataSetFromMatrix(countData = round(df),
colData = meta,
design= ~ condition + gender + age)
dds <- DESeq(dds, test = "LRT", reduced = ~ gender + age)
Would it be correct in saying that when I leave out "condition", I am testing to see if there is a significant impact from "condition" while controlling for age and gender? In other words, I am a bit confused on what should or should not be in the reduced formula.