The analyses were conducted using the statistical modeling program Mplus (Muthén and Muthén, 1998; 2005). The default estimator for the MIMIC models was a robust, weighted, least-squares estimator. Model fit indices included chi-square, comparative fit index (CFI), root mean square residuals (RMSR), and root mean square error approximation (RMSEA). We used the following cutoff values as indicators of good fit: CFI > 0.95; RMSR <0.05; RMSEA < 0.06, as suggested by Muthén and Muthén (2005). Sampling weights and design effects of the NSDUH were taken into account for parameter estimation as well as standard error and model fit calculations.