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, 겹침(overlap)은 0 s로 설정하였습니다.
5. 에포크 검증
각 참가자에 대해 10s의 완전한 에포크만 유지하였습니다. 기록 끝부분에 남은 10 s 미만의 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%) 하위 집합으로 분할되었습니다. 이후 각 참가자에게 속한 예측 에포크들에 대해 다수결 투표(majority voting)를 실시하여 피험자 수준의 예측값을 얻었습니다. 주요 평가는 동일한 참가자의 에포크가 두 하위 집합 모두에 나타나는 것을 방지하기 위해 피험자 수준의 그룹 훈련-테스트 분할을 사용했습니다. 피험자 수준의 예측은 각 참가자에게 속한 예측 에포크들 간의 다수결 투표를 통해 도출되었습니다. 데이터셋의 결측값을 점검하였으며, 결측 관측치는 해당되는 방식에 따라 제거하거나 대체(impute)하였습니다. 진단 레이블은 인코딩되었습니다. StandardScaler를 훈련 데이터에 맞춘(fit) 후, 이를 훈련 및 테스트 데이터셋 모두에 적용하였습니다.
8. 랜덤 포레스트 모델 개발
80개의 트리, 최대 트리 깊이 10, 최대 4개의 특성, 리프당 최소 샘플 수 5, 그리고 random state 42를 사용하여 랜덤 포레스트 분류기를 초기화했습니다. 이 분류기는 표준화된 훈련 데이터셋을 사용하여 학습되었습니다.
9. 모델 평가
훈련 및 테스트 데이터셋 모두에 대해 클래스 레이블을 예측하였습니다. 혼동 행렬(confusion matrix)을 구성하고, 분류 보고서와 함께 정확도, 정밀도, 재현율 및 F1-score를 산출하였습니다. 훈련 및 테스트 정확도를 기록하였습니다.
참가자별 평가를 위해, 특정 참가자의 모든 에포크가 훈련 세트 또는 테스트 세트 중 하나에만 배타적으로 할당되었는지 확인하십시오. Random Forest 분류기는 기본 분석과 동일한 하이퍼파라미터를 사용하여 훈련되었습니다.
10. 교차 검증
추가적인 모델 평가 절차로서, shuffle = True 및 random_state = 30을 사용하여 5-겹 층화 교차 검증(five-fold stratified cross-validation)을 수행하였다. 5개 폴드에 대한 평균 정확도와 표준 편차를 산출하였다. 이 분석은 기본 피험자 수준의 홀드아웃(holdout) 평가와는 별도로 고려되었다.
11. 특성 안정성 분석
0부터 9까지의 랜덤 시드를 사용하여 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 및 연령을 포함하는 축소된 특성 행렬을 구축하였다. 동일한 하이퍼파라미터를 사용하여 랜덤 포레스트 분류기를 재학습시켰다. 축소된 특성 세트를 사용하여 모델 학습, 테스트, 교차 검증 및 수신자 조작 특성(ROC) 분석을 반복 수행하였다.
16. ROC 분석
클래스 확률은 최적화된 랜덤 포레스트(Random Forest) 분류기에서 얻었습니다. 다중 클래스 ROC 곡선은 일대다(one-vs-all) 방식의대비-휴식 전략. 클래스별 곡선 아래 면적(AUC) 값과 평균 AUC를 산출하였다.
17. 성능 비교
전체 특징 모델의 성능을 엔트로피 특징 제거(ablation) 후 얻은 축소 특징 모델의 성능과 비교하였다.