Generalized estimating equations (GEEs) were used for analyses to account for the correlation between responses from members of the same family. Wald chi-square statistics (χ2), degrees of freedom (df), and p-values (p) from univariable models are reported comparing the frequency of the risky behavior across drinking groups. The risky behavior was the outcome variable and the drinking onset status group was the predictor variable.
