Eine schematische Übersicht des EEG-Präprozessing-, Merkmalsextraktions- und Klassifizierungsworkflows ist in Abbildung 1 dargestellt. Die in der Studie verwendeten Materialien und Software sind in der Tabelle der Materialien aufgelistet.

Abbildung 1: Ablauf der EEG-Vorverarbeitung, Merkmalsextraktion und Klassifizierungspipeline. Schematische Übersicht des Studienablaufs, einschließlich Überprüfung der EEG-Daten, Bandpassfilterung und Neubezug, Extraktion der Effektivwertmerkmale (RMS), der spektralen Leistungsdichte (PSD) und der Entropiemerkmale, Hinzufügen des Alters, Aufteilung in Trainings- und Testdatensatz, Entwicklung des Random-Forest-Modells, Analyse der Merkmalsstabilität, SHAP-Analyse, Cohen’s-d-Effektgrößenanalyse, integrierte Merkmalsauswahl, Ablation der Entropiemerkmale, erneutes Training des Modells und Leistungsbewertung. Bitte klicken Sie hier, um eine größere Version dieser Abbildung anzusehen.
1. Datenerfassung
Der öffentlich zugängliche Ruhezustands-EEG-Datensatz wurde aus dem OpenNeuro-Repository bezogen (Datensatz ds004504, Version 1.0.9)33. Die Teilnehmerbezeichnungen und demografischen Informationen wurden überprüft. Der Datensatz umfasste 88 Teilnehmer, darunter 36 Teilnehmer mit Alzheimer-Krankheit (AD), 23 mit frontotemporaler Demenz (FTD) und 29 gesunde Kontrollpersonen (HC).
2. Import der EEG-Aufzeichnung
Für jeden Teilnehmer wurde die EEG-Aufzeichnung mithilfe der MNE-Python-Bibliothek geladen. Jede EEG-Datei wurde überprüft, um sicherzustellen, dass sie zugänglich war, erfolgreich importiert werden konnte und eine gültige Probanden-ID enthielt. Aufzeichnungen mit fehlenden oder beschädigten Dateien wurden ausgeschlossen. Die ursprüngliche EEG-Abtastrate von 500 Hz wurde in diesem Schritt beibehalten. Die Signale wurden heruntergetastet, um den Rechenaufwand zu verringern, während alle klinisch relevanten EEG-Informationen für diese Studie erhalten blieben.
3. Bandpassfilterung
Die kontinuierlichen EEG-Aufzeichnungen wurden mit einem Bandpassfilter vierter Ordnung mit einer unteren Grenzfrequenz von 0,5 Hz und einer oberen Grenzfrequenz von 40 Hz gefiltert.
4. Gemeinsame Durchschnittsreferenzierung und Epochierung
Die gemeinsame Durchschnittsreferenzierung (CAR) wurde auf die gefilterten EEG-Aufzeichnungen angewendet. Für jede Zeitstichprobe wurde das durchschnittliche Signal über alle verfügbaren EEG-Kanäle berechnet und von jedem Kanal subtrahiert. Die referenzierten EEG-Signale wurden für die anschließende Analyse beibehalten. Jede referenzierte EEG-Aufzeichnung wurde mithilfe der Funktion make_fixed_length_epochs() in MNE-Python in Epochen fester Länge unterteilt. Dabei wurde eine Epochenlänge von 10 s mit 0 s Überlappung verwendet.
5. Epoche-Validierung
Für jeden Teilnehmer wurden nur vollständige 10-s-Epochen beibehalten. Alle verbleibenden EEG-Segmente, die am Ende einer Aufnahme kürzer als 10 s waren, wurden verworfen. Jede beibehaltene Epoche wurde anschließend als individuelle Probe für die Merkmalsextraktion verwendet.
6. Merkmalsaufbereitung
Obwohl Mini-Mental-State-Examination-(MMSE-)Daten im Quelldatensatz verfügbar waren, wurden sie von der vorliegenden Analyse ausgeschlossen. Die extrahierten EEG-Merkmale wurden mit dem demografischen Teilnehmerinformationen wie Alter und Diagnosebezeichnungen verknüpft, die aus der Metadatendatei participants.csv stammten. Als gemeinsamer Schlüssel diente die Teilnehmerkennung, und es wurden nur Teilnehmer beibehalten, für die sowohl EEG-Merkmalsdatensätze als auch entsprechende demografische Informationen vorlagen.
Die mehrdimensionale Merkmalsmatrix umfasste fünf Effektivwertmerkmale (RMS; delta_rms, theta_rms, alpha_rms, beta_rms und gamma_rms), fünf Leistungsdichtespektrum-Merkmale (PSD; delta_psd, theta_psd, alpha_psd, beta_psd und gamma_psd) sowie fünf Entropie-Merkmale (delta_entropy, theta_entropy, alpha_entropy, beta_entropy und gamma_entropy). Das Alter wurde als ergänzende demografische Variable hinzugefügt. Da sich die Altersverteilungen zwischen den diagnostischen Gruppen unterscheiden können, ließ sich eine demografische Störung nicht vollständig ausschließen. Die diagnostische Gruppe, definiert als AD, HC oder FTD, wurde als Zielkennung zugewiesen.
7. Daten-Vorverarbeitung und Aufteilung in Trainings- und Testdatensatz
Während der Vorverarbeitung wurden die Aufnahmen von 500 Hz auf 250 Hz heruntergerechnet, um den Rechenaufwand zu verringern, gleichzeitig aber die interessierenden EEG-Frequenzinformationen beizubehalten.
Der Datensatz wurde mithilfe einer gruppierten Aufteilungsstrategie auf Ebene der Probanden in einen Trainings- (80 %) und einen Testanteil (20 %) unterteilt. Probandenspezifische Vorhersagen wurden anschließend durch Mehrheitsabstimmung über die vorhergesagten Epochen jedes Teilnehmers ermittelt. Die primäre Bewertung erfolgte anhand einer probandenspezifischen, gruppierten Aufteilung in Trainings- und Testdatensatz, um sicherzustellen, dass Epochen desselben Teilnehmers nicht in beiden Teilmengen enthalten waren. Probandenspezifische Vorhersagen wurden durch Mehrheitsabstimmung über die den einzelnen Teilnehmern zugehörigen vorhergesagten Epochen gewonnen. Der Datensatz wurde auf fehlende Werte überprüft, und fehlende Beobachtungen wurden entsprechend entfernt oder imputiert. Diagnosekennungen wurden codiert. Ein StandardScaler wurde an den Trainingsdaten angepasst und anschließend sowohl auf die Trainings- als auch auf die Testdatensätze angewandt.
8. Entwicklung des Random-Forest-Modells
Ein Random-Forest-Klassifikator wurde mit 80 Bäumen, einer maximalen Baumtiefe von 10, maximal vier Merkmalen, einem Minimum von fünf Proben pro Blatt und einem Zufallszustand von 42 initialisiert. Der Klassifikator wurde mit dem standardisierten Trainingsdatensatz trainiert.
9. Modellbewertung
Es wurden Klassenbezeichnungen sowohl für die Trainings- als auch für die Testdatensätze vorhergesagt. Eine Konfusionsmatrix wurde erstellt, und Genauigkeit, Präzision, Sensitivität sowie der F1-Score wurden zusammen mit dem Klassifizierungsbericht berechnet. Die Genauigkeiten für das Training und das Testen wurden dokumentiert.
Stellen Sie für die teilnehmerweise Auswertung sicher, dass alle Epochen eines bestimmten Teilnehmers ausschließlich entweder der Trainings- oder der Testuntergruppe zugeordnet wurden. Der Random-Forest-Klassifikator wurde mit denselben Hyperparametern wie in der primären Analyse trainiert.
10. Kreuzvalidierung
Als zusätzliche Modellbewertungsmethode wurde eine fünffache stratifizierte Kreuzvalidierung mit shuffle = True und random_state = 30 durchgeführt. Die mittlere Genauigkeit und die Standardabweichung über die fünf Folds wurden berechnet. Diese Analyse wurde separat von der primären Holdout-Evaluierung auf Ebene der Probanden betrachtet.
11. Analyse der Merkmalsstabilität
Die Merkmalsstabilität wurde überprüft, indem das Random-Forest-Training 10-mal mit Zufallszahlen-Startwerten („random seeds“) von 0 bis 9 wiederholt wurde. Für jeden Durchlauf wurden die Testgenauigkeit und die Werte zur Merkmalswichtigkeit („feature-importance“) erfasst. Der Mittelwert und die Standardabweichung des Merkmalswichtigkeitswerts für jeden Prädiktor wurden über die 10 Durchläufe hinweg berechnet, und die Prädiktoren wurden entsprechend ihrer Stabilität ranggeordnet. Diese Analyse diente dazu, die Konsistenz der Merkmalsrangfolgen zu bewerten, und nicht dazu, die primäre, auf Ebene der Versuchspersonen durchgeführte Auswertung zu ersetzen.
12. Erklärbare KI-Analyse
Der SHAP TreeExplainer wurde auf das trainierte Random-Forest-Modell angewendet. SHAP-Werte wurden berechnet, um den Beitrag jedes Merkmals zu den Modellvorhersagen abzuschätzen. Der mittlere absolute SHAP-Wert wurde für jedes Merkmal ermittelt, und die Merkmale wurden entsprechend ihrer SHAP-Beiträge gerankt. Merkmale mit durchgängig niedrigen SHAP-Werten wurden identifiziert und mit den Ergebnissen der Merkmalsstabilitäts- und Cohens-d-Analysen verglichen. Merkmale mit durchgängig geringen Beiträgen wurden für die anschließende Ablation und erneute Modelltrainierung ausgewählt.
13. Statistische Effektgrößenanalyse
Cohens d wurde für jeden EEG-Biomarker für die Vergleiche AD gegenüber HC, AD gegenüber FTD sowie FTD gegenüber HC berechnet. Die Größe der Effekte wurde anhand von Schwellenwerten von 0,20 für einen kleinen Effekt, 0,50 für einen mittleren Effekt und 0,80 für einen großen Effekt interpretiert.
14. Integrierte Merkmalsauswahl
Die Ergebnisse der Random-Forest-Merkmalswichtigkeitsanalyse, der SHAP-Analyse und der Cohen’s-d-Analyse wurden verglichen. Prädiktoren, die durchgängig eine geringe Merkmalswichtigkeit, einen geringen SHAP-Beitrag und kleine Effektgrößen aufwiesen, wurden als Kandidaten für das Entfernen identifiziert.
15. Merkmalsablation
Entropiemerkmale wurden in einer explorativen Ablationsanalyse bewertet, und eine reduzierte Merkmalsmatrix, die RMS, PSD und Alter enthielt, wurde erstellt. Der Random-Forest-Klassifikator wurde mit denselben Hyperparametern erneut trainiert. Modelltraining, -test, Kreuzvalidierung und Receiver-Operating-Characteristic-(ROC)-Analyse wurden mit dem reduzierten Merkmalssatz wiederholt.
16. ROC-Analyse
Die Klassenwahrscheinlichkeiten wurden aus dem optimierten Random-Forest-Klassifikator erhalten. Multiklassen-ROC-Kurven wurden mithilfe eines One-gegenüber-Rest-Strategie. Klassenspezifische Flächen unter der Kurve (AUC)-Werte und die durchschnittliche AUC wurden berechnet.
17. Leistungsvergleich
Die Leistung des vollständigen Merkmalsmodells wurde mit der des reduzierten Merkmalsmodells verglichen, das nach der Entropie-Merkmalsablation erhalten wurde.