Een schematisch overzicht van de workflow voor EEG-preprocessing, kenmerkextractie en classificatie wordt getoond in Figuur 1. De materialen en software die in de studie zijn gebruikt, staan vermeld in de Tabel met Materialen.

Figuur 1: Workflow van de EEG-preprocessing, kenmerkextractie en classificatie-pijplijn. Schematisch overzicht van de workflow van de studie, inclusief EEG-gegevensverificatie, banddoorlaatfiltering en her-referentie, extractie van root mean square (RMS), power spectral density (PSD) en entropie-kenmerken, toevoeging van leeftijd, train-test splitsing, ontwikkeling van het Random Forest-model, kenmerkstabiliteitsanalyse, SHAP-analyse, Cohen’s d effectgrootte-analyse, geïntegreerde kenmerkselectie, entropie-kenmerk ablatie, hertraining van het model en prestatie-evaluatie. Klik hier om een grotere versie van deze figuur te bekijken.
1. Gegevensverwerving
De openbaar toegankelijke EEG-dataset in rusttoestand is verkregen uit de OpenNeuro-repository (dataset ds004504, versie 1.0.9)33. De labels van de deelnemers en de demografische informatie zijn geverifieerd. De dataset bestond uit 88 deelnemers, waaronder 36 deelnemers met de ziekte van Alzheimer (AD), 23 met frontotemporale dementie (FTD) en 29 gezonde controles (HC).
2. Import van EEG-registraties
Voor elke deelnemer werd de EEG-registratie geladen met behulp van de MNE-Python bibliotheek. Elk EEG-bestand werd gecontroleerd om te bevestigen dat het toegankelijk was, succesvol geïmporteerd kon worden en een geldige deelnemer-identificatie bevatte. Registraties met ontbrekende of beschadigde bestanden werden uitgesloten. De oorspronkelijke EEG-bemonsteringsfrequentie van 500 Hz werd in deze fase behouden. De signalen werden gedownsampled om de rekenlast te verminderen, terwijl alle klinisch relevante EEG-informatie voor deze studie behouden bleef.
3. Banddoorlaatfiltering
De continue EEG-registraties werden gefilterd met een vierde-orde banddoorlaatfilter met een lage afsnijfrequentie van 0.5 Hz en een hoge afsnijfrequentie van 40 Hz.
4. Gemiddelde referentiemeting en epoching
Gemiddelde referentiemeting (common average referencing, CAR) werd toegepast op de gefilterde EEG-opnames. Bij elk tijdstip werd het gemiddelde signaal over alle beschikbare EEG-kanalen berekend en van elk kanaal afgetrokken. De gerefenceerde EEG-signalen werden behouden voor verdere analyse. Elke gerefenceerde EEG-opname werd gesegmenteerd in epochs van vaste lengte met behulp van de functie make_fixed_length_epochs() in MNE-Python. Er werd een epoch-duur van 10 s met 0 s overlap gebruikt.
5. Validatie van epochs
Voor elke deelnemer werden alleen volledige epochs van 10 s behouden. Elk overgebleven EEG-segment korter dan 10 s aan het einde van een opname werd verworpen. Elke behouden epoch werd vervolgens gebruikt als een individueel monster voor kenmerkextractie.
6. Voorbereiding van kenmerken
Hoewel gegevens van de Mini-Mental State Examination (MMSE) beschikbaar waren in de brondata-set, werden deze uitgesloten van de huidige analyse. De geëxtraheerde EEG-kenmerken werden samengevoegd met demografische informatie van de deelnemers, waaronder leeftijd en diagnostische labels, verkregen uit het metadata-bestand participants.csv. De deelnemer-identificatie werd gebruikt als gemeenschappelijke sleutel, en alleen deelnemers met zowel EEG-kenmerkrecords als overeenkomstige demografische informatie werden behouden.
De multidomein-kenmerkenmatrix bevatte vijf root mean square (RMS)-kenmerken (delta_rms, theta_rms, alpha_rms, beta_rms en gamma_rms), vijf power spectral density (PSD)-kenmerken (delta_psd, theta_psd, alpha_psd, beta_psd en gamma_psd) en vijf entropiekenmerken (delta_entropy, theta_entropy, alpha_entropy, beta_entropy en gamma_entropy). Leeftijd werd meegenomen als aanvullende demografische variabele. Omdat leeftijdsverdelingen kunnen variëren tussen diagnostische groepen, kon demografische confounding niet volledig worden uitgesloten. De diagnostische groep, gedefinieerd als AD, HC of FTD, werd toegewezen alsHet doellabel.
7. Preprocessing van gegevens en train-test partitionering
Tijdens de preprocessing werden de opnames gedownsampled van 500 Hz naar 250 Hz om de computationele vereisten te verminderen, terwijl de relevante EEG-frequentie-informatie behouden bleef.
De dataset werd op subjectniveau opgesplitst in trainings- (80%) en testsets (20%) met behulp van een grouped splitting-strategie. Voorspellingen op subjectniveau werden vervolgens verkregen via majority voting over de voorspelde epochs die toebehoorden aan elke deelnemer. Voor de primaire evaluatie werd een subject-level grouped train-test split gebruikt om te voorkomen dat epochs van dezelfde deelnemer in beide subsets voorkwamen. Voorspellingen op subjectniveau werden verkregen via majority voting over de voorspelde epochs die toebehoorden aan elke deelnemer. De dataset werd gecontroleerd op ontbrekende waarden, waarbij ontbrekende observaties werden verwijderd of, indien van toepassing, geïmputeerd. Diagnostische labels werden gecodeerd. Een StandardScaler werd gefit op de trainingsdata en vervolgens toegepast op zowel de trainings- als de testdatasets.
8. Ontwikkeling van het Random Forest-model
Een Random Forest-classifier werd geïnitialiseerd met 80 bomen, een maximale boomdiepte van 10, maximaal vier kenmerken, minimaal vijf monsters per blad en een random state van 42. De classifier werd getraind met behulp van de gestandaardiseerde trainingsdataset.
9. Model-evaluatie
Klasselabels werden voorspeld voor zowel de trainings- als de testdatasets. Er werd een confusiematrix opgesteld, en de accuratesse, precisie, recall en F1-score werden berekend samen met het classificatierapport. De accuratesse van de training en het testen werd genoteerd.
Zorg er bij de evaluatie per deelnemer voor dat alle epochs van een bepaalde deelnemer uitsluitend werden toegewezen aan óf de trainings- óf de testsubset. De Random Forest-classifier werd getraind met dezelfde hyperparameters als in de primaire analyse.
10. Kruisvalidatie
Als aanvullende procedure voor model-evaluatie werd een gestratificeerde vijf-voudige kruisvalidatie uitgevoerd met shuffle = True en random_state = 30. De gemiddelde nauwkeurigheid en standaarddeviatie over de vijf folds werden berekend. Deze analyse werd apart beschouwd van de primaire holdout-evaluatie op subjectniveau.
11. Analyse van de stabiliteit van kenmerken
De stabiliteit van de kenmerken werd beoordeeld door de training van het Random Forest 10 keer te herhalen met willekeurige seeds van 0 tot en met 9. Voor elke run werden de testnauwkeurigheid en de scores voor het belang van de kenmerken geregistreerd. Het gemiddelde en de standaarddeviatie van de score voor het belang van de kenmerken voor elke voorspeller werden berekend over de 10 runs, en de voorspellers werden gerangschikt op basis van hun stabiliteit. Deze analyse werd gebruikt om de consistentie van de rangschikking van de kenmerken te beoordelen, in plaats van om de primaire evaluatie op subjectniveau te vervangen.
12. Analyse met uitlegbare kunstmatige intelligentie
SHAP TreeExplainer werd toegepast op het getrainde Random Forest-model. SHAP-waarden werden berekend om de bijdrage van elk kenmerk aan de voorspellingen van het model te schatten. Voor elk kenmerk werd de gemiddelde absolute SHAP-waarde berekend, waarna de kenmerken werden gerangschikt op basis van hun SHAP-bijdragen. Kenmerken met consistent lage SHAP-waarden werden geïdentificeerd en vergeleken met de resultaten van de analyses naar kenmerkstabiliteit en de d van Cohen. Kenmerken die consistent lage bijdragen vertoonden, werden geselecteerd voor daaropvolgende ablatie en hertraining van het model.
13. Statistische analyse van de effectgrootte
Cohen’s d werd berekend voor elke EEG-biomarker voor de vergelijkingen AD versus HC, AD versus FTD en FTD versus HC. De grootte van het effect werd geïnterpreteerd met drempelwaarden van 0,20 voor een klein effect, 0,50 voor een gemiddeld effect en 0,80 voor een groot effect.
14. Geïntegreerde kenmerkselectie
De resultaten van de Random Forest-analyse van kenmerkbelang, de SHAP-analyse en de Cohen’s d-analyse werden vergeleken. Predictoren die consistent een laag kenmerkbelang, een lage SHAP-bijdrage en kleine effectgroottes vertoonden, werden geïdentificeerd als kandidaten voor verwijdering.
15. Ablatie van kenmerken
Entropiekenmerken werden geëvalueerd in een exploratieve ablatieanalyse, waarna een gereduceerde kenmerkenmatrix bestaande uit RMS, PSD en leeftijd werd samengesteld. De Random Forest-classifier werd opnieuw getraind met behulp van dezelfde hyperparameters. De modeltraining, het testen, de kruisvalidatie en de receiver operating characteristic (ROC)-analyse werden herhaald met gebruik van de gereduceerde kenmerkenset.
16. ROC-analyse
Klassewaarschijnlijkheden werden verkregen uit de geoptimaliseerde Random Forest-classifier. Multiklasse ROC-curven werden gegenereerd met behulp van een one-versus-rest-strategie. De klasse-specifieke waarden voor het oppervlak onder de curve (AUC) en de gemiddelde AUC werden berekend.
17. Vergelijking van de prestaties
De prestaties van het volledige kenmerkenmodel werden vergeleken met die van het gereduceerde kenmerkenmodel dat was verkregen na entropie-kenmerkenablatie.