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

Abbildung 1: Ablaufschema der EEG-Vorverarbeitung, Merkmalsextraktion und Klassifizierung. Schematische Übersicht über den Studienablauf, einschließlich Überprüfung der EEG-Daten, Bandpassfilterung und Neubezug, Extraktion der Merkmale quadratischer Mittelwert (RMS), spektrale Leistungsdichte (PSD) und Entropie, Hinzufügen des Alters, Aufteilung in Trainings- und Testdatensatz, Entwicklung eines Random-Forest-Modells, Analyse der Merkmalstabilität, SHAP-Analyse, Effektgrößenanalyse nach Cohens d, integrierte Merkmalsauswahl, Abtragungsentropie-Analyse, erneutes Training des Modells und Leistungsbewertung. Bitte klicken Sie hier, um eine vergrößerte Version dieser Abbildung anzusehen.
1. Datenerfassung
Der öffentlich zugängliche Ruhezustand-EEG-Datensatz wurde aus dem OpenNeuro-Repository bezogen (Datensatz ds004504, Version 1.0.9)33. Teilnehmerbezeichnungen und demografische 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 Identifikationsnummer des Teilnehmers 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, dabei aber alle klinisch relevanten EEG-Informationen für diese Studie erhalten.
3. Bandpass-Filterung
Die kontinuierlichen EEG-Aufzeichnungen wurden mit einem Bandpass-Filter vierter Ordnung gefiltert, das eine untere Grenzfrequenz von 0,5 Hz und eine obere Grenzfrequenz von 40 Hz aufweist.
4. Gemeinsame Durchschnittsreferenzierung und Epochierung
Auf die gefilterten EEG-Aufzeichnungen wurde die gemeinsame Durchschnittsreferenzierung (CAR) 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. Merkmalsvorbereitung
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 diagnostischen Kennzeichnungen zusammengeführt, 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 enthielt fünf Root-Mean-Square-Merkmale (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 Verwechslung nicht vollständig ausschließen. Die diagnostische Gruppe, definiert als AD, HC oder FTD, wurde als Zielkennzeichnung festgelegt.
7. Daten-Vorverarbeitung und Aufteilung in Trainings- und Testdatensatz
Während der Vorverarbeitung wurden die Aufzeichnungen von 500 Hz auf 250 Hz heruntergetastet, um den Rechenaufwand zu verringern, während gleichzeitig die interessierenden EEG-Frequenzinformationen erhalten blieben.
Der Datensatz wurde auf Ebene der Versuchspersonen mithilfe einer gruppierten Aufteilungsstrategie in Trainings- (80 %) und Testteilmenge (20 %) unterteilt. Probandenbezogene Vorhersagen wurden anschließend durch Mehrheitsabstimmung über die den einzelnen Teilnehmern zugehörigen vorhergesagten Epochen gewonnen. Die primäre Bewertung erfolgte mittels einer gruppierten Aufteilung in Trainings- und Testmenge auf Ebene der Versuchspersonen, um sicherzustellen, dass Epochen ein und desselben Teilnehmers nicht in beiden Teilmengen enthalten waren. Probandenbezogene Vorhersagen wurden durch Mehrheitsabstimmung über die den einzelnen Teilnehmern zugehörigen vorhergesagten Epochen ermittelt. Der Datensatz wurde auf fehlende Werte überprüft, und fehlende Beobachtungen wurden gegebenenfalls entfernt oder durch Imputation ersetzt. Diagnosekategorien wurden kodiert. 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 Stichproben 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 (Accuracy), Präzision, Sensitivität (Recall) sowie der F1-Score wurden gemeinsam 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 dem Trainings- oder dem Testdatensatz 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-Evaluation auf Ebene der Versuchspersonen betrachtet.
11. Analyse der Merkmalsstabilität
Die Merkmalsstabilität wurde überprüft, indem das Random-Forest-Training 10-mal mit Zufallszahlen von 0 bis 9 als Startwerten wiederholt wurde. Für jeden Durchlauf wurden die Testgenauigkeit und die Werte für die Merkmalswichtigkeit aufgezeichnet. Der Mittelwert und die Standardabweichung des Wichtigkeitswerts jedes Prädiktors wurden über die 10 Durchläufe berechnet, und die Prädiktoren wurden entsprechend ihrer Stabilität gerankt. Diese Analyse diente dazu, die Konsistenz der Rangfolge der Merkmale zu bewerten, und nicht, um die primäre subjektorientierte 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 jeder Merkmalsgröße für die Modellvorhersagen abzuschätzen. Der mittlere absolute SHAP-Wert wurde für jedes Merkmal ermittelt, und die Merkmale wurden entsprechend ihrer SHAP-Beiträge ranggeordnet. Merkmale mit durchgängig niedrigen SHAP-Werten wurden identifiziert und mit den Ergebnissen der Merkmalsstabilitäts- und Cohen’s-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 Effektstärken 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. Vorhersagevariablen, die durchgängig eine geringe Merkmalswichtigkeit, einen geringen SHAP-Beitrag und kleine Effektstärken 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 unter Verwendung eines One-gegenüber-Reststrategie. Klassenspezifische Werte für die Fläche unter der Kurve (AUC) sowie die durchschnittliche AUC wurden berechnet.
17. Leistungsvergleich
Die Leistung des vollständigen Merkmalmodells wurde mit der des reduzierten Merkmalmodells verglichen, das nach der Entropie-Merkmals-Ablation erhalten wurde.