(57)【要約】【課題】医療現場において気分障害を定量的かつ簡易に測定することが可能な気分障害測定装置および方法を提供する。【解決手段】気分障害測定装置は、被験者の頭部上に配置された複数の電位センサを用いて異なる複数の部位で測定された脳電位信号を取得する脳電位信号取得手段と、取得された脳電位信号を周波数解析することにより、複数の部位のうち所定の一つの部位での予め決められた第１の周波数帯域の総パワー値と、複数の部位のうち他の所定の一つの部位での予め決められた、第１の周波数帯域よりも低い第２の周波数帯域の総パワー値とをそれぞれ算出する第１のパワー値算出手段と、第１の周波数帯域の総パワー値と第２の周波数帯域の総パワー値との差分の、第１の周波数帯域および第２の周波数帯域の総パワー値に対する比率を被験者の気分障害を判断するための指標として算出する第１の演算手段と、を備える。【選択図】図１
2019154789.tif 000002
【特許請求の範囲】
【請求項１】
被験者の頭部上に配置された複数の電位センサを用いて異なる複数の部位で測定された脳電位信号を取得する脳電位信号取得手段と、
前記取得された脳電位信号を周波数解析することにより、前記複数の部位のうち所定の一つの部位での予め決められた第１の周波数帯域の総パワー値と、前記複数の部位のうち他の所定の一つの部位での予め決められた、第１の周波数帯域よりも低い第２の周波数帯域の総パワー値とをそれぞれ算出する第１のパワー値算出手段と、
前記第１の周波数帯域の総パワー値と前記第２の周波数帯域の総パワー値との差分の、前記第１の周波数帯域および前記第２の周波数帯域の総パワー値に対する比率を被験者の気分障害を判断するための指標として算出する第１の演算手段と、を備えることを特徴とする気分障害測定装置。
【請求項２】
国際１０−２０法で規定される電極配置において、前記所定の一つの部位はＰ３であり、前記他の所定の一つの部位はＯ２であることを特徴とする、請求項１に記載の気分障害測定装置。
【請求項３】
前記第１の周波数帯域はβ帯であり、前記第２の周波数帯域はδ帯および／またはθ帯であることを特徴とする、請求項１または２に記載の気分障害測定装置。
【請求項４】
前記脳電位信号取得手段は、被験者の頭部上の互いに異なる３つの部位での脳電位信号を取得するよう構成されており、
前記取得された各脳電位信号から脳深部の活動に起因する特定の周波数帯域の時系列データをそれぞれ抽出し、該それぞれ抽出された時系列データの位相関係に基づいて前記３つの部位で測定された脳電位信号の相関関係を示す相関値を算出し、算出された相関値に基づいて被験者の気分障害を判断するための指標を算出する第２の演算手段を備えることを特徴とする、請求項１から３までのいずれか一項に記載の気分障害測定装置。
【請求項５】
国際１０−２０法で規定される電極配置において、前記３つの部位は、Ｐ３、Ｐ４およびＯ２であることを特徴とする、請求項４に記載の気分障害測定装置。
【請求項６】
前記特定の周波数帯域はβ帯域であることを特徴とする、請求項４または５に記載の気分障害測定装置。
【請求項７】
前記脳電位信号取得手段は、被験者の前側頭部に配置された１つの電位センサを用いて該前側頭部での脳電位信号を取得し、
前記取得された脳電位信号を周波数解析することにより、予め決められた第３の周波数帯域の総パワー値と、該第３の周波数帯域よりも高い第４の周波数帯域の総パワー値とをそれぞれ算出する第２のパワー値算出手段と、
前記第４の周波数帯域の総パワー値に対する前記第３の周波数帯域の総パワー値の比率を被験者の気分障害を判断するための指標として算出する第３の演算手段を備えることを特徴とする、請求項１から６までのいずれか一項に記載の気分障害測定装置。
【請求項８】
国際１０−２０法で規定される電極配置において、前記１つの電位センサが配置される部位はＦ７であることを特徴とする、請求項７に記載の気分障害測定装置。
【請求項９】
前記第３の周波数帯域はδ帯であり、前記第４の周波数帯域はθ帯および／またはβ帯であることを特徴とする、請求項８に記載の気分障害測定装置。
【請求項１０】
前記第１の演算手段で算出された前記指標、前記第２の演算手段で算出された前記指標および前記第３の演算手段で算出された前記指標に対して主成分分析を行って複数の主成分の主成分得点を算出する主成分得点算出手段と、
算出した主成分の主成分得点に基づいて開眼時と閉眼時を分類する開眼閉眼分類手段と、を備えることを特徴とする、請求項４から６までのいずれかを引用する請求項７から９までのいずれか一項に記載の気分障害測定装置。
【請求項１１】
被験者の頭部上に配置された複数の電位センサを用いて異なる複数の部位で測定された脳電位信号を取得するステップと、
前記取得された脳電位信号を周波数解析することにより、前記複数の部位のうち所定の一つの部位での予め決められた第１の周波数帯域の総パワー値と、前記複数の部位のうち他の所定の一つの部位での予め決められた、第１の周波数帯域よりも低い第２の周波数帯域の総パワー値とをそれぞれ算出するステップと、
前記第１の周波数帯域の総パワー値と前記第２の周波数帯域の総パワー値との差分の、前記第１の周波数帯域および前記第２の周波数帯域の総パワー値に対する比率を被験者の気分障害を判断するための指標として算出するステップと、を含むことを特徴とする気分障害測定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、頭部に設置した複数の電極により観測される脳電位の時系列データを解析することで被験者の気分障害の度合いを定量的に測定する装置および方法に関する。
【背景技術】
【０００２】
気分障害とは、慢性疼痛などの疾病などが原因で、常に気分が落ち込んだり、あるいは常に高まったりすることで日常生活に様々な支障をきたしてしまう心の病気である。医療現場では、うつ病等の気分障害を診断する際に、問診による診断が主である。また、末梢性疼痛に関しては、慢性的に強い痛みが生じている場合、 2点識別覚の閾値の増加といった中枢神経系を含む機能異常が報告されている。また、外部からの侵害刺激に対する実験痛を定量的に評価する熱的方法（Hardyの輻射熱法）があり、痛みを定量化する方法として優れていると言われている（非特許文献１）。
【０００３】
しかし、心因性の疼痛に関しては患者の主観的評価に大きく依存し、その客観化と定量化は非常に難しいとされており、定量的な診断方法の手段がない。
そこで、大規模な装置を使用することなく、簡易的に診断を行うことができ、気分障害の度合いを定量的に評価できる手法が求められる。
【０００４】
脳が悲しみや恐怖といった精神的ストレスを受けると、情動の中枢である大脳辺縁系に大きな影響が及ぶが、一次的な情動反応は、主として、扁桃体部分で生成されると考えられる。例えば、非特許文献２には、ペットを亡くした２０名の被験者に当該ペットに関連する言語的タスクを課してｆＭＲＩで活動部位を調査した結果、悲しみや忌避に対応して扁桃体を中心とした部分に優位な活性が観測されることが開示されている。この一次的な神経活動は、同じく大脳辺縁系の中にある視床下部の神経核に働きかけ、交換神経及び副交感神経からなる自律神経系を介して、二次的な生理反応（例えば血圧や、脈拍、瞳孔反応等）を惹き起こすことが知られている。また二次的な反応は前頭部にも及び、側坐核を中心とした「喜びの中枢」としての自己報酬系の活動を抑制することが知られている。更には、精神疾患により喜びの感情が低下した被験者では、側坐核の活動低下とともに、前帯状皮質の脳波（ＥＥＧ）のδ成分が増加することが知られている（非特許文献３）。また、同一の電極について定義される脳波のβ波とθ波のパワーの比率を用いて、気分障害の度合いの指標としている技術もある（非特許文献４）。これは、閉眼時、気分障害によりβ波の総パワー値がθ波に比べ、増加することを特徴としている。
【先行技術文献】
【非特許文献】
【０００５】
【非特許文献１】山本光璋 「痛みの測定」日本臨床麻酔学会誌 Vol.7 No.2, pp.146-156, Apr. 1987
【非特許文献２】Peter J. Freed, Ted K. Yanagihara, Joy Hirsch, and J. John Mann, “Neural mechanism of grief regulation”, Biol. Psychiatry, vol.66, no.1, pp.33-40, July 2009.
【非特許文献３】Jan Wacker, Daniel G. Dillon, Diego A. Pissagalli, “The role of the nucleus accumbens and rostral anterior cingulate cortex in anhedonia: Integration of resting EEG, fMRI, and Volumetric Techniques, Neuroimage, vol.46,  no.1,pp.327-337, May 2009.
【非特許文献４】Bachmann M1, Lass J, Suhhova A, Hinrikus H.,’ Spectral asymmetry and Higuchi's fractal dimension measures of depression electroencephalogram.’, Comput Math Methods Med. 2013;2013:251638. doi: 10.1155/2013/251638. Epub 2013 Oct 22.
【発明の概要】
【発明が解決しようとする課題】
【０００６】
気分障害における脳活動を捉えるためには、ＭＲＩやＰＥＴのような大規模な装置を必要とするが、医療現場においても患者にも経済的な負担が大きい。一方、脳電位（脳波）は大がかりな装置を必要としない。そこで、本発明は、医療現場において気分障害を定量的かつ簡易に測定することが可能な気分障害測定装置および方法を提供することを主目的とする。
【課題を解決するための手段】
【０００７】
上記の課題は以下の特徴を有する本発明によって解決される。すなわち、本発明の一態様としての気分障害測定装置は、被験者の頭部上に配置された複数の電位センサを用いて異なる複数の部位で測定された脳電位信号を取得する脳電位信号取得手段と、前記取得された脳電位信号を周波数解析することにより、前記複数の部位のうち所定の一つの部位での予め決められた第１の周波数帯域の総パワー値と、前記複数の部位のうち他の所定の一つの部位での予め決められた、第１の周波数帯域よりも低い第２の周波数帯域の総パワー値とをそれぞれ算出する第１のパワー値算出手段と、前記第１の周波数帯域の総パワー値と前記第２の周波数帯域の総パワー値との差分の、前記第１の周波数帯域および前記第２の周波数帯域の総パワー値に対する比率を被験者の気分障害を判断するための指標として算出する第１の演算手段と、を備える。
【０００８】
本発明の一態様としての気分障害測定装置において、国際１０−２０法で規定される電極配置において、前記所定の一つの部位はＰ３であり、前記他の所定の一つの部位はＯ２であることが好ましい。
【０００９】
本発明の一態様としての気分障害測定装置において、前記第１の周波数帯域はβ帯であり、前記第２の周波数帯域はδ帯および／またはθ帯であることが好ましい。
【００１０】
本発明の一態様としての気分障害測定装置において、前記脳電位信号取得手段は、被験者の頭部上の互いに異なる３つの部位での脳電位信号を取得するよう構成されており、前記取得された各脳電位信号から脳深部の活動に起因する特定の周波数帯域の時系列データをそれぞれ抽出し、該それぞれ抽出された時系列データの位相関係に基づいて前記３つの部位で測定された脳電位信号の相関関係を示す相関値を算出し、算出された相関値に基づいて被験者の気分障害を判断するための指標を算出する第２の演算手段を備えることが好ましい。
【００１１】
本発明の一態様としての気分障害測定装置において、国際１０−２０法で規定される電極配置において、前記３つの部位は、Ｐ３、Ｐ４およびＯ２であることが好ましい。
【００１２】
本発明の一態様としての気分障害測定装置において、前記特定の周波数帯域はβ帯域であることが好ましい。
【００１３】
本発明の一態様としての気分障害測定装置において、前記脳電位信号取得手段は、被験者の前側頭部に配置された１つの電位センサを用いて該前側頭部での脳電位信号を取得し、前記取得された脳電位信号を周波数解析することにより、予め決められた第３の周波数帯域の総パワー値と、該第３の周波数帯域よりも高い第４の周波数帯域の総パワー値とをそれぞれ算出する第２のパワー値算出手段と、前記第４の周波数帯域の総パワー値に対する前記第３の周波数帯域の総パワー値の比率を被験者の気分障害を判断するための指標として算出する第３の演算手段を備えることが好ましい。
【００１４】
本発明の一態様としての気分障害測定装置において、国際１０−２０法で規定される電極配置において、前記１つの電位センサが配置される部位はＦ７であることが好ましい。
【００１５】
本発明の一態様としての気分障害測定装置において、前記第３の周波数帯域はδ帯であり、前記第４の周波数帯域はθ帯および／またはβ帯であることが好ましい。
【００１６】
本発明の一態様としての気分障害測定装置において、前記第１の演算手段で算出された前記指標、前記第２の演算手段で算出された前記指標および前記第３の演算手段で算出された前記指標に対して主成分分析を行って複数の主成分の主成分得点を算出する主成分得点算出手段と、算出した主成分の主成分得点に基づいて開眼時と閉眼時を分類する開眼閉眼分類手段と、を備えることが好ましい。
【００１７】
本発明の一態様としての気分障害測定方法は、被験者の頭部上に配置された複数の電位センサを用いて異なる複数の部位で測定された脳電位信号を取得するステップと、前記取得された脳電位信号を周波数解析することにより、前記複数の部位のうち所定の一つの部位での予め決められた第１の周波数帯域の総パワー値と、前記複数の部位のうち他の所定の一つの部位での予め決められた、第１の周波数帯域よりも低い第２の周波数帯域の総パワー値とをそれぞれ算出するステップと、前記第１の周波数帯域の総パワー値と前記第２の周波数帯域の総パワー値との差分の、前記第１の周波数帯域および前記第２の周波数帯域の総パワー値に対する比率を被験者の気分障害を判断するための指標として算出するステップと、を含む。
【発明の効果】
【００１８】
本発明によれば、２つの電位センサのうち一方から取得した脳電位信号の第１の周波数帯域の総パワー値と、他方の電位センサから取得した脳電位信号の、第１の周波数帯域よりも低い第２の周波数帯域のパワー値との差分の、第１の周波数帯域および第２の周波数帯域の全パワー値に対する比率を気分障害を測定する際の指標としたことにより、簡単に医療現場において気分障害を簡易に測定することができ、しかも異なる皮質部位の活動の相互関係を反映して高い精度で気分障害を評価することが可能になる。
【図面の簡単な説明】
【００１９】
【図１】本発明の実施形態に係る気分障害測定装置の概略図である。
【図２】国際１０−２０法の電極配置を説明する図である。
【図３】第１実施形態の気分障害測定装置の処理部における処理ブロック（機能）を示す図である。
【図４】第１実施形態の気分障害測定方法の処理の流れを示すフローチャートである。
【図５】第２実施形態の気分障害測定装置の処理部における処理ブロック（機能）を示す図である。
【図６】第２実施形態の気分障害測定方法の処理の流れを示すフローチャートである。
【図７】第３実施形態の気分障害測定装置の処理部における処理ブロック（機能）を示す図である。
【図８】第３実施形態の気分障害測定方法の処理の流れを示すフローチャートである。
【図９】第４実施形態の気分障害測定装置の処理部における処理ブロック（機能）を示す図である。
【図１０】第４実施形態の気分障害測定方法の処理の流れを示すフローチャートである。
【図１１】実施例１の測定結果を示す図である。
【図１２】実施例２の測定結果を示す図である。
【図１３】実施例２において、ＮＸＰＲ値とＲＦ７値を２次元プロットした図である。
【図１４】実施例３における電位センサの配置を示した図である。
【図１５】実施例３において脳電位の３重相関値を算出するための処理ブロックを示す図である。
【図１６】実施例３において３重相関値Ｓｉを算出する処理の流れを示すフローチャートである。
【図１７】実施例３における３重相関表示部を示す図である。
【図１８】（ａ）、（ｂ）は、実施例３において、２つの遅延パラメータ（τ１、τ２）が形成する特徴空間上にプロットされた３重相関値分布の３次元表示を示す図である。
【図１９】３重相関値分布の３次元表示を上からみた図であって、３つの信号が同符号をとる領域を白で表し、３つの信号のいずれか１つの符号が異なる領域を黒で表した図である。
【図２０】実施例３において、３重相関値の標準偏差Ｓ値、ＳＤ値から作成したガウシアン分布を示す図である。
【図２１】実施例３において、ガウシアン分布の等高線の高さの値とＮＸＰＲ値とを２次元プロットした図である。
【図２２】実施例４における主成分分析結果を示す図である。
【発明を実施するための形態】
【００２０】
以下、図面を参照し本発明の実施の形態の気分障害測定装置および方法について説明する。
【００２１】
第１実施形態の気分障害測定装置は、被験者が気分障害状態にあるか否かを判別する際に用いられる第１の気分障害指標を生成する。図１は、第１実施形態の気分障害測定装置の装置構成例を概念的に示す図である。気分障害測定装置は、被験者の頭部の異なる部位に対応して少なくとも２つの電位センサ１０２が取り付けられて脳内神経活動に基づく電位を測定する頭部装着部１０４と、各電位センサ１０２と信号ケーブル１０６で接続された増幅器・帯域フィルタ１０８と、増幅器・帯域フィルタ１０８と信号ケーブル１１０で接続された解析用コンピュータ１１２とを備える。また、気分障害測定装置は、基準電位測定用の基準電位センサ１１４をさらに有する。頭部装着部１０４は、図２に示す国際１０−２０法の電極配置に従い各電位センサ１０２が配置されるように、被験者の頭部に装着する。この場合、頭部装着部１０４は、それぞれ扁桃体に対応する中心部（Ｃ３）、頭頂部（Ｐ３）、中側頭部（Ｔ３）および後側頭部（Ｔ５）並びに後頭部（Ｏ２）のうちの２つの部位に電位センサ１０２がそれぞれ配置されるように被験者へ装着するのが好ましい。頭部装着部１０４は、国際１０−２０法に基づくヘルメット型電位センサを用いて、選択的に２つの電位センサ１０２を使用することができ、この場合も、Ｃ４、Ｐ３、Ｔ３、Ｔ５およびＯ２のうちのいずれか２つの電極配置に対応する電位センサ１０２を使用するのが好ましい。なお、被験者頭部の所定の部位に電位センサ１０２を配置できる限り、頭部装着部１０４の形態に制限はない。測定用の少なくとも２つの電位センサ１０２および基準電位センサ１１４は無線通信機能を有し、同様に無線通信機能を有する解析用コンピュータ１１２へ、測定用の各電位センサ１０２と基準電位センサ１１４から得られる脳電位信号の差分を２つの脳電位信号として、無線で送信するようにしてもよい。
【００２２】
解析用コンピュータ１１２は、各種電算処理や算出を行うものであり、少なくとも１つのＣＰＵ（Ｃｅｎｔｒａｌ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）を含む処理部と、算出結果などを表示する表示部（ディスプレイ、プリンタなど）と、各種データやプログラムを格納する記憶部と、有線および／または無線通信を行う通信部とを有する。なお、解析用コンピュータ１１２は、増幅器・帯域フィルタの機能を有してもよく、この場合、増幅器・帯域フィルタ１０８は省略され得る。このような帯域フィルタの機能や請求項に記載した各種手段などは、解析用コンピュータ１１２において処理部のＣＰＵが所定のプログラムを実行することによって実現することができる。解析用コンピュータ１１２はＰＣ（Ｐｅｒｓｏｎａｌ  Ｃｏｍｐｕｔｅｒ）のような汎用コンピュータであってもよいし、専用コンピュータであってもよい。処理部は、ＣＰＵに代えてまたは加えて、特定用途向けの集積回路（ＡＳＩＣ）、ＤＳＰ（Ｄｉｇｉｔａｌ  Ｓｉｇｎａｌ  Ｐｒｏｃｅｓｓｏｒ）、ＧＰＵ（Ｇｒａｐｈｉｃｓ  Ｐｒｏｃｅｓｓｉｎｇ  Ｕｎｉｔ）等の少なくとも１つを含んでいてもよい。
【００２３】
第１実施形態の気分障害測定装置は、例えば処理部のＣＰＵが所定のプログラムを実行することによって得えられる機能部として、図３に示すように、脳電位信号取得手段としての脳電位信号取得部１１６と、第１のパワー値算出手段としての第１のパワー値算出部１１８と、第１の演算手段としての第１の演算部１２０とを有する。
【００２４】
第１実施形態の気分障害測定装置による気分障害測定方法の処理を図４のフローチャートに示すように、脳電位信号取得部１１６は、被験者の頭部の所定部位（例えばＰ３およびＯ２の電極位置）に配置された２つの電位センサ１０２から脳内神経活動に基づく脳電位信号を一定のサンプリング間隔で取得する（Ｓ１０１）。この後、第１のパワー値算出部１１８は、取得された脳電位信号を周波数解析することにより、上記２つの部位のうち一方の部位（例えばＰ３）での予め決められた第１の周波数帯域の総パワー値と、他方の部位（例えばＯ２）での予め決められた、第１の周波数帯域よりも低い第２の周波数帯域の総パワー値とをそれぞれ算出する（Ｓ１０２）。好ましくは、第１の周波数帯域はβ帯域であり、第２の周波数帯域はδ帯またはδ〜θ帯である。
【００２５】
そして、第１の演算部１２０は、下記式（１）に示すように第１の周波数帯域（例えばβ帯域）の総パワー値Ｐｉ（β）と第２の周波数帯域（例えばθ帯域）の総パワー値Ｐｊ（θ）との差分の、第１の周波数帯域および第２の周波数帯域の全パワー値に対する比率（ＮＸＰＲ：Ｎｏｒｍａｌｉｚｅｄ  Ｃｒｏｓｓ  Ｐｏｗｅｒ  Ｒａｔｉｏ）を算出する（Ｓ１０３）。第１実施形態では、このＮＸＰＲ値を第１の気分障害指標として用いることにより、気分障害の判定を行う（Ｓ１０４）。
2019154789.tif 000003
【００２６】
特に扁桃体付近および／または後頭部に配置された２つの電位センサ１０２のうち一方から取得した脳電位信号の第１の周波数帯域の総パワー値と、他方の電位センサから取得した脳電位信号の第２の周波数帯域のパワー値との差分の、第１の周波数帯域および第２の周波数帯域の全パワー値に対する比率（ＮＸＰＲ値）は、気分に何らかの不安や不快感など、気分障害があるときには、周波数帯域の高い第１の周波数帯域の総パワー値に比べて周波数帯域の低い第２の周波数帯域の総パワー値が増えることで負の値を取る傾向がある。上述した非特許文献４に記載されている方法では、同一の電位センサから測定されるβ波とθ波のパワー値の差を正規化し表示しているのに対し、本実施形態では、相互に異なる２つの電位センサ１０２からの脳電位信号を用い、脳電位信号別に異なる周波数帯（例えばβ波とθ波）のパワー値の差を正規化している。これにより、異なる皮質部位の活動の相互関係を反映した評価が可能になる。
【００２７】
第２実施形態の気分障害測定装置は、第１実施形態の気分障害装置に第２の演算手段を付加したものであり、装置のハードウェア構成は、３つの電位センサを要すること以外は第１実施形態とほぼ同じであり、例えば処理部のＣＰＵが所定のプログラムを実行することにより得られる機能部として、図５に示すように第２の演算手段として第２の演算部１２２を有する点で異なる。
【００２８】
第２実施形態の気分障害測定装置による気分障害測定方法の処理を図６のフローチャートに示すように、脳電位信号取得部１１６は、被験者の頭部上に所定部位に配置された３つの電位センサ１０２もしくはこれよりも多い電位センサのうち選択された３つの電位センサ１０２から脳内神経活動に基づく脳電位信号を一定のサンプリング間隔で取得する（Ｓ２０１）。好ましくは、３つの電位センサ１０２は国際１０−２０法で規定される電極配置において、Ｐ３，Ｐ４およびＯ２に配置されるものである。次いで第２の演算部１２２は、３つの電位が同符号の信号を選択する（Ｓ２０２）。電位の符号を判定する際の基準電位は、例えば皮質活動を直接反映しない耳朶が用いられうるが、増幅器の帯域フィルタで直流分は遮断されるので、実質的には、各々の電位センサごとの時間平均から見た正負の符号を判定することになる。なお当然のことながら、基準電位の取り方はこれらに限定されず、導電性ゴム電極も使用されうる。さらには、無線通信機能を有する測定用の３つの電位センサから得られる脳電位信号と、当該３つの電位センサの中央に配置される基準電位センサから得られる脳電位信号の差分を、３つの脳電位信号として無線で送信する構成とすることもできる。この場合は、前述のとおり、解析用コンピュータ１１２が帯域フィルタの機能を有する。
【００２９】
続いて第２の演算部１２２は、３重相関値を算出する（Ｓ２０３）。３重相関値は、３つの電位センサ１０２からのβ帯域、好ましくは２０Ｈｚ〜４０Ｈｚの周波帯域の電位信号をそれぞれＥＶＡ（ｔ）、ＥＶＢ（ｔ）、ＥＶＣ（ｔ）としたとき、１つの電極の電位信号に対し、τ１、τ２の時間ずれのある信号との積を使用する。以下に示す式（２）は３重相関値Ｓｔの１つの例示である。Ｔは３重相関値の演算対象時間であり、Δｔは各電位信号のデータサンプリング周期であり、Ｎは規格化するための定数であって、例えば３つの信号の積の計算回数である。
2019154789.tif 000004
【００３０】
３重相関値を算出することにより脳の深部の状態を観測することで、気分障害を伴う脳内活性をより的確に捉え、被験者の気分状態を推定することができる。３重相関値を算出した後、解析用コンピュータで所定の演算を行うことにより３重相関値の標準偏差Ｓ値（３重相関値の空間軸方向のばらつき）および３重相関値の時間軸方向の間隔の標準偏差ＳＤ値（３重相関値の時間軸方向のばらつき）という指標を算出し、当該指標により気分障害の識別判定を行うこともできる（Ｓ２０４）。
【００３１】
第３実施形態の気分障害測定装置は、第１実施形態の気分障害装置または第２実施形態の気分障害装置に第３の演算手段を付加したものであり、装置のハードウェア構成は第１実施形態とほぼ同じであり、例えば処理部のＣＰＵが所定のプログラムを実行することにより得られる機能部として、図７に示すように第２のパワー値算出手段としての第２のパワー値算出部１２４と、第３の演算手段として第３の演算部１２６とを有する点で異なる。
【００３２】
第３実施形態の気分障害測定装置による気分障害測定方法の処理を図８のフローチャートに示すように、脳電位信号取得部１１６は、被験者の頭部上に所定部位に配置された１つの電位センサ１０２から脳内神経活動に基づく脳電位信号を一定のサンプリング間隔で取得する（Ｓ３０１）。脳電位信号取得部１１６は、好ましくは、国際１０−２０法で規定される電極配置において、ｆ７に配置された電位センサ１０２から脳電位信号を取得する。次いで第２のパワー値算出部１２４は、取得された脳電位信号を周波数解析することにより、予め決められた第３の周波数帯域の総パワー値ＰＦ７_ｌｏｗと、該第３の周波数帯域よりも高い第４の周波数帯域の総パワー値ＰＦ７_ｈｉｇｈとをそれぞれ算出する（Ｓ３０２）。第３の周波数帯域は好ましくはδ帯、より好ましくは１．５６Ｈｚ〜３．１２Ｈｚであり、第４の周波数帯域はθ帯および／またはβ帯、より好ましくは４．８Ｈｚ〜１８．８Ｈｚである。第３の周波数帯域の総パワー値ＰＦ７_ｌｏｗと第４の周波数帯域の総パワー値ＰＦ７_ｈｉｇｈは下記式（３）、（４）をそれぞれ用いて求めることができる。
2019154789.tif 000005
ここで、ｋとｌ、ｍとｎは，パワースペクトルの周波数ビンであり、ｋ＜ｌ、ｍ＜ｎである。
【００３３】
そして、第３の演算部１２６は、例えば分解能１．５６Ｈｚで第４の周波数帯域の総パワー値ＰＦ７_ｈｉｇｈに対する第３の周波数帯域の総パワー値ＰＦ７_ｌｏｗの比率（ＲＦ７）を算出する（Ｓ３０３）。ＲＦ７は、下記式（５）を用いて求めることができる。第３実施形態では、この比率を第３の気分障害指標として用いることにより、気分障害の判定を行う（Ｓ３０４）。
2019154789.tif 000006
【００３４】
ストレス等の不快な気分を感じたとき、特に前側部（国際１０−２０法の電極配置のＦ７）の比較的低い周波数帯域におけるパワー値が増える傾向がある。したがって、本実施形態のように、予め決められた第３の周波数帯域の総パワー値ＰＦ７_ｌｏｗと、該第３の周波数帯域よりも高い第４の周波数帯域の総パワー値ＰＦ７_ｈｉｇｈとの比率（ＲＦ７）を気分障害指標とすることで、気分障害の判定を行うことができる。
【００３５】
第４実施形態は、上記に示した３の指標を用いて主成分分析を行い、第１主成分と第３主成分の主成分得点を用いて、閉眼時と開眼時とを分類するものである。具体的には、第４実施形態の気分障害測定装置は、第１〜第３実施形態の気分障害測定装置に主成分得点算出手段および開眼閉眼分類手段を付加したものであり、装置のハードウェア構成は第３実施形態とほぼ同じであり、例えば処理部のＣＰＵが所定のプログラムを実行して得られる機能部として、図９に示すように主成分得点算出手段としての主成分得点算出部１２８と開眼閉眼分類手段としての開眼閉眼分類部１３０とを有する点で異なる。
【００３６】
主成分得点算出部１２８は、第１の演算部１２０により得られる指標（ＮＸＰＲ値）、第２の演算部１２２により得られる指標（３重相関値）および第３の演算部１２６により得られる指標（ＲＦ７）に対して主成分分析を行って複数の主成分の主成分得点を算出するよう構成されている。開眼閉眼分類部１３０は、算出した主成分の例えば第１主成分と第３主成分の主成分得点に基づいて開眼時と閉眼時を分類するよう構成されている。
【００３７】
本実施形態によれば、気分障害の測定に際し被験者に開眼か閉眼かを指示しなくても、図１０に示すように、自動的に開眼か閉眼かが判定されるので、判定結果に基づき開眼時および閉眼時の気分障害評価をそれぞれ行うことが可能となる。
【実施例】
【００３８】
［実施例１］
実施例１では、第１実施形態の気分障害測定装置を用いて気分障害の定量的な評価を行った。被験者に視覚による気分障害を与えるため、被験者１７名に対して、気分が悪くなりそうな動画を見せ、開眼時の脳波を測定した。具体的な手順は以下の通りである。
（１）各被験者の安静閉眼（180秒）の脳波を測定
（２）各被験者の安静開眼（120秒）の脳波を測定
（３）各被験者に嫌悪コンテンツ（蛇の動画）を１２０秒間見せ、脳波を測定
（４）各被験者に目をつぶらせ、脳波の波形チェック（３０秒程度）
（５）各被験者に嫌悪コンテンツ（ごきぶりの動画）を２００秒間見せ、脳波を測定
（６）各被験者に目をつぶらせ、脳波の波形チェック（３０秒程度）
（７）各被験者の安静開眼（１２０秒）の脳波を測定
電位センサは国際１０−２０法電極配置に従うＰ３とＯ２にそれぞれ配置してこれらの２箇所の脳電位信号を取得し、Ｐ３の電位センサから取得した脳電位信号についてはβ帯域、Ｏ２の電位センサから取得した脳電位信号についてはθ帯域を用いて上記式（１）によりＮＸＰＲ値を算出して気分障害指標とした。
【００３９】
測定結果を図１１に示す。被験者１７人に、蛇、ごきぶりが写っている動画を見せて、先に定義したＮＸＰＲ値を算出した結果、安静開眼時には値が正を取るのに対し、蛇、ごきぶりの動画を見ているときには負の値をとった。下記表１では、気分障害は、ＮＸＰＲ値が正の値をとった場合に「気分障害無し」、負の値をとった場合に「気分障害有り」と評価した。また、今回のコンテンツにおいて、全評価数に対する、安静・開眼時の気分障害無し評価数および蛇、ごきぶり動画視覚時の気分障害有り評価数の合計の比率をＮＸＰＲ値を指標とする気分障害測定の判別率として算出したところ、下記表１に示すように７３％と高い判別率であった。このことにより、ＮＸＰＲ値が開眼時の気分障害の有効な指標であることが確認された。
（表１）
2019154789.tif 000007
【００４０】
［実施例２］
実施例２では、第３実施形態の気分障害測定装置を用いて気分障害の定量的な評価を行った。被験者に視覚による気分障害を与えるため、実施例１と同じ手順により被験者１７名に対して、気分が悪くなりそうな動画を見せ、開眼時の脳波を測定した。電位センサは国際１０−２０法電極配置に従うＦ７に配置してこれから脳電位信号を取得し、Ｆ７の電位センサから取得した脳電位信号の分解能１．５６Ｈｚで１．５６Ｈｚ〜３．１２Ｈｚの総パワー値と４．８Ｈｚ〜１８．８Ｈｚの総パワー値を上記式（３）、（４）を用いてそれぞれ算出し、１．５６Ｈｚ〜３．１２Ｈｚの総パワー値を４．８Ｈｚ〜１８．８Ｈｚの総パワー値で割った値（ＲＦ７）を気分障害指標とした。
【００４１】
測定結果を図１２に示す。実施例１と同様、被験者１７名に、へび動画、ごきぶり動画を見せたところ、図１２に示すように、蛇、ごきぶり動画ともに、安静開眼時に比べ、ＲＦ７指標値が大きくなることがわかった。したがって、ＲＦ７の値を気分障害指標として使用することができることが確認された。
【００４２】
次に、実施例１に示したＮＸＰＲ値とＲＦ７値の両者の値を図１３に示すように２次元プロットし、気分障害の評価を行った。気分障害がある場合、ＮＸＰＲ値は負の値をとり、ＲＦ７値は正の値をとる傾向にあることから、ＮＸＰＲ値とＲＦ７値の２次元分布において気分障害判定領域を第２象限とすることができる。今回のコンテンツにおいて、全評価数に対する、安静・開眼時の気分障害無し評価数および蛇、ごきぶり動画視覚時の気分障害有り評価数の合計の比率をＮＸＰＲ値およびＲＦ７値を指標とする気分障害測定の判別率として算出したところ、下記表２に示すように７５％と高い判別率であった。このことにより、ＮＸＰＲ値およびＲＦ７値が開眼時の気分障害の有効な指標であることが確認された。なお、図１３中、楕円状の判定線で囲んだ領域に不快に感じる動画が集中していることがわかる。この判定線は、例えばデータ点の分布を近似する２次元正規分布関数に閾値操作を施して得ることができ、ＮＸＰＲ値とＲＦ７値の２次元分布において第２象限を気分障害判定領域とするのに加えてあるいは代えて、このような判定線に基づく気分障害判定を行ってもよい。この判定線を用いた場合の判別率は８６％となり、有効な指標になることが分かる。
（表２）
2019154789.tif 000008
【００４３】
［実施例３］
実施例３では、第２実施形態の気分障害測定装置を用いて気分障害の定量的な評価を行った。音響による気分障害を測定するため、被験者１５名に対して、ストレス音楽とヒーリング音楽を閉眼状態で聞かせ、閉眼時の脳波を測定した。具体的な手順は以下の通りである。
（１）各被験者を１分間安静にさせる
（２）各被験者の安静開眼（１２０秒）の脳波を測定
（３）各被験者に閉眼状態でヒーリング音楽を３００秒間聞かせ脳波を測定
（４）閉眼状態で各被験者の脳波の波形チェック（３０秒程度）
（５）安静開眼状態で各被験者の脳波を測定（１２０秒）
（６）各被験者に閉眼状態でストレス音楽を３００秒間聞かせ脳波を測定
電位センサは国際１０−２０法電極配置に従うＰ３、Ｐ４およびＯ２にそれぞれ配置してこれらの３箇所の脳電位信号の時系列データを取得し、それぞれ抽出された時系列データの位相関係に基づいて３つの部位で測定された脳電位信号の相関関係を示す相関値を算出した。相関値は次のように求めることができる。
【００４４】
図１４に示すように、３つの電位センサＥＡ（１０２）、ＥＢ（１０２）、ＥＣ（１０２）が三角形の各頂点部分に配置され、別途設置される基準電位センサ１１４と、各電位センサとの差として、電位信号ＶＡ（ｔ）、ＶＢ（ｔ）、ＶＣ（ｔ）が計測される。なお前述の通り、Ｐ３、Ｐ４、Ｏ２の位置の３つの電位センサが使用されるのが好ましい。各脳電位信号は第２実施形態の気分障害測定装置の第２の演算部１２２によって処理される。図１５は第２実施形態の気分障害装置のより詳細な処理ブロックを示す図であり、例えば３ｃｈ増幅器・帯域フィルタと解析用コンピュータによって実現される。図１５に示すように、脳電位増幅器６０１によって増幅された信号はバンドパスフィルタ６０２によって、β波帯を主とする特定の周波数帯、例えば２０Ｈｚ〜４０Ｈｚの脳電位波形が抽出される。これは、脳深部に仮定したダイポール電位活動を定量的に解析し、脳深部の活動を評価するために行われる。ただし、帯域フィルタはβ波を主とする特定の周波数を検出できるものであり、図１５における周波数の値に限定されない。
【００４５】
次に、これら３つの信号による３重相関値の算出方法について示す。抽出された信号は３重相関値算出部６０３によって、図１６のフローチャートに示すように処理される。図１６は、ｉ秒からｉ＋１秒における３重相関値Ｓｉ（ｉ＝１、２、…、Ｔ）を算出する処理のフローチャートを示す。なお、ここで実施される処理は、趣旨を逸脱しない範囲において変更することができる。
【００４６】
前述の通り３つの信号が入力されるとサンプリング周期でデータが抽出され（Ｓ７０１）、それぞれの電位センサの電位ごとに標準偏差（σA、σB、σC）で割って規格化される（Ｓ７０２）。この規格化処理は１秒ごとに行うのが好ましいが、これに限定されない。
2019154789.tif 000009
なおバンドパスフィルタによる周波数抽出処理は、規格化処理の前後いずれかに行われる。また規格化処理の前には、ノイズ処理を行うのが好ましい。ノイズ処理は、例えば、１）±１００μＶ以上のセグメントを除く、２）フラットな電位（２５ｍｓｅｃ以上一定の電位だった場合）を除く、３）±１μＶ以内の電位が１秒以上続く場合は除く、という処理から構成される。
【００４７】
ここで、上記、３つの信号は、電位センサＥＡに対し、電位センサＥＢはτ１、電位センサＥＣはτ２の時間のずれがあるものとする。続いて、３つの信号の符号がすべて正（ＥＶＡ（ｔ）＞０、ＥＶＢ（ｔ）＞０、ＥＶＣ（ｔ）＞０）または、すべて負（ＥＶＡ（ｔ）＜０、ＥＶＢ（ｔ）＜０、ＥＶＣ（ｔ）＜０）の信号のみを計算対象とする処理をする（Ｓ７０３）。式（９）に示すように３重相関値は、時間ずれのある３つの電位信号の積を加算することで求められる（Ｓ７０４）。この処理は、ｔがｔ＝ｉ＋１秒となるまでΔｔ秒ずつずらして行われる（Ｓ７０６、Ｓ７０７）。なお図１６ではｔ＝ｉ秒からｉ＋１秒における３重相関値Ｓｉを算出していることからも分かる通り、全データ（Ｔ秒）について一度に計算するのではなく、所定時間ごとに、本実施例においては１秒ごとに３重相関値Ｓｉを求めＴ個の３重相関値の平均値を最終的には３重相関値とし、時間ずれτ１、τ２も１秒の中でΔｔ秒ずつずらして３重相関値を算出する。例えば、電位データサンプリング周波数をｆｓ（Ｈｚ）とすると、ｆｓ＝２００Ｈｚの場合はΔｔ＝１／ｆｓ＝０．００５秒ずつずらして、３つの電位信号の積を算出する。また、１秒ごとに３つの信号が正または負になった時の回数Ｎを求め（Ｓ７０５）、最後に割る（Ｓ７０８）。式（９）に１秒ごとの、３重相関値Ｓｉの計算式を示す。
2019154789.tif 000010
（ｉ＝１、２、…、T、τ１＝Δｔ、２Δｔ、…、１（秒）、τ２＝Δｔ、２Δｔ、…、１（秒））
【００４８】
このようにして、１秒ごとにＳｉを全データＴ秒まで計算する（Ｓ1、Ｓ2、・・・、ＳT）。Ｔ（秒）は好ましくは１０（秒）である。ただしＳｉは１秒ごとに算出されることに限定されない。τ１及びτ２の取りうる値はサンプリング周期の整数倍に等しい１秒以下の時間であるが、これらの値の大きさの最大値は１秒に限定されない。またサンプリング周期は０．００５秒に限定されない。なお３重相関値は、３つの信号の符号判定を行わずに、式９によって算出することもできる。
【００４９】
この結果は、図１７に示すように、３重相関表示部６０４によって２つの遅延パラメータ（τ１、τ２）が形成する特徴空間上にプロットすることで、２つの遅延パラメータ（τ１、τ２）が形成する特徴空間上にプロットされた３重相関値分布の疑似３次元表示をすることができる。図１８は参考として、α波帯域の３重相関値分布の疑似３次元表示を示すものであるが、相関を有しないデータの影響を排除するため、予め定められたｔの値、例えばｔ＝ｉ＋１、においてＥＶＡ（ｔ）、ＥＶＢ（ｔ−τ１）及びＥＶＣ（ｔ−τ２）のすべてが同符号であったＳｉ（τ１，τ２）のみをプロットしたものである。図１８（ａ）は３重相関値の空間的、時間的乱れが小さいものであり、（ｂ）は３重相関値の空間的、時間的乱れが比較的大きいものを例示したものである。
【００５０】
図１９は、図１８に例示したような３次元表示の図を上から見た図で、３つの波形が同符号をとる領域を白で表示し、３信号のどれか１つ符号が異なる領域を黒で表し、格子縞の乱れが生じていることが分かる。そこで、この乱れ（時間軸τ１、ｔ２方向のばらつき）を定量化するために、以下のように標準偏差ＳＤを定義する。
【００５１】
図１９に示す白い四角形の領域は、隣接する白い四角形の領域と、縦横方向にそれぞれ間隔を有する。その間隔を図１９に示すように、ｄｘｉ（ｉ＝１，２、…、ｍ）、ｄｙｊ（ｊ＝１，２、…、ｎ）とする。次いで式（１０）、式（１１）に示すように、ｍ個のｄｘｉの標準偏差Ｓｔｄ＿ｄｘとｎ個のｄｙｊの標準偏差Ｓｔｄ＿ｄｙを算出し、式（１２）に示すように２つの標準偏差の平均値を指標値ＳＤとする。
2019154789.tif 000011
【００５２】
また、式９で算出された３重相関値の標準偏差（空間軸方向のばらつき）Ｓは次式（１３）より求めることができる。
2019154789.tif 000012
【００５３】
そして、このようにして求めたＳ値、ＳＤ値からガウシアン分布を作成すると図２０に示すようになる。図２０に示すように、楕円の中心から距離が遠いところにヒーリング音楽が分布していることがわかる。
【００５４】
次に，上記ガウシアン分布により作成された等高線の高さの値とＮＸＰＲ値を用いて評価を行うため、図２１の散布図に示すようにストレス音楽、安静併願、ヒーリング音楽の３つの群をプロットした。
【００５５】
安静閉眼とヒーリング音楽を定常状態とし、ストレス音楽との線形判別分析を行うと、等高線からの値Ｔ値（標高値）とＮＸＰＲ値において、式（１４）の判別関数で７３％の判別率で分類できる。そこで、これを気分障害指標として使用し、値が大きいほど気分障害の度合いが大きいということができる。
2019154789.tif 000013
【００５６】
下の表３に示すように、全評価数に対する、安静・閉眼とヒーリング音楽視聴時の気分障害無し評価数およびストレス音楽視聴時の気分障害有り評価数の合計の比率をＮＸＰＲ値を指標とする気分障害測定の判別率として算出したところ、６０％の判別率（閾値をＮＸＰＲ＝０とする）であったが、３重相関値に基づく指標を併用することで判別率が１３％向上したことから、２つの指標を用いることで閉眼時の気分障害を測定する際の有効な気分障害指標になり得ることが分かる。
（表３）
2019154789.tif 000014
【００５７】
［実施例４］
実験で行った閉眼時と開眼時のすべてのコンテンツを用いて主成分分析を行い、第１主成分と第３主成分の主成分得点を２次元にプロットした結果、開眼時と閉眼時を分類することが可能となった。図２２に示すように、線形判別分析等により得られる直線（境界）より上が開眼時コンテンツ、直線より下が閉眼時コンテンツである。このことより、開眼と閉眼を分離することが可能となり、それぞれの測定条件に特化した処理を行うことが可能になる。例えば、被検者が開眼状態であると判定された場合には、実施例１に記載したＮＸＰＲ値に基づく評価が、また、閉眼と判定された場合には、実施例３に記載した３重相関値とＮＸＰＲ値とを組み合わせた指標に基づく評価が有効な評価手段となり得る。
【符号の説明】
【００５８】
１０２  電位センサ
１０４  頭部装着部
１０８  増幅器・帯域フィルタ
１１２  解析用コンピュータ
１１４  基準電位センサ
１１６  脳電位信号取得部
１１８  第１のパワー値算出部
１２０  第１の演算部
１２２  第２の演算部
１２４  第２のパワー値算出部
１２６  第３の演算部
１２８  主成分得点算出部
１３０  開眼閉眼分類部
【図１】
2019154789.tif 000015
【図２】
2019154789.tif 000016
【図３】
2019154789.tif 000017
【図４】
2019154789.tif 000018
【図５】
2019154789.tif 000019
【図６】
2019154789.tif 000020
【図７】
2019154789.tif 000021
【図８】
2019154789.tif 000022
【図９】
2019154789.tif 000023
【図１０】
2019154789.tif 000024
【図１１】
2019154789.tif 000025
【図１２】
2019154789.tif 000026
【図１３】
2019154789.tif 000027
【図１４】
2019154789.tif 000028
【図１５】
2019154789.tif 000029
【図１６】
2019154789.tif 000030
【図１７】
2019154789.tif 000031
【図１８】
2019154789.tif 000032
【図１９】
2019154789.tif 000033
【図２０】
2019154789.tif 000034
【図２１】
2019154789.tif 000035
【図２２】
2019154789.tif 000036