Functional Connectivity-guided spectral band selection for Motor Imagery Brain-Computer Interfaces (MI-BCIs). Ranks frequency bands using phase-based connectivity (wPLI, PLV, PLI) across sensorimotor channels to identify subject-specific discriminative bands, reducing CSP pipeline dimensionality while maintaining classification accuracy.
Instrucciones de origen · Vista previa de solo lectura
name
fc-guided-band-selection-bci
description
Functional Connectivity-guided spectral band selection for Motor Imagery Brain-Computer Interfaces (MI-BCIs). Ranks frequency bands using phase-based connectivity (wPLI, PLV, PLI) across sensorimotor channels to identify subject-specific discriminative bands, reducing CSP pipeline dimensionality while maintaining classification accuracy.
["FC-guided band selection","functional connectivity BCI","MI-BCI spectral selection","wPLI PLV band selection","filter bank CSP band selection","phase connectivity motor imagery","subject-specific frequency bands","hemispheric coupling BCI"]
This skill implements Functional Connectivity (FC)-guided spectral band selection for Motor Imagery Brain-Computer Interfaces (MI-BCIs), based on the methodology described in "Functional Connectivity-Guided Band Selection for Motor Imagery Brain-Computer Interfaces" (arXiv: 2605.00746, Araújo do Carmo & Nagarajan, 2026).
The core insight: rather than using predefined frequency sub-bands in Filter Bank CSP (FBCSP), this method identifies the most discriminative spectral bands by calculating phase-based functional connectivity across sensorimotor channels, ranking bands by the effect size of their hemispheric coupling differences, and pruning to the top-K bands for feature extraction and classification.
Problem Statement
CSP Performance Dependency on Spectral Range
The Common Spatial Pattern (CSP) algorithm is the cornerstone of MI-BCI decoding, but its performance depends critically on the spectral range of the input EEG data:
Individual variability: User-specific neural rhythms (μ and β bands) vary significantly across individuals in both center frequency and bandwidth
Heuristic limitations: Standard FBCSP uses predefined, evenly-spaced frequency sub-bands (e.g., 4-8, 8-12, 12-16 Hz) that are not selected using subject-specific physiological criteria
Dimensionality waste: Using all bands in a filter bank increases computational cost and may include non-discriminative frequency ranges, introducing noise and overfitting
Volume conduction artifacts: Traditional coherence measures are contaminated by zero-lag volume conduction, obscuring true neural connectivity
Key Questions Addressed
Which frequency bands carry the most discriminative information for a given subject?
What is the minimum number of bands (K*) required to maintain performance within a 2% equivalence zone of the full filter bank baseline?
Which phase-based connectivity metric (wPLI, PLV, or PLI) provides the best trade-off between dimensionality reduction and inter-session robustness?
Methodology
1. Filter Bank Design
The EEG signal is decomposed into a filter bank of 9 contiguous sub-bands spanning 4–40 Hz:
Band Index
Frequency Range (Hz)
Physiological Interpretation
1
4–8
Theta
2
8–12
μ (mu) — sensorimotor rhythm
3
12–16
Low β (beta)
4
16–20
Low-mid β
5
20–24
Mid β
6
24–28
Mid-high β
7
28–32
High β
8
32–36
High β / low γ
9
36–40
Low γ
Each band is extracted using a zero-phase bandpass filter (e.g., 4th-order Butterworth, forward-backward via scipy.signal.filtfilt) to avoid phase distortion.
2. Phase-Based Connectivity Metrics
Three complementary phase-based connectivity metrics are computed for each band:
Weighted Phase Lag Index (wPLI)
Measures the asymmetry of the imaginary component of the cross-spectrum, weighted by magnitude:
wPLI = |E[|Im{X}| · sgn(Im{X})]| / E[|Im{X}|]
Advantage: Robust to volume conduction and common source artifacts
Mitigates volume conduction by weighting non-zero-lag components
More stable across recording sessions with varying electrode setups
Better generalization to unseen data
Weaknesses:
May require more bands (higher K*) to reach equivalence zone
Computationally slightly more expensive than PLV
May miss genuine zero-lag functional connectivity
Practical Recommendation
Use Case
Recommended Metric
Maximum compression, single-session
PLV
Cross-session / longitudinal studies
wPLI
Quick screening / resource-constrained
PLI
Clinical / real-world BCI deployment
wPLI
Research / maximum accuracy
Compare all three
Usage Examples
Basic Usage
from your_module import FCGuidedBandSelector, FBCSPWithFCSelection
# Initialize with PLV for aggressive dimensionality reduction
pipeline = FBCSPWithFCSelection(sfreq=250, metric='plv', n_csp_components=4)
# Fit on training data (shape: n_trials x n_channels x n_timepoints)
pipeline.fit(X_train, y_train, k=4) # Select top 4 bands# Predict
predictions = pipeline.predict(X_test)
# Access selected bandsprint(f"Selected bands: {pipeline.selector.selected_bands_}")
print(f"Band rankings: {pipeline.selector.band_rankings_}")
Finding Optimal K*
# Evaluate all K values
selector = FCGuidedBandSelector(sfreq=250, metric='wpli')
# ... (compute connectivity and rank bands as above) ...# Get accuracies for each K (requires cross-validation loop)
accuracies_k = {}
for k inrange(1, 9):
acc = cross_validate_with_k(X, y, k, ranked_indices)
accuracies_k[k] = acc
# Find optimal K*
baseline_acc = cross_validate_with_k(X, y, 9, ranked_indices)
k_star = selector.find_optimal_k(ranked_indices, accuracies_k, baseline_acc)
print(f"Optimal K*: {k_star} (baseline accuracy: {baseline_acc:.3f})")
Working with MNE-Python Datasets
import mne
from mne.datasets import eegbci
# Load BCI Competition IV-2a data (via MNE)# Note: Use the official dataset download for Competition IV-2a
raw = mne.io.read_raw_edf('A01T.edf', preload=True)
raw.filter(0.5, 100)
raw.notch_filter([50])
# Extract epochs
events = mne.find_events(raw, stim_channel='STI 014')
epochs = mne.Epochs(raw, events, event_id={
'left_hand': 1, 'right_hand': 2,
'foot': 3, 'tongue': 4
}, tmin=0, tmax=4, baseline=None, preload=True)
# Pick sensorimotor channels
epochs.pick_channels(['C3', 'CP3', 'C4', 'CP4'])
# Convert to numpy array
X = epochs.get_data()
y = epochs.events[:, 2]
# Run FC-guided band selection
pipeline = FBCSPWithFCSelection(sfreq=250, metric='wpli')
pipeline.fit(X, y, k=4)
Configuration Parameters
Parameter
Type
Default
Description
sfreq
float
250
EEG sampling frequency (Hz)
metric
str
'wpli'
Connectivity metric: 'wpli', 'plv', 'pli'
band_edges
list
9 bands 4–40 Hz
Custom filter bank definitions
n_csp_components
int
4
Number of CSP spatial filters per band
k
int
None (use all)
Number of top bands to select
equivalence_tolerance
float
0.02
2% accuracy tolerance for K* selection
channel_pairs
list
4 pairs
Sensorimotor channel indices for FC
Limitations and Considerations
Static vs. Adaptive: This method uses static band selection (computed once per subject). For adaptive BCIs, consider re-ranking bands online as neural patterns shift.
Channel Selection: The 4-channel sensorimotor selection is optimal for hand MI. Foot and tongue imagery may require additional channel pairs (e.g., Cz, FCz).
Filter Design: The 4 Hz band width is a design choice. Narrower bands increase frequency resolution but may reduce signal-to-noise ratio.
Effect Size Sensitivity: Cohen's d assumes approximately normal distributions. For small trial counts, consider permutation-based effect size estimation.
Multi-class Extension: The current methodology is demonstrated on 2-class MI. For 4-class (Competition IV-2a), apply the ranking per class pair and aggregate via voting or mean effect size.
Computational Cost: Connectivity computation scales as O(n_bands × n_channel_pairs × n_trials × n_timepoints). For large datasets, consider subsampling trials or using parallel computation.
Related Skills
eeg-brain-connectivity-bci: General EEG-based brain connectivity analysis for BCI applications
fc-guided-band-selection-mi-bci: Variant focusing on motor imagery-specific implementations
hermes-brain-connectivity: Comprehensive brain connectivity toolkit with multiple modalities
References
Araújo do Carmo, N., & Nagarajan, A. (2026). Functional Connectivity-Guided Band Selection for Motor Imagery Brain-Computer Interfaces. arXiv:2605.00746 [q-bio.NC, eess.SP].
Blankertz, B., et al. (2008). Optimizing spatial filters for robust EEG single-trial analysis. IEEE Signal Processing Magazine, 25(1), 41-56.
Vinck, M., et al. (2011). An improved index of phase-synchronization for electrophysiological data in the presence of volume-conduction, noise and sample-size bias. NeuroImage, 55(4), 1548-1565.
Lachaux, J.-P., et al. (1999). Measuring phase synchrony in brain signals. Human Brain Mapping, 8(4), 194-208.
Stam, C. J., et al. (2007). Phase lag index: Assessment of functional connectivity from multi channel EEG and MEG with diminished bias from common sources. Human Brain Mapping, 28(11), 1178-1193.
Ang, K. K., et al. (2008). Filter bank common spatial pattern (FBCSP) in brain-computer interface. IJCNN, 2390-2397.
