De voorgestelde methode segmenteert EEG-signalen strategisch in het frequentiedomein om de spectrale resolutie te verbeteren. Vervolgens worden de vastgelegde EEG-signalen opgesplitst in 10 segmenten om de temporele resolutie en consistentie in feature-berekening te verbeteren. De voorgestelde methode implementeert een op RF, XGBoost en SVM gebaseerde classifier op afgeleide EEG-kenmerken, waarmee hij bijdraagt aan de ontwikkeling van robuuste, schaalbare methoden voor EEG-analyse. Het RF-model met bepaalde hyperparameterwijzigingen is het best presterende van de andere standaardmodellen gebleken. Daarom wordt een hybride, gefuseerde, gestapelde metalearner-model voorgesteld dat RF en XG-Boost combineert om de voorspellingsnauwkeurigheid verder te verbeteren.
De dataset die bij deze methode werd gebruikt, werd verzameld volgens institutionele ethische normen. De datasetproviders kregen toestemming van deelnemers voorafgaand aan de gegevensverzameling.
1. Datasetspecificatie
Om de robuustheid van het systeem te bewijzen, zijn twee datasets gebruikt. De eerste set gegevens is verzameld uit het EEG-datarepository20, bestaande uit 88 deelnemers in de rusttoestand, met gesloten ogen, waarvan 36 werden gediagnosticeerd met AD, 23 van hen FTD en 29 gezond. EEG-opnames werden in rusttoestand en ogen gesloten gemaakt met een 19-kanaals systeem na de internationale 10-20 montage. In de dataset die in Figuur 2 wordt getoond, werden alle signalen gesampled op 500 Hz en geleverd in BIDS-formaat met voorlopige voorbewerking door de dataset-bijdragers.
De tweede set gegevens wordt verzameld uit een externe openbare repository21 met 35 deelnemers in rusttoestand. Van de 13 deelnemers zijn gediagnosticeerd met de ziekte van Alzheimer (AD), 7 zijn patiënten met milde cognitieve beperking (MCI) en 15 zijn gezonde ouderen. Alleen de resting-state baseline segmenten van de tweede dataset werden gebruikt om consistentie te behouden met de eerste dataset20. De tweede set data21 vereiste een volledige pre-processing pipeline, zoals vermeld in Figuur 3.
2. Pre-processing
Bestaande studies zijn sterk afhankelijk van complexe pre-processing pipelines, handmatige artefactcorrectiestappen of ICA-ASR om spierbewegingen te verwijderen, wat de reproduceerbaarheid in routinematige klinische workflows beperkt. Om deze beperkingen aan te pakken, richt de voorgestelde pijplijn zich op een gestroomlijnde, uitsluitend EEG-benadering die de noodzaak van computationeel intensieve procedures voor het verwijderen van artefacten elimineert en in plaats daarvan de nadruk legt op gecontroleerde filtering, epochsegmentatie en frequentie-specifieke feature-computation.
- EEG-bestandssortering
Voor de preprocessing werden alle EEG-bestanden automatisch gesorteerd om een consistentie van de volgorde per deelnemer en sessie te behouden. Dit zorgde ervoor dat feature-extractie, epoching en label-uitlijning uniform waren over alle onderwerpen en voorkwam fouten in bestandsafstemming. Deze stap wordt weergegeven in Figuur 3 als "Sort EEG-bestand." Alle voorbewerking gebeurt met behulp van een standaard open-source EEG-verwerkingstoolkit.
- Bandpassfiltering en herreferentie
Een Butterworth-banddoorlaatfilter werd toegepast in het bereik van 0,5-40 Hz om relevante EEG-componenten te behouden en drift en hoogfrequente ruis te dempen. Signalen werden opnieuw verwezen naar het gemiddelde van alle EEG-kanalen, waardoor ruimtelijke bias werd geminimaliseerd en de signaal-ruisverhouding werd verbeterd.
- Nyquist-frequentieverschuiving
Om de resolutie van laagfrequente componenten te verbeteren, werd een Nyquist-shift-verbetering toegepast door de bemonsteringsfrequentie te halveren van 500Hz naar 250Hz. Na bovenstaande stap wordt het door de Nyquist-modificatie teruggebracht tot 125Hz. Omdat de Delta-, Theta-, Alpha-, Beta- en Gamma-frequenties langzame golven zijn variërend van 0,5-40 Hz, verschuift dit de effectieve Nyquist-limiet en verbetert het de helderheid van laagfrequente oscillaties die relevant zijn voor dementie. Deze aanpassing hielp om langzame hersengolven beter vast te leggen. Vergelijking (1) toont het medicatieresultaat van Nyquist duidelijk.
F-sampling= 500 Hz
Nbemonstering= Fbemonstering/2
Nnyquist= Nbemonstering/2 (1)
Het is belangrijk op te merken dat Nnyquist niet de operationele bemonsteringsfrequentie is, maar simpelweg de maximaal representatieve frequentie , en daarom de epochsegmentatie niet beïnvloedt. Het is slechts een theoretische limiet die wordt gebruikt om te garanderen dat de frequentie van belang (0,5-40 Hz) volledig behouden blijft.
- Epoch-bemonsteringsfrequentie naar 10 segmenten
De EEG-opnames zijn verdeeld in 10 segmenten van gelijke lengte. Deze stap is essentieel om de niet-stationaire activiteit van het EEG-signaal in de loop van de tijd vast te leggen. De korte, uniforme segmenten maken stabiele berekening van bandspecifieke kenmerken mogelijk en verrijken trainingsvoorbeelden voor machine-learning classificatie. Daardoor worden de 250 monsters per seconde 250 × 10 = 2500 monsters. Tot slot worden per onderwerp 10 epochs van gelijke lengte genomen voor epochgewijze analysestrategieën. Deze analysestrategie genereerde 880 steekproeven voor 88 daadwerkelijke deelnemers.
De noodzaak voor epoch-gewijze analyse is om veranderingen in hersensignalen minutieus te monitoren, zodat het model voldoende temporele data krijgt om te leren en overfitting te voorkomen. Het belang van deze aanpak is het waarborgen van fijnmazigere EEG-epochgegevens, waardoor analyses mogelijk worden die een hogere temporele resolutie vereisen zonder de totale featurewaarden in gevaar te brengen. Door het genereren van 10-seconden epochs ondersteunt de methodologie verbeterde modellering van EEG-dynamiek en sluit deze aan bij standaard epoche-duren in EEG-analyseprotocollen.
3. Feature-extractie
De voorbewerkte EEG-signalen worden doorgegeven aan feature-extractiemethoden. EEG-signalen worden banddoorlaatgefilterd in 5 standaard frequentiebanden: Delta (0,5-4 Hz), Theta (4-8 Hz), Alpha (8-13 Hz), Beta (13-25 Hz) en Gamma (25-40 Hz).
Voor elk gefilterd signaal werd de Root Mean Square (RMS)-waarde berekend over alle EEG-kanalen. De wiskundige formulering van RMS wordt gegeven in Vergelijking (1):
RMS =
(1)
Hier is xi de EEG-signaalamplitude bij het i-de tijdmonster. N is het totale aantal samples in het signaalsegment.
RMS werd gekozen als primaire kenmerk vanwege het vermogen om de energie van oscillerende activiteit in elke frequentieband te kwantificeren. De ziekte van Alzheimer wordt vaak geassocieerd met verhoogde delta- en theta-activiteit en verminderde alfa- en bèta-activiteit. Fronto-temporale dementie (FTD) kan duidelijke patronen vertonen over deze banden20,21. Omdat PSD wiskundig overbodig is ten opzichte van RMS, wordt het niet als een functie beschouwd. Om een compacte dataset te behouden, worden de overige kenmerken, zoals Hjorth en Entropie, uitgesloten.
4. Hybride fusie
Geëxtraheerde tijdsdomein-RMS-functies werden gecombineerd met frequentiedomeinkenmerken ontvangen uit banddoorlaatfiltering door het een hybride ruimte voor modelinvoer te maken. De details van de compilatie worden weergegeven in Figuur 3. De uiteindelijke genormaliseerde dataset is georganiseerd in een tabelformaat met rijen die deelnemers vertegenwoordigen en kolommen die tijd-frequentiebandkenmerken en demografische kenmerken zoals leeftijd, geslacht en groep weergeven. De pijplijn die in Figuur 2 wordt getoond, haalt efficiënt EEG-kenmerken uit die relevant zijn voor Alzheimeronderzoek. De geëxtraheerde kenmerken bieden inzicht in neurale activiteit binnen sleutelfrequentiebanden en kunnen worden gebruikt voor machine learning en statistische analyse.
5. Selectie van functies
Om de prestaties van het model te verbeteren en de dimensie van kenmerken te verminderen, werd een tweerichtingsmethode op basis van Analysis of Variance (ANOVA) toegepast. De ANOVA-gebaseerde selectie werd toegepast met Groep en leeftijd als onafhankelijke factoren, en RMS-frequentiebanden als afhankelijke kenmerken. Deze analyse evalueerde de effecten van groep-, leeftijds- en leeftijdsgroepinteractie op elk kenmerk. Kenmerken met p-waarden < 0,05 voor ten minste één factor werden geselecteerd voor verdere classificatie. De gedetailleerde resultaten op basis van ANOVA worden weergegeven in het resultatengedeelte. Deze kenmerken zijn de meest informatieve variabelen voor classificatie, wat hielp de nauwkeurigheid en generalisatie van het uiteindelijke voorspellende model te verbeteren. Als gevolg van de ANOVA-test werden alleen de meest informatieve functies geselecteerd voor modelcreatie. En de extra niet-numerieke attributen, Gender en participant_id, werden uitgesloten van de analyse. De doelvariabele Groep werd label-gecodeerd voor classificatie.
6. Modelbeschrijving
Om de ziekte in 3 klassen te classificeren, namelijk Alzheimer, Control en frontale dementie, werden drie begeleide machine learning-modellen gebouwd, zoals XG-Boost, Random Forest, Support Vector Machine SVM en een gestapeld model. Alle genoemde modellen gebruikten 70% van de trainingsdataset en 30% van de testdataset. Elk model werd geselecteerd vanwege zijn bewezen prestaties in de analyse van gezondheidszorggegevens en zijn vermogen om niet-lineaire en hoog-dimensionale featureruimtes te verwerken. Evaluatie van het model werd uitgevoerd op de validatiedataset om de generaliseerbaarheid van het model te beoordelen.
De modelprestaties werden geëvalueerd met behulp van een verwarringsmatrix, bestaande uit true positive (TP), false positive (FP), true negative (TN) en false negative (FN). Prestatie-indicatoren omvatten nauwkeurigheid en een classificatierapport met precisie, terugroeping en F1-scores voor elke klasse. De volgende vergelijking illustreert prestatie-metrics.
(2)
(3)
(4)
(5)
1. XG-Boost classifier
Deze methode implementeerde een multi-klasse classificatiemodel met het XGBoost-algoritme om EEG-afgeleide cognitieve toestanden te voorspellen uit geëxtraheerde kenmerken, zoals weergegeven in output 1 van Figuur 3. Output 1 bevat RMS-functies en de MMSE-score. Bij het bouwen van het XGBoost-model wordt de MMSE-score opgenomen om het niet-lineaire gedrag van de RMS-functies te testen. Het wordt gebruikt als benchmark en basis voor de systeemontwikkeling. MMSE werd in latere fasen uitgesloten om het doel van een volledig geautomatiseerd, therapeutisch onafhankelijk systeem te bereiken.
XG-Boost is een gradient boosting-framework dat geoptimaliseerd is voor efficiëntie en nauwkeurigheid, waardoor het zeer geschikt is voor het verwerken van gestructureerde data in classificatietaken. De weinige hyperparameters worden afgestemd op de EEG-dataset met technieken zoals crossvalidatie en grid search. De hieronder gedefinieerde sleutelparameterwaarden zijn gekozen om overfitting via regularisatie te beheersen en om de complexiteit van de EEG-kenmerken vast te leggen zonder te complex te zijn.
De belangrijkste parameters van het model zijn: (1) Maximale boomdiepte = 8, wat de diepte van beslissingsbomen beperkt om overfitting te verminderen en generalisatie te verbeteren. Omdat een ondiepe boom onnodige informatie zou missen, omdat een diepere boom kan overpassen. (2) L2-regularisatie (λ = 10) voegt een straf toe voor grote coëfficiënten om overfitting te minimaliseren door de modelcomplexiteit te beheersen. (2) L1-regularisatie (α = 5) wordt geïntroduceerd om extra sparsiteit in het model te geven en zo de interpreteerbaarheid en robuustheid te verbeteren. Het aantal Schatters (n_estimators=8) beperkte het aantal boosting-iteraties om de rekenefficiëntie te behouden terwijl de prestaties in balans worden gebracht.
2. Willekeurige Bosclassificator
Random Forest is een ensembleclassifier die meerdere beslissingsbomen bouwt op willekeurige deelverzamelingen van data en features, waarbij hun output wordt geaggregeerd voor de uiteindelijke voorspelling. Het biedt robuustheid tegen ruis, verwerkt niet-lineaire data goed en vermindert overfitting door het gemiddelde van diverse boomvoorspellingen te nemen. In dit werk werd het Random Forest-model geconfigureerd met sleutelparameters door gebruik te maken van de trial-and-error-methode. Dergelijke parameters worden hieronder uitgelegd:
1) n_estimators=100
Deze wordt gekozen om het aantal beslissingsbomen in het bos te beperken. Hoe meer bomen, hoe beter de generalisatie en hoe stabieler de voorspellingen.
2) max_depth=10
Het geeft de maximale diepte van elke boom aan. Als de boom een lage diepte heeft, is er een kans op onderfitting. Hoewel hoge diepte details kan vastleggen en de trainingsgegevens kan onthouden, kan het overfit veroorzaken. Daarom is het kiezen van 10 een ideale positie, vooral bij ruisachtige of complexe signalen zoals EEG.
3) random_state=40
Het fixt het willekeurige getallengeneratiezaad dat wordt gebruikt voor Bootstrap-sampling en Tree-constructie. Dit zorgt voor reproduceerbaarheid
4) n_jobs=-1
Door de waarde van n_jobs op -1 te zetten, maakt het model efficiënt gebruik van alle CPU-kernen.
7. Support Vector Machine (SVM)
Support Vector Machine is een marge-gebaseerde classifier die bekend staat om zijn lineaire prestaties in binaire en multi-klasse classificatie22. Dit model gebruikte een radiale basisfunctie-kernel vanwege het vermogen om niet-lineaire EEG-signalen vast te leggen. Het model vertoonde echter suboptimale prestaties, waarschijnlijk door het gebrek aan feature-scaling en de beperkte datasetgrootte. Hoewel SVM theoretisch krachtig is, kan de gevoeligheid voor hyperparameterafstemming en datadistributie de lagere nauwkeurigheid verklaren die wordt waargenomen in vergelijking met boomgebaseerde modellen.
8. Voorgesteld HY-fusiemodel
Hoewel het Random Forest-model een hoge nauwkeurigheid behaalde, werd om generaliseerbaarheid en schaalbaarheid te waarborgen een gestapeld model gebouwd met Random Forest en XG-Boost. De output wordt als input ingevoerd in een logistisch regressiemodel. En logistische regressie fungeert als een metalearner. Het architectuurdiagram van het model is weergegeven in Figuur 4.
Naast deeerste set invoer 20 nam het gestapelde model ook invoer van de tweede set data23. De dataset van Mendely21 is vooraf verwerkt tot een modelbegrijpelijke vorm door RMS te extraheren en te epocheren tot 10 s. De resultaten die elk model verkregen, worden weergegeven in de resultatensectie.