A Compact Stance-Indexed Anterior-Posterior COP Representation for Parkinson's Disease Classification from Plantar VGRF
Organizations: Department of Computer Science and Engineering, University of Rajshahi, Rajshahi, Bangladesh
Abstract
Parkinson's disease alters gait and bilateral coordination, but machine-learning performance also depends on how continuous gait signals are represented. This study investigates whether preserving anterior-posterior center-of-pressure (AP-COP) information at fixed locations across normalized stance provides a compact and informative representation of plantar-force gait signals. Bilateral vertical ground reaction force recordings from 165 participants in the Gait in Parkinson's Disease Database were evaluated using repeated fully nested participant-level cross-validation. We propose AP-COP10, comprising AP-COP position and bilateral asymmetry across five stance windows. AP-COP10 achieved an AUC of 0.894 and outperformed three harmonized literature-derived COP representations under the same evaluation pipeline. The complementary 25 non-AP-COP descriptors alone achieved an AUC of 0.856, while the complete 35-feature representation achieved 0.908. Removing AP-COP10 from the complete representation produced a statistically supported loss in discrimination, whereas adding the complementary descriptors to AP-COP10 yielded only a small, unsupported improvement. Feature competition indicated that the most informative stance-indexed descriptors were concentrated in early and early-mid stance, while source-study holdout and sensor-perturbation analyses supported the robustness of the representation. These findings indicate that stance-indexed AP-COP retains discriminative information that is not readily recovered by broader engineered gait descriptors, supporting compact and interpretable representations for machine-learning analysis of pathological gait.
Figures & tables
| Study | Group | Age (y) | Female | Weight (kg) | Speed (m/s) | |
|---|---|---|---|---|---|---|
| Ga | Control | 18 | 69 [67–77] | 44% | 70 [64–82] | 1.2 [1.1–1.3] |
| Ga | PD | 29 | 71 [68–78] | 31% | 74 [62–84] | 1.0 [0.8–1.1] |
| Ju | Control | 25 | 65 [60–68] | 52% | 72 [60–80] | 1.2 [1.1–1.3] |
| Ju | PD | 29 | 68 [64–74] | 45% | 68 [65–80] | 1.1 [0.9–1.2] |
| Si | Control | 29 | 57 [53–62] | 38% | 74 [67–83] | 1.3 [1.1–1.4] |
| Si | PD | 35 | 63 [56–68] | 37% | 74 [66–80] | 1.1 [1.0–1.2] |
| Family | Information retained | |
|---|---|---|
| Temporal rhythm | 6 | Step/stride timing, stance and swing variability, step-to-step irregularity. |
| Dynamic support | 4 | Double-support organization and stance-related support structure. |
| Bilateral asymmetry | 4 | Normalized left–right differences in temporal and loading summaries. |
| Broader COP rollover | 6 | AP range, path length, median speed, speed variability, monotonicity, early slope. |
| Vertical loading | 5 | Normalized peak force, impulse, loading and unloading behavior, waveform variability. |
| Stance-indexed AP-COP | 10 | Bilateral AP-COP position and asymmetry at five fixed normalized-stance windows. |
| Representation | Candidates | AUC (95% CI) | BAcc | Sens. | Spec. | Brier | AUC (AP-COP10 comparator) (95% CI) | |
|---|---|---|---|---|---|---|---|---|
| Complete engineered multidomain (Full35) | 35 | 0.908 (0.862–0.948) | 0.817 | 0.828 | 0.806 | 0.125 | ( to 0.015) | 0.361 S |
| Proposed stance-indexed AP-COP (AP-COP10) | 10 | 0.894 (0.843–0.937) | 0.828 | 0.796 | 0.861 | 0.135 | – | – |
| Complementary non-AP-COP engineered (Non-AP-COP25) | 25 | 0.856 (0.796–0.910) | 0.751 | 0.753 | 0.750 | 0.154 | 0.038 ( to 0.090) | 0.134 S |
| Alam-derived COP4 | 4 | 0.809 (0.742–0.871) | 0.705 | 0.688 | 0.722 | 0.180 | 0.085 (0.037–0.135) | 0.0012 H |
| Jin-derived COP7 | 7 | 0.784 (0.714–0.852) | 0.724 | 0.753 | 0.694 | 0.190 | 0.110 (0.050–0.171) | 0.0012 H |
| Tong-derived COP18 | 18 | 0.746 (0.665–0.821) | 0.700 | 0.677 | 0.722 | 0.206 | 0.148 (0.084–0.213) | 0.0012 H |
| Model | Ga held out | Ju held out | Si held out | Mean |
|---|---|---|---|---|
| Proposed AP-COP10 | 0.837 | 0.937 | 0.946 | 0.907 |
| Jin gait-only a (SVM-RBF) | 0.845 | 0.500 | 0.833 | 0.726 |
| Descriptor | Family | Retention |
|---|---|---|
| Asymmetry 5–15% | AP-COP | 0.99 |
| Position 25–35% | AP-COP | 0.99 |
| Position 5–15% | AP-COP | 0.90 |
| Asymmetry 25–35% | AP-COP | 0.40 |
| Position 65–75% | AP-COP | 0.23 |
| Position 45–55% | AP-COP | 0.15 |
| ID | Feature identifier | Window | Role |
|---|---|---|---|
| AP1 | cop_landing_position | 5–15% | Bilateral AP-COP position |
| AP2 | cop_pos_25_35 | 25–35% | Bilateral AP-COP position |
| AP3 | cop_midstance_position | 45–55% | Bilateral AP-COP position |
| AP4 | cop_pos_65_75 | 65–75% | Bilateral AP-COP position |
| AP5 | cop_terminal_position | 85–95% | Bilateral AP-COP position |
| AP6 | cop_landing_asymmetry | 5–15% | Bilateral AP-COP asymmetry |
| ID | Feature identifier | Family | Information represented |
|---|---|---|---|
| Temporal/rhythm | |||
| 1 | step_interval_median | Temporal/rhythm | Median alternating-foot step interval. |
| 2 | step_interval_cv | Temporal/rhythm | Step-interval variability. |
| 3 | step_interval_mad_ratio | Temporal/rhythm | Step-interval MAD relative to central scale. |
| 4 | stride_interval_cv | Temporal/rhythm | Same-foot stride-interval variability. |
| 5 | stance_duration_cv | Temporal/rhythm | Stance-duration variability. |
| Set | Reproduced descriptors and harmonization decisions | |
|---|---|---|
| Alam4 | 4 | Mean/SD ML-COP and mean/SD AP-COP from the left-foot reference. A bilateral eight-candidate extension was additionally evaluated as a sensitivity analysis. |
| Jin7 | 7 | Left/right AP-COP SD, left/right ML-COP SD, left/right 2-D COP path length, and bilateral COP-path symmetry. Forefoot-to-rearfoot loading was excluded because the controlled comparison was restricted to COP-specific descriptors; Jin7 is therefore not the complete source-paper model. |
| Tong18 | 18 | Nine descriptor types per foot: RMS ML-COP, RMS AP-COP, RMS total COP magnitude, RMS 2-D COP velocity, RMS acceleration, RMS jerk, COP path efficiency, ML-COP sample entropy, and AP-COP sample entropy. CSIP/CISP were excluded because their reconstruction was under-specified; source-specific persistent-homology processing, oversampling, and SVM classification were not reproduced. Sample entropy used and SD. |
| Family | Mean selected/fold | Folds represented |
|---|---|---|
| Stance-indexed AP-COP | 3.73 | 100% |
| Dynamic support | 0.93 | 93% |
| Temporal/rhythm | 0.72 | 72% |
| Broader COP rollover | 0.39 | 37% |
| Vertical loading | 0.09 | 9% |
| Other bilateral asymmetry | 0.06 | 6% |
| Feature | Family | Retention |
|---|---|---|
| cop_landing_asymmetry | AP-COP asymmetry | 0.99 |
| cop_pos_25_35 | AP-COP position | 0.99 |
| cop_landing_position | AP-COP position | 0.90 |
| cop_asym_25_35 | AP-COP asymmetry | 0.40 |
| cop_pos_65_75 | AP-COP position | 0.23 |
| cop_midstance_position | AP-COP position | 0.15 |
| Representation | Candidates | AUC (95% CI) | BAcc | Sens. | Spec. |
|---|---|---|---|---|---|
| Full35 | 35 | 0.908 (0.862–0.948) | 0.817 | 0.828 | 0.806 |
| AP-COP10 | 10 | 0.894 (0.843–0.937) | 0.828 | 0.796 | 0.861 |
| Non-AP-COP25 | 25 | 0.856 (0.796–0.910) | 0.751 | 0.753 | 0.750 |
| Feature | Outer-fold retention |
|---|---|
| Step-interval MAD ratio | 0.98 |
| Double-support proportion | 0.96 |
| COP AP range | 0.94 |
| COP early slope | 0.81 |
| Peak-force asymmetry | 0.72 |
| Normalized loading rate | 0.66 |
| Analysis | Setting | AUC | 95% CI |
|---|---|---|---|
| Budget | 1 | 0.802 | 0.734–0.866 |
| Budget | 2 | 0.865 | 0.810–0.914 |
| Budget | 3 | 0.885 | 0.835–0.931 |
| Budget | 4 | 0.901 | 0.855–0.941 |
| Budget | 5 | 0.906 | 0.861–0.947 |
| Budget | 6 | 0.908 | 0.862–0.947 |
| Classifier | AUC | BAcc | Sens. | Spec. |
|---|---|---|---|---|
| Histogram gradient boosting | 0.913 | 0.851 | 0.882 | 0.819 |
| Logistic regression | 0.908 | 0.817 | 0.828 | 0.806 |
| RBF-SVM | 0.899 | 0.803 | 0.828 | 0.778 |
| Random forest | 0.898 | 0.817 | 0.871 | 0.764 |
| Condition | AUC | BAcc | AUC loss |
|---|---|---|---|
| Clean | 0.908 | 0.817 | ref. |
| 5% gain SD | 0.907 | 0.812 | 0.001 |
| 10% gain SD | 0.902 | 0.814 | 0.006 |
| 20% gain SD | 0.881 | 0.789 | 0.027 |
| One-sensor dropout | 0.858 | 0.787 | 0.050 |
| Two-sensor dropout | 0.802 | 0.738 | 0.106 |
| Analysis | Representation | AUC | Relevant paired comparison |
|---|---|---|---|
| Asymmetry-definition sensitivity | |||
| Primary relative asymmetry (AP-COP10) | 0.894 | Reference. | |
| Absolute bilateral AP-COP difference | 0.896 | AUC vs primary (95% CI to 0.008; ). | |
| AP-range-normalized bilateral difference | 0.894 | AUC vs primary (95% CI to 0.008; ). | |
| Representation attribution | |||
| Five AP-COP positions only | 0.857 | Position-only reference. | |
| Rep. | Pooled | Ga | Ju | Si |
|---|---|---|---|---|
| Full35 | 0.872 (0.813–0.921) | 0.807 (0.670–0.920) | 0.932 (0.859–0.983) | 0.867 (0.767–0.948) |
| AP-COP10 | 0.907 (0.860–0.947) | 0.837 (0.713–0.939) | 0.937 (0.865–0.988) | 0.946 (0.890–0.986) |
| Setting | Value |
|---|---|
| Participants | 165 (93 PD, 72 controls) |
| Outer evaluation | 20 repetitions 5 folds |
| Outer stratification | Study Diagnosis |
| Stability subsamples | 100 per outer fold |
| Stability threshold | 60% nonzero-selection frequency |
| Redundancy threshold | (Spearman) |