I want to determine if there are genes that are differentially expressed between males and females for each treatment (none, low, and high). For gender my reference level is set as female and for treatment my reference level is set as none. Samples and script below.
#design matrix >dds = DESeqDataSetFromMatrix(countData=countdata, colData=coldata, design=~gender+treatment+gender:treatment) #set reference levels >dds$treatment=relevel(dds$treatment, ref="none") >dds$gender=relevel(dds$gender, ref="female") >dds = DESeq(dds) >res=results(dds) >resultsNames(dds)  "Intercept" "gender_male_vs_female" "treatment_low_vs_none" "treatment_high_vs_none" "gendermale.treatmentlow" "gendermale.treatmenthigh"
1) So is my design set up correctly for the question I want to answer, which genes are differentially expressed between males and females under a specific treatment?
2) And does "gendermale.treatmentlow" look for DEGs between males and females for low treatments, whereas "gendermale.treatmenthigh" looks for DEGS between males and females for high treatment?
3) Why is there no "gendermale.treatmentnone"?