August 2007
Volume 48, Issue 8
Free
Clinical and Epidemiologic Research  |   August 2007
An Interval-Censored Model for Predicting Myopic Regression after Laser In Situ Keratomileusis
Author Affiliations
  • Yun-I Chen
    From the Department of Ophthalmology, Taipei Hospital, Taipei, Taiwan; the
    Institute of Preventive Medicine, the
    Department of Ophthalmology, National Taiwan University, Taipei, Taiwan; and the
  • Kuo-Liong Chien
    Institute of Preventive Medicine, the
  • I-Jong Wang
    Department of Ophthalmology, National Taiwan University, Taipei, Taiwan; and the
  • Amy Ming-Fang Yen
    Institute of Preventive Medicine, the
  • Li-Sheng Chen
    Institute of Preventive Medicine, the
  • Pi-Jung Lin
    Universal Eye Center, Taipei, Taiwan.
  • Tony Hsiu-Hsi Chen
    Institute of Preventive Medicine, the
    Centre of Biostatistics Consultation, and the
    Division of Biostatistics, Graduate Institute of Epidemiology, College of Public Health, and the
Investigative Ophthalmology & Visual Science August 2007, Vol.48, 3516-3523. doi:https://doi.org/10.1167/iovs.06-1044
  • Views
  • PDF
  • Share
  • Tools
    • Alerts
      ×
      This feature is available to authenticated users only.
      Sign In or Create an Account ×
    • Get Citation

      Yun-I Chen, Kuo-Liong Chien, I-Jong Wang, Amy Ming-Fang Yen, Li-Sheng Chen, Pi-Jung Lin, Tony Hsiu-Hsi Chen; An Interval-Censored Model for Predicting Myopic Regression after Laser In Situ Keratomileusis. Invest. Ophthalmol. Vis. Sci. 2007;48(8):3516-3523. https://doi.org/10.1167/iovs.06-1044.

      Download citation file:


      © ARVO (1962-2015); The Authors (2016-present)

      ×
  • Supplements
Abstract

purpose. A time-varying statistical model was proposed to predict the risk of regression toward myopia after laser in situ keratomileusis (LASIK) and to identify significant predictors within a time frame.

methods. A total of 615 eyes of 311 patients derived from a retrospective cohort who underwent LASIK in 2003 were analyzed. Refraction outcomes were recorded at 1 day, 1 week, and 1, 3, 6, 9, and 12 months or longer after LASIK. A cross-validated design was used, to split data into trained (n = 308) and validated (n = 307) data sets. These data sets were used in an interval-censored model to predict the probability of regression toward myopia and to assess the predictors including demographic features and preoperative and postoperative variables.

results. Myopia regression was observed in 164 (26.7%) of 615 eyes during the follow-up period of 12 months or longer after LASIK. Significant predictors for myopia regression after LASIK included preoperative manifest spherical equivalent (P = <0.0001), mean preoperative central corneal curvature (P = 0.001), size of optic zone (P = 0.0043), undercorrection (P = 0.04), and age (P = 0.0734). The risk of regression toward myopia after LASIK increased rapidly within 1 month, slowed down between 1 and 6 months, and became steady after 6 months, regardless of risk group. The risk of myopia regression up to 6 months after LASIK was 21% in average-risk eyes (based on all eyes).

conclusions. The proposed interval-censored model was useful not only for predicting the probability of myopia regression after LASIK but also for identifying the evolution of patients within low, moderate, and high-risk groups.

As both prevalence and severity of myopia have been noted over the past two decades in Taiwanese schoolchildren, 1 treatment with LASIK (laser in situ keratomileusis) has increasingly gained attention. However, the incidence of regression toward myopia (hereafter, myopia regression) after LASIK varies across studies, ranging from 5.5% to 27.7%. 2 3 4 5 6 7 Most of the previous studies in which investigators identified clinical correlates responsible for refractive outcome after LASIK were focused only on cases involving retreatment, 2 4 8 9 10 rather than on all patients with different degrees of regression toward myopia, from mild to severe. However, patients’ tolerance of residual myopia and their ophthalmologists’ threshold for determining the need for retreatment play vital roles in their decisions for retreatments. 4 As retreatment cases may represent only a small proportion of patients who had severe myopia regression after LASIK, this selection bias may lead to the failure of completely identifying clinical correlates related to myopia regression and may underestimate the risk of regression. 
In addition to including all cases with myopia regression, it is also valuable to throw light on whether and when myopia regression may occur after LASIK. However, it is very difficult, from a practical aspect, to know the exact time of onset of myopia regression after LASIK between two follow-up visits. An interval-censored model, 11 considering the degree of myopia regression occurring during a specific time interval after LASIK, was proposed to tackle this problem. The purposes of this study, conducted in all identifiable cases of myopia regression after LASIK, were therefore (1) to report the risk of myopia regression within a time frame; (2) to identify significant predictors; and (3) to develop a time-varying predictive model for estimating the probability of myopia regression by different time intervals of follow-up based on (1) and (2). 
Methods
Study Subjects
Subjects were enrolled from a retrospective cohort consisting of 372 patients with 735 eyes undergoing LASIK for myopia and myopic astigmatism by one experienced surgeon (PJL) with an excimer laser (EC-5000; Nidek Co., Ltd., Gamagori, Japan) at a franchise of ophthalmology clinic (Universal Eye Center, Taipei, Taiwan) between June and December 2003. The study was conducted according to the tenets of the Declaration of Helsinki of the World Medical Association regarding scientific research on human subjects. 
The criteria for LASIK eligibility consisted of (1) patients over 18 years of age with myopia or myopic astigmatism and with stable refraction for at least 1 year; (2) no autoimmune diseases (e.g., lupus, rheumatoid arthritis) or uncontrolled diabetes; (3) no degenerative, neurotrophic corneal disease or scarring; (4) no current pregnancy or nursing; (5) no HIV (human immunodeficiency virus) or hepatitis C virus infection; (6) no previous ocular trauma or surgery within 6 months of the LASIK procedure. 
Nasally oriented hinged flaps were created with a microkeratome (MK-2000; Nidek Co., Ltd.) with a flap diameter ranging from 8.5 to 9.0 mm and a thickness of 130 and 160 μm. The eyes were treated within an optic zone of 5.3 to 6.0 mm plus an additional transition zone, 1 mm wider than the optic zone. The standard laser ablation nomogram was used. A pulse repetition rate of 30 Hz and a laser energy pulse between 105 and 135 mJ were the laser parameters. The goal of our nomogram was emmetropia in all eyes. The targeted correction was calculated according to a customized clinical nomogram based on the Nidek photorefractive keratectomy (PRK) computer algorithm. The PRK nomogram was customized based on the prior surgeon’s experience with LASIK and the Nidek EC-5000 excimer laser. The targeted cylinder correction (minus cylinder format) entered into the computer was equal to the baseline cylinder by manifest refraction. The target spherical correction (minus-cylinder format) was equal to (greater than −10 D) or was increased by 2% (−6 to −10 D) or 5% (lower than −6 D) by the amount of manifest spherical component (minus cylinder format), but was further reduced by 1% of the cylinder component if the target cylinder was more than −3 D. 
After the ablation, the corneal flap was irrigated with a physiologic saline solution and then repositioned in the stromal bed with a cannula. A dripping-wet sponge was brushed gently across the cornea, lifting the excess fluid from the interface. The flap was allowed to adhere for 2 minutes. The slit beam of the laser was used to inspect possible wrinkles and debris. Plastic shields were placed after the surgery. The surgical site was inspected again with a standard slit-lamp biomicroscope to ensure proper flap position after blinking. 
The postoperative regimen included a topical dexamethasone/neomycin solution (Decaneomycin; Union, Taiwan) applied four times a day for the first week followed by the administration of 0.1% fluorometholone (Fluorometholone; Oasis, Taiwan) four times a day with tapering dose until the end of the month. Preservative-free artificial tears in polyvinyl alcohol 14 mg/mL (Refresh; Allergan, Irvine, CA) were used during the postoperative course, with the date of withdrawal suggested by the ophthalmologist. 
After exclusion of eyes with previous ocular surgeries (n = 4), those undergoing surgery with complications (n = 7), those lacking complete preoperative measurement records (n = 23), those having postoperative visits for <1 month (n = 78), and those having postoperative complications such as SOS (sands of the Sahara) syndrome (n = 8), 615 eyes (311 subjects) met the inclusion criteria and were used in the analysis. 
Variables and Outcome Measurement
Preoperative variables consisted of age, sex, manifest refraction measured by a phoropter with fogging, nomogram-adjusted refraction, and cycloplegic refraction. Refraction was recorded as sphere, cylinder, axis, and spherical equivalent (KR-7100; Topcon, Tokyo, Japan); central corneal thickness (CCT; Corneo-Gage Plus; Sonogage, Cleveland, OH); preoperative central corneal curvature (Mean K, KR 7100; Topcon); intraocular pressure (IOP); and Schirmer test II result (basic tear reflex). Also recorded were the variables calculated by the standard laser ablation nomograms, including the diameters of the optic zone (Size_OZ) and transition zone (Size_TZ, 1 mm wider than the optic zone), the ablation depths of the optic and transition zones, and corneal flap diameter and depth. Residual bed thickness (RBT) was calculated by subtracting the flap thickness and optic zone thickness from the CCT. 
Patients were examined at subsequent time points: 1 day, 1 week, and 1, 3, 6, 9, and 12 months or longer after the surgery. Manifest refraction results were recorded at each visit. 
Myopia regression was defined as a residual myopia of −1.0 D or greater, detected at the last visit, and a 0.5-D or greater shift toward myopia during follow-up visits. The eyes were otherwise defined as nonregressive. Undercorrection was defined as a residual refraction error of −1 D or greater detected 1 week after LASIK. 
To evaluate the influence of preoperative manifest refraction on postoperative refraction change, all eyes we stratified into three groups, in accordance with a previous study 6 : group 1 (low myopia), baseline refraction below −6.0 D; group 2 (moderate myopia), baseline refraction from −6.0 D to less than −10 D; and group 3 (high myopia), baseline refraction greater than or equal to −10 D. 
Most of the eyes were treated with an optic zone of 6.0 mm or above (352 eyes), in light of two previous studies, 12 13 and therefore the eyes were stratified into another two groups for assessment of the refraction change after LASIK: group 1, eyes with optic zone greater than or equal to 6.0 mm; group 2, eyes with optic zone <6.0 mm. 
To build up a predictive model for myopia regression, a cross-validation design was used with random allocation of one eye to the trained data set and the other eye to the validated data set for those individuals with both eyes used in the study (304 subjects; 608 eyes). Seven other individuals with only one eye used were randomly allocated to the trained data set or the validated data set. Thus, there were 308 eyes used in the trained model and 307 eyes in the validated model. 
The trained data set was first used to estimate regression coefficients 1 (Table 2) , which, together with the values of predictors from the validated data, were further used to predict occurrence of myopia regression, which was compared with the observed regression, to assess the difference between the observed and the predicted, yielding sensitivity and 1 − specificity on which analysis of the receiver operating (ROC) curve was based. In principle, the smaller the difference the better internal validity the model was considered to show, by ROC curve (see below). 
Preprocessing and Selecting Highly Correlated Predictors for Myopia Regression
To avoid high correlation between similar clinical correlates—the so-called multicollinearity problem—we checked the predictors with high correlation with each other and preselected one among them before performing the following multivariate analysis. We first dealt with the relationships of the variables with a linear function (e.g., MSE is a linear function of Msph and Mcyl (MSE = Msph + ½ Mcyl) (abbreviations defined in the next paragraph) by retaining only one of the variables. In the case of this example, we retained only MSE and dropped Msph and Mcyl. Similarly, diameter of the optic zone (Size_OZ) was retained and the diameter of the transition zone (Size_TZ) was dropped. 
Pearson correlation analysis was further applied, to assess the extent of the correlations between any two independent variables. Three groups were identified in the light of ranking correlation: correlation coefficients greater than 0.6 (group 1, including manifest axis [Maxis], nomogram-adjusted axis [Nomo_axis] and cycloplegic adjusted axis [Caxis]); correlation coefficients greater than 0.8 (group 2, including MSE, Msph, nomogram-adjusted spherical equivalent [Nomo-SE], nomogram-adjusted sphere [Nomo_sph], cycloplegic spherical equivalent [CSE], cycloplegic spherical equivalent [Csph], depth of optic zone [Depth_OZ] and depth of transition zone [Depth_TZ]), and correlation coefficients equal to 0.9 (group 3, including CCT and RBT). 
According to the principle of reducing superfluous variables due to multicollinearity (high correlation between two variables) and taking into account clinical significance, we selected one variable from each study group: MSE for group 1, Maxis for group 2, and RBT for group 3. These three selected variables, together with the remaining variables—age, sex, preoperative central corneal curvature (Mean K), diameter of optic zone (Size_OZ), flap thickness and diameter, Schirmer test II result, intraocular pressure (IOP), residual bed thickness (RBT), and undercorrection—formed a constellation of predictors for myopia regression included in the following multivariable model. 
Statistical Analysis
Student’s t-test and ANOVA were used to assess the differences of relevant clinical correlates across two or three groups. Differences were considered statistically significant at P < 0.05. To detect multicollinearity, we applied the Pearson product moment correlation analysis to examine the correlation between independent variables and to classify three groups in the light of correlation coefficients. 
To guarantee that there was no multicollinearity problem among the variables in the final model, a multiple linear regression between independent variables and the final manifest refraction was performed to check whether the variance inflation factor (VIF), one of the indicators for detecting multicollinearity, was less than 10. 
Because myopia regression was detected during follow-up visits without knowledge of the exact occurrence time, we adopted an interval-censored model underpinning the concept of the generalized model with complementary log-log transformation of the probability of myopia regression. A detailed description of the model has been provided by Collett. 11  
The model can be expressed as  
\[\mathrm{log}{[}{-}\mathrm{log}(1{-}{\pi}_{ij}){]}{=}{\gamma}_{j}{+}\mathrm{SCORE}_{i}\]
 
\[{\pi}_{ij}{=}1{-}\mathrm{exp}{[}{-}\mathrm{exp}({\gamma}_{j}{+}\mathrm{SCORE}_{i}){]},\]
where π ij is the probability that the ith subject will have myopia regression in the jth period and γ j is the regression coefficient of the jth follow-up period denoted by γ1, γ2, … , γ5, corresponding to five intervals between follow-up visits: 1 week to 1 month, 1 month to 3 months, 3 months to 6 months, 6 months to 9 months, 9 months to 12 months or longer after surgery. Note that since the interval between examinations was not equivalent, the variable “period” was treated as a categorical variable. The risks for the development of myopia regression in each period (γ j ) were also estimated. 
Score i represents the magnitude of the risk of having myopia regression based on clinical correlates, for the ith subject. It can be written as follows  
\[\mathrm{Score}{=}{\alpha}{+}{\Sigma}{\beta}X\]
where α is the intercept of the predictive model, X is a set of predictors, and β is the corresponding regression coefficients that can be trained by the application of the interval-censored model using equations 1 and 2with the trained data set. 
To assess the adequacy of the predictive model, the receiver operating characteristic (ROC) curve was further used by plotting the sensitivity (y-axis) against 1 − specificity (x-axis), both of which were calculated by cross-tabulating the predicted with the observed, using the trained and validated data sets, respectively, given a series of cutoff points based on scores in equation 2 . The larger the area under the ROC, the better the internal validity of the model. Whether the lower limit of the 95% confidence interval (CI) for the area under the ROC curve was larger than 50% was used to assess statistical significance for the adequacy of the predictive model. Note that the calculations were performed for both the trained and validated data sets, but the result from the validated model is more important than that from the trained model, because the former was not involved in training regression coefficients. All statistical analyses were performed using the SAS 9.1 program for Windows. 
Results
Baseline Characteristics
Table 1presents the descriptive results of preoperative characteristics of patients. In the 311 patients, there was a preponderance of women (n = 242; 77.8%), and the mean age was 31.1 ± 6.2 (SD) years (interquartile range: 27, 35 years). The mean manifest spherical equivalent refraction was −7.08 ± 2.49 (SD) D (interquartile range: −8.88, −5.13 D). The mean ablation depth of the optic zone was 99.27 ± 26.95 μm (interquartile range: 82.2, 120.6 μm) with a mean optic zone diameter of 5.81 ± 0.34 mm (interquartile range: 5.6, 6.0 mm). 
Refraction Change after LASIK
Myopia regression after LASIK was observed in 164 (26.7%) of 615 eyes during a 12-month or longer follow-up period. Among the 81 eyes with undercorrection, 47 (58%) eyes were susceptible to regression. The preoperative mean manifest spherical equivalent refraction in undercorrected eyes was −8.88 ± 2.56 D and was −6.81 ± 2.37 D in adequately corrected eyes. 
Greater preoperative myopia led to a higher likelihood of having myopia regression after LASIK. Statistically significant differences were found 1 day and 1, 3, 6, and 12 months after surgery across three groups (low myopia with baseline refraction below −6.0 D, moderate myopia with baseline refraction from −6.0 to −10.0 D, high myopia with baseline refraction greater than −10 D; Fig. 1 ). 
Eyes with smaller optic zones had a higher likelihood of having a greater amount of regression toward myopia after LASIK. Statistically significant differences between eyes with a small optic zone (diameter <6.0 mm) and those with a large optic zone (diameter ≥6.0 mm) were found 1, 6, 9, and 12 months after surgery (Fig. 2)
Time-Varying Predicted Probability for Myopia Regression with Predictors
Table 2shows the results of multivariate analysis in an interval-censored model. Significant factors identified in the multivariate analysis included manifest spherical equivalent (β = −0.3570, P = <0.0001), mean preoperative corneal curvature (β = −0.1807, P = 0.0009), diameter of optic zone (β = −0.6593, P = 0.0043), undercorrection (β = 0.5530, P = 0.0396), and age (β = 0.0297, P = 0.0734). The risks of myopia regression in the five periods (γ1, γ2, … , γ5) are also shown in Table 2
Using our interval-censored model, Table 3shows the predicted probability of myopia regression in each period and cumulative probability up to a certain period using three selected eyes, representing groups with low, moderate, or high myopia. Taking the first case as an example, the score is calculated as:  
\[\mathrm{intercept}{+}{\beta}_{1}\ {\times}\ \mathrm{MSE}{+}{\beta}_{2}{\times}\mathrm{Mean\ K}{+}{\beta}_{3}{\times}\mathrm{size_OZ}{+}{\beta}_{4}\ {\times}\mathrm{undercorrection}{+}{\beta}5{\times}\mathrm{age}{=}5.9856{+}({-}0.3570){\cdot}({-}13.625){+}({-}0.1807){\cdot}44.24{+}({-}0.6593){\cdot}5.5{+}0.5530{\cdot}0{+}0.0297{\cdot}23{=}{-}0.0875,\]
where β12, … , β5 are the coefficients for each period in Table 2 . Calculation of equation 1shows that the predicted probability (π1) of myopia regression in the first period (from 1 week to 1 month after LASIK) is  
\[{\pi}_{1}{=}1{-}\mathrm{exp}{\{}{-}{[}\mathrm{exp}({\gamma}_{1}{+}\mathrm{SCORE}){]}{\}}{=}0.6002,\]
where γ1 is the regression coefficient of the first follow-up period (Table 3) . The predicted probability of other periods (π2, π3, π4, and π5) for the same patient can be calculated in a similar manner. Therefore, the cumulative probability of myopia regression until period 5 is equivalent to  
\[{\pi}(\mathrm{recurrence\ in\ period}\ 1)\ {+}{\pi}(\mathrm{recurrence\ in\ period\ 2\ and\ no\ recurrence\ in\ period\ 1})\ {+}{\ldots}{+}{\pi}(\mathrm{recurrence\ in\ period\ 5},\]
 
\[\mathrm{no\ recurrence\ in\ period\ 1,\ .\ .\ ,\ and\ period\ 4})\]
 
\[{=}{\pi}_{1}{+}{\pi}_{2}(1{-}{\pi}_{1}){+}{\pi}_{3}(1{-}{\pi}_{1})(1{-}{\pi}_{2})\]
 
\[{+}{\pi}_{4}(1{-}{\pi}_{1})(1{-}{\pi}_{2})(1{-}{\pi}_{3}){+}{\pi}_{5}(1{-}{\pi}_{1})(1{-}{\pi}_{2})\]
 
\[{\times}(1{-}{\pi}_{3})(1{-}{\pi}_{4}){=}1{-}{[}(1{-}{\pi}_{1})(1{-}{\pi}_{2}){\times}(1{-}{\pi}_{3})(1{-}{\pi}_{4})(1{-}{\pi}_{5}){]}{=}1{-}{[}(1{-}0.6002)(1{-}0.2975)\]
 
\[{\times}(1{-}0.5092)(1{-}0.3211)(1{-}0.3126){]}{=}0.9357.\]
(Table 3) . The eyes with a higher score had a greater probability of myopia regression. The predicted probabilities of average risk based on the mean of all eyes are also shown in Table 3 . Figure 3shows the cumulative probability up to a a certain period in the three groups. The risk of myopia regression increased rapidly within a month, slowed between 1 and 6 months, and became steady after 6 months, regardless of risk group. The risks of myopia regression up to 6 months after LASIK were 86%, 42%, and 2% for high, moderate, and low risks, respectively. The corresponding risk estimate for the average risk based on all eyes was 21%. 
Model Validation of Predictive Model
The validity of this model was tested by cross-validation using the ROC curve. The area under the curve of the interval-censored model was 82.96% (95% CI: 77.77%–88.16%) in the trained group and 78.43% (95% CI: 73.04%–83.82%) in the validated group (Fig. 4) . The lower limit of the 95% CI of the area under the ROC were significantly higher than 50% in the validated group, which suggests that the predicted validity of our model was good. 
Discussion
Variation of Regression and Retreatment Rate
The complete ascertainment of myopia regression after LASIK, together with the application of an interval-censored model suggests that the risk of myopia regression after LASIK is highly dependent on the definition and inclusion criteria of myopia regression and the reported follow-up time after LASIK. This dependency on inclusion criteria accounts for why myopia regression after LASIK reported in previous studies has ranged from 5.5% to 27.7%. 2 3 4 5 6 Lyle et al. 2 reported 114 eyes re-treated for regression with an overall regression rate of 5.5% in a large case study (114/2065 eyes with preoperative myopia of −6.11 ± 2.35 D). The lower rate may be attributable to the exclusion of the cases of regression that did not undergo retreatment. The definition of regression in Lyle et al. 2 was a residual myopia of greater than −0.5 and a 0.25-D or greater shift toward myopia between follow-ups except in the cases of undercorrection. Chayet et al. 3 found a regression rate of 28% (13/47 eyes were retreated for regression) in patients with mean preoperative myopia of −14.02 ± 5.3 D based on the definition of regression as a 0.25-D or greater myopic shift occurred between follow-up visits. If we take different operative definitions of myopia regression in previous studies 2 3 5 6 14 15 as a reference, our regression rate defined by a residual myopia of −1.0 D or greater and myopic shift by 0.5 D or greater during follow-up visits is 26.7%, close to the 28.01% predicted value reported in Table 3 . The incidence of myopia regression found in our study is probably more comprehensive than that reported in the previous studies, because we analyzed the evolution of all cases undergoing LASIK by one surgeon instead of retreatment cases only. 
The use of an interval-censored model provides a comparison of the predicted time of onset of myopia regression against the reported time frame available in the literature. The myopia regression rate after LASIK has been reported within 6 months after surgery in several studies, 3 8 16 17 18 but only one study included follow-up to 1 to 2 years after surgery. 2 In cases of high preoperative myopia, regression has been observed to peak within 6 months after the initial LASIK surgery and then to stabilize afterward. 3 16 This finding was consistent with our estimated results based on the average-risk group. Cumulative risk for myopia regression increased from 10% within a month to 21% within 6 months after LASIK but then became stable (Table 3) . In eyes with preoperative myopia of −14 D, Chayet et al. 3 reported a regression rate of 27.7% within 3 months after surgery that then stabilized between 3 and 6 months. Their result was higher than the 14% reported in our study (Table 3) , because patients treated with LASIK in their study had a higher grade of myopia than did those in our study (−14 D vs. −7.08 D). Perez-Santonja et al. 16 reported a regression of 0.53 D between 1 and 3 months after surgery in eyes with preoperative mean myopia of −13 D, but no significant regression after 3 months after surgery. Based on Figure 2in Perez-Santonja et al., 16 the regression rate when −1 D was used as a criterion 6 months after surgery was ∼19.56%, which was close to the 21% in our study (Table 3)
Predictors of Myopia Regression
Factors reported to be associated with refractive outcome after LASIK vary from study to study, including preoperative refractive error, 2 3 4 6 7 16 19 20 21 22 corneal curvature, 2 3 14 16 corneal thickness, 3 6 15 20 flap thickness, 23 24 25 ablation depth, 2 optic zone size, 6 12 chronic dry eye, 5 age, 4 surgeon, 4 IOP, 2 16 and humidity. 2 Using an interval-censored model, we not only identified clinical correlates related to myopia regression but also predicted time-varying risk for myopia regression. We considered multicollinearity between variables and chose the appropriate ones, given the orthogonal properties between two of each. Significant predictors identified in multivariate analysis of the interval-censored model included manifest refraction, preoperative keratometry reading, size of optic zone, and undercorrection. Our significant predictors were consistent with those identified in the previous findings. Greater preoperative refraction 2 3 4 6 7 16 19 20 21 22 increased the probability of myopia regression or retreatment; size of optic zone 6 12 16 was associated with refractive outcome after LASIK; flatter preoperative keratometry was associated with a more myopic outcome after LASIK 16 ; and preoperative keratometry <43.5 D predicted greater undercorrection (greater residual myopia shift at 3 months after surgery, according to the definition in Rao et al., which may contain some regressed eyes, according to our definition). 14  
Different from previous studies on retreatment after LASIK, which did not regard initial undercorrection cases as an independent variable, 2 3 26 our study treated undercorrection as one of our independent variables in multivariate analysis instead and found that undercorrection led to a higher probability of myopia regression. However, a high association between undercorrection and preoperative MSE is probable. This probability raises a concern as to whether both, if included in the same model, may cause colinearity. However, as tested by the VIF, a method for detecting colinearity, no apparent colinearity problem was revealed and thus undercorrection and preoperative MSE should be thought of as two independent predictors of myopia regression. 
Usefulness of the Predictive Model
In most previous studies, the assessment of the overall effects of multiple factors on outcome after LASIK have been addressed with only one or two factors, and very few studies have developed a model based on a series of predictors. Hu et al. 4 developed a mathematical predictive formula for retreatment, with multiple variables including age, preoperative cycloplegic sphere, and surgeon. Our study had two characteristics that were different from those in Hu et al. 4 : In their study, the risk of myopia regression was not reported in a longitudinal time frame, and only retreatment cases were included. A logistic regression model that treats retreatment or not as a binary variable (yes/no) predicted an 18% rate of retreatment, which is slightly lower than the predicted 28% rate a year after LASIK— the time when myopia regression reached a plateau. The higher rate in our study compared with that in Hu et al. was partly due to the inclusion of myopia regression in our study rather than only retreatment as in Hu et al. and also was partly due to different types of predictive models in our study (interval-censored model) and their study (logistic regression model). Another disparity is that the “surgeon factor” was taken into account in their study in which patients were treated by two surgeons who had variations in threshold for retreatment, whereas only one ophthalmologist was involved in our study. Our model focused on each eye rather than a subject (person) with cross-validated design by randomly selecting one eye from each subject to build up the predictive model. This method may eliminate the problem of correlation between both eyes of a patient. 
An interval-censored model for predicting the risk of myopia regression with follow-up time after LASIK has never been developed before, to the best of our knowledge. The novel idea of applying an interval-censored model is to predict when and what sort of subjects may be at risk of myopia regression, corresponding to the evolution of refraction after LASIK over time. When the exact time point of the onset of myopia regression is not discernible, analysis with this model to examine occurrence between two visits, given a time frame, will be more accurate in identifying related risk factors and in predicting the probability of myopia regression after LASIK, which will be helpful in counseling patients during follow-up. 
There are two limitations to our study. Since our outcomes were based on changes in refractive error, as measured by the optometer, the result in determining regression could also be influenced by the accommodation of the subject. Based on the biomechanical factors observed to accompany myopia regression, 3 15 16 the definition of myopia regression would have been more valid if the measurement of the postoperative Mean K and CCT had been collected together. Second, one of our variables in the model is undercorrection, which is identified after surgery, and thus may limit the model’s application in patient counseling before LASIK but is still useful for predicting the evolution of refraction during follow-up after LASIK. 
In conclusion, we propose a useful time-varying predictive model to predict myopia regression after LASIK for myopia and myopic astigmatism. Information obtained from the predictive model enables a clinician to identify who has a high risk of myopia regression and to estimate the probability within the follow-up time frame more precisely. 
 
Table 1.
 
Basic Characteristics of Preoperative Variables
Table 1.
 
Basic Characteristics of Preoperative Variables
Variable Mean SD Lower Quartile Upper Quartile
Age (y) 31.1 6.2 27 35
Manifest refraction
 Msph (D) −6.61 2.40 −8.50 −4.75
 Mcyl (D) −0.95 0.76 −1.25 −0.50
 Maxis (deg) 104.88 79.99 15.00 170.00
 MSE (D) −7.08 2.49 −8.88 −5.13
Cycloplegic refraction
 Csph (D) −6.39 2.84 −8.00 −4.50
 Ccyl (D) −1.13 3.11 −1.25 −0.50
 Caxis (deg) 96.95 74.02 11.00 170.00
 CSE (D) −6.95 3.24 −8.75 −4.88
Nomogram-adjusted refraction
 Nomo_sph (D) −6.73 2.39 −8.50 −5.00
 Nomo_cyl (D) −0.93 0.83 −1.25 −0.50
 Nomo_axis (deg) 102.61 182.04 5.00 170.00
 NomoSE (D) −7.19 2.50 −9.00 −5.25
MeanK (D) 43.67 1.95 42.83 44.63
Flap_thick (μm) 131.71 6.96 130.00 130.00
Flap_diameter (mm) 8.97 0.17 9.00 9.00
Size_OZ (mm) 5.81 0.34 5.60 6.00
Size_TZ (mm) 6.83 0.37 6.70 7.00
Depth_OZ (μm) 99.27 26.95 82.20 120.60
Depth_TZ (μm) 17.15 6.69 13.50 20.40
Central corneal thickness (CCT) (μm) 546.75 55.81 529.00 571.00
RBT (μm) 315.76 60.93 287.00 344.20
IOP (mm Hg) 15.53 6 13.00 17.00
Schirmer test (mm) 20.35 10.31 14.00 26.00
Figure 1.
 
The comparisons of mean spherical equivalent across three groups stratified by the degree of myopia over, with significance (*P < 0.05) assessed by analysis of variance (ANOVA).
Figure 1.
 
The comparisons of mean spherical equivalent across three groups stratified by the degree of myopia over, with significance (*P < 0.05) assessed by analysis of variance (ANOVA).
Figure 2.
 
The comparisons of mean spherical equivalent between the two groups stratified by diameter of optic zone, with significant result (*P < 0.05) assessed by independent t-test.
Figure 2.
 
The comparisons of mean spherical equivalent between the two groups stratified by diameter of optic zone, with significant result (*P < 0.05) assessed by independent t-test.
Table 2.
 
Estimated Results of Multivariate Analysis in an Interval-Censored Model
Table 2.
 
Estimated Results of Multivariate Analysis in an Interval-Censored Model
Variable Regression Coefficient SE P (Wald Test)
Intercept 5.9856 2.8238 0.0340
Period (γ) 0.0219
 Period 1 0
 Period 2 −0.9541 0.3496
 Period 3 −0.2530 0.3280
 Period 4 −0.8615 0.5998
 Period 5 −0.8945 0.6016
MSE (β1) −0.3570 0.0505 <0.0001
Mean K (β2) −0.1807 0.0544 0.0009
Size_OZ (β3) −0.6593 0.2311 0.0043
Undercorrection (β4) 0.5530 0.2688 0.0396
Age (β5) 0.0297 0.0166 0.0734
Table 3.
 
Variables and Time-Varying Predicted Probabilities in Three Selected Cases in an Interval-Censored Model
Table 3.
 
Variables and Time-Varying Predicted Probabilities in Three Selected Cases in an Interval-Censored Model
Eye MSE MEANK Size_OZ Under correction* Age Score Period, † Probability, ‡ Cumulative Probability, §
1 −13.625 44.24 5.5 0 23 −0.0875 1 0.6002 0.6002
2 0.2975 0.7191
3 0.5092 0.8621
4 0.3211 0.9064
5 0.3126 0.9357
2 −8.25 42.75 5.5 0 35 −1.3807 1 0.2223 0.2223
2 0.0923 0.2941
3 0.1773 0.4193
4 0.1008 0.4778
5 0.0977 0.5288
3 −2.75 44.69 6.5 0 23 −4.7105 1 0.0091 0.0091
2 0.0035 0.0124
3 0.0071 0.0194
4 0.0039 0.0232
5 0.0017 0.0248
All eyes (mean) −7.08 43.67 5.81 0.13171 31 −2.2094 1 0.1040 0.1040
2 0.0414 0.1411
3 0.0817 0.2113
4 0.0453 0.2470
5 0.0439 0.2801
Figure 3.
 
Cumulative risk of myopia regression for the overall mean of all eyes and three selected cases in an interval-censored model.
Figure 3.
 
Cumulative risk of myopia regression for the overall mean of all eyes and three selected cases in an interval-censored model.
Figure 4.
 
Area under ROC curves and 95% CI, based on the trained and the validation data sets, in an interval-censored model.
Figure 4.
 
Area under ROC curves and 95% CI, based on the trained and the validation data sets, in an interval-censored model.
LinLL, ShihYF, HsiaoCK, ChenCJ. Prevalence of myopia in Taiwanese schoolchildren: 1983 to 2000. Ann Acad Med Singapore. 2004;33:27–33.
LyleWA, JinGJ. Retreatment after initial laser in situ keratomileusis. J Cataract Refract Surg. 2000;26:650–659. [CrossRef] [PubMed]
ChayetAS, AssilKK, MontesM, et al. Regression and its mechanisms after laser in situ keratomileusis in moderate and high myopia. Ophthalmology. 1998;105:1194–1199. [CrossRef] [PubMed]
HuDJ, FederRS, BastiS, et al. Predictive formula for calculating the probability of LASIK enhancement. J Cataract Refract Surg. 2004;30:363–368. [CrossRef] [PubMed]
AlbietzJM, LentonLM, McLennanSG. Chronic dry eye and regression after laser in situ keratomileusis for myopia. J Cataract Refract Surg. 2004;30:675–684. [CrossRef] [PubMed]
LianJ, ZhangQ, YeW, ZhouD, WangK. An analysis of regression after laser in situ keratomileusis for treatment of myopia. Zhonghua Yan Ke Za Zhi. 2002;38:363–366. [PubMed]
CondonPI, MulhernM, FulcherT, Foley-NolanA, O’KeefeM. Laser intrastromal keratomileusis for high myopia and myopic astigmatism. Br J Ophthalmol. 1997;81:199–206. [CrossRef] [PubMed]
PatelNP, ClinchTE, WeisJR, AhnC, LunderganMK, HeidenreichK. Comparison of visual results in initial and re-treatment laser in situ keratomileusis procedures for myopia and astigmatism. Am J Ophthalmol. 2000;130:1–11. [CrossRef] [PubMed]
OzdamarA, ArasC, BahceciogluH, SenerB. Secondary laser in situ keratomileusis 1 year after primary LASIK for high myopia. J Cataract Refract Surg. 1999;25:383–388. [CrossRef] [PubMed]
GimbelHV, PennoEE, van WestenbruggeJA, FerensowiczM, FurlongMT. Incidence and management of intraoperative and early postoperative complications in 1000 consecutive laser in situ keratomileusis cases. Ophthalmology. 1998;105:1839–1847. [CrossRef] [PubMed]
CollettD. Modelling Survival Data in Medical Research. 2003; 2nd ed. 237–251.Taylor & Francis Bristol, UK.
DitzenK, HandzelA, PiegerS. Laser in situ keratomileusis nomogram development. J Refract Surg. 1999;15:S197–S201. [PubMed]
FelthamMH, WolfeRJ. Some variables to consider to avoid the need for LASIK surgical enhancements. Clin Exp Optom. 2000;83:76–81. [CrossRef] [PubMed]
RaoSK, ChengAC, FanDS, LeungAT, LamDS. Effect of preoperative keratometry on refractive outcomes after laser in situ keratomileusis. J Cataract Refract Surg. 2001;27:297–302. [CrossRef] [PubMed]
PanQ, GuYS, WangJ, et al. Differences between regressive eyes and non-regressive eyes after LASIK for myopia in the time course of corneal changes assessed with the Orbscan. Ophthalmologica. 2004;218:96–101. [CrossRef] [PubMed]
Perez-SantonjaJJ, BellotJ, ClaramonteP, IsmailMM, AlioJL. Laser in situ keratomileusis to correct high myopia. J Cataract Refract Surg. 1997;23:372–385. [CrossRef] [PubMed]
PallikarisIG, SiganosDS. Excimer laser in situ keratomileusis and photorefractive keratectomy for correction of high myopia. J Refract Corneal Surg. 1994;10:498–510. [PubMed]
ZaldivarR, DavidorfJM, ShultzMC, OscherowS. Laser in situ keratomileusis for low myopia and astigmatism with a scanning spot excimer laser. J Refract Surg. 1997;13:614–619. [PubMed]
PallikarisIG, SiganosDS. Laser in situ keratomileusis to treat myopia: early experience. J Cataract Refract Surg. 1997;23:39–49. [CrossRef] [PubMed]
MagallanesR, ShahS, ZadokD, et al. Stability after laser in situ keratomileusis in moderately and extremely myopic eyes. J Cataract Refract Surg. 2001;27:1007–1012. [CrossRef] [PubMed]
ZhengL, ZhuQ, SunY, YeM. The complications of excimer laser in-situ keratomileusis for myopia. Zhonghua Yan Ke Za Zhi. 1999;35:363–365. [PubMed]
FarahSG, AzarDT, GurdalC, WongJ. Laser in situ keratomileusis: literature review of a developing technique. J Cataract Refract Surg. 1998;24:989–1006. [CrossRef] [PubMed]
No authors listed. Induced astigmatism following laser in situ keratomileusis for myopia with a free cap. J Refract Surg. 2000;16:375–379. [PubMed]
EleftheriadisH, PrandiB, Diaz-RatoA, MorcilloM, SabaterJB. The effect of flap thickness on the visual and refractive outcome of myopic laser in situ keratomileusis. Eye. 2005;19:1290–1296. [CrossRef] [PubMed]
WilsonSE, MohanRR, HongJW, LeeJS, ChoiR, MohanRR. The wound healing response after laser in situ keratomileusis and photorefractive keratectomy: elusive control of biological variability and effect on custom laser vision correction. Arch Ophthalmol. 2001;119:889–896. [CrossRef] [PubMed]
Perez-SantonjaJJ, AyalaMJ, SaklaHF, Ruiz-MorenoJM, AlioJL. Retreatment after laser in situ keratomileusis. Ophthalmology. 1999;106:21–28. [CrossRef] [PubMed]
Figure 1.
 
The comparisons of mean spherical equivalent across three groups stratified by the degree of myopia over, with significance (*P < 0.05) assessed by analysis of variance (ANOVA).
Figure 1.
 
The comparisons of mean spherical equivalent across three groups stratified by the degree of myopia over, with significance (*P < 0.05) assessed by analysis of variance (ANOVA).
Figure 2.
 
The comparisons of mean spherical equivalent between the two groups stratified by diameter of optic zone, with significant result (*P < 0.05) assessed by independent t-test.
Figure 2.
 
The comparisons of mean spherical equivalent between the two groups stratified by diameter of optic zone, with significant result (*P < 0.05) assessed by independent t-test.
Figure 3.
 
Cumulative risk of myopia regression for the overall mean of all eyes and three selected cases in an interval-censored model.
Figure 3.
 
Cumulative risk of myopia regression for the overall mean of all eyes and three selected cases in an interval-censored model.
Figure 4.
 
Area under ROC curves and 95% CI, based on the trained and the validation data sets, in an interval-censored model.
Figure 4.
 
Area under ROC curves and 95% CI, based on the trained and the validation data sets, in an interval-censored model.
Table 1.
 
Basic Characteristics of Preoperative Variables
Table 1.
 
Basic Characteristics of Preoperative Variables
Variable Mean SD Lower Quartile Upper Quartile
Age (y) 31.1 6.2 27 35
Manifest refraction
 Msph (D) −6.61 2.40 −8.50 −4.75
 Mcyl (D) −0.95 0.76 −1.25 −0.50
 Maxis (deg) 104.88 79.99 15.00 170.00
 MSE (D) −7.08 2.49 −8.88 −5.13
Cycloplegic refraction
 Csph (D) −6.39 2.84 −8.00 −4.50
 Ccyl (D) −1.13 3.11 −1.25 −0.50
 Caxis (deg) 96.95 74.02 11.00 170.00
 CSE (D) −6.95 3.24 −8.75 −4.88
Nomogram-adjusted refraction
 Nomo_sph (D) −6.73 2.39 −8.50 −5.00
 Nomo_cyl (D) −0.93 0.83 −1.25 −0.50
 Nomo_axis (deg) 102.61 182.04 5.00 170.00
 NomoSE (D) −7.19 2.50 −9.00 −5.25
MeanK (D) 43.67 1.95 42.83 44.63
Flap_thick (μm) 131.71 6.96 130.00 130.00
Flap_diameter (mm) 8.97 0.17 9.00 9.00
Size_OZ (mm) 5.81 0.34 5.60 6.00
Size_TZ (mm) 6.83 0.37 6.70 7.00
Depth_OZ (μm) 99.27 26.95 82.20 120.60
Depth_TZ (μm) 17.15 6.69 13.50 20.40
Central corneal thickness (CCT) (μm) 546.75 55.81 529.00 571.00
RBT (μm) 315.76 60.93 287.00 344.20
IOP (mm Hg) 15.53 6 13.00 17.00
Schirmer test (mm) 20.35 10.31 14.00 26.00
Table 2.
 
Estimated Results of Multivariate Analysis in an Interval-Censored Model
Table 2.
 
Estimated Results of Multivariate Analysis in an Interval-Censored Model
Variable Regression Coefficient SE P (Wald Test)
Intercept 5.9856 2.8238 0.0340
Period (γ) 0.0219
 Period 1 0
 Period 2 −0.9541 0.3496
 Period 3 −0.2530 0.3280
 Period 4 −0.8615 0.5998
 Period 5 −0.8945 0.6016
MSE (β1) −0.3570 0.0505 <0.0001
Mean K (β2) −0.1807 0.0544 0.0009
Size_OZ (β3) −0.6593 0.2311 0.0043
Undercorrection (β4) 0.5530 0.2688 0.0396
Age (β5) 0.0297 0.0166 0.0734
Table 3.
 
Variables and Time-Varying Predicted Probabilities in Three Selected Cases in an Interval-Censored Model
Table 3.
 
Variables and Time-Varying Predicted Probabilities in Three Selected Cases in an Interval-Censored Model
Eye MSE MEANK Size_OZ Under correction* Age Score Period, † Probability, ‡ Cumulative Probability, §
1 −13.625 44.24 5.5 0 23 −0.0875 1 0.6002 0.6002
2 0.2975 0.7191
3 0.5092 0.8621
4 0.3211 0.9064
5 0.3126 0.9357
2 −8.25 42.75 5.5 0 35 −1.3807 1 0.2223 0.2223
2 0.0923 0.2941
3 0.1773 0.4193
4 0.1008 0.4778
5 0.0977 0.5288
3 −2.75 44.69 6.5 0 23 −4.7105 1 0.0091 0.0091
2 0.0035 0.0124
3 0.0071 0.0194
4 0.0039 0.0232
5 0.0017 0.0248
All eyes (mean) −7.08 43.67 5.81 0.13171 31 −2.2094 1 0.1040 0.1040
2 0.0414 0.1411
3 0.0817 0.2113
4 0.0453 0.2470
5 0.0439 0.2801
×
×

This PDF is available to Subscribers Only

Sign in or purchase a subscription to access this content. ×

You must be signed into an individual account to use this feature.

×