Models for dichotomous outcomes were performed as relative risk regressions using a quasi-likelihood approach in SAS PROC GENMOD (SAS Institute, 2008). The predicted outcome was p (the probability of being a drinker) while the residual variance was allowed to be heteroscedastic, estimated as p/(1 − p). Predicted outcomes were constrained to lie in the interval (0,1).
