I'm modeling something that exclusively uses dummy variables, both the dependent and independent variables.

I'm using logit. Doesn't seem as though there is any obvious preference for logit or probit as far as I can tell, though if someone thinks otherwise I'd love to hear it.

However, I'm wondering what makes most sense in terms of predictive margins. dydx(*) default? dydx at means? Something else?

For what it's worth, the coefficients don't change much based upon my selection. But i still want my decision to be defensible in theory.