As all treatment arms were designed to deliver a similar average aerobic intensity, we planned to pool a priori all completers throughout the study to analyze the effects of aerobic exercise with higher statistical power. We used Wilcoxon Signed Ranks tests or paired t-tests to compare final vs. baseline outcomes. When a significant change in outcomes was observed, we used regression models to assess and adjust for the effect of different settings and training methods, calendar year, and change in levodopa equivalent. We also used Pearson correlations and regression models to quantify associations of changes in outcomes with changes in aerobic fitness.