Comparative Analysis of Corneal Higher-Order Aberrations after Laser-Assisted In Situ Keratomileusis, Photorefractive Keratectomy, and Small Incision Lenticule Extraction with Correlations to Change in Myopic Q-Value and Spherical Equivalent with and without Astigmatism

Background: This retrospective chart review compared the higher-order aberrations (HOAs) among photorefractive keratectomy (PRK), laser-assisted in situ keratomileusis (LASIK), and small incision lenticule extraction (SMILE) alongside changes in spherical equivalent (SEQ) and corneal shape (Q-value). Methods: Analyzing 371 myopic eyes, including 154 LASIK, 173 PRK, and 44 SMILE cases, Pentacam imaging was utilized pre-operatively and at one-year post-operative visits. Results: All procedures resulted in 100% of patients achieving an uncorrected distance visual acuity (UDVA) of 20/40 or better, with 87% of LASIK and PRK, and 91% of SMILE patients having 20/20 or better. Significant increases in HOAs were observed across all procedures (p < 0.05), correlating positively with SEQ and Q-value changes (LASIK (0.686, p < 0.05), followed by PRK (0.4503, p < 0.05), and SMILE (0.386, p < 0.05)). Vertical coma and spherical aberration (SA) were the primary factors for heightened aberration magnitude among the procedures (p < 0.05), with the largest contribution in SMILE, which is likely attributed to the centration at the corneal apex. Notably, PRK showed insignificant changes in vertical coma (−0.197 µm ± 0.0168 to −0.192 µm ± 0.0198, p = 0.78), with an increase in oblique trefoil (p < 0.05). Conclusions: These findings underscore differences in HOAs among PRK, LASIK, and SMILE, helping to guide clinicians.


Introduction
By 2050, it is estimated that nearly 4.8 billion individuals, equivalent to approximately 50% of the world's population, will be affected by myopia [1].In response to this growing challenge, the ophthalmic field has witnessed transformative advancements with the LASIK: Flaps were created to be 8.7-9 mm in diameter and 100 µm thick.They were completed using the FS200 (Alcon Laboratories, Inc., Fort Worth, TX, USA) femtosecond laser.The WaveLight EX500 Excimer Laser System (Alcon Laboratories, Inc., Fort Worth, TX, USA) was the excimer laser used for corneal ablation.Ablation zones were 6.5 mm with a transition zone of up to 9.0 mm.Post-operative treatment included topical ofloxacin or moxifloxacin (0.3 or 0.5%, respectively) and prednisolone acetate 1%.The steroid was to be used every hour for the remainder of the day following the procedure.The steroid and antibiotic were then used q.i.d. for 1 week.
PRK: The corneal epithelium was removed after 20 s application of 18% alcohol.This was followed by ablation using WaveLight EX500 Excimer Laser System (Alcon Laboratories, Inc., Fort Worth, TX, USA) with an ablation a zone of 6.5 mm central ablation and up to a 9 mm transitional zone.If the ablation depth went beyond 65 µm, then Mitomycin C 0.02% was given for approximately 20 s.Post-operatively, an Acuvue Oasys (Johnson & Johnson Vision Care, Inc., Jacksonville, FL, USA) bandage contact was applied for five to seven days, until epithelialization was complete.Topical moxifloxacin 0.5% was prescribed for four drops daily for a total of seven days, and prednisolone acetate was prescribed at four drops daily for one month.After one month, fluorometholone 0.1% drops were applied three times per day for three weeks and then tapered: b.i.d. for another three weeks and then q.d. for three weeks.SMILE: The VisuMax 500 kHz femtosecond laser (Carl Zeiss Meditec, Jena, Germany) was used to create the lenticule, which was preset to be either 6.0 mm with a 0.5 mm transition zone for toric treatments or 6.5 mm for spherical treatment.This was completed by using a laser energy of 125 nJ, with cap thickness and diameter inputted at 120 µm and 7.5 mm, respectively.A 3.5 mm superior incision was used at a 90-degree meridian with a side-cut angle of 90 degrees.The measurements for lenticular spot separation, side cut, cap cut, and cap angle were 3.7 µm, 2.0 µm, 3.8 µm, and 2.0 µm, respectively.Post-operative drops included topical ofloxacin 0.3% or moxifloxacin 0.5% and prednisolone acetate 1%.The antibiotic was dosed q.i.d. for 1 week, and the steroid was tapered: q.i.d. for one week, b.i.d. for one week, and q.d. for two weeks.

Study Measurements
Measurements included root-mean-square (RMS) values for the full thickness of the cornea.The primary outcome variables were visual metrics one year post-surgery, including uncorrected distance visual acuity (UDVA), corneal HOAs, manifest refraction, changes in corrected distance visual acuity (CDVA), SEQ, and safety/efficacy indices.The HOAs analyzed, which were characterized by their respective Zernicke coefficients, included horizontal trefoil (Z3,3), horizontal coma (Z3,1), vertical coma (Z3,−1), oblique trefoil (Z3,−3), and spherical aberration (SA) (Z04,4), which were based on a 6 mm optical zone.Differences between pre-operative and post-operative HOA values were analyzed for statistical significance.Additionally, ratios of change in HOAs to change in SEQ (∆HOA/∆SEQ), along with changes in Q-value (∆Q), and SA (∆SA) were measured and examined.

Statistical Methods and Data Analysis
The statistical analysis, along with the generation of graphs and tables, were conducted using Microsoft Excel 2016 (Microsoft Inc., Redmond, VA, USA) and IBM SPSS Statistics for Windows, Version 29.0 (IBM Corp., Armonk, NY, USA).Basic descriptive statistical computations, including means and standard deviations were carried out to summarize the data.The Shapiro-Wilk test was employed to assess the normality of data distributions, guiding the choice between parametric and non-parametric tests.For quantitative demographic data that conformed to normal distribution and homogeneity of variances, one-way ANOVA and independent samples t-tests were utilized for comparisons across the three surgical procedures.Qualitative or categorical data were analyzed using Chi-squared tests to evaluate the distribution of frequencies across different categories.
In instances where data did not meet the assumptions of normality, non-parametric tests were applied.The Wilcoxon rank-sum test was utilized to compare median values of non-normally distributed variables between two independent groups.For multiple pairwise comparisons, the Bonferroni adjustment was employed to correct for the increased risk of Type I errors.
Pre-operative data, post-operative data, and their associated delta values were assessed through generalized estimating equations (GEEs) to account for inter-eye variability.Statistical significance was determined to be p < 0.05.

Pre-Operative Characteristics
The levels of myopic sphere and SEQ correction were greater for SMILE patients, being recorded at −5.185 D ± 0.350 and −5.406 D ± 0.346, respectively (both p < 0.001), while a reduced pre-operative cylinder was observed, at −0.443 D ± 0.072 (p = 0.011) compared to PRK and LASIK.

Visual Outcomes
Across all procedures, 100% of patients achieved a UDVA of 20/40 or better.Specifically, 87% of LASIK and PRK patients, along with 91% of SMILE patients, attained a UDVA of 20/20 or better (Figure A1).CDVA analysis revealed a loss of one Snellen line in 9% of LASIK, 22% of PRK, and 16% of SMILE patients.Conversely, 24% of LASIK, 16% of PRK, and 25% of SMILE patients experienced an improvement of one Snellen line in CDVA.
Notably, no patient across all procedures had a change of two or more Snellen lines in CDVA (Figure A3).
Safety indices (post-operative BCVA/pre-operative BCVA) for LASIK, PRK, and SMILE were 1.027, 1.003, and 1.031, respectively.Efficacy indices (post-operative UCVA/preoperative BCVA) were recorded at 0.993 for LASIK, 0.974 for PRK, and 0.968 for SMILE.Regarding procedure accuracy, 93% of LASIK, 95% of PRK, and 86% of SMILE patients were within 0.5 D of the intended SEQ.When the range was extended to 1.0 D, these figures jumped to 99% for LASIK, 100% for PRK, and 100% for SMILE.
Table 2. Mean values of total HOAs and subtypes pre-and post-procedure.This table presents the mean values of total HOAs and their five subtypes pre-and post-LASIK, PRK, and SMILE procedures.The table also includes standard error values for both pre-test and post-test measurements, along with p-values to assess the significance of changes between pre-and post-procedure values.PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations.Pairwise comparisons indicated significant pre-test differences between SMILE and PRK and between PRK and LASIK for total HOAs (p-value < 0.05).Post-test, LASIK differed significantly from both PRK and SMILE.Similarly, vertical coma showed significant pre-test differences in LASIK vs. PRK and PRK vs. SMILE, with all post-test comparisons (LASIK vs. PRK, LASIK vs. SMILE, PRK vs. SMILE) also revealing significant differences (p-values < 0.05) (Figure 1).(LASIK vs. PRK, LASIK vs. SMILE, PRK vs. SMILE) also revealing significant differences (p-values < 0.05) (Figure 1).There were statistically significant differences in the median changes in total HOAs, spherical aberration, and vertical coma across all procedures (p < 0.05).Significant differences in total HOAs and vertical coma changes were identified post-SMILE compared to both LASIK and PRK (p < 0.05).Additionally, significant changes in spherical aberration and oblique trefoil were noted between LASIK and PRK (p < 0.05).SMILE (0.2880 µm ± 0.1607) resulted in a median increase in total HOAs that was greater than those of both LASIK (0.2070 µm ± 0.1974, effect size: 0.335) and PRK (0.1680 µm ± 1.3092, effect size: 0.410).For spherical aberration, PRK showed a median increase of 0.1460 µm ± 0.1666 compared to that of LASIK (0.1160 µm ± 0.1665, effect size: 0.160).The median change in vertical coma with SMILE (0.2240 µm ± 0.2403) was greater than both LASIK (effect size: −0.4786) and PRK (effect size: −0.559).Moreover, the median change in oblique trefoil with PRK (−0.0318 µm ± 0.1802) was greater than with LASIK (effect size: −0.194) (Figure 2).There were statistically significant differences in the median changes in total HOAs, spherical aberration, and vertical coma across all procedures (p < 0.05).Significant differences in total HOAs and vertical coma changes were identified post-SMILE compared to both LASIK and PRK (p < 0.05).Additionally, significant changes in spherical aberration and oblique trefoil were noted between LASIK and PRK (p < 0.05).SMILE (0.2880 µm ± 0.1607) resulted in a median increase in total HOAs that was greater than those of both LASIK (0.2070 µm ± 0.1974, effect size: 0.335) and PRK (0.1680 µm ± 1.3092, effect size: 0.410).For spherical aberration, PRK showed a median increase of 0.1460 µm ± 0.1666 compared to that of LASIK (0.1160 µm ± 0.1665, effect size: 0.160).The median change in vertical coma with SMILE (0.2240 µm ± 0.2403) was greater than both LASIK (effect size: −0.4786) and PRK (effect size: −0.559).Moreover, the median change in oblique trefoil with PRK (−0.0318 µm ± 0.1802) was greater than with LASIK (effect size: −0.194) (Figure 2).

Correlation Analysis
Analysis indicated non-normal distribution for ∆ SEQ and ∆ total HOAs in LASIK and PRK, with p-values of <0.05 for LASIK and PRK.In contrast, for SMILE, ∆ SEQ and ∆ total HOAs showed normal distribution, with p-values of 0.11 and 0.67, respectively.In the LASIK cohort, ∆ total HOAs showed a Spearman's correlation of 0.529 with ∆ SEQ (p = <0.05,95% CI: 0.033 to 0.056).For PRK, the correlation between ∆ total HOAs and ∆ SEQ was 0.373 (p < 0.05, 95% CI: 0.034 to 0.070).Similarly in the SMILE cohort, ∆ total HOAs exhibited a Pearson's correlation of 0.570 with ∆ SEQ (p < 0.05, 95% CI: 0.033 to 0.085) (Figure 3).Across the three procedures, significant positive correlations were consistently observed between changes in SA and Q-value.Δ SA versus Δ Q in LASIK displayed the strongest Spearman's correlation of 0.686 (p < 0.05, 95% CI: 0.592 to 0.762).In PRK, a moderate positive Spearman's correlation of 0.4503 (p < 0.05, 95% CI: 0.323 to 0.562) was noted, while SMILE showed a relatively weaker yet significant correlation of 0.386 (p < 0.05, 95% CI: 0.101 to 0.613).A comparison between the correlations of LASIK and PRK indicated a significant difference in correlation strengths between Δ SA and Δ Q when comparing LASIK to both PRK (z = 3.17, p < 0.05) and SMILE (z = 2.46, p < 0.05).However, no significant difference was observed between PRK and SMILE (z = 0.446, p = 0.66) (Figure 4).

Discussion
Variations in pre-operative sphere, cylinder, and SEQ within the SMILE cohort can be attributed to the capability of SMILE to correct up to 10.0 D of myopia, compared to LASIK or PRK, which correct up to 8.0 D and 9.75 D, respectively [5,6].While LASIK and PRK procedures can correct up to 6 D of astigmatism, SMILE can only correct up to 3 D [7], accounting for the lower average cylinder observed in pre-operative SMILE patients (Table 1).Similarly, PRK showed significant differences in pre-operative pachymetry, keratometry measurements, total HOAs, and vertical coma, as individuals with a pachymetry of <500 µm were excluded from LASIK and SMILE, leading to a selection of PRK for individuals with thinner corneas.Across the three procedures, significant positive correlations were consistently observed between changes in SA and Q-value.∆ SA versus ∆ Q in LASIK displayed the strongest Spearman's correlation of 0.686 (p < 0.05, 95% CI: 0.592 to 0.762).In PRK, a moderate positive Spearman's correlation of 0.4503 (p < 0.05, 95% CI: 0.323 to 0.562) was noted, while SMILE showed a relatively weaker yet significant correlation of 0.386 (p < 0.05, 95% CI: 0.101 to 0.613).A comparison between the correlations of LASIK and PRK indicated a significant difference in correlation strengths between ∆ SA and ∆ Q when comparing LASIK to both PRK (z = 3.17, p < 0.05) and SMILE (z = 2.46, p < 0.05).However, no significant difference was observed between PRK and SMILE (z = 0.446, p = 0.66) (Figure 4).Across the three procedures, significant positive correlations were consistently ob served between changes in SA and Q-value.Δ SA versus Δ Q in LASIK displayed the strongest Spearman's correlation of 0.686 (p < 0.05, 95% CI: 0.592 to 0.762).In PRK, a mod erate positive Spearman's correlation of 0.4503 (p < 0.05, 95% CI: 0.323 to 0.562) was noted while SMILE showed a relatively weaker yet significant correlation of 0.386 (p < 0.05, 95% CI: 0.101 to 0.613).A comparison between the correlations of LASIK and PRK indicated a significant difference in correlation strengths between Δ SA and Δ Q when comparing LASIK to both PRK (z = 3.17, p < 0.05) and SMILE (z = 2.46, p < 0.05).However, no signifi cant difference was observed between PRK and SMILE (z = 0.446, p = 0.66) (Figure 4).

Discussion
Variations in pre-operative sphere, cylinder, and SEQ within the SMILE cohort can be attributed to the capability of SMILE to correct up to 10.0 D of myopia, compared to LASIK or PRK, which correct up to 8.0 D and 9.75 D, respectively [5,6].While LASIK and PRK procedures can correct up to 6 D of astigmatism, SMILE can only correct up to 3 D [7], accounting for the lower average cylinder observed in pre-operative SMILE patients (Table 1).Similarly, PRK showed significant differences in pre-operative pachymetry, ker atometry measurements, total HOAs, and vertical coma, as individuals with a pachymetry of <500 µm were excluded from LASIK and SMILE, leading to a selection of PRK for indi viduals with thinner corneas.

Discussion
Variations in pre-operative sphere, cylinder, and SEQ within the SMILE cohort can be attributed to the capability of SMILE to correct up to 10.0 D of myopia, compared to LASIK or PRK, which correct up to 8.0 D and 9.75 D, respectively [5,6].While LASIK and PRK procedures can correct up to 6 D of astigmatism, SMILE can only correct up to 3 D [7], accounting for the lower average cylinder observed in pre-operative SMILE patients (Table 1).Similarly, PRK showed significant differences in pre-operative pachymetry, keratometry measurements, total HOAs, and vertical coma, as individuals with a pachymetry of <500 µm were excluded from LASIK and SMILE, leading to a selection of PRK for individuals with thinner corneas.
All procedures resulted in increased total HOAs, aligning with established outcomes for corneal refractive surgery due to corneal shape modification [8,9].The increase seen in both post-LASIK and SMILE vertical coma may be linked to the inherent structural manipulation of the procedures, which could introduce more vertical asymmetry than PRK's surface ablation technique [10].Effect size analysis further revealed that SMILE, like its impact on total HOAs, had a more pronounced effect on vertical coma than LASIK and PRK.This finding is supported by Chen X et al. [10], who reported higher vertical coma after SMILE versus wavefront-guided femtosecond LASIK.These differences might be associated with SMILE's unique centration technique, which employs the corneal vertex as a reference point, diverging from the pupillary center used in LASIK and PRK, aligning with findings from the existing literature [6,[10][11][12][13][14].In LASIK and PRK, variations in pupil size relative to treatment zones result in distinct HOA inductions.Further context to these findings is provided in Table A1, which outlines the rates of HOA induction found in previous studies.
The changes in SA and Q-value following LASIK, PRK, and SMILE could be explained by the central corneal flattening inherent in these procedures [15].While this study observed variations within individual pre-and post-test groups, the overall lack of significant differences between the pre-surgical groups or the post-surgical groups in SA across all procedures may indicate a broadly similar impact on this type of aberration.Interestingly, an assessment of the variations in post-vs.pre-surgical change in SA among the procedures via effect size analysis revealed nuanced differences, particularly between LASIK and PRK.PRK showed a median absolute change in SA of 1.54 times greater compared to LASIK, suggesting that while the general trend in aberration increase is similar, the extent can vary depending on the specific procedure.Miraftab et al. [16] found no significant differences in spherical aberration changes between LASIK and PRK.However, studies by Jahadi Hosseini et al. [17] and Russo et al. [18] noted a greater induction of spherical aberration in LASIK compared to PRK.LASIK's strong correlation of ∆SA/∆Q suggests that altering the Q-value might have a greater impact on minimizing spherical aberration induction with LASIK compared to PRK or SMILE.In contrast, rates of change in SA for a given Q for SMILE and PRK were not significantly different, suggesting that there is likely no optimal procedure among these two which would induce SA at a lower rate for a given change in Q-value [19,20].
The PRK group also demonstrated a greater induction of oblique trefoil, further characterized by a rank-biserial correlation of −0.19 and a 12-fold difference in relation to LASIK.However, despite the statistical significance of PRK's larger impact on oblique trefoil, the clinical implications are potentially limited due to the minor influence of oblique trefoil on overall retinal image quality [16].Biscevic et al. [21] also found that PRK induced greater trefoil changes compared to LASIK.Nonetheless, the unique impact of PRK on vertical coma and oblique trefoil may highlight a specific interaction between the technique and the aberrations.
One of the main limitations of this study is the disparity in sample sizes between the SMILE cohort and the LASIK/PRK cohorts.However, SMILE is a newer procedure at our institution, and we do not believe the uneven group sizes affected the results related to higher-order aberrations, mirroring the current literature [14,16].Additionally, it is important to emphasize that all of our patients were emmetropic, with satisfactory visual outcomes, and none of them required enhancement, reflecting that visual disturbances were due to HOAs not lower-order aberrations such as myopia or astigmatism.Some may argue that there is a selection bias, as the SMILE and PRK cohorts were not similar pre-operatively, but this is the apparent nature of patient selection, as those who typically undergo SMILE have a higher SEQ and sphere with a lower cylinder, while those who undergo PRK typically have thinner corneas [22].These parameters were determined jointly by patients and clinicians rather than being randomly assigned.While it could be argued that one eye per patient should have been randomly selected rather than using both eyes, the GEE formula was employed to rectify inter-eye variability.

Conclusions
Each of the three procedures exhibited a significant rise in total HOAs, with a positive correlation observed between SEQ, Q-value, and the magnitude of induced HOAs.Vertical coma and spherical aberration emerged as the primary contributors to the overall increase in aberration magnitude, with SMILE showing the most pronounced effect, potentially attributed to its unique corneal apex fixation compared to the central axis alignment in LASIK and PRK.Notably, PRK displayed a noteworthy stability in vertical coma, contrasting with the other procedures, and demonstrated a distinct elevation in oblique trefoil, although the clinical implications of the latter remain uncertain.Table A1.Comparative analysis of HOA induction after refractive procedures.This table presents a comprehensive comparison of the induction of higherorder aberrations (HOAs) based on our study findings and other published reports.The included studies utilized diverse imaging modalities to evaluate the impact of different refractive procedures on ocular aberrations.Variations in the imaging modalities, number of eyes examined, procedures performed, and follow-up duration are recorded across the studies.Results that were not statistically significant are denoted as "NS", while "NR" indicates data that were not reported in the original publication.PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations.

Figure 1 .
Figure 1.Visualization of cumulative HOA induction for each surgical procedure.The colored bars differentiate the surgical procedures: LASIK is in red, PRK is in gray, and SMILE is in black.Each aberration type has its own column, showcasing pre-operative values on the left and postoperative values on the right for each procedure.Error bars signify the standard error of the measurements.PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations.Statistical significance is indicated by an asterisk (*) above the bars, with p < 0.05.

Figure 1 .
Figure 1.Visualization of cumulative HOA induction for each surgical procedure.The colored bars differentiate the surgical procedures: LASIK is in red, PRK is in gray, and SMILE is in black.Each aberration type has its own column, showcasing pre-operative values on the left and post-operative values on the right for each procedure.Error bars signify the standard error of the measurements.PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations.Statistical significance is indicated by an asterisk (*) above the bars, with p < 0.05.

Figure 2 .
Figure 2. Differential HOA progression patterns post-refractive procedures.A composite analysis of total HOAs (A), spherical aberration (B), vertical coma (C), and oblique trefoil (D) showcases the variability in vision damage changes after LASIK, PRK, and SMILE surgeries.Medians are emphasized with white lines across the colored boxes-red for LASIK, blue for PRK, and black for SMILE.Asterisk (*) = statistical significance, Diamonds = individual data points that are beyond the whiskers of the box plots, PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations.

Figure 2 .
Figure 2. Differential HOA progression patterns post-refractive procedures.A composite analysis of total HOAs (A), spherical aberration (B), vertical coma (C), and oblique trefoil (D) showcases the variability in vision damage changes after LASIK, PRK, and SMILE surgeries.Medians are emphasized with white lines across the colored boxes-red for LASIK, blue for PRK, and black for SMILE.Asterisk (*) = statistical significance, Diamonds = individual data points that are beyond the whiskers of the box plots, PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations.

Figure 3 .
Figure 3. Correlation and regression analysis of Δ HOAs vs. Δ SEQ.This figure consists of three plots representing the correlation and regression analysis for the change in Δ HOAs against the change in Δ SEQ across LASIK (A), PRK (B), and SMILE (C).Each scatter plot with a linear regression line shows the equation of the regression line and the corresponding R 2 value.PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations, SEQ = spherical equivalent.

Figure 4 .
Figure 4. Correlation and regression analysis of Δ SA vs. Δ Q.This figure comprises three plots representing the correlation and regression analysis for the change in Δ SA against the change in Q-value (Δ Q) across LASIK (A), PRK (B), and SMILE (C).Each scatter plot with a linear regression line shows the equation of the regression line and the corresponding R 2 value.PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations, SA = spherical aberration.

Figure 3 .
Figure 3. Correlation and regression analysis of ∆ HOAs vs. ∆ SEQ.This figure consists of three plots representing the correlation and regression analysis for the change in ∆ HOAs against the change in ∆ SEQ across LASIK (A), PRK (B), and SMILE (C).Each scatter plot with a linear regression line shows the equation of the regression line and the corresponding R 2 value.PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations, SEQ = spherical equivalent.

Figure 3 .
Figure 3. Correlation and regression analysis of Δ HOAs vs. Δ SEQ.This figure consists of three plots representing the correlation and regression analysis for the change in Δ HOAs against the change in Δ SEQ across LASIK (A), PRK (B), and SMILE (C).Each scatter plot with a linear regression line shows the equation of the regression line and the corresponding R 2 value.PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations, SEQ = spherical equivalent.

Figure 4 .
Figure 4. Correlation and regression analysis of Δ SA vs. Δ Q.This figure comprises three plots representing the correlation and regression analysis for the change in Δ SA against the change in Q-value (Δ Q) across LASIK (A), PRK (B), and SMILE (C).Each scatter plot with a linear regression line shows the equation of the regression line and the corresponding R 2 value.PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations, SA = spherical aberration.

Figure 4 .
Figure 4. Correlation and regression analysis of ∆ SA vs. ∆ Q.This figure comprises three plots representing the correlation and regression analysis for the change in ∆ SA against the change in Q-value (∆ Q) across LASIK (A), PRK (B), and SMILE (C).Each scatter plot with a linear regression line shows the equation of the regression line and the corresponding R 2 value.PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations, SA = spherical aberration.

Figure A9 .
Figure A9.Comparison of total HOAs pre-and post-operatively.This figure visualizes the changes in total HOAs for each surgical procedure, represented by bars (pre-operative: diagonal stripe patterns, post-operative: solid colors).Red is used for LASIK, blue for PRK, and black for

Figure A9 .
Figure A9.Comparison of total HOAs pre-and post-operatively.This figure visualizes the changes in total HOAs for each surgical procedure, represented by bars (pre-operative: diagonal stripe patterns, post-operative: solid colors).Red is used for LASIK, blue for PRK, and black for

Figure A9 .
Figure A9.Comparison of total HOAs pre-and post-operatively.This figure visualizes the changes in total HOAs for each surgical procedure, represented by bars (pre-operative: diagonal stripe patterns, post-operative: solid colors).Red is used for LASIK, blue for PRK, and black for

Figure A9 .
Figure A9.Comparison of total HOAs pre-and post-operatively.This figure visualizes the changes in total HOAs for each surgical procedure, represented by bars (pre-operative: diagonal stripe patterns, post-operative: solid colors).Red is used for LASIK, blue for PRK, and black for Figure A9.Comparison of total HOAs pre-and post-operatively.This figure visualizes the changes in total HOAs for each surgical procedure, represented by bars (pre-operative: diagonal stripe patterns, post-operative: solid colors).Red is used for LASIK, blue for PRK, and black for SMILE.Error bars indicate standard error, and asterisks (*) denote statistically significant differences (p < 0.05).PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations.J. Clin.Med.2024, 13, x FOR PEER REVIEW 13 of 1

Figure A10 .
Figure A10.Comparison of spherical aberration (Z4, 0) pre-and post-operatively.This figure illustrates the changes in spherical aberration, with pre-operative measures in diagonal stripe patterns and post-operative in solid colors.Red indicates LASIK, blue PRK, and black SMILE.Error bars represent standard error, and asterisks (*) signify statistically significant differences (p < 0.05) PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations.

Figure A10 .
Figure A10.Comparison of spherical aberration (Z4, 0) pre-and post-operatively.This figure illustrates the changes in spherical aberration, with pre-operative measures in diagonal stripe patterns and post-operative in solid colors.Red indicates LASIK, blue PRK, and black SMILE.Error bars represent standard error, and asterisks (*) signify statistically significant differences (p < 0.05).PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations.

Figure A10 .
Figure A10.Comparison of spherical aberration (Z4, 0) pre-and post-operatively.This figure illustrates the changes in spherical aberration, with pre-operative measures in diagonal stripe patterns and post-operative in solid colors.Red indicates LASIK, blue PRK, and black SMILE.Error bars represent standard error, and asterisks (*) signify statistically significant differences (p < 0.05).PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations.

Figure A11 .
Figure A11.Comparison of vertical coma (Z3, −1) pre-and post-operatively.This figure showcases the changes in vertical coma for each procedure, with pre-operative values in diagonal stripes and post-operative in solid colors.Red denotes LASIK, blue PRK, and black SMILE.Error bars denote standard error, and asterisks (*) highlight statistically significant differences (p < 0.05).PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations.

Figure A12 .
Figure A12.Comparison in horizontal coma (Z3, 0) pre-and post-operatively.This figure presents the changes in horizontal coma, with pre-operative bars in diagonal stripes and post-operative in solid colors.Colors represent different procedures: red for LASIK, blue for PRK, and black for SMILE.Error bars indicate standard error.PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations.

Figure A12 .
Figure A12.Comparison in horizontal coma (Z3, 0) pre-and post-operatively.This figure presents the changes in horizontal coma, with pre-operative bars in diagonal stripes and post-operative in solid colors.Colors represent different procedures: red for LASIK, blue for PRK, and black for SMILE.Error bars indicate standard error.PRK = photorefractive keratectomy, LASIK = laserassisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations.

Figure A13 .
Figure A13.Comparison of horizontal trefoil (Z3, 3) pre-and post-operatively.This figure depicts the changes in horizontal trefoil for each procedure, with pre-operative values in diagonal stripes and post-operative in solid colors.Red signifies LASIK, blue PRK, and black SMILE.Error bars denote standard error.PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations.

Figure A13 . 1 Figure A14 .
Figure A13.Comparison of horizontal trefoil (Z3, 3) pre-and post-operatively.This figure depicts the changes in horizontal trefoil for each procedure, with pre-operative values in diagonal stripes and post-operative in solid colors.Red signifies LASIK, blue PRK, and black SMILE.Error bars denote standard error.PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations.J. Clin.Med.2024, 13, x FOR PEER REVIEW 15 of 1

Figure A14 .
Figure A14.Comparison of oblique trefoil (Z3, −3) pre-and post-operatively.This figure displays the changes in oblique trefoil, with pre-operative values in diagonal stripes and post-operative in solid colors.The color scheme includes red for LASIK, blue for PRK, and black for SMILE.Error bars represent standard error, and asterisks (*) indicate statistically significant changes (p < 0.05).PRK = photorefractive keratectomy, LASIK = laser-assisted in situ keratomileusis, SMILE = small incision lenticule extraction, HOAs = higher-order aberrations.