All data analysis and statistics were performed in MATLAB R2016 (MathWorks, Natick, MA, USA). First, we removed outliers using methods described by Iglewicz and Hoaglin.
24 In brief, we calculated a
z-score based on the median and the median absolute deviation, with values greater than 3.5 indicating an outlier. After excluding poor cannulations and outliers, we obtained the following numbers of animals per group: young sham (n = 5), young OVX (n = 6), middle-aged sham (n = 5), and middle-aged OVX (n = 5). Similar to an earlier approach, a log transformation was applied to outflow facility and ocular compliance, and all datasets were then assessed for normality using an Anderson–Darling test. All fitting parameters (
Cr, β, ϕ
r, and γ)
are detailed in the
Table and were analyzed using a two-way ANOVA, with age (young vs. middle-aged) and menopausal status (sham vs. OVX) as the two main factors. Based on established statistical methodology,
25,26 when the interaction effect between age and menopausal status was not significant in the two-way ANOVA, we report only the main effects of age and of menopause on our outcome measures. Thus, we compare only young versus middle-aged or sham versus OVX. A difference was considered significant at
P < 0.05. Comparisons between groups are represented as the percent difference and 95% confidence interval (CI). Data are presented as box-and-whisker plots (lines represent the median, and bars are the minimum to maximum).