$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
Die vorgeschlagene Methode segmentiert EEG-Signale strategisch in den Frequenzbereich, um die spektrale Auflösung zu verbessern. Anschließend werden die erfassten EEG-Signale in 10 Segmente aufgeteilt, um die zeitliche Auflösung und Konsistenz in der Merkmalsberechnung zu verbessern. Die vorgeschlagene Methode implementiert einen RF-, XGBoost- und SVM-basierten Klassifikator auf abgeleiteten EEG-Merkmalen und trägt so zur Entwicklung robuster, skalierbarer Methoden zur EEG-Analyse bei. Das RF-Modell mit bestimmten Hyperparameteränderungen hat sich als das leistungsstärkste unter den Standardmodellen erwiesen. Daher wird ein hybrides, fusioniertes, gestapeltes Metalearner-Modell vorgeschlagen, das RF und XG-Boost kombiniert, um die Vorhersagegenauigkeit weiter zu verbessern.
Der bei dieser Methode verwendete Datensatz wurde gemäß institutionellen ethischen Standards erhoben. Die Datensatzanbieter haben vor der Datenerhebung die Zustimmung der Teilnehmer eingeholt.
1. Datensatzspezifikation
Um die Robustheit des Systems zu beweisen, wurden zwei Datensätze verwendet. Der erste Datensatz stammt aus dem EEG-Datenarchiv20 und umfasst 88 Teilnehmer im Ruhezustand mit geschlossenen Augen, wobei 36 als AD diagnostiziert wurden, 23 davon FTD und 29 gesund. EEG-Aufnahmen wurden im Ruhezustand mit geschlossenen Augen und einem 19-Kanal-System nach der internationalen 10-20-Montage aufgenommen. Im in Abbildung 2 gezeigten Datensatz wurden alle Signale mit 500 Hz abgetastet und im BIDS-Format mit vorläufiger Vorverarbeitung durch die Datensatz-Beitragenden bereitgestellt.
Der zweite Datensatz wird aus einem externen öffentlichen Repository21 mit 35 Teilnehmern im Ruhezustand gesammelt. Von den 13 Teilnehmern wurde Alzheimer (AD) diagnostiziert, 7 Patienten mit leichter kognitiver Beeinträchtigung (MCI) und 15 sind gesunde ältere Menschen. Nur die Ruhezustands-Baseline-Segmente des zweiten Datensatzes wurden verwendet, um die Konsistenz mit dem ersten Datensatz20 zu gewährleisten. Der zweite Datensatz21 erforderte eine vollständige Vorverarbeitungspipeline, wie in Abbildung 3 erwähnt.
2. Vorbearbeitung
Bestehende Studien stützen sich stark auf komplexe Vorverarbeitungspipelines, manuelle Artefaktkorrekturschritte oder ICA-ASR, um Muskelbewegungen zu entfernen, die die Reproduzierbarkeit in routinemäßigen klinischen Arbeitsabläufen einschränken. Um diese Einschränkungen zu beheben, konzentriert sich die vorgeschlagene Pipeline auf einen stromlinienförmigen, ausschließlich EEG-Ansatz, der die Notwendigkeit rechenintensiver Artefaktentfernungsverfahren eliminiert und stattdessen kontrollierte Filterung, Epochensegmentierung und frequenzspezifische Merkmalsberechnung betont.
- EEG-Dateisortierung
Vor der Vorverarbeitung wurden alle EEG-Dateien automatisch sortiert, um eine konsistente Reihenfolge nach Teilnehmer und Sitzungen zu gewährleisten. Dies stellte sicher, dass Feature-Extraktion, Epoching und Label-Ausrichtung bei allen Probanden einheitlich waren und Dateifehler verhinderten. Dieser Schritt ist in Abbildung 3 als "EEG-Datei sortieren" dargestellt. Alle Vorverarbeitung erfolgt mit einem standardisierten Open-Source-EEG-Verarbeitungstoolkit.
- Bandpassfilterung und Neureferenzierung
Ein Butterworth-Bandpassfilter wurde im Bereich von 0,5–40 Hz angelegt, um relevante EEG-Komponenten zu erhalten und gleichzeitig Drift und hochfrequentes Rauschen zu dämpfen. Die Signale wurden auf den Durchschnitt aller EEG-Kanäle umbezogen, wodurch die räumliche Vorspannung minimiert und das Signal-Rausch-Verhältnis verbessert wurde.
- Nyquist-Frequenzverschiebung
Um die Auflösung von Niederfrequenzkomponenten zu verbessern, wurde eine Nyquist-Verschiebungsverstärkung angewendet, indem die Abtastfrequenz von 500 Hz auf 250 Hz halbiert wurde. Nach dem obigen Schritt reduziert die Nyquist-Modifikation sie auf 125 Hz. Da die Delta-, Theta-, Alpha-, Beta- und Gamma-Frequenzen langsame Wellen von 0,5 bis 40 Hz sind, verschiebt sich dadurch die effektive Nyquist-Grenze und verbessert die Klarheit von niederfrequenten Schwingungen, die für Demenz relevant sind. Diese Modifikation half, langsame Gehirnwellen besser einzufangen. Gleichung (1) zeigt das Ergebnis der Nyquist-Medikation eindeutig.
F-Abtastung = 500 Hz
NAbtastung=F-Abtastung/2
NNyquist= NStichproben/2 (1)
Es ist wichtig zu beachten, dassN-Nyquist nicht die Betriebsabtastrate ist, sondern lediglich die maximal darstellbare Frequenz , und daher die Epochensegmentierung nicht beeinflusst. Es ist nur eine theoretische Grenze, die dazu dient, sicherzustellen, dass der relevante Frequenzgehalt (0,5–40 Hz) vollständig erhalten bleibt.
- Epochen-Sampling-Frequenz auf 10 Segmente
Die EEG-Aufnahmen sind in 10 gleich lange Segmente unterteilt. Dieser Schritt ist unerlässlich, um die nichtstationäre Aktivität des EEG-Signals über die Zeit zu erfassen. Die kurzen, einheitlichen Segmente ermöglichen eine stabile Berechnung bandspezifischer Merkmale und bereichern Trainingsbeispiele für die Klassifikation des maschinellen Lernens. Daher werden aus 250 Proben pro Sekunde 250 × 10 = 2500 Proben. Abschließend werden pro Fach 10 gleich lange Epochen für epochenweise Analysestrategien gewählt. Diese Analysestrategie erzeugte 880 Stichproben für 88 tatsächliche Teilnehmer.
Die Notwendigkeit einer epochenweisen Analyse besteht darin, Veränderungen der Gehirnsignale minutiös zu überwachen, wodurch dem Modell eine ausreichende Menge zeitlicher Daten zur Verfügung gestellt wird, um zu lernen und Überanpassungen zu verhindern. Die Bedeutung dieses Ansatzes besteht darin, feinkörnigere EEG-Epochendaten zu gewährleisten und Analysen zu ermöglichen, die eine höhere zeitliche Auflösung erfordern, ohne die Gesamtwerte der Merkmale zu beeinträchtigen. Durch die Erzeugung von 10-Sekunden-Epochen unterstützt die Methodik eine verbesserte Modellierung der EEG-Dynamik und entspricht den Standard-Epochendauern in EEG-Analyseprotokollen.
3. Merkmalsextraktion
Die vorverarbeiteten EEG-Signale werden an Merkmalsextraktionsmethoden weitergeleitet. EEG-Signale werden bandpassgefiltert in 5 Standardfrequenzbänder: Delta (0,5–4 Hz), Theta (4–8 Hz), Alpha (8–13 Hz), Beta (13–25 Hz) und Gamma (25–40 Hz).
Für jedes gefilterte Signal wurde der Root Mean Square (RMS)-Wert über alle EEG-Kanäle berechnet. Die mathematische Formulierung des RMS ist in Gleichung (1) gegeben:
RMS =
(1)
Hier ist xi die EEG-Signalamplitude an der i-ten Zeitprobe. N ist die Gesamtzahl der Samples im Signalsegment.
RMS wurde als Hauptmerkmal ausgewählt, da es die Energie der oszillatorischen Aktivität in jedem Frequenzband quantifizieren kann. Die Alzheimer-Krankheit ist oft mit erhöhter Delta- und Theta-Aktivität sowie mit reduzierter Alpha- und Beta-Aktivität verbunden. Fronto-temporale Demenz (FTD) kann über diese Bänder hinweg deutliche Muster zeigen(20,21). Da PSD mathematisch redundant zu RMS ist, gilt es nicht als Feature. Um einen kompakten Datensatz aufrechtzuerhalten, werden die übrigen Merkmale wie Hjorth und Entropie ausgeschlossen.
4. Hybridfusion
Extrahierte Zeitdomänen-RMS-Funktionen wurden mit Frequenzdomänenmerkmalen kombiniert, die aus der Bandpassfilterung erhalten wurden, indem sie einen hybriden Raum für die Modelleingabe bildeten. Die Details der Kompilierung sind in Abbildung 3 dargestellt. Der endgültige normalisierte Datensatz ist in einem tabellarischen Format organisiert, mit Zeilen, die Teilnehmer repräsentieren, und Spalten, die Zeit-Frequenzband-Merkmale sowie demografische Merkmale wie Alter, Geschlecht und Gruppe darstellen. Die in Abbildung 2 gezeigte Pipeline extrahiert effizient EEG-Merkmale, die für die Alzheimer-Forschung relevant sind. Die extrahierten Merkmale liefern Einblicke in neuronale Aktivitäten innerhalb wichtiger Frequenzbänder und können für maschinelles Lernen und statistische Analysen genutzt werden.
5. Funktionsauswahl
Um die Leistung des Modells zu verbessern und die Merkmalsdimensionalität zu reduzieren, wurde ein auf Varianz basierender 2-Wege-Ansatz (ANOVA) angewandt. Die ANOVA-basierte Auswahl wurde mit Gruppe und Alter als unabhängige Faktoren und RMS-Frequenzbändern als abhängige Merkmale angewendet. Diese Analyse bewertete die Auswirkungen von Gruppen-, Alters- und Altersgruppen-Interaktionen auf jedes Merkmal. Merkmale mit p-Werten < 0,05 für mindestens einen Faktor wurden für eine weitere Klassifizierung ausgewählt. Die detaillierten ANOVA-basierten Ergebnisse sind im Ergebnisabschnitt dargestellt. Diese Merkmale sind die informativsten Variablen für die Klassifikation, was dazu beigetragen hat, die Genauigkeit und Verallgemeinerung des endgültigen prädiktiven Modells zu verbessern. Als Ergebnis des ANOVA-Tests wurden nur die informativsten Funktionen für die Modellerstellung ausgewählt. Und die zusätzlichen nicht-numerischen Attribute, Geschlecht und participant_id, wurden aus der Analyse ausgeschlossen. Die Zielvariable Gruppe wurde für die Klassifikation labelkodiert.
6. Modellbeschreibung
Um die Krankheit in drei Klassen zu klassifizieren, nämlich Alzheimer, Kontroll- und Frontaldemenz, wurden drei überwachte Machine-Learning-Modelle entwickelt, wie XG-Boost, Random Forest, Support Vector Machine SVM und ein Stacked-Modell. Alle genannten Modelle verwendeten 70 % des Trainingsdatensatzes und 30 % des Testdatensatzes. Jedes Modell wurde aufgrund seiner nachgewiesenen Leistung in der Analyse von Gesundheitsdaten und seiner Fähigkeit, nichtlineare und hochdimensionale Merkmalsräume zu verarbeiten, ausgewählt. Die Bewertung des Modells erfolgte am Validierungsdatensatz, um die Generalisierbarkeit des Modells zu bewerten.
Die Modellleistung wurde mittels einer Verwirrungsmatrix bewertet, die aus wahr positiv (TP), falsch positiv (FP), echtnegativ (TN) und falsch negativ (FN) bestand. Leistungskennzahlen umfassten Genauigkeit und einen Klassifikationsbericht, der Genauigkeit, Rückruf und F1-Ergebnisse für jede Klasse detailliert darstellte. Die folgende Gleichung veranschaulicht Leistungskennzahlen.
(2)
(3)
(4)
(5)
1. XG-Boost-Klassifikator
Diese Methode implementierte ein Multiklassifikationsmodell mit dem XGBoost-Algorithmus, um EEG-abgeleitete kognitive Zustände aus extrahierten Merkmalen vorherzusagen, wie in Ausgabe 1 von Abbildung 3 gezeigt. Ausgabe 1 enthält RMS-Funktionen und den MMSE-Score. Beim Bau des XGBoost-Modells wird der MMSE-Score einbezogen, um das nichtlineare Verhalten der RMS-Funktionen zu testen. Sie dient als Maßstab und Basis für die Systementwicklung. MMSE wurde in späteren Phasen ausgeschlossen, um das Ziel eines vollautomatisierten, klinikerunabhängigen Systems zu erreichen.
XG-Boost ist ein Gradient-Boosting-Framework, das auf Effizienz und Genauigkeit optimiert ist und es gut geeignet macht, strukturierte Daten in Klassifikationsaufgaben zu verarbeiten. Die wenigen Hyperparameter werden auf dem EEG-Datensatz mit Techniken wie Crossvalidierung und Rastersuche abgestimmt. Die unten definierten Schlüsselparameterwerte werden gewählt, um Overfitting durch Regularisierung zu steuern und die Komplexität der EEG-Merkmale zu erfassen, ohne zu komplex zu sein.
Die Schlüsselparameter des Modells sind: (1) Maximale Baumtiefe = 8, was die Tiefe der Entscheidungsbäume begrenzt, um Überanpassung zu reduzieren und die Verallgemeinerung zu verbessern. Denn ein flacher Baum würde unnötige Informationen verpassen, da ein tieferer Baum überpassen könnte. (2) L2-Regularisierung (λ = 10) fügt eine Strafe für große Koeffizienten hinzu, um Überanpassungen durch Kontrolle der Modellkomplexität zu minimieren. (2) L1-Regularisierung (α = 5) wird eingeführt, um zusätzliche Sparsamkeit im Modell zu gewährleisten und so die Interpretierbarkeit und Robustheit zu verbessern. Die Anzahl der Schätzer (n_estimators=8) begrenzte die Anzahl der Boosting-Iterationen, um die Recheneffizienz aufrechtzuerhalten und gleichzeitig die Leistung auszubalancieren.
2. Random Forest-Klassifikator
Random Forest ist ein Ensemble-Klassifikator, der mehrere Entscheidungsbäume auf zufälligen Teilmengen von Daten und Features erstellt und deren Ausgaben für die endgültige Vorhersage aggregiert. Es bietet Robustheit gegenüber Rauschen, verarbeitet nichtlineare Daten gut und reduziert Overfitting durch die Mittelung verschiedener Baumvorhersagen. In dieser Arbeit wurde das Random Forest-Modell mit Schlüsselparametern mittels der Trial-and-Error-Methode konfiguriert. Solche Parameter werden im Folgenden erklärt:
1) n_estimators=100
Sie wird gewählt, um die Anzahl der Entscheidungsbäume im Wald zu begrenzen. Je mehr Bäume, desto besser die Verallgemeinerung und desto stabiler sind die Vorhersagen.
2) max_depth=10
Es gibt die maximale Tiefe jedes Baumes an. Wenn der Baum eine geringe Tiefe hat, besteht die Möglichkeit eines Underfittings. Während hohe Tiefe Details erfassen und die Trainingsdaten auswendig lernen kann, kann sie überfitten. Daher ist die Wahl von 10 ein Sweet Spot, besonders für verrauschte oder komplexe Signale wie das EEG.
3) random_state=40
Es behebt den Seed zur Zufallszahlengenerierung, der für Bootstrap-Sampling und Baumbau verwendet wird. Dies gewährleistet die Reproduzierbarkeit
4) n_jobs=-1
Indem der Wert von n_jobs auf -1 gesetzt wird, nutzt das Modell alle CPU-Kerne effizient.
7. Support Vector Machine (SVM)
Die Support Vector Machine ist ein marginbasierter Klassifikator, der für seine lineare Leistung in binärer und multiklassiger Klassifikation22 bekannt ist. Dieses Modell verwendete einen radialen Basisfunktionskern, da er nichtlineare EEG-Signale erfassen kann. Das Modell zeigte jedoch eine suboptimale Leistung, vermutlich aufgrund fehlender Feature-Skalierung und der begrenzten Datensatzgröße. Obwohl SVM theoretisch leistungsfähig ist, könnte seine Empfindlichkeit gegenüber Hyperparameter-Tuning und Datenverteilung die geringere Genauigkeit im Vergleich zu baumbasierten Modellen erklären.
8. Vorgeschlagenes HY-Fusionsmodell
Obwohl das Random Forest-Modell eine hohe Genauigkeit erzielte, wurde zur Gewährleistung von Generalisierbarkeit und Skalierbarkeit ein gestapeltes Modell mit Random Forest und XG-Boost erstellt. Die Ausgabe wird als Eingabe in ein logistisches Regressionsmodell eingespeist. Und logistische Regression wirkt als Metalerner. Das Architekturdiagramm des Modells ist in Abbildung 4 dargestellt.
Neben demersten Satz von Eingabe20 nahm das gestapelte Modell auch Eingaben aus dem zweiten Datensatz23. Der Datensatz aus Mendely21 wird vorverarbeitet und in eine modellverständliche Form verarbeitet, indem RMS extrahiert und auf 10 Sekunden epochiert werden. Die von jedem Modell erhaltenen Ergebnisse werden im Ergebnisabschnitt dargestellt.