A total of 1097 P dc data points were collected. On average, each subject provided two sets (SD = 2; range, 1–16) of P dc measurements of both eyes, where 10% of subjects provided five sets or more; 21% of 2 to 4 sets, and the remaining 69% only one set of data. Fifty-nine observations with central corneal staining greater than 1 were excluded. Based on an a priori criterion, 11 P dc data below the lower limit and 30 beyond the upper limit were excluded from analysis. The final data set pooled into statistical analysis contained 997 observations from 276 subjects (42% Asians, 58% non-Asians; 49% females, 51% males). The P dc measurement time from awakening ranged from 3 minutes to 17 hours.
The descriptive statistics of hourly
P dc are summarized in
Table 1. The highest mean
P dc was 0.136 [0.123, 0.150] nm/sec, observed within 1 hour of awakening. The hourly mean
P dc with 95% confidence interval (CI) is plotted in
Figure 1. Results from post hoc
t-test with Tukey adjustments showed that mean
P dc at 1 hour was significantly higher than all subsequent hourly means (
P < 0.05). There was no significant
P dc change after 2 hours (
P > 0.05) from awakening. A scatterplot of
P dc with a nonparametric LOESS curve in
Figure 2 characterized the
P dc dynamics: a fast drop of
P dc from awakening up till nearly 2 hours, which can be regarded as a recovery period, followed by a plateau.
The results from mixed-effect modeling are presented in
Table 2.
P dc values were transformed to natural logarithm to approximate normal distribution. Model 1 investigated the
P dc recovery period from awakening to 2 hours. Model 2 examined the
P dc plateau period. Covariate candidates included age, sex, ethnicity, awakening time, and a quadratic term of awakening time. Other covariates considered the interaction variables of age and sex, age and ethnicity, sex and ethnicity, awakening time and age, awakening time and sex, and awakening time and ethnicity. Due to the long duration of the study, a variable “Days
x machine,” representing the number of days from the first use of each machine, was added in the models to adjust for the effect of potential instrument and observer bias.
Model 1 indicated a significant decline of
P dc with length of awake time (
P = 0.000). Model 2 showed no significant evidence of association between
P dc and awake time after 2 hours, confirming the observation of plateau. Both models showed that
P dc in Asians was significantly higher than in non-Asians and the difference was relatively constant (
Fig. 3).
P dc in males was higher than in females during the first 2 hours with marginal significance (
P = 0.098). Such sex difference was not observed after 2 hours. Age was positively associated with
P dc during the whole 17-hour period, indicating a higher
P dc (i.e., weaker epithelial barrier function) with age in the study age group of 18 to 38 years. None of the interaction terms between ethnicity, sex, age, and
T A showed statistical significance. Both models suggested a need to adjust for the duration of instrument usage of one of the fluorotrons (
P = 0.003 in Model 1,
P = 0.016 in Model 2).