Hello I have the following experimental design of RNASEQ experiment: 12 samples, 6 are males 6 are females, and two conditions of treatment (Before and After) according to this table:

ID(subject)    sex    treat
1              M      B
1              M      A
2              M      B
2              M      A
3              M      B
3              M      A
4              F      B
4              F      A
5              F      B
5              F      A
6              F      B
6              F      A

I found that there is a large gene expression variation among subjects from the same sex (ID). My assumption, based on biological observations, is that there should be genes that are deferentially expressed according to the treatment in the males but are not deferentially expressed in the females. What would be the right software and model design to find those DEG.

I tried DESeq2 but it did not accept the model design ~ID+sex+treat+treat:sex because the ID is limited to one sex so it is not a full model.

Thank you