Differences between groups on all outcomes will be assessed by ANCOVAs for uni-dimensional questionnaires, and MANCOVAs for multi-dimensional questionnaires. To control for baseline differences in the primary outcome, results will present the treatment effect as the estimated mean group difference in change scores from baseline to post-intervention controlling for baseline demographic differences on parental education and number of children in the family, at a conventional p value of .05, the effect sizes using Cohen's d, and the estimated effect at 95% confidence intervals (Altman, 2005). Change scores will be calculated by simple subtraction, such as that a positive change score represents greater improvement in the desired therapeutic direction. The estimated marginal means will be reported as means scores adjusted for the covariates used (Austin, Manca, Zwarrenstein, Juurlink, & Stanbrook, 2010). Significant effects will be further explored using ANCOVAs, and pairwise comparisons determined any significant differences between conditions. Bonferroni adjustments will be made based on the number of pairwise comparisons conducted.