A linear mixed modeling (LMM) growth curve framework was used to analyze the data. Analyses were run using the lme4 package (Bates & Maechler, 2010) in R version 2.11.1 (R Development Core Team, 2010). Percentages of missing data for variables in the study ranged from 2 to 8%. Missing data were multiply imputed (van Buuren & Groothuis-Oudshoorn, 2009). The equation, expressed in HLM form, for the primary analysis examining the effects of early caregiving and young adult romantic quality on the intercepts and slopes of anxious-depressed symptoms included one Level 1 and two Level 2 equations: Level 1: Anxious-Depressed Symptomsij = ß1i + ß2itij + eijLevel 2: ß1i = ß1 + ß3Genderi + ß4Early Caregivingi + ß5 Romantic Qualityi + ß6 (Early Caregivingi * Romantic Qualityi) + b1iLevel 2: ß2i = ß2 + ß7Early Caregivingi + ß8Romantic Qualityi + ß9(Early Caregivingi * Romantic Qualityi )+ b2i where ß1 is the intercept (the mean anxious-depressed symptom score at age 23), ß2 is the slope for the linear time term, ß3 is the effect of gender on the intercept, ß4 is the
