目前有一份问卷数据，里面有大学毕业生户籍（农村/城市），性别，家庭背景，在校成绩，学校批次，毕业后的收入。

以毕业生收入为因变量y, 户籍（农村/城市)，性别，家庭背景，在校成绩，学校批次等为自变量x, 随着x的不同，而建立不同的多元回归模型，发现性别（female=1, male=0)这个变量当在各个回归模型中（无论其他自变量怎么变）都是显著的（系数为负，女性收入少）。导师说想看一下把男性和女性毕业生分开来后的对比结果。

除了用 reg income gender x2 x3 x4 if gender==1
                reg income gender x2 x3 x4 if gender==0  以外，还有没有其他对比男性女性的方法？

我想到了anova,
比如 anova income gender gpa school       (这意味着控制gpa,school后看一下男性和女性在收入上的对比是否显著）

大家还有没有别的关于男性女性的对比方法？ 或者在以上两个方法（reg 和anova）中有没有什么建议？

感谢！