My ultimate goal is for each locality, see how males and females are differentially expressed during each separate life stage.
Example: In Chicago (locality), how is my plant at the protonema stage, differentially expressed in males compared to females?
When I construct my object dds, I was wondering how can I make my dataset include only 1 locality but also include Sex, and life stages?
I am using the following so far:
dds <- DESeqDataSetFromMatrix(countData = countData, colData = colData, design = ~ Sex + Stage + Locality)
I have 3 replicates for each different locality, as well for each different sex, and life stage.