Data were processed in SAS v9. Binary and multinomial logistic regression was conducted in STATA to examine the association between each alcohol measure (dependent variables: ever drinking, timing of first alcohol use, and AUD) and religious service attendance (independent variables: C-RA and A-RA). A robust variance estimator was used to account for familial clustering.
