Three separate Cholesky decomposition models were run, one for each environmental variable. Full details of the model fitting process, including saturated model fit is included in the Appendix (Tables 3, 4, 5). Cholesky decompositions were transformed into correlated factors solutions for ease of interpretation.Table 3Model fitting results for the maternal negativity Cholesky decompositionModel−2LL (DF)AIC χ 2 (df) p 1. Saturated26496.63 (10485)5526.632. Cholesky with scalars on all variables27058.22 (10854)5350.22 2a. Cholesky with scalars on depression and physical aggression 27062.47 (10858) 5346.47 4.24 (4) 0.37 Best fitting model is highlighted in bold Table 4Model fitting results for the paternal negativity Cholesky decompositionModel−2LL (DF)AIC χ 2 (df) p 1. Saturated26316.88 (10371)5574.882. Cholesky with scalars on all variables26907.67 (10740)5422.67 2a. Cholesky with scalars on depression and physical aggression 26907.68 (10744) 5419.68 5.01 (4) 0.29 Best fitting model is highlighted in bold Table 5Model fitting results for the negative life events Cholesky decompositionModel−2LL (DF)AIC χ 2 (df) p 1. Saturated26563.08 (10602)5359.082. Cholesky with scalars on all variables27148.15 (10971)5206.15 2a. Cholesky with scalars on depression and physical aggression 27153.43 (10975) 5203.43 5.28 (4) 0.26 Best fitting model is highlighted in bold