本研究では、アルツハイマー病、前頭側頭型認知症、および健常対照群を区別するための脳波バイオマーカーについて調査しています。二乗平均平方根、パワースペクトル密度、およびエントロピーの特徴量を用いたランダムフォレストモデルを、説明可能な人工知能(XAI)および効果量分析を用いて解析しました。影響の少ないエントロピー特徴量を除去することで分類性能が向上し、テスト精度94%を達成しました。
本研究では、アルツハイマー病、前頭側頭型認知症、および健常対照群を区別するための脳波バイオマーカーについて調査しています。二乗平均平方根、パワースペクトル密度、およびエントロピーの特徴量を用いたランダムフォレストモデルを、説明可能な人工知能(XAI)および効果量分析を用いて解析しました。影響の少ないエントロピー特徴量を除去することで分類性能が向上し、テスト精度94%を達成しました。
脳波検査(EEG)は、神経変性疾患における脳活動を特性づけるための非侵襲的なアプローチとして研究されてきました。本研究では、マルチドメインEEGバイオマーカーが、アルツハイマー病(AD)、前頭側頭型認知症(FTD)、および健常対照群(HC)を区別できるかどうかを評価しました。公開データセットには、AD患者36名、FTD患者23名、および健常対照群29名の計88名から得られた、閉眼安静時のEEG記録が含まれていました。機械学習による分類のために、マルチドメインEEGバイオマーカーと人口統計学的特徴が抽出されました。EEG由来の特徴量には、二乗平均平方根(RMS)、パワースペクトル密度(PSD)、およびエントロピー指標が含まれています。統合された特徴量セットを用いてランダムフォレスト分類器を学習させ、説明可能な人工知能を用いて解釈を行いました。特徴量の寄与度を評価するために、特徴量の安定性とCohen's d効果量の分析を行い、続いてアブレーション分析を実施しました。その結果、エントロピーベースの特徴量が分類への寄与度が最も低いことが示されました。完全なマルチドメインモデルでは、学習精度96%、テスト精度92%を達成しました。エントロピー特徴量を削除した後、テスト精度は94%に向上しました。これらの結果は、ランダムフォレストの枠組みにおける、選択されたマルチドメインEEGおよび人口統計学的特徴の分類性能を記述したものです。
アルツハイマー病(AD)は、認知機能および運動機能の低下を伴う進行性の神経疾患です。その有病率は加齢とともに大幅に増加し、認知症は高齢者の罹患率および死亡率の主要な原因となっています1,2,3,4,5。診断および疾患分類は依然として困難であるため、認知症に伴う神経生理学的変化を特徴付ける手段として、脳信号の解析が研究されています。
脳波検査(EEG)は、脳の電気活動を記録するための非侵襲的な手法であり、アルツハイマー病およびその他の神経疾患において研究されてきた6,7,8,9,10,11。EEG信号はニューロン活動によって生成される電位を反映しており、頭皮上の複数の部位から非侵襲的に記録される12,13,14。脳活動は一般的に、delta、theta、alpha、beta、およびgammaの周波数帯域を用いて特徴づけられる15。アルツハイマー病は、EEGの徐波化、信号複雑性の低下、および同期性の変化に関連している16。スペクトル解析では、健常対照群と比較して、AD患者においてdeltaおよびtheta活動の増加と、alphaおよびbeta活動の減少が報告されている17,18,19。これまでのEEG分類研究では、ADと健常対照群(HC)、および前頭側頭型認知症(FTD)とHCを区別する高い性能が報告されている。対照的に、ADとFTDの識別は一般的により困難である10,11。
機械学習の手法は、神経疾患に関連するEEG信号のパターンを特定するために用いられてきた20,21。したがって、特徴抽出はEEGに基づいた分類における重要な構成要素となる。先行研究では、スペクトル、時間、複雑性、および機能的結合性の特徴が、機械学習モデルと組み合わせて調査されてきた3,15,22,23。スペクトル特徴は、ADで報告されているアルファ活動の低下やシータ活動の増加など、脳活動における周波数依存性の変化を定量化できるため、特に重要である23,24,25,26。
EEGベースの複数の研究において、特徴量、分類器、および検証戦略のさまざまな組み合わせが評価されてきました。コネクティビティに基づく解析はAD、FTD、およびHCを区別するために用いられてきましたが、報告されている分類性能は、比較対象や特徴量選択の手順によって異なります27。また、フーリエ変換ベースの特徴量についても、AD、軽度認知障害(MCI)、および健康対照群を区別するために評価されていますが、前処理の不足や、検証戦略に関連した過学習の可能性などが限界として挙げられています28。非線形EEG複雑度指標は、AD分類において識別能を示すことが明らかになっており、一部では高い精度、感度、特異度が達成されていますが、MCIへの適用や疾患ステージの識別における有用性は依然として限定的です21。
時間-周波数EEG表現に基づくディープラーニング手法も、AD、MCI、およびHCの多クラス分類に向けて研究されています29。高い分類性能が報告されていますが、小規模または限定的なコホート、クラス不均衡、疾患重症度の不均一性、外部検証の不足、およびモデルの解釈性の低下などの限界があります29。検証戦略も重要な方法論的検討事項です。参加者独立アプローチにより、同一参加者のエポックがトレーニングセットとテストセットの両方に含まれる場合、従来のエポックレベルの交差検証では性能を過大評価する可能性があることが示されています30。この限界は、異なる検証スキームを使用している研究間で分類結果を比較する場合に特に重要となります。
複雑性と接続性の指標がサポートベクターマシン分類器を用いて評価されており、位相ベースの接続性とマルチスケールエントロピーがさまざまなレベルの判別性能を示している31,32。フーリエ変換およびウェーブレット変換に基づく特徴量を用いた別の研究では、AD、MCI、およびHCの比較において性能にばらつきがあることが報告されており、一般的に多クラス分類は二値分類よりも困難なままである14。さらに、エントロピーおよび同期に基づくアプローチでもADと対照群の判別が実証されているが、対照群の特性が観察される効果量に影響を与える可能性がある2。
本研究では、人口統計学的変数としての年齢に加え、EEG周波数帯域から抽出した平方根平均二乗(RMS)、パワースペクトル密度(PSD)、およびエントロピー指標からなるマルチドメインの特徴量セットを評価します。RMSは各周波数帯域内の信号振幅を表し、PSDは周波数領域の電力を表します。エントロピーは信号の複雑性の尺度として導入されています。これらの特徴量はEEG活動の相補的な記述を提供し、ランダムフォレスト分類器を用いて評価されます。さらに、抽出した特徴量の相対的な寄与度を評価するため、特徴量の重要度、特徴量の安定性、説明可能な人工知能(XAI)、および効果量分析について検討します。本研究では、機械学習を用いて、スペクトル、時間、および複雑性に基づく指標を含むマルチドメインEEG特徴量が、アルツハイマー病(AD)、前頭側頭型認知症(FTD)、および健常対照群(HC)を識別する上でどの程度有効であるかを調査します。また、PSD、RMS、およびエントロピー指標を用いた際に、どのEEG周波数帯域の特徴量が分類に最も強く寄与するかを検討します。さらに、PSDやRMSを含む電力ベースの特徴量と、エントロピーベースの特徴量の寄与を比較し、分類モデルに対する年齢の寄与を評価します。
EEGの前処理、特徴抽出、および分類ワークフローの概略図を図1に示します。本研究で使用した材料およびソフトウェアは、材料表に記載されています。

図 1: EEGの前処理、特徴抽出、および分類パイプラインのワークフロー。EEGデータの検証、バンドパスフィルタリングとリリファレンス、二乗平均平方根(RMS)、パワースペクトル密度(PSD)、およびエントロピー特徴の抽出、年齢の追加、学習・テストデータの分割、ランダムフォレストモデルの開発、特徴安定性分析、SHAP分析、Cohen’s d 効果量分析、統合的な特徴選択、エントロピー特徴のアブレーション、モデルの再学習、および性能評価を含む、本研究のワークフローの概略図。こちらのリンクをクリックして、この図の拡大版を表示してください。
1. データ収集
公開されている安静時EEGデータセットをOpenNeuroリポジトリ(データセットds004504、バージョン1.0.9)から取得した33。被験者のラベルおよび人口統計学的情報を確認した。本データセットは、アルツハイマー病(AD)患者36名、前頭側頭型認知症(FTD)患者23名、および健常対照群(HC)29名の計88名の被験者で構成されていた。
2. EEG記録のインポート
各参加者のEEG記録をMNE-Pythonライブラリを用いてロードした。各EEGファイルについて、アクセス可能であること、正常にインポートできること、および有効な参加者識別子が含まれていることを確認した。ファイルが欠落しているか破損している記録は除外した。この段階では、元のEEGサンプリング周波数である500 Hzを維持した。その後、本研究において臨床的に関連のあるすべての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メタデータファイルから取得した年齢や診断ラベルを含む参加者のデモグラフィック情報と統合した。参加者IDを共通キーとして使用し、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を適合させ、その後、トレーニングデータセットとテストデータセットの両方に適用しました。
8. ランダムフォレストモデルの開発
決定木数80本、最大木の深さ10、最大特徴量数4、葉あたりの最小サンプル数5、および乱数シード(random state)42としてランダムフォレスト分類器を初期化した。この分類器を、標準化したトレーニングデータセットを用いて学習させた。
9. モデルの評価
トレーニングセットおよびテストセットの両方に対してクラスラベルの予測を行った。混同行列を作成し、分類レポートとともに、正解率(accuracy)、適合率(precision)、再現率(recall)、およびF1スコアを算出した。トレーニングおよびテストの正解率を記録した。
参加者ごとの評価では、特定の参加者のすべてのエポックが、トレーニングセットまたはテストセットのいずれか一方にのみ排他的に割り当てられていることを確認してください。ランダムフォレスト分類器は、主要分析と同じハイパーパラメータを用いてトレーニングされました。
10. 交差検証
追加のモデル評価手順として、shuffle = Trueおよびrandom_state = 30を用いて、5分割層化交差検証を実施した。5つのフォールドにおける平均正解率と標準偏差を算出した。この解析は、主要な被験者レベルのホールドアウト評価とは別に検討した。
11. 特徴量の安定性解析
0から9までの乱数シードを用いてRandom Forestのトレーニングを10回繰り返し、特徴量の安定性を評価した。各ランについて、テスト精度と特徴量重要度スコアを記録した。10回のランにおける各予測因子の特徴量重要度スコアの平均値と標準偏差を算出し、その安定性に基づいて予測因子をランク付けした。この解析は、主要な被験者レベルの評価に代わるものではなく、特徴量のランキングの整合性を評価するために用いられた。
12. 説明可能な人工知能(XAI)解析
学習済みランダムフォレストモデルにSHAP TreeExplainerを適用した。各特徴量がモデルの予測に与える寄与度を推定するため、SHAP値を算出した。各特徴量について平均絶対SHAP値を算出し、SHAP寄与度に基づいて特徴量の順位付けを行った。一貫してSHAP値が低い特徴量を特定し、特徴量の安定性およびCohen's d解析の結果と比較した。寄与度が一貫して低い特徴量を、その後のアブレーション解析およびモデルの再学習の対象として選択した。
13. 統計学的効果量の分析
ADに対する各EEGバイオマーカーについて、Cohen's dを算出した。 対比 健常対照群, アルツハイマー病 対照 前頭側頭型認知症、および前頭側頭型認知症 対比 HCとの比較。効果量の大きさは、小効果を0.20、中効果を0.50、大効果を0.80とする閾値を用いて解釈した。
14. 統合的な特徴量選択
ランダムフォレストの特徴量重要度解析、SHAP解析、およびCohen’s d解析の結果を比較した。特徴量重要度が低く、SHAP寄与度が低く、かつ効果量が小さいことが一貫して認められた予測因子を、削除の候補として特定した。
15. 特徴量のアブレーション
エントロピー特徴量を探索的なアブレーション解析で評価し、RMS、PSD、および年齢を含む削減特徴量行列を構築した。同一のハイパーパラメータを用いて、ランダムフォレスト分類器を再学習させた。削減された特徴量セットを用いて、モデルの学習、テスト、交差検証、および受信者動作特性(ROC)解析を繰り返した。
16. ROC解析
クラス確率は、最適化されたランダムフォレスト分類器から得られた。マルチクラスROC曲線は、1つのー対比-休止戦略。クラス別の曲線下面積(AUC)値および平均AUCを算出した。
17. 性能比較
完全な特徴量モデルの性能を、エントロピー特徴量アブレーション後に得られた削減特徴量モデルの性能と比較した。
Random Forestモデルは、トレーニング精度96%、テスト精度92%を達成した。マクロ平均F1スコアは92%であり、加重平均F1スコアも同様に92%であったことから、3つのクラス全体で同等の性能を示した。分類性能はアルツハイマー型認知症(AD)クラスで最も高く、一方で前頭側頭型認知症(FTD)クラスは相対的に低い性能を示した(表1に提示)。
分類しきい値におけるクラス分離を評価するために、受信者動作特性(ROC)分析を用いた。図2に示すように、3つのクラスすべてでROC曲線下面積(AUC)が0.98を超えており、AD、FTD、および健康対照(HC)群の間で高い分離能があることが示された。
Random Forest分類モデルに対する特徴量の寄与度を評価するため、TreeExplainerを用いたSHAP解析を実施した。図3の要約バープロットは、電力スペクトル密度(PSD)、二乗平均平方根(RMS)、エントロピー、および人口統計学的特徴量の相対的な寄与度を示している。SHAP値の絶対値が高い特徴量ほどモデルの予測への寄与が強く、一方でエントロピー関連の特徴量の寄与は比較的低かった。
指定されたランダムシードを用いてRandom Forestのトレーニングを繰り返し、特徴量の安定性解析を行った。各ランにおいてテスト精度と特徴量重要度スコアを記録し、特徴量重要度の平均値と標準偏差を算出した。この解析において、モデルは平均分類精度95%(標準偏差0.17%)を達成した。これらの繰り返しランは、主に特徴量ランキングの一貫性を評価するために用いられ、主要な被験者レベルのホールドアウト評価とは別に検討された。表2に示す通り、特徴量安定性解析ではAgeが高くランク付けされた。また、theta_psdおよびtheta_rmsの特徴量も高いランクにあり、一方でalpha_psd、alpha_rms、beta_psd、およびbeta_rmsは中程度の特徴量重要度を示した。
群間差の大きさを評価するために、Cohen's dを算出し、次のように示した。 表3アルファ波の減衰(alpha_rms, d = 0.97; alpha_psd, d = 0.87)およびシータ波活動の増加(theta_rms, d = 0.72)に、AD患者と健康対照者の間で報告された最大の違いが見られた。これらの知見は、認知症のEEG研究で報告されている皮質緩徐化(cortical slowing)と一致していた。34エントロピーに基づく特徴量は、比較的小さい効果量を示しました。ADとFTDの比較では、主に小から中程度の効果量が得られました。
AD群において、アルファ波帯域の特徴量は大きな効果量を示した。対-HCおよびFTD-対比-健常対照群(HC)との比較。AD(アルツハイマー病)について 対比 HCにおいて、報告された効果量はalpha_psdで−0.87、alpha_rmsで−0.97であった。FTDについては 対比 HCにおいて、アルファ帯域の減衰は−0.73であった。theta_psdという特徴量は、ADに対して0.61の効果量を示した。 対比 HCであり、一方theta_rmsの効果量は0.72であった。ADではより小さな差が観察された。 対 前頭側頭型認知症(FTD)
ベータ波帯域の特徴は、主に小〜中程度の群間効果量を示しました。beta_psdの効果量は、AD versus HCで−0.35、AD versus FTDで−0.41であり、一方でbeta_rmsは主に小さな効果を示しました。ガンマ波帯域の特徴は、一部の比較において中程度の効果を示しました。gamma_psdの効果量は、FTD versus HCで0.47、AD versus HCで0.20であり、一方でgamma_rmsはFTD versus HCで0.50の効果量を示しました。
Delta PSDおよびRMSは、ADに対して小さい正の効果量を示した。 対比 HC(それぞれ0.29および0.34)およびFTD 対比 HC(それぞれ0.25および0.25)。ADの効果量 対 FTDはほぼゼロであった。デルタエントロピーは、3つの比較すべてにおいて小さな効果量(0.08、0.20、および−0.11)を示した。
特徴量の安定性、SHAP、および効果量の解析により、エントロピー特徴量の寄与度は比較的低いことが示された。その結果、エントロピー特徴量を除外し、RMS、PSD、および年齢を保持したアブレーション解析を実施した。特徴量アブレーション後、ランダムフォレストモデルは訓練精度96%、テスト精度94%を達成した。この結果は、評価した特徴量選択フレームワーク内における探索的な比較として解釈された。AD、HC、およびFTDクラスの適合率、再現率、F1スコア、およびサポートを表4に示す。
エントロピー特徴量のアブレーション前後の性能を表5に比較して示す。エントロピーを除去した後、テスト精度は92%から94%に向上した。クラスごとの変化は表6に提示している。FTDについては、適合率が4%向上し、F1スコアが2%向上した。HC群の再現率は4%向上したが、適合率は2%低下した。特徴量アブレーション後、マクロ平均および加重平均の性能指標は約1%~2%向上した。

図2アルツハイマー病、前頭側頭型認知症、および健常対照群の受信者動作特性曲線。 1対多(one-vs-all)方式を用いて作成された多クラス受信者動作特性(ROC)曲線対アルツハイマー病(クラスA)、前頭側頭型認知症(クラスF)、および健康対照群(クラスC)に対する静止戦略。各クラスについて曲線下面積(AUC)を算出した。 この図の拡大版を表示するには、こちらをクリックしてください。

図3学習済みランダムフォレストモデルのSHAPサマリー解析。 Random Forest分類モデルに対する、EEG由来の特徴量および人口統計学的特徴量の相対的な寄与度を示す要約棒グラフ。二乗平均平方根(RMS)、パワースペクトル密度(PSD)、エントロピー、および年齢の特徴量の寄与度を順位付けするために、平均絶対SHAP値を用いた。 こちらの図の拡大版を表示するには、ここをクリックしてください。
| 適合率 | 再現率 | F1スコア | サポート数 | |
| アルツハイマー病 | 0.92 | 0.95 | 0.94 | 579 |
| コントロール | 0.93 | 0.92 | 0.93 | 480 |
| 前頭側頭型認知症 | 0.93 | 0.87 | 0.9 | 314 |
| マクロ平均 | 0.93 | 0.92 | 0.92 | 1373 |
| 加重平均 | 0.93 | 0.92 | 0.92 | 1373 |
表1:テストデータの分類レポート。テストデータセットにおけるアルツハイマー病、健康対照群、前頭側頭型認知症の各クラスの適合率(Precision)、再現率(Recall)、F1スコア(F1-score)、およびサポート(Support)を、マクロ平均および加重平均の性能指標とともに示す。
| 特徴量 | 平均値 | 標準偏差 |
| Age | 0.1966 | 0.0017 |
| theta_psd | 0.0936 | 0.0027 |
| theta_rms | 0.0911 | 0.0023 |
| alpha_psd | 0.0824 | 0.0022 |
| alpha_rms | 0.0824 | 0.0024 |
| beta_psd | 0.0638 | 0.0013 |
| beta_rms | 0.0629 | 0.0012 |
| gamma_psd | 0.0534 | 0.001 |
| gamma_rms | 0.0527 | 0.0014 |
| delta_psd | 0.0403 | 0.0007 |
| delta_rms | 0.0359 | 0.001 |
| gamma_entropy | 0.0326 | 0.0008 |
| alpha_entropy | 0.0325 | 0.0006 |
| beta_entropy | 0.0283 | 0.0009 |
| theta_entropy | 0.0258 | 0.0004 |
| delta_entropy | 0.0255 | 0.0005 |
表 2: Random Forest モデルの特徴量の安定性ランキング。0から9までのランダムシードを用いて Random Forest の学習を繰り返すことで得られた、特徴量の重要度の平均値および対応する変動性。特徴量は、実行回数全体における安定性に基づいてランク付けされた。
| 特徴 | アルツハイマー病 対 HC | アルツハイマー病 対 前頭側頭型認知症 | 前頭側頭型認知症 対 HC |
| デルタPSD | 0.293 | -0.069 | 0.253 |
| シータ波PSD | 0.612 | 0.23 | 0.355 |
| α-PSD | -0.869 | -0.219 | -0.64 |
| ベータPSD | -0.352 | -0.408 | 0.201 |
| ガンマPSD | 0.201 | -0.326 | 0.468 |
| デルタエントロピー | 0.081 | 0.199 | -0.111 |
| シータエントロピー | 0.208 | 0.033 | 0.163 |
| アルファエントロピー | -0.061 | 0.071 | -0.137 |
| ベータエントロピー | -0.37 | -0.042 | -0.332 |
| ガンマエントロピー | -0.338 | 0.047 | -0.391 |
| デルタRMS | 0.342 | 0.041 | 0.251 |
| シータ波RMS | 0.718 | 0.354 | 0.318 |
| Alpha RMS | -0.977 | -0.205 | -0.73 |
| ベータRMS | -0.494 | -0.365 | 0.015 |
| ガンマRMS | 0.187 | -0.33 | 0.5 |
表3:診断グループ間におけるCohenのd効果量の比較。 アルツハイマー病におけるEEGバイオマーカーのCohen's d値 対 健常対照群、アルツハイマー病 対比 前頭側頭型認知症、および前頭側頭型認知症 対比 健常対照群との比較。
| 適合率(Precision) | 再現率(Recall) | F1スコア | サポート数 | |
| アルツハイマー病 | 0.94 | 0.95 | 0.94 | 579 |
| 対照群 | 0.91 | 0.96 | 0.94 | 480 |
| 前頭側頭型認知症 | 0.97 | 0.89 | 0.92 | 314 |
| マクロ平均 | 0.94 | 0.93 | 0.93 | 1373 |
| 加重平均 | 0.94 | 0.94 | 0.94 | 1373 |
表4:特徴量アブレーション後の分類レポート。エントロピー特徴量の削除およびRandom Forestモデルの再学習後における、アルツハイマー病、健常対照群、および前頭側頭型認知症クラスの適合率、再現率、F1スコア、およびサポート。
| 指標 | アブレーション前 | アブレーション後 | 変更 |
| 学習精度 | 0.96 | 0.96 | 0%の変化 |
| 精度の検証 | 0.92 | 0.94 | 2%の向上 |
| マクロ精度 | 0.93 | 0.94 | 1%の改善 |
| マクロ想起 | 0.92 | 0.93 | 1%の改善 |
| マクロF1スコア | 0.92 | 0.93 | 1%の向上 |
| 加重適合率 | 0.93 | 0.94 | 1%の向上 |
| 重み付き再現率 | 0.92 | 0.94 | 2%の向上 |
| 加重F1スコア | 0.92 | 0.94 | 2%の向上 |
表 5:特徴量アブレーション前後の性能比較エントロピー特徴量を削除した後の削減特徴量セットと、完全なマルチドメイン特徴量セットにおけるモデル性能指標の比較。
データの利用可能性:
モデルの開発および分析に使用された処理済みCSVデータセットは、付録ファイル1として提供されています。前処理コードは、Pythonノートブック(https://colab.research.google.com/drive/1PwidMZCbbIHdthDhat5OEZa7ya6f0-aC?usp=sharing)で利用可能であり、モデル構築および分析コードは、https://colab.research.google.com/drive/1q8wYJjQpofAqZGYxz9xHKB0emGF4YGwM?usp=sharing で利用可能です。また、提供されている前処理コードを使用して、処理済みCSVデータセットを再生成することも可能です。
| クラス | 指標 | 前 | 後 | 変更 |
| アルツハイマー病 | 精度 | 0.92 | 0.94 | 0.02 |
| 再現率 | 0.95 | 0.95 | 変更なし | |
| F1スコア | 0.94 | 0.94 | 変更なし | |
| 対照 | 精度 | 0.93 | 0.91 | −2% |
| 再現率 | 0.92 | 0.96 | 0.04 | |
| F1スコア | 0.93 | 0.94 | 0.01 | |
| 前頭側頭型認知症 | 精度 | 0.93 | 0.97 | 0.04 |
| 再現率 | 0.87 | 0.89 | 0.02 | |
| F1スコア | 0.9 | 0.92 | 0.02 |
表6:特徴量アブレーション前後のクラス別性能変化。エントロピー特徴量を削除した後の、アルツハイマー病、健康対照群、および前頭側頭型認知症クラスにおける適合率、再現率、およびF1スコアの変化。
補足ファイル 1:モデル開発および解析に使用された処理済みエポックレベル EEG 特徴データセット。こちらをクリックしてファイルをダウンロードしてください。
ランダムフォレスト分類器を用い、電力スペクトル密度(PSD)、二乗平均平方根(RMS)、およびエントロピー指標を含む多ドメイン特徴量を使用して、脳波(EEG)サンプルをアルツハイマー病(AD)、前頭側頭型認知症(FTD)、および健常対照(HC)群に分類した。補完的な人口統計学的変数として年齢を組み込んだ。診断グループ間で年齢分布が異なる可能性があるため、人口統計学的な混同変数の影響を完全に排除することはできず、報告された性能は人口統計学的情報と電気生理学的情報の複合的な寄与を反映している。本解析では、多ドメインEEG特徴量がAD、FTD、およびHC群を識別できる能力を検証し、デルタ波、シータ波、アルファ波、ベータ波、およびガンマ波帯域の特徴量が分類にどのように寄与しているかを評価した35。特徴量の安定性と効果量の解析により、PSDやRMSを含む電力ベースの特徴量が、エントロピーベースの特徴量よりも強く寄与していることが示された。抽出されたEEGバイオマーカーは、リズムの徐波化およびアルファ活動の低下に関連するパターンも反映していた。
先行研究では、異なる機械学習アプローチを用いて、エントロピー、PSD、およびRMSに基づくEEG特徴量が評価されています2,14,27,28,36,37。エントロピーはEEGの複雑性と不規則性を特徴づけるために用いられ、PSDはスペクトル特性を定量化するために用いられてきました2,14,27,28,36。また、RMSは周波数帯域の特徴量として評価されています37。同一の基礎データセットを用いた先行研究では、1被験者抜き交差検証によるPSDベースの二値分類が適用されていました33。それに対し、本解析では、マルチクラス分類フレームワーク内において、アルファ波、ベータ波、ガンマ波、シータ波、およびデルタ波の各帯域におけるRMS、エントロピー、PSDを組み合わせました。同じデータセットを用い、1被験者抜き二値分類戦略で評価された軽量EEGトランスフォーマーでは、83.28%の精度が報告されています。分類タスクおよび検証手順が異なるため、本結果との直接的な数値比較は慎重に解釈されるべきです38
ADとHC、およびFTDとHCを区別する上で、アルファ波帯の特徴が最大の効果量を示し、またシータ波活動の増加もADとHCの識別に寄与しました。ベータ波帯の特徴はグループ間での効果が小さかった一方、エントロピーベースの特徴は識別能が比較的限定的でした。ガンマ波帯の特徴は、FTD versus HCを含むいくつかの比較において中程度の効果を示しました。先行研究では、ADおよび関連する診断比較において、アルファパワーの低下とシータ活動の変化が報告されています23,25,39。エントロピーについても、神経変性疾患における潜在的に有用なバイオマーカーとして記述されており36、その寄与はデータセット、特徴定義、および解析アプローチによって異なる可能性が示唆されています。
これまでのEEG分類研究では、サポートベクターマシンが頻繁に使用されてきた3,21,24,37。本解析で使用したランダムフォレスト分類器は、特徴量の重要度の推定が可能であり、本研究で用いた構造化特徴量セットに適用可能であった40,41。エントロピー特徴量の除去(アブレーション)後、学習精度96%およびテスト精度94%が得られた。エントロピー特徴量の除去によってモデルの性能が向上したことは、安定性、SHAP、および効果量解析においてこれらの特徴量の相対的な寄与度が低いことが観察された結果を支持している。ディープラーニングやトランスフォーマーベースのアプローチもEEG分類に適用されているが、こうしたアプローチは一般的により大規模なデータセットを必要とし、より高次元または系列的な表現向けに設計されている27,39,42,43,44,45。
EEG信号はノイズが多く非定常であるため、構造化された特徴表現を分析する際には、アンサンブル決定木ベースのモデルが堅牢性を提供できます。本研究では、Random Forestモデルへの相補的な入力として、RMS、PSD、エントロピー、および年齢を評価しました。RMSは周波数帯域内の信号振幅を、PSDはスペクトル密度を、エントロピーは信号の複雑性を表しています。解析の結果、評価したフレームワーク内では、PSDとRMSがエントロピーよりも分類に強く寄与することが示されました。年齢もモデルに寄与しましたが、その導入により人口統計学的な混同変数が生じる可能性がありました。今後の研究では、汎用性と分類性能をさらに評価するために、より大規模なデータセットを用いて、追加の人口統計学的変数や、代替となる機械学習またはディープラーニングの手法を検討することが考えられます。
本研究で報告された内容に影響を及ぼす可能性のある、既知の競合する利益、個人的関係、金銭的関係、またはその他の活動は申告されていません。
本研究で使用した公開EEGデータセット、および研究目的で当該データセットを利用可能にした研究者の貢献に感謝いたします。公的機関、営利企業、または非営利団体からの特定の資金提供は受けていません。
| 名前 | 会社 | カタログ番号 | コメント |
|---|---|---|---|
| Google Colab | [環境の確認]Colabランタイム | 解析ノートブックの実行に使用されるクラウドベースの計算環境 | |
| MNE-Python | MNE-Pythonプロジェクト | 1.12.1 | EEGの前処理、リファレンスの設定、フィルタリング、およびエポック分割に使用されます。 |
| NumPy | NumPyプロジェクト | (2.0.2) | 数値計算および配列ベースの計算に使用されます |
| OpenNeuro EEGデータセット | OpenNeuro | ds004504, (1.0.9) | AD、FTD、およびHCの記録を含む、公開利用可能な安静時EEGデータセット |
| Python | Python Software Foundation | [3.12.13 (main, Mar 4 2026, 09:23:07) [GCC 11.4.0]] | データ処理およびモデル開発に使用したプログラミング環境 |
| scikit-learn | scikit-learnプロジェクト | 1.6.1 | 前処理、グループ化データ分割、ランダムフォレスト分類、および性能評価に使用される |
| SciPy | SciPyプロジェクト | (1.16.3) | 信号処理および統計計算に使用される |
| SHAP(Shapley Additive Explanations) | SHAPプロジェクト | 0.52.0 | TreeExplainerを用いた説明可能な人工知能分析に使用 |