Statistical analysis was performed using PASW 18.0 for Windows (IBM) and the R System for Statistical Computing, version 2.15.2 (provided in the public domain,
http://www.r-project.org). Results are demonstrated as mean and SD or frequency and percentage if appropriate for all subgroups. Comparisons between groups and variables were performed using
t-test or Mann-Whitney
U test for quantitative variables and the χ
2 test for qualitative variables. Linear mixed model analysis
13 was used to analyze the rate of change in mean RNFL thickness during follow-up. The overall time trend and the mean differences among the groups (healthy, glaucoma without progression, and glaucoma with progression) were modeled as fixed effects. Individual distributions of the intercepts of the patients were modeled as normally distributed random effects. In addition, a fixed interaction term between the time trend and the groups was included in the linear mixed model. This strategy resulted in individual trend lines for each of the three groups (
Fig. 2). By definition of the model, group-specific rates of change in mean RNFL thickness could be quantified by the estimated slopes of the trend lines. In the following, the slopes are denoted by
βi,
i = 1, 2, 3. To investigate whether the rates of change in mean RNFL thickness differed among the groups, differences between the group-specific rates were tested against zero. This was done by using general linear hypothesis tests.
14 The method by Bretz et al.
14 automatically adjusted
P values for multiple comparisons. A significance level of α = 0.05 was used for all statistical hypothesis tests.