# The art of nomograms

- Samuel Arba Mosquera
^{1, 2, 3}Email author, - Diego de Ortueta
^{4}and - Shwetabh Verma
^{2, 5, 6, 7}

**Received: **3 August 2017

**Accepted: **8 January 2018

**Published: **25 January 2018

## Abstract

### Background

To retrospectively analyse strategies for adjusting refractive surgery plans with reference to the preoperative manifest refraction.

### Methods

We constructed seven nomograms based on the refractive outcomes (sphere, cylinder, axis [SCA]) of 150 consecutive eyes treated with laser in situ keratomileusis for myopic astigmatism. We limited the initial data to the SCA of the manifest refraction. All nomograms were based on the strategy: if for x diopters (D) of attempted metric, y D is achieved; we can reverse this sentence and state for achieving y D of change in the metric, x D will be planned. The effects of the use of plus or minus astigmatism notation, spherical equivalent, sphere, principal meridians notation, cardinal and oblique astigmatism, and astigmatic axis were incorporated.

### Results

All nomograms detected subtle differences in the spherical component (*p* < 0.0001). Nomograms 5 and 7 (using power vectors) and 6 (considering axis shifts) detected significant astigmatic differences (nomogram 5, *p* < 0.001; nomogram 6, *p* < 0.05; nomogram 7, *p* < 0.005 for cardinal astigmatism, *p* = 0.1 for oblique astigmatism). We observed mild clinically relevant differences (~ 0.5 D) in sphere or astigmatism among the nomograms; differences of ~ 0.25 D in the proposals for sphere or cylinder were not uncommon. All nomograms suggested minor improvements versus actual observed outcomes, with no clinically relevant differences among them.

### Conclusions

All nomograms anticipated minor improvements versus actual observed outcomes without clinically relevant differences among them. The minimal uncertainties in determining the manifest refraction (~ 0.6 D) are the major limitation to improving the accuracy of refractive surgery nomograms.

## Keywords

## Background

Nomograms have been used from the beginning of refractive surgery. In 1998, Yang et al. [1] evaluated a commercially available neural network program for calculating photorefractive keratectomy treatment nomograms and concluded that neural networks offer a potential means of segmenting and refining treatment nomograms to account for patient demographics, preoperative examinations, surgeon style, and equipment bias.

The primary reasons for developing nomograms in refractive surgery are to avoid surgical retreatments [2] and achieve an optimum refractive result. The need for having nomograms is catalysed by variability. To achieve optimum results, the laser systems available in the markets can be optimized based on their technical specifications, but the results of refractive surgery also depend on the subtle differences among lasers, surgeons, operation room environment and the patient population and demographics. Thus, to compensate for these differences and obtain optimum results, different values are used by surgeons when compared to the actual patient data. These adjustments are determined using Nomograms. Nomograms are designed after precisely analysing a range of patient data based on several factors contributing variability. The precision and success of a nomogram rests on several factors such as valid algorithm design, accurate patient data and correct grouping of the data. In addition to this, nomograms must be constantly updated, as lasers, surgeons, and techniques change over time.

Feltham and Wolfe, who retrospectively analysed the effects of ablation size, refractive errors, patient age, and corneal curvature on the retreatments, found that the older the patient and the larger the refractive error, the greater the risk of not achieving a residual refractive error of ±0.50 diopter (D) 3 months postoperatively and that corneas that are steeper preoperatively have a greater risk of retreatment.

Most reported nomograms have provided accurate results for myopic spheres; however, for hyperopic spheres and astigmatic outcomes the results are less predictable. Moniz and Fernandes [3] analysed a nomogram for treating astigmatism with laser in situ keratomileusis (LASIK). Similarly, Alpins and Goggin [4] provided a method to analyse refractive outcomes of astigmatic refractive surgery. The Alpins methodology uses three principal vectors and the various ratios among them provide an aggregate analysis of astigmatic change with parallel indices for spherical correction. A comparative analysis using arithmetic and vectorial means and necessary nomogram adjustments for refining spherical and astigmatic treatments can also be derived. These advanced techniques, together with their suitability for statistical analysis, comprehensively address the outcome analysis requirements of the entire cornea and the refractive correction for examining successful refractive surgery outcomes.

However, nomograms are only useful if the effects of their proposed surgical plans can be positively verified. Gailitis [5] compared the outcomes between two different excimer laser platforms using optimized nomograms and reported good results for both platforms, but superior results were obtained with the platform using the more recent nomogram.

Nomograms can consider all kinds of variables for analysis. A recent nomogram proposed the coupling effects between preoperative high-order aberrations and refractive outcomes [6]. The authors found that patient satisfaction was slightly higher than that of patients who underwent previous laser refractive surgery at the same clinic and concluded that the advanced nomogram increased treatment accuracy regarding the uncorrected visual acuity (UDVA) and the mean postoperative refraction and reduced the rate of hyperopic overcorrection compared with earlier studies. The need for retreatment procedures decreased, and patient satisfaction was high.

There is no single approach for adjusting the surgical parameters based on retrospective analyses of previous outcomes. Arnalich-Montiel et al. [7] examined four systematic strategies and one intuitive approach for adjusting the ablation sphere in myopic wavefront LASIK with reference to preoperative manifest refraction. Surprisingly, they found that the postoperative manifest refraction spherical equivalent varied lesser when non-systematic, intuitive adjustments to the ablation sphere were used. There was a strong trend toward reduced variability in the results in patients with a larger wavefront diameter. The authors concluded that back-calculation to model results with different pre-treatment ablation adjustment strategies might be useful to eliminate unpromising new approaches before clinical trials.

A nomogram can be designed based on several factors. One can analyse the factors as variables that are used to design an equation that describes an empirical data precisely. Looking at this process purely mathematically, the more variables that are included in the analysis the better the fit becomes. However, there is no ideal criterion to judge the relevance of the analysed factors on the refractive outcomes. These factors are selected mostly based on scientific studies, common sense and sometimes even a feeling. Although before including other factors, reaching an optimum basis for the refractive surgery design is imperative. In the current study, we aim at analysing the basic criteria that can be used to develop a nomogram, the Sphere Cylinder and axis (SCA) component of the manifest refraction. We analyse seven systematic strategies for adjusting the surgical plan in refractive surgery with reference to the preoperative manifest refraction (SCA), with an aim to find an optimum fundamental nomogram that can be further optimized with several other popular factors affecting refractive outcomes.

## Methods

One hundred and fifty consecutive eyes of 75 patients that had been treated with the Amaris (SCHWIND eye-tech-solutions, Kleinostheim, Germany) “aberration neutral” (Aberration-Free™) aspheric ablation profiles were analysed retrospectively.

The inclusion criteria were a bilateral surgery on the same day targeted for emmetropia, preoperative bilateral corrected distance visual acuity (CDVA) ≥ 20/25 (logarithm of the minimum angle of resolution [logMAR] ≤ 0.1), and no signs of amblyopia.

The 6-month follow-up data were available for all 150 eyes (100%). The mean preoperative manifest defocus refraction was − 3.60 ± 1.54 diopters (D) (range, − 7.50 to − 1.25 D) and the mean preoperative manifest astigmatism was 0.79 ± 0.74 D (range, 0.00 to 4.00 D). We measured the corneal topography [8] for all eyes and obtained the corneal wavefront aberrations [9, 10] up to the 7th Zernike order (36 terms) (Keratron-Scout, OPTIKON2000, Rome, Italy), manifest refraction, UDVA, and CDVA. The measurements were performed preoperatively and at 1, 3, and 6 months postoperatively.

All ablations were non-wavefront-guided but were based on aspheric [11] aberration-neutral profiles (and not on the profiles proposed by Munnerlyn [12]) to balance the induction of spherical aberration [13, 14] (prolateness optimization [15]). This approach included a multidynamic aspheric transition zone, aberration and focus shift compensation due to ablation, pseudomatrix-based spot positioning, and enhanced compensation for the loss of efficiency [16]; all were based on theoretical equations validated with ablation models and clinical evaluations.

We used a 6.3-mm central, fully corrected optical zone (OZ) for myopia and a 7.0-mm OZ for high astigmatism, together with a variable transition size that was provided automatically by the laser depending on the planned refractive correction (range, 6.5–9.2 mm). The ablation was performed using the AMARIS excimer laser (SCHWIND eye-tech-solutions, Kleinostheim, Germany), which is a flying-spot laser system that uses a real ablative spot volume locally considered through a self-constructing algorithm that controls for the local repetition rates to minimize the thermal load of the treatment [17].

The AMARIS laser system works at a repetition rate of 750 Hz and produces a beam (size, 0.54 mm) (full-width-at-half-maximum) and a super-Gaussian spot profile [18, 19]. High-speed eye tracking (pupil and limbus tracker with cyclotorsional tracking [20]) with a 1050 Hz acquisition rate is accomplished with a 3-ms latency time [21].

We based our nomogram analyses on three manifest refraction values (sphere, cylinder, axis [SCA]), and for simplicity we ignored other known factors of refractive deviation such as treatment duration [22], in that increased dehydration leads to different laser tissue interaction attributes; OZ [23–26], in that larger zones provide more natural corneal shapes but smaller zones save ablated tissue; age [2, 24], in that corneal water content and accommodation decrease with age; treated eye [27, 28], in the sense of bilateral symmetry for astigmatism; keratometric values [2, 27], in that greater loss of efficiency affects ablation of steeper corneas; wavefront refraction [7, 29–31]; coupling effects [28]; effects of wave aberration on the manifest refraction [24, 32]; and neural adaptation [25, 33].

For all nomogram construction, we used linear correlation as a nomogram proposal: if for x D of attempted metric, y D is achieved; we can reverse this sentence and state that for achieving y D of change in the metric, x D will be planned.

### Nomogram 1: sphere and negative cylinder

In this approach, we correlated the achieved spherical change (considering astigmatism as the negative convention) with the attempted change and the magnitude of the achieved cylindrical change (vectorial analysis with astigmatism as the negative convention) with the attempted change (all at the corneal plane where the ablation is performed).

### Nomogram 2: sphere and positive cylinder

We correlated the achieved spherical change (considering astigmatism as the positive convention) with the attempted change and the magnitude of the achieved cylindrical change (vectorial analysis with astigmatism as the positive convention) with the attempted change (all at the corneal plane where the ablation is performed).

### Nomogram 3: spherical equivalent and astigmatism

We correlated the achieved spherical equivalent change with the attempted change, and the absolute magnitude of the achieved cylindrical change (vectorial analysis) with the attempted change (all at the corneal plane where the ablation is performed).

### Nomogram 4: principal meridians

_{1}and Cyl

_{2}are the meridional powers of the principal meridians. We correlated the achieved meridional power change (vectorial analysis) with the attempted power change for the principal meridians at the corneal plane where the ablation is performed.

### Nomogram 5: spherical equivalent and cardinal and oblique astigmatism

_{+}, M, C

_{x}].

_{+}the cardinal astigmatism, and C

_{x}the oblique astigmatism. The three components, respectively, represent the power of a Jackson crossed-cylinder with axes at 0 degree and 90 degrees, the spherical equivalent power, and the power of a Jackson crossed-cylinder with axes at 45 degrees and 135 degrees. We correlated the achieved spherical equivalent change with the attempted change, and the achieved cylindrical change (cardinal and oblique components analysis) with the attempted change (all at the corneal plane where the ablation is performed).

### Nomogram 6: spherical equivalent and astigmatism considering residual cyclotorsion

We correlated the achieved spherical equivalent change with the attempted change, and the achieved cylindrical change (cardinal and oblique components analysis) with the attempted change (all at the corneal plane where the ablation is performed). Moreover, for astigmatism we compared the achieved astigmatic axis to the attempted axis.

The factor that minimizes residual astigmatism is the one that minimizes the square root in the formula above, thus, the factor that minimizes the quadratic function is inside the square root.

As π/4 is much larger than typical cyclotorsion errors, it can be assumed that for practical purposes F is always positive.

### Nomogram 7: power vectors analysis

We correlated the achieved spherical equivalent change to the attempted change, and the achieved cylindrical change (separate cardinal and oblique components analyses) to the attempted change (all at the corneal plane where the ablation is performed).

### Statistical analyses

We assessed the statistical significance of the nomogram correction compared to the preoperative refraction using paired Student’s t-tests.

We assessed the statistical significance of the correlations using the Student’s t-tests. The coefficient of determination (r^{2}) was used, and the significance of the correlations has been evaluated considering a metric distributed approximately as t with N-2 degrees of freedom where N is the size of the sample.

The level of statistical significance was *p* < 0.05.

## Results

### Refractive outcomes

*p*< 0.0001) and mean residual astigmatism 0.16 ± 0.21 D (range, 0.00 to 0.75 D;

*p*< 0.001). In 85% of eyes (

*n*= 127), the spherical equivalent was within ±0.50 D of emmetropia, and in 97% of eyes (

*n*= 145), the astigmatism was within ±0.50 D of emmetropia (Table 1).

Comparison of refractive outcomes 6 months after surgery for all 150 eyes

Pre-op (Mean ± SD) | Post-op (Mean ± SD) |
| |
---|---|---|---|

Spherical equivalent (D) | −3.60 ± 1.54 | − 0.08 ± 0.36 | < 0.0001* |

Cylinder (D) | 0.79 ± 0.74 | 0.16 ± 0.21 | < 0.001* |

Predictability within ±0.50 D (%) | – | 85% for Seq 97% for Ast | – |

Predictability within ±1.00 D (%) | – | 97% for Seq 100% for Ast | – |

### Nomogram 1: sphere and negative cylinder

*p*< 0.0001), and the magnitude of the achieved cylindrical change (vectorial analysis with astigmatism as a negative convention) was correlated with the attempted one (

*p*< 0.0001) (Fig. 1).

The nomogram correction differed significantly from the preoperative refraction for sphere (*p* < 0.0001) with a mean nomogram correction of − 0.23 ± 0.39 D (range, − 1.49 D to + 0.76 D; median − 0.25 D) but not for astigmatism (*p* = 0.3).

### Nomogram 2: sphere and positive cylinder

*p*< 0.0001), and the magnitude of the achieved cylindrical change (vectorial analysis with astigmatism as a positive convention) was correlated with the attempted change (

*p*< 0.0001) (Fig. 2).

The nomogram correction differed significantly from the preoperative refraction for sphere (*p* < 0.0001) with a mean nomogram correction of − 0.23 ± 0.36 D (range, − 1.46 D to + 0.76 D; median − 0.25 D) but not for astigmatism (*p* = 0.3).

### Nomogram 3: spherical equivalent and astigmatism

*p*< 0.0001), and the absolute magnitude of the achieved cylindrical change (vectorial analysis) was correlated with the attempted change (

*p*< 0.0001) (Fig. 3).

The nomogram correction differed significantly from the preoperative refraction for sphere (*p* < 0.0001), with a mean nomogram correction of − 0.23 ± 0.36 D (range, − 1.47 D to + 0.76 D; median − 0.25 D) but not for astigmatism (*p* = 0.3).

### Nomogram 4: principal meridians

*p*< 0.0001) (Fig. 4), and the nomogram correction differed significantly from the preoperative refraction (

*p*< 0.0001), with a mean nomogram correction of − 0.23 ± 0.37 D (range, − 1.49 D to + 0.76 D; median − 0.25 D).

### Nomogram 5: spherical equivalent and cardinal and oblique astigmatism

*p*< 0.0001), and the achieved cylindrical change (cardinal and oblique components analysis) was correlated with the attempted change (

*p*< 0.0001) (Fig. 5).

The nomogram correction differed significantly from the preoperative refraction for sphere (*p* < 0.0001), with a mean nomogram correction of − 0.23 ± 0.36 D (range, − 1.47 D to + 0.76 D; median, − 0.25 D) and for astigmatism (*p* < 0.001), with a mean nomogram correction of + 0.02 ± 0.09 D (range, − 0.29 D to + 0.37 D; median, 0.00 D).

### Nomogram 6: spherical equivalent and astigmatism considering residual cyclotorsion

The achieved spherical equivalent change was correlated with the attempted change (*p* < 0.0001), and the absolute magnitude of the achieved cylindrical change (vectorial analysis) was correlated with the attempted change (*p* < 0.0001) (Fig. 3).

The nomogram correction differed significantly from the preoperative refraction for sphere (*p* < 0.0001) with a mean nomogram correction of − 0.23 ± 0.36 D (range, − 1.47 D to + 0.76 D; median − 0.25 D) and for astigmatism (*p* < 0.05), with a mean nomogram correction of − 0.02 ± 0.21 D (range, − 0.75 D to + 0.66 D; median 0.00 D).

A refraction of − 2.50 D sphere with − 1.50 D of astigmatism at 63 degrees should be planned as − 2.76 D sphere with − 1.49 D astigmatism at 63 degrees.

### Nomogram 7: power vectors analysis

*p*< 0.0001), and the achieved cylindrical change (separate cardinal and oblique components analyses) was correlated with the attempted change (

*p*< 0.0001 for cardinal astigmatism;

*p*< 0.0001 for oblique astigmatism) (Fig. 6).

The nomogram correction differed significantly from the preoperative refraction for sphere (*p* < 0.0001) with a mean nomogram correction of − 0.23 ± 0.36 D (range, − 1.47 D to + 0.76 D; median − 0.25 D) and for cardinal astigmatism (*p* < 0.005) with a mean nomogram correction of + 0.03 ± 0.10 D (range, − 0.29 D to + 0.37 D; median 0.00 D) but not for oblique astigmatism (*p* = 0.1).

### Comparison of the different nomogram proposals

Comparison of the different nomogram proposals

Attempted | Nom 1: Sph + negCyl | Nom 2: Sph + posCyl | Nom 3: SEq + Ast | Nom 4: Cyl | Nom 5: SEq + (C | Nom 6: Cyclotorsion | Nom 7: SEq + C | Max. internomogram difference | ||||||||||||||||||
---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|---|

Sphere (D) | Cylinder (D) | Axis (deg) | Sphere (D) | Cylinder (D) | Axis (deg) | Sphere (D) | Cylinder (D) | Axis (deg) | Sphere (D) | Cylinder (D) | Axis (deg) | Sphere (D) | Cylinder (D) | Axis (deg) | Sphere (D) | Cylinder (D) | Axis (deg) | Sphere (D) | Cylinder (D) | Axis (deg) | Sphere (D) | Cylinder (D) | Axis (deg) | Sphere (D) | Cylinder (D) | Axis (deg) |

−3.00 | 0.00 | 0 | − 3.25 | 0.02 | 0 | −3.29 | 0.02 | 0 | −3.28 | 0.02 | 0 | −3.27 | 0.00 | 0 | −3.26 | −0.02 | 0 | −3.28 | 0.02 | 0 | −3.23 | − 0.08 | 98 | 0.06 | 0.09 | 98 |

0.00 | −2.00 | 15 | −0.41 | −2.05 | 15 | −0.27 | − 2.04 | 15 | − 0.34 | − 2.04 | 15 | − 0.38 | − 1.93 | 15 | − 0.34 | − 2.03 | 15 | −0.35 | − 2.01 | 15 | − 0.39 | − 1.93 | 14 | 0.14 | 0.12 | 1 |

−3.00 | − 2.00 | 150 | − 3.25 | − 2.04 | 150 | − 3.18 | − 2.03 | 150 | −3.21 | − 2.03 | 150 | − 3.27 | − 1.92 | 150 | − 3.21 | − 2.02 | 150 | −3.22 | −2.00 | 150 | −3.28 | − 1.88 | 150 | 0.10 | 0.15 | 0 |

−1.50 | −1.00 | 45 | −1.83 | − 1.01 | 45 | − 1.78 | −1.01 | 45 | −1.81 | − 1.01 | 45 | − 1.82 | −0.96 | 45 | − 1.80 | −1.02 | 45 | −1.82 | − 1.00 | 45 | − 1.86 | − 0.91 | 47 | 0.08 | 0.11 | 2 |

−4.75 | −3.25 | 120 | −4.90 | −3.31 | 120 | −4.82 | −3.30 | 120 | −4.84 | −3.31 | 120 | −4.95 | −3.11 | 120 | −4.86 | −3.26 | 120 | −4.86 | −3.26 | 120 | −4.91 | −3.15 | 118 | 0.13 | 0.20 | 2 |

4.75 | −3.25 | 75 | 4.11 | −3.36 | 75 | 4.43 | −3.33 | 75 | 4.27 | −3.35 | 75 | 4.21 | −3.14 | 75 | 4.25 | −3.29 | 75 | 4.25 | −3.30 | 75 | 4.26 | −3.32 | 77 | 0.32 | 0.22 | 2 |

2.00 | −2.00 | 90 | 1.49 | −2.05 | 90 | 1.68 | −2.04 | 90 | 1.58 | −2.05 | 90 | 1.55 | −1.93 | 90 | 1.57 | −2.03 | 90 | 1.57 | −2.02 | 90 | 1.62 | −2.12 | 90 | 0.19 | 0.19 | 0 |

3.00 | −6.00 | 165 | 2.44 | −6.20 | 165 | 2.86 | −6.14 | 165 | 2.67 | −6.17 | 165 | 2.52 | −5.78 | 165 | 2.62 | −6.04 | 165 | 2.63 | −6.08 | 165 | 2.58 | −5.96 | 166 | 0.41 | 0.42 | 1 |

### Adverse events

No adverse events or complications were observed intraoperatively or postoperatively. No patient needed or requested retreatment of either eye.

## Discussion

In the current study, we used non-customized aberration-neutral profiles, i.e., the ablations were optimized to induce no change in wavefront aberration within the OZ other than sphere and cylinder components, leaving all existing higher order aberrations (HOAs) unchanged because the CDVA was unaffected by the pre-existing aberrations [25]. Thus, to compensate for the induced aberrations observed with other types of profile definitions [26], several sources of aberration might be considered; some of those sources of aberration are related to the loss of efficiency of the laser ablation for non-normal incidence [34–40].

The aim of the current study was to evaluate the differences among seven systematic strategies for adjusting the surgical refractive plan in reference to the preoperative manifest refraction, based on retrospective analysis of 150 eyes treated with the AMARIS system that used an Aberration-Free ablation profile. The advantage of the Aberration-Free ablation profile is that it aims to be neutral for HOAs, leaving the visual print of the patient as it was preoperatively with the best spectacle correction.

All patients were followed for 6 months postoperatively, and, although no nomogram adjustments were applied, no patient required retreatment. A study with a longer follow-up and more eyes would provide more proof of stable outcomes, even though the refractive spherical and astigmatic results were stable after 3 months.

The average residual defocus was about − 0.1 D and the residual cylinder about 0.2 D, with 85% of the eyes within 0.50 D and 97% within 1.0 D of emmetropia. This might be specific to the SCHWIND laser, which is calibrated for all treatments to achieve − 0.2 D of residual myopia.

Although this small series of treated eyes does not allow for definitive conclusions to be drawn or evidence-based statements, our preliminary results are promising. We compared seven different approaches for constructing a nomogram for planned refraction for laser refractive surgery. For all approaches, we limited the initial data to the SCA values of the manifest refraction. All nomogram proposals were derived based on the strategy: if for x D of attempted metric, y D is achieved; we can reverse this sentence and state that for achieving y D of change in the metric, x D will be planned.

Six nomograms considered only the SC values independent of the astigmatism axis, and nomogram 7 also considered the astigmatism axis.

Nomograms 1 and 2 were identical except for the fact that they used opposite cylinder conventions (negative versus positive). These approaches have the disadvantage of depending on the astigmatic sign, which provide different results.

Nomogram 3 is the most established method for analysing refractive outcomes, and we think that together with the intuitive approach, nomogram 3 is the most extended approach for constructing nomograms for laser refractive surgery. This approach is independent of the cylinder convention, but is based on enhancement of the spherical equivalent and astigmatism regression lines without considering the astigmatic axis. The problem is that this approach considers a perfect result as that in which the attempted correction is − 1 D cylinder at 0 degrees, and the achieved correction is − 1 D cylinder at 90 degrees. That means that it may fail for datasets with only low astigmatic values, small datasets, or datasets in which large rotations of the astigmatic angle are achieved. Nomogram 6 is a modification of nomogram 3, which further considers the effects of cyclotorsion on residual astigmatism [28], i.e., it targets reducing postoperative residual astigmatism instead of enhancing the achieved astigmatism change.

Nomogram 4 uses the astigmatism value only indirectly, since it is based on correlation of the changes in the meridional optical powers. This approach is independent from the cylinder convention, but it is based on the enhancement of the meridional regression line without considering the astigmatic axis. It can be useful for small datasets since the statistical power is enhanced using only one correlation with 2 N points, instead of two correlations with N points each.

Nomograms 5 and 7 use SCA values expressed as power vectors, whereas in nomogram 5 the astigmatism nomogram is based on analysis of the cardinal and oblique astigmatism components together, and in nomogram 7 the astigmatism nomogram is based on separate analyses of the cardinal and oblique astigmatism components. These approaches are independent of the used cylinder convention and based on enhancement of the regression lines considering the astigmatic axis.

Only nomograms 5 and 7, which used power vectors, and nomogram 6, which considered the effects of cyclotorsion on residual astigmatism [28], detected significant differences for astigmatism (nomogram 5, *p* < 0.001; nomogram 6, *p* < 0.05; and nomogram 7, *p* < 0.005 for cardinal astigmatism, but *p* = 0.1 for oblique astigmatism). This means that only by considering SCA, with emphasis on the axis, subtle astigmatic differences can be detected.

Comparing the different nomogram proposals, we observed mild significant differences of almost 0.5 D in sphere or astigmatism among the different nomogram proposals; differences of about 0.25 D in the proposals for sphere or cylinder were not uncommon.

Ditzen et al. [41] compared the correction values entered into the laser with the achieved change in refraction for these eyes and incorporated the effect of OZ size and patient age. Scatterplots comparing the laser settings to achieved postoperative refractions showed a clear 20% trend toward overcorrection. This trend increased with patient age and OZ diameter.

Sheludchenko and Fadeykina [42] studied the results of excimer laser correction for mixed astigmatism using bitoric nomograms for cylinder and standard monotoric procedures and found that bitoric ablations for astigmatic corrections resulted in fewer retreatments compared to standard monotoric procedures.

Anderson et al. [43] analysed the attempted versus achieved changes in refraction. Factors such as age, corneal thickness (pachymetry), preoperative spherical equivalent refraction, preoperative cylinder, and OZ were studied to evaluate their roles in predicting the refractive outcome 6 months after LASIK. The preoperative spherical equivalent refraction and OZ size were strong predictors of the 6-month refractive outcomes, that is, with the preoperative spherical equivalent refraction and OZ size, the postoperative refraction of the patient could be strongly estimated. However, there was no correlation between the postoperative refraction and the age, preoperative cylinder, or surgeon.

Caster et al. [44] evaluated the refractive results of conventional (non-wavefront) LASIK for treating myopia and myopic astigmatism using the Alcon LADARVision 4000 excimer laser system and nomogram adjustment techniques. They found that the outcomes improved substantially throughout the development of an accurate nomogram because of continually updated regression analysis of previous refractive results.

Zaldivar et al. [45] compared five nomogram refinements accounting for accommodation, use of a 7.0-mm OZ and a 9.5-mm transition zone, a targeted mean flap diameter of 10.5-mm, sequential interruption of the laser ablation, and cleaning of the interface. They concluded that the outcomes support the observation that five surgical and technical modifications of the hyperopic LASIK procedure resulted in excellent visual quality and refractive outcomes and a low regression rate.

Mrochen et al. [46] stated that nomograms are efficient tools for improving the predictability of refractive procedures by using statistical methods to analyse the preoperative and postoperative refractive data. The authors found that using individual nomograms significantly improved the predictability of the refractive outcome. However, theoretical investigation showed that homogeneous data distribution within cohorts is a key factor for predictable nomogram calculations. Nomograms are helpful for improving refractive outcomes. However, they are currently limited to about 90% within ±0.5 D of emmetropia.

The limitations of the current study included short follow-up and no clinical study groups for which the proposed nomogram adjustments have been planned. Long-term follow-up of these eyes will help determine the stability of these results.

Refined models to analyse the refractive outcomes and derivation of adjusted nomograms may increase the accuracy of the results; however, the high validation cost and the analysis of sufficiently large datasets to convey sufficient statistical power must be examined together with the potential benefits of nomograms.

Artificial intelligence applied to nomograms may also contribute to auto-updating of the proposals for the surgical plans, e.g., weighting the refractive outcomes based on the time they were planned (i.e., recent treatments weigh more than older ones).

We have not clinically validated the proposed nomograms, since on the one hand the results observed in this setting were already highly accurate (range residual spherical equivalent between − 1.12 D and + 0.75 D, and residual astigmatism below 0.75 D), and on the other hand no patient requested an enhancement.

Comparison of the virtual refractive outcomes for all nomograms

Post-op SEq (D) (Mean ± SD) (Range) | Post-op Cyl (D) (Mean ± SD) (Range) |
| |
---|---|---|---|

Actual outcomes | −0.08 ± 0.36 − 1.12 to + 0.75 | 0.16 ± 0.21 0.00 to 0.75 | – |

Nomogram 1: Sph + negCyl | 0.00 ± 0.35 − 1.20 to + 0.91 | 0.12 ± 0.17 0.00 to 0.75 | < 0.0001 |

Nomogram 2: Sph + posCyl | 0.00 ± 0.35 − 1.19 to + 0.95 | 0.12 ± 0.17 0.00 to 0.75 | < 0.0001 |

Nomogram 3: SEq + Ast | 0.00 ± 0.35 − 1.19 to + 0.93 | 0.12 ± 0.17 0.00 to 0.75 | < 0.0001 |

Nomogram 4: Cyl | 0.00 ± 0.35 − 1.20 to + 0.93 | 0.13 ± 0.17 0.00 to 0.78 | < 0.0001 |

Nomogram 5: SEq + (C | 0.00 ± 0.35 − 1.19 to + 0.93 | 0.14 ± 0.16 0.02 to 0.79 | < 0.0001 |

Nomogram 6: Cyclotorsion | 0.00 ± 0.35 − 1.19 to + 0.93 | 0.14 ± 0.16 0.00 to 0.74 | < 0.0001 |

Nomogram 7: SEq + C | 0.00 ± 0.35 − 1.20 to + 0.91 | 0.19 ± 0.17 0.05 to 0.80 | < 0.0005 |

One fact we must consider is that there are a number of sources of uncertainty (due to fluctuations in accommodation, uncertainty in point of final refraction, errors in working distance, uncertainties in the power of the trial lenses, errors in trial lens vertex distance) in the quoted values for manifest refraction [47]. Some of these can be analysed statistically, using the International Organization for Standardization (ISO) guidelines and can be used to estimate an uncertainty in the final recorded values. These uncertainties should be regarded as minimal because other factors such as unwanted accommodation are involved, which will affect the final outcome. The analysis showed that the uncertainty, which provides a 95% confidence level, would then be 0.6 D in refractive measurement, which is much larger than the 0.25 D steps in which manifest refraction used to be measured.

We based our analysis on simple power vectors [48] but not more comprehensive power matrices [49] that represent dioptric power in its full character because power vectors (with three components) are useful for adding, subtracting, and averaging powers of thin lens systems.

## Conclusion

The use of linear regression analysis to derive percentage and bias adjustments to attempted sphere and cylinder in laser vision correction has been well described in the past. Our analysis further suggests that for all the tested nomograms, minor improvements could be anticipated versus the actual observed outcomes, however without clinically relevant differences among them. The minimal uncertainties in determining the manifest refraction (~ 0.6 D) are the major limitation to improving the accuracy of refractive surgery nomograms; hence, more emphasis shall be placed on the estimation of manifest refraction with higher precision, for improving the treatment planning and accuracy in postoperative outcomes.

## Declarations

### Funding

None

### Authors’ contributions

DdO: Conception and design, critical revision of the article, final approval of the article, data collection; SAM: Conception and design, analysis and interpretation, writing the article, final approval of the article, statistical expertise, literature search; SV: Analysis and interpretation, writing the article, final approval of the article, literature search. All authors read and approved the final manuscript.

### Competing interests

S Arba-Mosquera and Shwetabh Verma are employees of SCHWIND eye-tech-solutions GmbH. D de Ortueta is a consultant for SCHWIND eye-tech-solutions GmbH. Neither author received any financial support associated with the conduct of this study.

**Open Access**This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.

## Authors’ Affiliations

## References

- Yang SH, Van Gelder RN, Pepose JS. Neural network computer program to determine photorefractive keratectomy nomograms. J Cataract Refract Surg. 1998;24:917–24.View ArticlePubMedGoogle Scholar
- Feltham MH, Wolfe RJ. Some variables to consider to avoid the need for LASIK surgical enhancements. Clin Exp Optom. 2000;83:76–81.View ArticlePubMedGoogle Scholar
- Moniz N, Fernandes ST. Nomogram for treatment of astigmatism with laser in situ keratomileusis. J Refract Surg. 2002;18(Suppl 3):S323–6.PubMedGoogle Scholar
- Alpins NA, Goggin M. Practical astigmatism analysis for refractive outcomes in cataract and refractive surgery. Surv Ophthalmol. 2004;49:109–22.View ArticlePubMedGoogle Scholar
- Gailitis RP. Comparison of LASIK outcomes with the Alcon LADARVision4000 and the VISX STAR S2 excimer lasers using optimized nomograms. J Refract Surg. 2005;21:683–90.PubMedGoogle Scholar
- Lapid-Gortzak R, van der Linden JW, van der Meulen IJ, Nieuwendaal CP. Advanced personalized nomogram for myopic laser surgery: first 100 eyes. J Cataract Refract Surg. 2008;34:1881–5.View ArticlePubMedGoogle Scholar
- Arnalich-Montiel F, Wilson CM, Morton SJ, Allan BD. Back-calculation to model strategies for pretreatment adjustment of the ablation sphere in myopic wavefront laser in situ keratomileusis. J Cataract Refract Surg. 2009;35:1174–80.View ArticlePubMedGoogle Scholar
- Mattioli R, Tripoli NK. Corneal geometry reconstruction with the Keratron videokeratographer. Optom Vis Sci. 1997;74:881–94.View ArticlePubMedGoogle Scholar
- Salmon TO. Corneal contribution to the wavefront aberration of the eye. PhD Dissertation. 1999.Google Scholar
- Mrochen M, Jankov M, Bueeler M, Seiler T. Correlation between corneal and total wavefront aberrations in myopic eyes. J Refract Surg. 2003;19:104–12.PubMedGoogle Scholar
- Mrochen M, Büeler M. Aspheric optics: physical fundamentals. Ophthalmologe. 2008;105:224–33.View ArticlePubMedGoogle Scholar
- Munnerlyn CR, Koons SJ, Marshall J. Photorefractive keratectomy: a technique for laser refractive surgery. J Cataract Refract Surg. 1988;14:46–52.View ArticlePubMedGoogle Scholar
- Yoon G, MacRae S, Williams DR, Cox IG. Causes of spherical aberration induced by laser refractive surgery. J Cataract Refract Surg. 2005;31:127–35.View ArticlePubMedGoogle Scholar
- Hersh PS, Fry K, Blaker JW. Spherical aberration after laser in situ keratomileusis and photorefractive keratectomy. Clinical results and theoretical models of etiology. J Cataract Refract Surg. 2003;29:2096–104.View ArticlePubMedGoogle Scholar
- Calossi A. Corneal asphericity and spherical aberration. J Refract Surg. 2007;23:505–14.PubMedGoogle Scholar
- Arba-Mosquera S, de Ortueta D. Geometrical analysis of the loss of ablation efficiency at non-normal incidence. Opt Express. 2008;16:3877–95.View ArticlePubMedGoogle Scholar
- Bende T, Seiler T, Wollensak J. Side effects in excimer corneal surgery. Corneal thermal gradients. Graefe Arch Ophthalmol. 1988;226:277–80.Google Scholar
- Huang D, Arif M. Spot size and quality of scanning laser correction of higher-order wavefront aberrations. J Cataract Refract Surg. 2002;28:407–16.View ArticlePubMedGoogle Scholar
- Guirao A, Williams DR, MacRae SM. Effect of beam size on the expected benefit of customized laser refractive surgery. J Refract Surg. 2003;19:15–23.PubMedGoogle Scholar
- Arba-Mosquera S, Merayo-Lloves J, de Ortueta D. Clinical effects of pure cyclotorsional errors during refractive surgery. Invest Ophthalmol Vis Sci. 2008;49:4828–36.View ArticlePubMedGoogle Scholar
- Bueeler M, Mrochen M. Simulation of eye-tracker latency, spot size, and ablation pulse depth on the correction of higher order wavefront aberrations with scanning spot laser systems. J Refract Surg. 2005;21:28–36.PubMedGoogle Scholar
- Wirbelauer C, Aurich H, Pham DT. Online optical coherence pachymetry to evaluate intraoperative ablation parameters in LASIK. Graefes Arch Clin Exp Ophthalmol. 2007;245:775–81.View ArticlePubMedGoogle Scholar
- Camellin M, Arba Mosquera S. Aspheric Optical Zones: The Effective Optical Zone with the SCHWIND AMARIS. J Refract Surg. 2011;27:135–46.Google Scholar
- Steinert RF, Fynn-Thompson N. Relationship between preoperative aberrations and postoperative refractive error in enhancement of previous laser in situ keratomileusis with the LADARVision system. J Cataract Refract Surg. 2008;34:1267–72.View ArticlePubMedGoogle Scholar
- Chen L, Artal P, Gutierrez D, Williams DR. Neural compensation for the best aberration correction. J Vis. 2007;7:9.1–9.View ArticleGoogle Scholar
- Marcos S, Cano D, Barbero S. Increase in corneal asphericity after standard laser in situ keratomileusis for myopia is not inherent to the Munnerlyn algorithm. J Refract Surg. 2003;19:S592–6.Google Scholar
- Arbelaez MC, Vidal C, Arba-Mosquera S. Bilateral Symmetry before and Six Months after Aberration-Free™ Correction with the SCHWIND AMARIS TotalTech Laser: Clinical Outcomes. J Optom. 2010;3:20–8.Google Scholar
- De Ortueta D, Haecker C. Laser in situ keratomileusis for mixed astigmatism using a modified formula for bitoric ablation. Eur J Ophthalmol. 2008;18:869–76.View ArticlePubMedGoogle Scholar
- Cheng X, Bradley A, Thibos LN. Predicting subjective judgment of best focus with objective image quality metrics. J Vis. 2004;4:310–21.View ArticlePubMedGoogle Scholar
- Marsack JD, Thibos LN, Applegate RA. Metrics of optical quality derived from wave aberrations predict visual performance. J Vis. 2004;4:322–8.View ArticlePubMedGoogle Scholar
- Watson AB, Ahumada AJ Jr. Predicting visual acuity from wavefront aberrations. J Vis. 2008;8(4):17.1–19.Google Scholar
- Subbaram MV, MacRae SM. Customized LASIK treatment for myopia based on preoperative manifest refraction and higher order aberrometry: the Rochester nomogram. J Refract Surg. 2007;23:435–41.PubMedGoogle Scholar
- Artal P, Chen L, Fernández EJ, Singer B, Manzanera S, Williams DR. Neural compensation for the eye's optical aberrations. J Vis. 2004;4:281–7.View ArticlePubMedGoogle Scholar
- Mrochen M, Seiler T. Influence of corneal curvature on calculation of ablation patterns used in photorefractive laser surgery. J Refract Surg. 2001;17:S584–7.PubMedGoogle Scholar
- Jiménez JR, Anera RG, Jiménez del Barco L, Hita E. Effect on laser-ablation algorithms of reflection losses and nonnormal incidence on the anterior cornea. Appl Phys Lett. 2002;81:1521–3.View ArticleGoogle Scholar
- Jiménez JR, Anera RG, Jiménez del Barco L, Hita E, Pérez-Ocón F. Correlation factor for ablation algorithms used in corneal refractive surgery with gaussian-profile beams. Opt Express. 2005;13:336–43.Google Scholar
- Dorronsoro C, Cano D, Merayo-Lloves J, Marcos S. Experiments on PMMA models to predict the impact of corneal refractive surgery on corneal shape. Opt Express. 2006;14:6142–56.View ArticlePubMedGoogle Scholar
- Dorronsoro C, Marcos S. Experiments on PMMA model to predict the impact of corneal refractive surgery on corneal shape: reply. Opt Express. 2007;15:7245–6.View ArticlePubMedGoogle Scholar
- Kwon Y, Choi M, Bott S. Impact of ablation efficiency reduction on post-surgery corneal asphericity: simulation of the laser refractive surgery with a flying spot laser beam. Opt Express. 2008;16:11808–21.View ArticlePubMedGoogle Scholar
- Kwon Y, Bott S. Postsurgery corneal asphericity and spherical aberration due to ablation efficiency reduction and corneal remodelling in refractive surgeries. Eye (Lond). 2009;23:1845–50.View ArticleGoogle Scholar
- Ditzen K, Handzel A, Pieger S. Laser in situ keratomileusis nomogram development. J Refract Surg. 1999;15(Suppl 2):S197–201.PubMedGoogle Scholar
- Sheludchenko VM, Fadeykina T. Comparative results between standard and bitoric nomograms for astigmatism correction. J Refract Surg. 2001;17(Suppl 2):S238–41.PubMedGoogle Scholar
- Anderson NJ, Hardten DR, Davis EA, Schneider TL, Samuelson TW, Lindstrom RL. Nomogram considerations with the Technolas 217A for treatment of myopia. J Refract Surg. 2003;19:654–60.PubMedGoogle Scholar
- Caster AI, Hoff JL, Ruiz R. Nomogram adjustment of laser in situ keratomileusis for myopia and myopic astigmatism with the Alcon LADARVision system. J Refract Surg. 2004;20:364–70.PubMedGoogle Scholar
- Zaldivar R, Oscherow S, Bains HS. Five techniques for improving outcomes of hyperopic LASIK. J Refract Surg. 2005;21(Suppl 5):S628–32.PubMedGoogle Scholar
- Mrochen M, Hafezi F, Iseli HP, Löffler J, Seiler T. Nomograms for the improvement of refractive outcomes. Ophthalmologe. 2006;103:331–8.View ArticlePubMedGoogle Scholar
- Smith G. Refraction and visual acuity measurements: what are their measurement uncertainties? Clin Exp Optom. 2006;89:66–72.View ArticlePubMedGoogle Scholar
- Thibos LN, Horner D. Power vector analysis of the optical outcome of refractive surgery. J Cataract Refract Surg. 2001;27:80–5.View ArticlePubMedGoogle Scholar
- Harris WF. Power vectors versus power matrices, and the mathematical nature of dioptric power. Optom Vis Sci. 2007;84:1060–3.View ArticlePubMedGoogle Scholar