> model1 <- lme(FOD_FODs ~ 0 + Condition_FODs*Age_FODs, random=~1|subject_FODs, + na.action=na.omit, + method = "ML", data=indvar_FODs,control = lmeControl(opt = "optim")) > ano <- anova(model1) > round(ano, 3) numDF denDF F-value p-value Condition_FODs 4 3139 0.347 0.846 Age_FODs 15 3139 0.084 1.000 Condition_FODs:Age_FODs 45 3139 0.046 1.000 > eta_squared(anova(model1), partial=TRUE) # Effect Size for ANOVA Parameter | Eta2 (partial) | 95% CI ------------------------------------------------------- Condition_FODs | 4.42e-04 | [0.00, 1.00] Age_FODs | 4.02e-04 | [0.00, 1.00] Condition_FODs:Age_FODs | 6.52e-04 | [0.00, 1.00] - One-sided CIs: upper bound fixed at [1.00]. > library(MuMIn) > r.squaredGLMM(model1) R2m R2c [1,] 0.001171854 0.001178491 Can you please report these results in APA style?

F(4, 3139) = 0.34, p = 0.846, η2 = 0.000442.
F(15, 3139) = 0.08, p = 1.000, η2 = 0.000402.
F(45, 3139) = 0.05, p = 1.000, η2 = 0.000652.

R2m = 0.00117
R2c = 0.00118