Data were analyzed by applying a cohort sequential design, and analyses were conducted using Mplus version 7. All analyses incorporated stratification and survey weights. Survey weights were incorporated to compensate for differences in selection probabilities of cases within sampling units, differential nonresponse rates, or chance fluctuation of the sample from the population. Missing data were handled with maximum likelihood estimation.
