Increase from baseline [0.03, 95 CI (-0.12, 0.19), F(1,83)=0.16, p=.69]. Number of Symbolic Play Types–The number of symbolic play types revealed zero inflation [F(1,84)=4.48, p=.04]. Consequently, we modeled it using a hurdle model. The hurdle model did not reveal group differences in the change from below the hurdle to over the hurdle [0.14, 95 CI (-.84, 1.14), F(1,82)=0.57, p=.45], nor did it reveal group differences in the change of the observed level of symbolic play from baseline to the end of the study for participants who were within the measurement range [0.07, 95 CI (-0.56, 0.69), F(1,82)=0.84, p=.36]. In addition, we tested for Vasoactive Intestinal Peptide (human, rat, mouse, rabbit, canine, porcine) msds overall changes over time as the interactions were not significant. There were no overall changes for the participants below the hurdle [0.97, 95 CI (-0.61, 2.57), F(1,82)=1.49, p=.22], and the mean score for participants above the mean did not change [-0.04, 95 CI (-1.03, 0.95), F(1,82)=0.01, p=. 93]. At the follow-up time point, there was no difference between the groups in the number of order Avermectin B1a children that moved from out of the measurement range into the measurement range [0.14, 95 CI (-.0.84, 1.13), F(1,82)=0.08, p=.77]; however, there was a significant increase over time of the proportion of children in the measurement range [0.84, 95 CI (0.10, 1.59), F(1,82)=5.10, p=.02]. Although there was no overall difference between groups in the increase of the scores for children who were within the range of the measurement [0.07, 95 CI (-0.56, 0.69), F(1,82)=.05, p=.82], some results revealed increased scores between baseline and follow-up [0.35, 95 CI (-.06, 0.76), F(1,82)=2.85, p=.09]. Highest Play Level Achieved–Highest play level achieved did not show significant zero inflation [F(1,84)=0.39, p=.53], but the skew of the data led to a best-fit model of a Poisson distribution (AIC=833) rather than a normal distribution (AIC=839). Using the Poisson model, there was a significant treatment by time interaction for the highest play level achieved [0.36, 95 CI (0.11, 0.61), F(1,84)=9.07, p<.01], such that the JASPER group increased more than the PEI group in highest play level achieved, although the effect size was small (Cohen's f2=.11).Author Manuscript Author Manuscript Author Manuscript Author ManuscriptJ Consult Clin Psychol. Author manuscript; available in PMC 2016 June 01.Kasari et al.PageAt follow-up, the difference between the treatment groups was no longer significant [0.10, 95 CI (-.02, 0.22), F(1,83)=2.91, p=.09], and there were no overall sustained treatment gains from baseline across the groups [0.07, 95 CI (-0.02, 0.16), F(1,83)=2.54, p=.11]. Reynell Receptive Language--This variable did not show zero inflation [F(1,84)=0.00, p=.98]. The skew of the data led to a best-fit model of a Poisson distribution (AIC=1137) rather than a normal distribution (AIC=1793). There was no significant treatment effect on the rate of change [-0.11, 95 CI (-0.31, 0.08), F(1,84)=1.35, p=.25], but there was an overall increase of receptive language over time across both groups [1.23, 95 CI (1.09, 1.38), F(1,84)=278.59, p<.01]. This same pattern was obtained at follow-up, with an overall significant increase from baseline [0.62, 95 CI (0.54, 0.69), F(1,84)=278.59, p<.01], but no differences were found between the treatment groups in the degree of change [-0.06, 95 CI (-0.16, 0.04), F(1,84)=1.35, p=.25]. Both groups increased in receptive language by nearly 17 months over the 9-month study. Reynell Expressi.Increase from baseline [0.03, 95 CI (-0.12, 0.19), F(1,83)=0.16, p=.69]. Number of Symbolic Play Types--The number of symbolic play types revealed zero inflation [F(1,84)=4.48, p=.04]. Consequently, we modeled it using a hurdle model. The hurdle model did not reveal group differences in the change from below the hurdle to over the hurdle [0.14, 95 CI (-.84, 1.14), F(1,82)=0.57, p=.45], nor did it reveal group differences in the change of the observed level of symbolic play from baseline to the end of the study for participants who were within the measurement range [0.07, 95 CI (-0.56, 0.69), F(1,82)=0.84, p=.36]. In addition, we tested for overall changes over time as the interactions were not significant. There were no overall changes for the participants below the hurdle [0.97, 95 CI (-0.61, 2.57), F(1,82)=1.49, p=.22], and the mean score for participants above the mean did not change [-0.04, 95 CI (-1.03, 0.95), F(1,82)=0.01, p=. 93]. At the follow-up time point, there was no difference between the groups in the number of children that moved from out of the measurement range into the measurement range [0.14, 95 CI (-.0.84, 1.13), F(1,82)=0.08, p=.77]; however, there was a significant increase over time of the proportion of children in the measurement range [0.84, 95 CI (0.10, 1.59), F(1,82)=5.10, p=.02]. Although there was no overall difference between groups in the increase of the scores for children who were within the range of the measurement [0.07, 95 CI (-0.56, 0.69), F(1,82)=.05, p=.82], some results revealed increased scores between baseline and follow-up [0.35, 95 CI (-.06, 0.76), F(1,82)=2.85, p=.09]. Highest Play Level Achieved--Highest play level achieved did not show significant zero inflation [F(1,84)=0.39, p=.53], but the skew of the data led to a best-fit model of a Poisson distribution (AIC=833) rather than a normal distribution (AIC=839). Using the Poisson model, there was a significant treatment by time interaction for the highest play level achieved [0.36, 95 CI (0.11, 0.61), F(1,84)=9.07, p<.01], such that the JASPER group increased more than the PEI group in highest play level achieved, although the effect size was small (Cohen's f2=.11).Author Manuscript Author Manuscript Author Manuscript Author ManuscriptJ Consult Clin Psychol. Author manuscript; available in PMC 2016 June 01.Kasari et al.PageAt follow-up, the difference between the treatment groups was no longer significant [0.10, 95 CI (-.02, 0.22), F(1,83)=2.91, p=.09], and there were no overall sustained treatment gains from baseline across the groups [0.07, 95 CI (-0.02, 0.16), F(1,83)=2.54, p=.11]. Reynell Receptive Language--This variable did not show zero inflation [F(1,84)=0.00, p=.98]. The skew of the data led to a best-fit model of a Poisson distribution (AIC=1137) rather than a normal distribution (AIC=1793). There was no significant treatment effect on the rate of change [-0.11, 95 CI (-0.31, 0.08), F(1,84)=1.35, p=.25], but there was an overall increase of receptive language over time across both groups [1.23, 95 CI (1.09, 1.38), F(1,84)=278.59, p<.01]. This same pattern was obtained at follow-up, with an overall significant increase from baseline [0.62, 95 CI (0.54, 0.69), F(1,84)=278.59, p<.01], but no differences were found between the treatment groups in the degree of change [-0.06, 95 CI (-0.16, 0.04), F(1,84)=1.35, p=.25]. Both groups increased in receptive language by nearly 17 months over the 9-month study. Reynell Expressi.