EEG 전처리, 특징 추출 및 분류 워크플로의 개략적인 개요가 그림 1에 제시되어 있습니다. 연구에 사용된 재료 및 소프트웨어는 재료 표에 나열되어 있습니다.

그림 1: EEG 전처리, 특징 추출 및 분류 파이프라인의 워크플로. EEG 데이터 검증, 밴드패스 필터링 및 재참조, 제곱평균제곱근(RMS), 전력 스펙트럼 밀도(PSD) 및 엔트로피 특징 추출, 연령 추가, 훈련-테스트 데이터 분할, Random Forest 모델 개발, 특징 안정성 분석, SHAP 분석, Cohen’s d 효과 크기 분석, 통합 특징 선택, 엔트로피 특징 제거(ablation), 모델 재학습 및 성능 평가를 포함한 연구 워크플로의 개략도. 이 그림의 더 큰 버전을 보려면 여기를 클릭하십시오.
1. 데이터 획득
공개적으로 이용 가능한 휴지기 EEG 데이터셋을 OpenNeuro 저장소(데이터셋 ds004504, 버전 1.0.9)에서 획득하였다33. 참가자 라벨과 인구통계학적 정보를 확인하였다. 해당 데이터셋은 알츠하이머병(AD) 환자 36명, 전두측두엽 치매(FTD) 환자 23명, 건강한 대조군(HC) 29명을 포함하여 총 88명의 참가자로 구성되었다.
2. EEG 기록 가져오기
각 참가자에 대해 MNE-Python 라이브러리를 사용하여 EEG 기록을 로드하였습니다. 각 EEG 파일의 접근 가능 여부, 성공적인 가져오기 가능 여부 및 유효한 참가자 식별자 포함 여부를 확인하였습니다. 파일이 누락되었거나 손상된 기록은 제외하였습니다. 이 단계에서는 500 Hz의 원래 EEG 샘플링 주파수를 유지하였습니다. 본 연구에 임상적으로 유의미한 모든 EEG 정보를 보존하면서 계산 부하를 줄이기 위해 신호를 다운샘플링하였습니다.
3. 대역 통과 필터링
연속 EEG 기록은 저역 차단 주파수 0.5 Hz 및 고역 차단 주파수 40 Hz의 4차 대역 통과 필터를 사용하여 필터링되었습니다.
4. 공통 평균 참조 및 에포킹
필터링된 EEG 기록에 공통 평균 참조(CAR)를 적용했습니다. 각 시간 샘플에서 사용 가능한 모든 EEG 채널의 평균 신호를 계산하여 각 채널에서 뺐습니다. 참조된 EEG 신호를 이후 분석을 위해 유지했습니다. 각 참조된 EEG 기록은 MNE-Python의 make_fixed_length_epochs() 함수를 사용하여 고정 길이의 에포크로 분할되었습니다. 에포크 길이는 10 s, 중첩 시간은 0 s로 설정했습니다.
5. 에포크 검증
각 참가자에 대해, 10s의 완전한 에포크만 유지하였습니다. 기록 끝부분에 남은 10s 미만의 EEG 세그먼트는 모두 폐기하였습니다. 이렇게 유지된 각 에포크는 이후 특징 추출을 위한 개별 샘플로 사용되었습니다.
6. 특성 준비
소스 데이터셋에 Mini-Mental State Examination (MMSE) 데이터가 포함되어 있었으나, 본 분석에서는 제외하였습니다. 추출된 EEG 특성을 participants.csv 메타데이터 파일에서 얻은 연령 및 진단 레이블을 포함한 참가자 인구통계학적 정보와 병합하였습니다. 참가자 식별자를 공통 키로 사용하였으며, EEG 특성 기록과 그에 해당하는 인구통계학적 정보가 모두 있는 참가자만 유지하였습니다.
다중 도메인 특성 행렬에는 5가지 제곱 평균 제곱근(RMS) 특성(delta_rms, theta_rms, alpha_rms, beta_rms 및 gamma_rms), 5가지 전력 스펙트럼 밀도(PSD) 특성(delta_psd, theta_psd, alpha_psd, beta_psd 및 gamma_psd), 그리고 5가지 엔트로피 특성(delta_entropy, theta_entropy, alpha_entropy, beta_entropy 및 gamma_entropy)이 포함되었습니다. 연령은 보완적인 인구통계학적 변수로 포함되었습니다. 진단 그룹 간에 연령 분포가 다를 수 있으므로, 인구통계학적 혼란 변수의 영향을 완전히 배제할 수는 없었습니다. AD, HC 또는 FTD로 정의된 진단 그룹을 타겟 라벨로 지정하였습니다.
7. 데이터 전처리 및 훈련-테스트 데이터 분할
전처리 단계에서는 관심 있는 EEG 주파수 정보를 유지하면서 계산 요구 사항을 줄이기 위해, 기록된 데이터를 500 Hz에서 250 Hz로 다운샘플링했습니다.
데이터셋은 그룹 분할 전략을 사용하여 피험자 수준에서 훈련(80%) 및 테스트(20%) 하위 집합으로 분할되었습니다. 이후 각 참가자에게 속한 예측 에포크들에 대해 다수결 투표를 실시하여 피험자 수준의 예측값을 얻었습니다. 주요 평가는 동일한 참가자의 에포크가 두 하위 집합 모두에 나타나는 것을 방지하기 위해 피험자 수준의 그룹 훈련-테스트 분할을 사용하였습니다. 피험자 수준의 예측값은 각 참가자에게 속한 예측 에포크들에 대한 다수결 투표를 통해 얻었습니다. 데이터셋의 결측값을 점검하였으며, 결측 관측치는 해당되는 경우 제거하거나 대체하였습니다. 진단 레이블은 인코딩되었습니다. StandardScaler를 훈련 데이터에 맞춘(fitted) 후 훈련 및 테스트 데이터셋 모두에 적용하였습니다.
8. 랜덤 포레스트 모델 개발
80개의 트리, 최대 트리 깊이 10, 최대 4개의 특성, 리프당 최소 5개의 샘플, 그리고 random state 42를 사용하여 랜덤 포레스트(Random Forest) 분류기를 초기화했습니다. 이 분류기는 표준화된 훈련 데이터셋을 사용하여 학습되었습니다.
9. 모델 평가
훈련 데이터셋과 테스트 데이터셋 모두에 대해 클래스 레이블을 예측했습니다. 혼동 행렬을 생성하고, 분류 보고서와 함께 정확도, 정밀도, 재현율 및 F1-스코어를 계산했습니다. 훈련 및 테스트 정확도를 기록했습니다.
참가자별 평가를 위해, 특정 참가자의 모든 에포크(epoch)가 훈련 세트 또는 테스트 세트 중 하나에만 배타적으로 할당되었는지 확인하십시오. Random Forest 분류기는 기본 분석과 동일한 하이퍼파라미터를 사용하여 훈련되었습니다.
10. 교차 검증
추가적인 모델 평가 절차로서, shuffle = True 및 random_state = 30을 사용하여 5-겹 층화 교차 검증(five-fold stratified cross-validation)을 수행하였다. 5개 폴드 전체의 평균 정확도와 표준 편차를 계산하였다. 이 분석은 1차 피험자 수준 홀드아웃 평가와는 별개로 고려되었다.
11. 특성 안정성 분석
0부터 9까지의 무작위 시드(random seed)를 사용하여 Random Forest 학습을 10회 반복함으로써 특성 안정성을 평가하였다. 각 실행마다 테스트 정확도와 특성 중요도 점수를 기록하였다. 10회 실행에 걸쳐 각 예측 변수의 특성 중요도 점수의 평균과 표준 편차를 계산하였으며, 예측 변수들을 안정성에 따라 순위 매겼다. 이 분석은 일차적인 피험자 수준의 평가를 대체하기보다는 특성 순위의 일관성을 평가하기 위해 사용되었다.
12. 설명 가능한 인공지능 분석
훈련된 Random Forest 모델에 SHAP TreeExplainer를 적용하였습니다. 각 특성이 모델의 예측에 기여하는 정도를 추정하기 위해 SHAP 값을 계산하였습니다. 각 특성에 대해 평균 절대 SHAP 값을 산출하였으며, SHAP 기여도에 따라 특성의 순위를 매겼습니다. 일관되게 낮은 SHAP 값을 보이는 특성들을 식별하여 이를 특성 안정성 및 Cohen’s d 분석 결과와 비교하였습니다. 일관되게 낮은 기여도를 보이는 특성들을 선정하여 이후의 절제 연구(ablation) 및 모델 재훈련에 활용하였습니다.
13. 통계적 효과 크기 분석
AD에 대한 각 EEG 바이오마커의 Cohen's d를 산출하였다. 대비 건강 대조군, 알츠하이머병 대비 전두측엽 치매 및 전두측엽 치매 대조 건강 대조군(HC)과의 비교. 효과 크기의 정도는 소효과 0.20, 중효과 0.50, 대효과 0.80의 임계값을 사용하여 해석하였다.
14. 통합 특성 선택
Random Forest 특성 중요도 분석, SHAP 분석 및 Cohen’s d 분석 결과를 비교하였습니다. 특성 중요도가 지속적으로 낮고, SHAP 기여도가 낮으며, 효과 크기가 작은 예측 변수들을 제거 대상으로 선정하였습니다.
15. 특성 제거 분석
탐색적 제거 분석을 통해 엔트로피 특성을 평가하였으며, RMS, PSD 및 연령을 포함하는 축소된 특성 행렬을 구축하였습니다. 동일한 하이퍼파라미터를 사용하여 랜덤 포레스트(Random Forest) 분류기를 재학습시켰습니다. 축소된 특성 세트를 사용하여 모델 학습, 테스트, 교차 검증 및 수신자 조작 특성(ROC) 분석을 반복 수행하였습니다.
16. ROC 분석
최적화된 Random Forest 분류기로부터 클래스 확률을 얻었습니다. 다중 클래스 ROC 곡선은 one-대비-휴식 전략. 클래스별 곡선 아래 면적(AUC) 값과 평균 AUC를 산출하였다.
17. 성능 비교
전체 특징 모델의 성능을 엔트로피 특징 제거(entropy-feature ablation) 후 얻은 축소 특징 모델의 성능과 비교하였습니다.