Thanks for your replies.
What I am trying to calculate is if the mean of a dummy
variable is different across the categories of a separate
categorical variable. So if the mean of a dummy variable
(e.g. let's say 1=has university degree, 0=does not have
university degree) is significantly different across a
nominal variable like religious affiliation which has five
possible values. If I had just two categories in the
religious affiliation variable, I could just prtest
university, by(religion). Since I have multiple categories,
however, this becomes impossible.
If my DV was continuous, I could do an anova and with
post-hoc estimations figure out where the significant
differences between categories were. However, because my DV
is not continuous, I have been told an anova here is not
appropriate, hence my confusion. Perhaps I am just being
pedantic?
I would really appreciate your opinion now that I have
fully explained myself!