The mean IOPs before and after CCT correction were 15.557 ± 3.740 mm Hg and 15.708 ± 3.690 mm Hg, respectively, when measured with GAT, and 15.564 ± 3.520 mm Hg and 15.715 ± 3.350 mm Hg, respectively, when measured with RT. For the Corvis ST, the mean corrected IOP was 16.105 ± 4.060 mm Hg (
Table 1). The mean CCTs were 556 ± 38 μm for the ultrasound pachymetry, 543 ± 52 μm for the Corvis ST, and 547 ± 54 μm for the Pentacam, and there was no significant difference among the CCT measurements (
P > 0.05, Wilcoxon paired test).
All values in the groups and subgroups were not normally distributed; thus, for statistical analysis, nonparametric tests were used. The exception was the linear Pearson coefficient, which was selected for cluster analysis because we built our model on linear dependency rather than on outlier values.
In a single IOP values analysis, the highest IOP was typically measured with the Corvis ST (48.5% of all measurements), and the lowest IOP was measured with the Icare RT (44.3% of all measurements) (
Fig. 1). In the Wilcoxon paired test analysis, no significant differences were observed between the IOP values measured with GAT and Icare RT (
P > 0.05). There was a significant difference between Corvis and GAT (Wilcoxon test,
P < 0.001) and between Corvis and RT (Wilcoxon test,
P < 0.03). Using the Spearman correlation test, the relation between the IOP and the CCT was checked; in addition, the CCT was corrected for GAT and RT. In this analysis, the CCT significantly affected the IOP values before correction (for GAT,
r = 0.319 and
P < 0.0001, and for RT,
r = 0.349 and
P < 0.0001). After CCT correction, the values showed negative, weak, but still significant correlations (for GAT,
r = −0.193 and
P = 0.03, and for RT,
r = −0.186 and
P = 0.05) (
Table 2,
Fig. 2). Bland-Altman plots were constructed to show agreement between GAT and RT, GAT and Corvis, and RT and Corvis before and after CCT correction. For the GAT-RT comparison, the mean of the difference was 0.0 and the 95% limit of agreement was +4.5 mm Hg before the CCT; these values were unchanged after CCT correction. For the GAT-Corvis comparison, the mean of the difference was −0.5 before CCT with 95% limit of agreement of +4.4 mm Hg; these values were −0.5 and +5.5 mm Hg, respectively, after CCT correction. For the RT-Corvis comparison, the mean of the difference was −0.5 and the 95% limit of agreement was +5.4 mm Hg before CCT correction; and after CCT correction, −0.5 and +6.3 mm Hg, respectively (
Fig. 3).
The corneal functional properties were analyzed using UHS ST. For each cornea, the following parameters were calculated: deformation maximum amplitude (Def. Amp. Max.), the maximum amplitude of deflection at the corneal apex (highest concavity) (mm); first A time, the time from the start until the first applanation (ms); first A length, the cord length of the first applanation (mm); A1 velocity, corneal speed during the first applanation moment (m/s); second A time, the time from the start until the second applanation (ms); second A length, the cord length of the second applanation (mm); A2 velocity, corneal speed during the second applanation moment (m/s); and HC time, the time from starting until the highest concavity is reached (ms) (
Fig. 4). These parameters were compared among normal, keratoconic, glaucomatous, and swollen corneas (CCT between 570 and 700 μm). There was a significant difference among the A1 and A2 times and velocities and the Def. Amp. Max. when comparing glaucomatous versus bullous corneas, keratoconic versus glaucomatous corneas, healthy versus glaucomatous corneas, keratoconic versus bullous corneas, and keratoconic versus healthy corneas (
Table 3).
An analysis of the correlations showed strong, significant dependence among the GAT, RT, and Corvis IOP values before the CCT correction and the corneal properties expressed as the A1 and A2 times, A1 and A2 velocities, and Def. Amp. Max. After CCT correction, there was medium or weak dependence between these values. The correlation coefficients are presented in
Table 4. In the cluster analysis using the Pearson correlation coefficient, the variables measured with Corvis were grouped in four supraclusters, where A1 time and A2 velocity were in one supracluster together with IOP (the strongest relation). The shortest distance from the IOP cluster was observed in the case of A1 time, A2 velocity, peak distance, and radius (
Fig. 5).
Finally, to evaluate the agreement between the two new IOP measurement methods (Corvis and RT), the three sites' difference was compared between each value. The mean difference in the IOP values was 0.01 mm Hg between the RT and the GAT and 0.44 mm Hg between the Corvis and the Icare. When analyzing individual patients' measurements, the difference was more than 3 mm Hg in 33 cases (17.1%) when comparing the GAT and the Icare, and in 45 cases (23.4%) when comparing the GAT and the Corvis (
Fig. 6).