Substance abuse outcomes were analyzed using generalized linear mixed models (GLMMs). For each of the 13 weeks in the TLFB, participants were scored either as drug and alcohol free or as having used substances (see Figure 2). GLMM is appropriate for the present case because it allows for full nesting of scores and can accommodate binary data (Breslow & Clayton, 1993). Wolfinger and O’Connell’s (1993) refinement of the Breslow and Clayton (1993) approach was used.