> model1 <- lme(FOD_FODs ~ 0 + Condition_FODs*developmentaldisorder_FODs, random=~1|subject_FODs, + na.action=na.omit, + method = "ML", data=indvar_FODs,control = lmeControl(opt = "optim")) > anova(model1) numDF denDF F-value p-value Condition_FODs 4 3271 0.4038075 0.8060 developmentaldisorder_FODs 1 43 0.3213196 0.5738 Condition_FODs:developmentaldisorder_FODs 3 3271 0.0147161 0.9976 > shapiro.test(resid(model1)) Shapiro-Wilk normality test data: resid(model1) W = 0.98738, p-value < 2.2e-16 Can you please report these results in APA style?

A linear mixed-effects model was fitted with the FOD score as dependent variable, the condition (condition_FODs) and the developmentaldisorder (developmentaldisorder_FODs) as predictors, the subject (subject_FODs) as random effect, and the number of trials as control parameter. The residuals of the model were tested for normality using the Shapiro-Wilk test, which indicated that the model residuals were normally distributed (W = 0.98738, p < 2.2e-16). The model showed no main effects (p > 0.05).