パワースペクトル密度（PSD）とは、EEG（脳波）信号を分解し、それぞれの速度（周波数）が全体の記録にどれだけのエネルギーをもたらしているかを教えてくれるツールです。PSDプロットが読めるようになれば、脳のリズムの楽譜を読んでいるようなものであり、どのテンポが優勢で、どのテンポが背景に消えていくかを示すチャートを読み解くことができます。
EEGにおけるパワースペクトル密度（PSD）とは？
生のEEG波形は、各楽器が異なる周波数の脳波を奏でるオーケストラによる和音のようなものだと考えてみてください。耳で聴くだけでは、それぞれの楽器がどれだけ寄与しているかを判断することはできません。
PSDは、各周波数における信号の振幅の2乗を計算し、ある記録時間の長さにわたってその測定値を平均化することでこの問題を解決します。その結果得られるのは、一方の軸に周波数、もう一方の軸にパワーをとった曲線であり、脳の電気活動が遅いリズムと速いリズムにどのように分布しているかを正確に示します。
PSDは記録期間全体を1本の曲線に圧縮するため、ブレンドが秒単位でどのように変化するかを追跡するのではなく、その時間枠における周波数の全体的なバランスを教えてくれます。個人の脳の状態は急速に変化する可能性がありますが、標準的なPSDプロットはそのような変化を平均化してしまうため、この点は極めて重要です。
医療グレードのヘッドセットとコンシューマー向けEEG headsetsを比較した研究者らは、この平均化にウェルチの修正ペリオドグラムと呼ばれる特定の手法を用いました。これは記録をオーバーラップするセグメントに分割し、それらのスペクトルを平均化して、よりクリーンでノイズの少ないPSD推定値を生成する手法です。
このような信号処理ステップは、ランダムな変動ではなく基礎となる脳活動を反映したPSD曲線を生成するための標準的なプラクティスです。
パワースペクトル密度の段階的な計算手順
PSDの計算には、通常、単一のコマンドではなく、信号品質および信号処理に関する一連の判断が必要です。まず、記録に欠損データ、不適切なチャンネル、顕著なアーチファクトがないか確認されます。その後、アナリストはセグメントを選択し、必要に応じて低周波のトレンドを除去またはモデル化し、窓関数を適用し、データを周波数領域に変換して、結果をスケーリングします。
得られたスペクトルは、元の時系列データと併せて確認する必要があり、該当する場合は体動、電極の外れ、眼球運動、筋活動、または実験イベントのアノテーションも併せて確認すべきです。この比較により、本物のリズミカルなパターンと、たまたまスペクトルパワーを生み出した一時的な乱れとを区別しやすくなります。
実践的なワークフローは、以下のステージに整理できます。
サンプリングレート、チャンネルラベル、基準電極、記録時間を検証する。
使用可能なエポックを特定し、アーチファクトや除外されたセグメントに対処する。
セグメント長、オーバーラップ、窓関数、周波数範囲を選択する。
PSDを推定してスケーリングし、ピークと帯域制限パワーを確認する。
結果を再現および比較できるよう、パラメータを記録する。
これらの手順により、計算の追跡可能性が確保されます。また、視覚的に見栄えの良いスペクトルが、それを生み出した前処理や推定の選択と無関係であるかのように扱われるのを防ぐこともできます。
一般的な手法：ペリオドグラムとウェルチ法
ペリオドグラムは、選択したセグメントに対するフーリエ変換の二乗絶対値を計算します。これは単純明快であり、セグメントが長ければ細かな公称周波数間隔を得ることができますが、その推定値は周波数ビンごとに大きく変動する可能性があります。このばらつきは、ノイズの多い生体記録において特に顕著です。
ウェルチ法は、信号を重複するセグメントに分割し、各セグメントに窓関数を適用してそれぞれペリオドグラムを計算し、その結果を平均化します。平均化により、通常はより滑らかで分散の低い推定値が得られますが、セグメントが短くなることで周波数分解能が低下する場合があります。最適なバランスは、分析において狭いピーク、広い帯域、あるいは安定した平均活動のいずれを重視するかによって異なります。
この違いは単に見た目だけの問題ではありません。スペクトルが滑らかになると広範なパターンの比較が容易になりますが、過度に平均化すると、短時間の変化や局所的な変化が見えなくなる可能性があります。時間とともに変化するEEGの場合、記録全体にわたる単一のPSDよりも、時間周波数アプローチの方が適切な場合があります。
健康な安静時EEGスペクトルの特徴
健康な成人に座ってもらい、リラックスした状態でEEGを記録すると、得られるPSD曲線には識別可能な特徴が現れます。パワーは最も低い周波数で最も高くなり、周波数が増加するにつれて着実に低下して、1/fのような減少傾向を示します。
この下降曲線は前述の非周期的（aperiodic）成分を反映しており、その急峻さはスペクトル指数（βと表記されることもあります）と呼ばれる単一の数値で表されます。減少が急なほど負の指数は大きくなり（絶対値が大きくなり）、減少がなだらかなほどパワーは周波数全体により均等に分布していることを意味します。
この背景の減少傾向の上に、閉眼時の安静な脳は通常、8〜12 Hz付近に明確な隆起（アルファピーク）を示します。この隆起は通常、後頭部で最も強くなります。同じ人が目を開けると、そのアルファピークは目に見えて縮小します。
医療グレードのEEGシステムとコンシューマー向けヘッドバンドを直接比較したRattiらの研究では、まさにこのパターンが記録されました。テストしたすべてのシステムにおいて、データ内の最大のリズム的特徴は、閉眼時と開眼時の条件に対して予測どおりに反応する後頭部アルファリズムでした。
その同じ比較研究では、より実践的な疑問も投げかけられました：同じ人物を別々の機会に2回記録した場合、PSD曲線は同じように見えるだろうか？
その答えは、機器の品質に大きく依存していました。医療グレードのシステムでは、訪問間で一貫したFp1電極のPSD（Fp1は前頭部に位置し、テストされた4つのシステムすべてに共通する唯一の電極でした）が得られたのに対し、コンシューマー向けシステムではテスト・再テストセッション間で最も高い相対的変動が見られました。このことは、安定したPSD曲線が単なる脳の特性にとどまらず、信号がどれほど慎重に捕捉されたかの特性でもあることを示しています。
麻酔が脳のスペクトルシグネチャーをどのように変化させるか
静的なPSDは、脳の状態が劇的に変化するときに特に多くの情報をもたらします。その最も直接的な例が麻酔です。Colomboらによる2019年の研究では、キセノン、プロポフォール、ケタミンの3種類の薬剤による麻酔の前および麻酔中の健康なボランティアにおいて、安静時EEGが記録されました。
各セッションの後、研究者らは無反応の期間中に意識的体験の記憶があったかどうかを尋ね、意識があったと報告された状態と、意識が存在しなかったか思い出すことができなかったとみられる状態とを区別しました。
スペクトル指数はこの区別を正確に追跡しました。意識的体験の喪失に関連する薬剤であるキセノンとプロポフォールはいずれも、各被験者内の負のスペクトル指数を大幅に減少させました。つまり、PSD曲線が顕著に急峻になったのです。1〜40 Hzの広範囲にわたるこの傾きの急峻化は、単なる無反応ではなく、意識的体験の欠如を示していると考えられます。
ケタミンは異なる様相を示しました。被験者はケタミン投与下で行動学的には無反応であったにもかかわらず、通常はある程度の意識的体験を保持しており、全体的なPSDの減衰は通常の覚醒状態と同様に見えました。唯一の例外は、特に20〜40 Hzの範囲における曲線の平坦化であり、これは他の2つの薬剤と比較したケタミン特有の作用機序と一致する特徴です。
おそらく最も説得力のある証拠は、このスペクトル指数と、摂動複雑性指数（PCI）と呼ばれる独立して得られた別の測定値との間の相関関係でした。PCIは、脳を刺激してその結果生じる反応の複雑さを測定する別の手法から生成されます。これら2つの全く異なる測定値の間に見られる高い相関関係は、スペクトル指数が偶然のパターンではなく、意識そのものの存在に結びついた真のマーカーであるという考えを支持しています。
ADHDがEEGパワースペクトルをどのように変化させるか
麻酔下で変化するのと同じ非周期的な減少傾向は、ADHDを持つ子どもでも異なることが明らかになっています。ADHDのある子どもと定型発達の同年代の子どもを比較した3〜7歳児の研究では、ピークアルファ周波数や帯域パワー比などの伝統的な指標とともに、スペクトルの傾きとオフセット（非周期曲線の高さ）が測定されました。
投薬歴のないADHDの子どもは、定型発達の子どもと比較して、アルファパワーが高く、非周期オフセットが大きく、スペクトルの傾きが急であることが示されました。これら2つの状態は重症度やメカニズムの点で明らかに比較できるものではありませんが、麻酔の研究で見られたのと同じ方向の変化を反映しているため、この傾きが急になるパターンは注目に値します。
スペクトルの傾きは、ADHD研究で長年使用されてきた従来のシータ/ベータ比とも相関しており、傾きが重複する情報を捉えつつ、よりクリーンで包括的なマーカーを提供できる可能性が示唆されました。
さらに別の詳細がこの発見の重要性を裏付けています。中枢神経刺激薬による治療を受けたことのあるADHDの子どもは、薬の有効成分の効果が切れているはずの24時間の休薬期間後であっても、定型発達グループと同等の傾きとオフセットを示しました。
この研究ではこれを確定した治療効果ではなく示唆として位置づけていますが、ADHDにおける変化したスペクトルパターンが固定的なものではない可能性を提起しています。スペクトルの傾きは、過去の研究において実行機能や脳の興奮性活動と抑制性活動のバランスとも関連づけられているため、著者らは、変化した傾きがシータ/ベータ比計算の統計的アーチファクトではなく、基礎となる状態に関する意味のある何かを反映している可能性があると示唆しています。
グループ | スペクトルの傾き | アルファパワー | オフセット |
|---|---|---|---|
定型発達 | 正常な傾き | 正常なパワー | 正常なオフセット |
ADHD（未投薬） | より急な傾き | より高いパワー | より大きなオフセット |
ADHD（刺激薬治療歴あり） | 定型発達と同等 | 報告なし | 定型発達と同等 |
パワースペクトル密度を用いたてんかん発作の分類
静的PSDが有用であることが証明されている領域は、意識と注意だけではありません。てんかんの分類は、より自動化された計算科学的な例を提供しています。
あるアプローチでは、てんかん患者のEEG記録を、研究者が「パワースペクトル密度エネルギーダイアグラム」と呼ぶもの（本質的にはPSDの視覚的表現）に変換し、これらの画像を画像状データのパターン認識に適した機械学習モデルの一種であるディープ畳み込みニューラルネットワークに入力しました。
モデルは、EEGセグメントを次の4つのカテゴリに分類するようにトレーニングされました。
発作間欠期（発作と発作の間）
2つの異なる発作前期間（発作開始の30分前と10分前）
発作状態そのもの
公開てんかんデータセットを用いたこの手法は、4つのカテゴリすべてにおいて90%を超える平均分類精度を達成しました。この結果は、信号の瞬間ごとの変化を追跡する必要なく、アルゴリズムがこれらの臨床的に重要な脳状態を確実に区別するのに十分な特徴的情報が、静的PSDの全体的な形状に含まれていることを裏付けています。
パワースペクトル密度の研究応用
パワースペクトル密度曲線の形状は、認知機能の低下から集中した瞑想まで、幅広い脳状態の調査に使用されてきました。以下の研究では、これらのスペクトルシグネチャーを捕捉するためにEmotivの技術が使用されました。
安静時PSDによる認知機能低下の識別
安静時EEGの研究において、パワースペクトル密度が健康な個人と認知機能障害を持つ個人を区別できるかどうかが調査されました。14チャンネルのEmotiv Epoc X headsetを使用し、研究者らは健常対照群、軽度認知障害（MCI）患者、および認知症患者における2分間の閉眼安静状態のPSDを分析しました。
著者らは、EEGが認知症患者と健常対照群を効果的に区別したものの、MCI群と健常群の間には有意差を示さなかったと報告しています。これは、静的PSDが進んだ認知機能低下の特定には期待が持てる一方で、初期段階の検出には課題があることを示唆しており、認知機能障害に関する研究でも指摘されています。
PSD測定における湿式電極と乾式電極の比較
EEGの取得に使用される電極のタイプは、得られるパワースペクトル密度に影響を与える可能性があります。ある調査では、2つのEmotivデバイス（湿式のEpoc Xと乾式のInsight）を使用して、湿式電極と乾式電極から得られたPSDを直接比較しました。記録は、音声刺激タスクの前後に側頭部および前頭前部の電極から取得されました。
電極タイプの比較において、刺激後のPSDの減少率は乾式電極よりも湿式電極の方が高いことが観察され、PSDベースのプロトコルを設計する研究者にとって実践的な検討事項が浮き彫りになりました。
マントラ詠唱中のPSDの変化
パワースペクトル密度は、瞑想実践中の脳活動の調査にも使用されています。熟練者を対象とした研究では、Emotiv Flex 32チャンネルデバイスで記録されたEEGに対する「OM」マントラ詠唱の効果が調査されました。
著者らはウェルチ法を用いてPSDを計算し、詠唱中にアルファ波、ガンマ波、ベータ波、デルタ波を含む複数の周波数帯域でパワーが有意に増加することを発見しました。マントラ詠唱に関する研究は、これらのPSDの変化が、実践中におけるより深いリラクゼーションと集中状態を示唆していると報告しています。
カフェインがEEGパワースペクトルに与える影響
一般的な興奮剤が脳のスペクトルプロファイルに与える影響が、Emotiv Flexシステムを使用したパイロット研究で捉えられました。研究者らは、参加者が適量のカフェインを摂取する前後に中心部電極からEEGを記録しました。パワースペクトル密度の分析により、絶対アルファパワーの抑制とベータパワーの増強が有意に増加していることが明らかになりました。
カフェインに関するパイロット研究の著者らは、真の振動変化であることが確認されたこれらのPSD変化が、ウェアラブルセンサーで検出可能な皮質覚醒の上昇を示していると報告しています。
これらの研究は、認知、覚醒、集中した注意に関連する脳活動の変化を捉えるためのツールとしてのパワースペクトル密度の汎用性を示しています。これらの知見は、個別の状態における一貫したスペクトルパターンを示唆していますが、それらは特定された、しばしば小規模なサンプル集団および制御された条件下から得られたものです。この研究は、PSDだけで何らかの状態を診断したり個人の状態を予測したりできることを立証するものではありませんが、この測定値が一連の有意義な神経変化に対して敏感であることを示唆しています。
PSDを解釈する際の一般的な落とし穴
PSD曲線の信頼性は、その取得元である記録の信頼性に依存します。いくつかの実用的な問題により、根本的な脳機能とは無関係な形で曲線の形状が歪む可能性があります。
アーチファクトは、最も直接的な懸念事項です。瞬きや筋肉の動きは独自の電気信号を生成し、特に前頭部位の頭皮電極によって拾われます。
コンシューマー向けEEGシステムと医療グレードEEGシステムの比較では、コンシューマー向けデバイスの方がこの種の混入を受けやすいことが判明しました。特に、あるコンシューマー向けヘッドバンドではスペクトル全体にわたって広帯域のパワー増加が見られ、これは本物の神経リズムではなく筋肉や瞬きのアーチファクトと一致する特徴でした。したがって、PSDを計算する前にデータを慎重にクリーニングし、アーチファクトの少ないセグメントを選択することは、オプションではなく必須のステップです。
容積伝導は、EEG研究全般において広く議論されているもう1つの限界点です。電場は頭皮に達する前に組織や頭蓋骨を通って広がるため、脳内の単一の活動源が近くにある複数の電極に同時に記録される可能性があり、特定のPSD測定値の空間的精度が曖昧になります。
基準電極の選択も同様の状況をもたらします。EEGは本質的に差動測定であり、すべての電極の信号は、基準点（全電極の平均、耳の後ろの乳様突起、またはCzのような頭皮中央の位置など）に対して相対的に計算されます。その基準点を変更すると、チャンネル間でパワーが再配分され、PSD曲線の絶対的な形状が変わる可能性があります。
したがって、どの基準電極が使用されたかを把握することは、異なる論文やデータセット間でPSDの知見を比較する際に確認すべき標準的な情報となります。
最後に、テスト・再テスト信頼性があります。被験者と機器を一定に保ったとしても、別々の機会に繰り返し記録を行うと、PSD曲線にある程度の自然なばらつきが生じる可能性があります。
医療グレードのシステムでは2回の訪問間でより一貫したFp1のPSDが得られたのに対し、テストされたコンシューマー向けシステムでは最大の相対的ばらつきが示されました。これは、単一の被験者で観察された小さなPSDの差異を過剰に解釈することに対する実践的な注意点であり、特に記録機器やプロトコルが厳密に管理されていない場合にはなおさらです。
パワースペクトル密度を計算するためのツールとソフトウェア
PSDは、科学計算プログラミングライブラリ、専用のEEG解析パッケージ、または一般的な信号処理ソフトウェアで計算できます。不可欠な要件は、サンプリングされたデータへのアクセス、ならびにサンプリングレート、セグメンテーション、窓関数、トレンド除去、オーバーラップ、スケーリング、周波数制限の透明性のある制御です。グラフィカルな出力だけでは、推定値がどのように生成されたかを立証するのに十分ではありません。
再現可能な実装では、スペクトルとともに前処理および推定パラメータが記録されます。また、チャンネル基準、単位、除外されたエポック、および計算前に適用された変換も保持されます。これらの詳細は、結果をグループ比較で使用する場合や、記録セッション間で統合する場合に特に重要です。
ワークフローによって数値配列と診断プロットの両方が生成されると、結果の監査が容易になります。有用なチェック項目には、元の波形、アーチファクトのマークされた区間、選択された周波数範囲、合理的なセグメント選択間での推定値の比較などが含まれます。このようなチェックによって解釈が自動化されるわけではありませんが、隠れた前提を検出しやすくなります。
結論
パワースペクトル密度プロットは脳のリズムのスコア（楽譜）として機能し、その背景にある減少の急峻さは、大きく異なる精神状態の間でも極めて一貫したマーカーとして浮かび上がります。
特定の麻酔薬によって意識を失った場合でも、子どもがADHDに関連する注意困難を示した場合でも、発作が始まった場合でも、この全体的な形状は予測可能な形で変化し、共有される測定可能なシグネチャーを通じてこれらの状態を結びつけます。このパターンは、脳の安静時の電気的プロファイルが、特定の周波数における孤立したピークだけでなく、その機能状態に関する意味のある横断的な手がかりを保持していることを示唆しています。
PSDは記録枠にわたる静的な平均値しか提供しませんが、意識のある期間と意識のない期間を確実に区別したり、アルゴリズムが高い精度で発作を特定できるようにしたりするその能力は、神経科学における実用的な価値を証明しています。これらの相関関係はこれらの状態の原因を証明するものではありませんが、意識、注意の違い、または発作のリスクを測るのに役立つ、曲線の傾きという分かりやすい単一の要約数値を提供します。
研究者や臨床医にとって、そのシンプルな測定値は脳のダイナミクスへの強力な入り口を提供し、抽象的なリズムと現実世界の差異との間を橋渡しします。
参考文献
Ratti, E., Waninger, S., Berka, C., Ruffini, G., & Verma, A. (2017). Comparison of Medical and Consumer Wireless EEG Systems for Use in Clinical Trials. Frontiers in human neuroscience, 11, 398. https://doi.org/10.3389/fnhum.2017.00398
Colombo, M. A., Napolitani, M., Boly, M., Gosseries, O., Casarotto, S., Rosanova, M., ... & Sarasso, S. (2019). The spectral exponent of the resting EEG indexes the presence of consciousness during unresponsiveness induced by propofol, xenon, and ketamine. NeuroImage, 189, 631-644. https://doi.org/10.1016/j.neuroimage.2019.01.024
Robertson, M. M., Furlong, S., Voytek, B., Donoghue, T., Boettiger, C. A., & Sheridan, M. A. (2019). EEG power spectral slope differs by ADHD status and stimulant medication exposure in early childhood. Journal of neurophysiology, 122(6), 2427-2437. https://doi.org/10.1152/jn.00388.2019
Gao, Y., Gao, B., Chen, Q., Liu, J., & Zhang, Y. (2020). Deep convolutional neural network-based epileptic electroencephalogram (EEG) signal classification. Front. Neurol. 11, 375 (2020). https://doi.org/10.3389/fneur.2020.00375
よくある質問（FAQ）
EEGにおけるパワースペクトル密度（PSD）とは何ですか？
PSDは、各周波数におけるEEG信号の振幅の2乗を計算し、これを記録期間にわたって平均化する手法です。得られた曲線は、遅い脳リズムと速い脳リズムにエネルギーがどのように分布しているかを示し、複雑な生の波形を分離して解き明かします。
健康な安静時のPSD曲線はどのような形状をしていますか？
パワーは低周波数で最も高く、周波数とともに着実に低下し、1/fのような背景を形成します。閉眼時には後頭部に8〜12 Hz付近の明確なアルファピークが現れるのが典型的で、開眼すると減少します。
非周期成分は振動性ピークとどのように異なりますか？
非周期的（1/f様）成分は、明確な隆起を形成することなく、周波数の増加とともにパワーが滑らかに低下する成分です。アルファリズムなどの振動性ピークは、特定の周波数でこの背景の上に突出し、脳機能に関する異なる情報を担っています。
麻酔は脳のスペクトルシグネチャーをどのように変化させますか？
キセノンやプロポフォールのように意識的体験の喪失をもたらす薬剤は、1〜40 HzにわたってPSD曲線をより急峻にします。無反応であっても意識的体験を維持するケタミンは、20〜40 Hzの範囲の平坦化を除いて、覚醒時と同様の曲線を示します。
ADHDの子どもにはどのようなPSDの違いが見られますか？
ADHDの子どもは、定型発達の同年代の子どもと比較して、より高いアルファパワー、より大きな非周期的オフセット、およびより急なスペクトルの傾きを示します。刺激薬による治療後は、休薬期間を設けても、傾きとオフセットが定型発達グループに類似することがあります。
PSDはどのようにしててんかん発作の状態を分類できますか？
静的なPSD曲線を画像に変換し、ディープニューラルネットワークに入力して、発作間欠期、発作前期、発作状態を区別することができます。このアプローチは少なくとも1つの研究で高い分類精度を達成しており、全体的なスペクトル形状に十分な識別情報が含まれていることが確認されました。
なぜPSDは静的なスナップショットと見なされるのですか？
PSDは記録期間全体を1本の曲線に圧縮し、秒単位の変化を平均化します。時間とともに混合がどのように変化するかを追跡するのではなく、その時間枠中の全体的な周波数バランスを明らかにします。
PSDを解釈する際の一般的な落とし穴は何ですか？
筋活動や眼球運動のアーチファクトが信号を汚染し、特にコンシューマー向けデバイスにおいて曲線を歪める可能性があります。容積伝導や基準電極の選択もPSDの形状に影響を与え、安定した機器であっても、自然なテスト・再テストのばらつきにより小さな差異には意味がない場合があります。
PSDと振幅スペクトルの違いは何ですか？
PSDは周波数あたりの二乗振幅を表し、パワーの計算に適しています。一方、振幅スペクトルは周波数成分の大きさ（絶対値）を表し、異なる単位を使用します。
PSDのピークが大きいことは、常に脳活動が強いことを意味しますか？
いいえ。ピークが大きいことは、生理学的活動だけでなく、アーチファクト、基準電極の影響、記録条件、または解析パラメータを反映している場合もあります。
Emotivは、アクセスしやすいEEGおよび脳データツールを通じて神経科学研究の進展を支援するニューロテクノロジー分野のリーダーです。
免責事項：当ウェブサイトで提供される情報は教育および情報提供のみを目的としており、医師の診断や健康のアドバイスを意図したものではありません。このコンテンツには誤りが含まれている可能性があり、人生を左右するような健康、医療、またはライフスタイルに関する選択を行う際の根拠とすべきではありません。病状や治療に関する疑問については、常に医師またはその他の資格を有する医療提供者の助言を求めてください。
クリスティアン・ブルゴス
