The myopia progression was negatively associated with RPRE, that is, the more hyperopic the a
M and a
P180, and the flatter the a
J0, the faster the myopia progression (
Fig. 3), respectively, but only a
J0 was associated with axial elongation (
Fig. 4).
Table 2 shows the changes in coefficients of determination, and the detailed statistical results for the hierarchical regressions are listed in
Supplementary Tables S1 and
S2. After controlling for the baseline refraction (model 2), the RPRE was independent of ΔM and ΔAL. Home scene parameters were then added as a covariate in the regression models. The introduction of the normality-transformed dioptric volume variable explained an additional 9% (
P = 0.03) and 8% (
P = 0.03) of variation in ΔM for the a
M and a
P180 models, respectively. The corresponding results for ΔAL were 11% (
P = 0.02) for both the a
M and a
P180 models. Furthermore, the introduction of the normality-transformed standard deviation of scene defocus variable explained an additional 10% (
P = 0.02) of the variation in ΔM for both the a
M and a
P180 models. The corresponding results for ΔAL were 21% (
P = 0.001), 18% (
P = 0.001), 19% (
P < 0.01), and 20% (
P < 0.001) for a
M, a
J0, a
P90, and a
P180 models, respectively.
The myopia progression in children was different in high versus low scene defocus dispersion (i.e., SD
D) and steep versus flat peripheral refraction (a
M, a
J0, and a
P180), by median split. (
Fig. 5). When the participants were equally divided into two groups according to their SD
D, the RPRE of participants with low SD
D was significantly associated with ∆M with improved correlation coefficients (low SD
D: a
M vs. ∆M:
r = −0.58,
P < 0.01; a
J0 vs. ∆M:
r = −0.50,
P = 0.01; a
P180 vs. ∆M:
r = −0.62,
P = 0.001), while those with high SD
D were independent of ∆M (high SD
D: a
M vs. ∆M:
r = −0.14,
P = 0.49; a
J0 vs. ∆M:
r = −0.27,
P = 0.20; a
P180 vs. ∆M:
r = −0.21,
P = 0.32).