JP 4978860 B2 2012.7.18(57)【特許請求の範囲】【請求項１】
観測対象者の脳活動を時系列的に観測し、得られた観測データに基づいて前記観測対象者の行動を予測する行動予測方法において、
前記観測データを前記観測対象者の脳内の複数の電流源における時系列データに写像するフィルタを設定するステップと、設定したフィルタによって写像される時系列データの入力に対して前記観測対象者が特定の行動を起こす確率を出力する確率モデルを設定する
ステップと、新たに取得した観測データを前記フィルタにより前記電流源における時系列
データに変換するステップと、変換した時系列データを設定した確率モデルに入力することにより前記確率を算出するステップとを有し、算出した確率に基づいて前記観測対象者

10の行動を予測することを特徴とする行動予測方法。
【請求項２】
予測すべき行動の数だけ前記フィルタを設定し、各フィルタの前記確率への寄与度を算出するステップを更に有することを特徴とする請求項１に記載の行動予測方法。
【請求項３】
前記確率モデルは、学習すべきパラメータを含み、前記観測対象者が前記特定の行動を実行した際に観測される観測データに基づいて前記パラメータを学習するステップを更に有することを特徴とする請求項１又は請求項２に記載の行動予測方法。
【請求項４】
前記パラメータをベイズ推定法を用いて学習するステップを有することを特徴とする請

20

(2)

JP 4978860 B2 2012.7.18

求項３に記載の行動予測方法。
【請求項５】
前記確率モデルは、前記時系列データの入力に応じたスカラー値を出力する判別関数と、該判別関数の出力値に応じて所定値域の値を出力するロジスティック関数とにより記述されることを特徴とする請求項１乃至請求項４の何れか１つに記載の行動予測方法。
【請求項６】
前記確率モデルは、スパースロジスティック回帰モデルであることを特徴とする請求項
１乃至請求項５の何れか１つに記載の行動予測方法。
【請求項７】
観測対象者の脳活動を時系列的に観測し、得られた観測データに基づいて前記観測対象

10

者の行動を予測する行動予測装置において、
前記観測データを前記観測対象者の脳内の複数の電流源における時系列データに写像するフィルタを設定する手段と、設定したフィルタによって写像される時系列データの入力に対して前記観測対象者が特定の行動を起こす確率を出力する確率モデルを設定する手段と、新たに取得した観測データを前記フィルタにより前記電流源における時系列データに変換する手段と、変換した時系列データを設定した確率モデルに入力することにより前記確率を算出する手段とを備え、算出した確率に基づいて前記観測対象者の行動を予測する
ようにしてあることを特徴とする行動予測装置。
【請求項８】
予測すべき行動の数だけ前記フィルタを設定し、各フィルタの前記確率への寄与度を算20出する手段を更に備えることを特徴とする請求項７に記載の行動予測装置。
【請求項９】
前記確率モデルは、学習すべきパラメータを含み、前記観測対象者が前記特定の行動を実行した際に観測される観測データに基づいて前記パラメータを学習する手段を更に備えることを特徴とする請求項７又は請求項８に記載の行動予測装置。
【請求項１０】
前記パラメータをベイズ推定法を用いて学習するようにしてあることを特徴とする請求
項９に記載の行動予測装置。
【請求項１１】
前記確率モデルは、前記時系列データの入力に応じたスカラー値を出力する判別関数と30、該判別関数の出力値に応じて所定値域の値を出力するロジスティック関数とにより記述されることを特徴とする請求項７乃至請求項１０の何れか１つに記載の行動予測装置。
【請求項１２】
前記確率モデルは、スパースロジスティック回帰モデルであることを特徴とする請求項
７乃至請求項１１の何れか１つに記載の行動予測装置。
【発明の詳細な説明】【技術分野】【０００１】
本発明は、脳活動の情報から行動を予測する行動予測方法及び行動予測装置に関する。
【背景技術】

40【０００２】
近年の生体計測技術の進歩はめざましく、従来では測定が困難で誤差も大きかった脳から発生する微弱電場（脳波）及び微弱磁場（脳磁波）の計測精度が年々向上している。
【０００３】
脳活動を計測する目的としては、脳のどの部位が感覚、運動などの認知に関わるかを知ることが挙げられる。そのため、脳内の活動を高い空間分解能、また場合によっては高い
時間分解能で観測することが求められている。非侵襲的に脳活動を計測する計測手段として、神経細胞脳電気活動由来の電気現象を計測するＥＥＧ（electroencephalogram）及びＭＥＧ（magnetoencephalogram）、神経細胞発火後の二次的な現象として起こる脳血流又は代謝の変化を計測するｆＭＲＩ（functional magnetic resonance imaging）、ＮＩＲ

50

(3)

JP 4978860 B2 2012.7.18

Ｓ（Near Infrared Spectroscopy）、ＰＥＴ（positron emission tomography）等が知られている。
【特許文献１】国際公開第０３／０５７０３５号パンフレット【発明の開示】【発明が解決しようとする課題】【０００４】
しかしながら、前述したＥＥＧ、ＭＥＧ計測は、時間分解能に優れており、神経活動に直接由来する電気現象を非侵襲的に観測することができるという利点を有している。しかしながら、その一方で、ＥＥＧ、ＭＥＧ計測は、計測地点が脳の内部ではなく頭蓋外部に配置したセンサにより計測しており、また、神経細胞の電気活動が物理法則に従って拡散10した電場又は磁場を計測しているため、ＥＥＧ、ＭＥＧ計測のみから脳内の電気活動を推定するには空間分解能に関して原理的な制限が存在する。
【０００５】
特に、外部から予期しない刺激が入力された場合、センサ空間における特徴量に対する外乱の影響を予測することは困難であるという問題点を有している。
【０００６】
本発明は斯かる事情に鑑みてなされたものであり、観測データを観測対象者の脳内の複数の電流源における時系列データに写像するフィルタを設定し、設定したフィルタによって写像される時系列データの入力に対して観測対象者が特定の行動を起こす確率を出力する確率モデルを設定し、新たに取得した観測データを前記フィルタにより電流源における

20

時系列データに変換し、変換した時系列データに基づいて特定の行動を起こす確率を算出
する構成とすることにより、外部から予期せぬ外乱が入力された場合であっても、行動予測を高精度に実施することができる行動予測方法及び行動予測装置を提供することを目的とする。
【課題を解決するための手段】【０００７】
第１発明に係る行動予測方法は、観測対象者の脳活動を時系列的に観測し、得られた観測データに基づいて前記観測対象者の行動を予測する行動予測方法において、前記観測データを前記観測対象者の脳内の複数の電流源における時系列データに写像するフィルタを設定するステップと、設定したフィルタによって写像される時系列データの入力に対して30

前記観測対象者が特定の行動を起こす確率を出力する確率モデルを設定するステップと、新たに取得した観測データを前記フィルタにより前記電流源における時系列データに変換するステップと、変換した時系列データを設定した確率モデルに入力することにより前記確率を算出するステップとを有し、算出した確率に基づいて前記観測対象者の行動を予測することを特徴とする。
【０００８】
第２発明に係る行動予測方法は、予測すべき行動の数だけ前記フィルタを設定し、各フィルタの前記確率への寄与度を算出するステップを更に有することを特徴とする。
【０００９】
第３発明に係る行動予測方法は、前記確率モデルは、学習すべきパラメータを含み、前40

記観測対象者が前記特定の行動を実行した際に観測される観測データに基づいて前記パラメータを学習するステップを更に有することを特徴とする。
【００１０】
第４発明に係る行動予測方法は、前記パラメータをベイズ推定法を用いて学習するステップを有することを特徴とする。
【００１１】
第５発明に係る行動予測方法は、前記確率モデルは、前記時系列データの入力に応じたスカラー値を出力する判別関数と、該判別関数の出力値に応じて所定値域の値を出力する
ロジスティック関数とにより記述されることを特徴とする。
【００１２】50

(4)JP 4978860 B2 2012.7.18

第６発明に係る行動予測方法は、前記確率モデルは、スパースロジスティック回帰モデ
ルであることを特徴とする。
【００１３】
第７発明に係る行動予測装置は、観測対象者の脳活動を時系列的に観測し、得られた観測データに基づいて前記観測対象者の行動を予測する行動予測装置において、前記観測データを前記観測対象者の脳内の複数の電流源における時系列データに写像するフィルタを設定する手段と、設定したフィルタによって写像される時系列データの入力に対して前記
観測対象者が特定の行動を起こす確率を出力する確率モデルを設定する手段と、新たに取得した観測データを前記フィルタにより前記電流源における時系列データに変換する手段と、変換した時系列データを設定した確率モデルに入力することにより前記確率を算出す

10る手段とを備え、算出した確率に基づいて前記観測対象者の行動を予測するようにしてあることを特徴とする。
【００１４】
第８発明に係る行動予測装置は、予測すべき行動の数だけ前記フィルタを設定し、各フィルタの前記確率への寄与度を算出する手段を更に備えることを特徴とする。
【００１５】
第９発明に係る行動予測装置は、前記確率モデルは、学習すべきパラメータを含み、前記観測対象者が前記特定の行動を実行した際に観測される観測データに基づいて前記パラメータを学習する手段を更に備えることを特徴とする。
【００１６】

20

第１０発明に係る行動予測装置は、前記パラメータをベイズ推定法を用いて学習するよ
うにしてあることを特徴とする。
【００１７】
第１１発明に係る行動予測装置は、前記確率モデルは、前記時系列データの入力に応じたスカラー値を出力する判別関数と、該判別関数の出力値に応じて所定値域の値を出力す
るロジスティック関数とにより記述されることを特徴とする。
【００１８】
第１２発明に係る行動予測装置は、前記確率モデルは、スパースロジスティック回帰モ
デルであることを特徴とする。
【００１９】

30

第１発明及び第７発明にあっては、観測データを観測対象者の脳内の複数の電流源における時系列データに写像するフィルタを設定し、設定したフィルタによって写像される時
系列データの入力に対して観測対象者が特定の行動を起こす確率を出力する確率モデルを設定し、新たに取得した観測データを前記フィルタにより電流源における時系列データに変換し、変換した時系列データに基づいて特定の行動を起こす確率を算出するようにしているため、ＥＥＧやＭＥＧなどのセンサ空間における特徴量を利用した従来の予測手法と異なり、脳内電流源を推定した上で行動予測する本願手法では、外乱などのセンサノイズの影響が最小限に抑えられる。
【００２０】
第２発明及び第８発明にあっては、予測すべき行動の数だけフィルタを設定し、各フィ40ルタの確率への寄与度を算出するようにしているため、観測データ自体がある種のカテゴ
リに分けられる場合に当該モデルを使用することができる。
【００２１】
第３発明及び第９発明にあっては、設定する確率モデルは、学習すべきパラメータを含
み、観測対象者が特定の行動を実行した際に観測される観測データに基づいてパラメータを学習するようにしているため、事前学習によって得られるパラメータを用いて観測対象者が実行する行動の予測が可能となる。
【００２２】
第４発明及び第１０発明にあっては、パラメータをベイズ推定法を用いて学習する。パ
ラメータの次元が高次である場合にはモンテカルロ法などの数値解析法は実用的ではない

50

(5)

JP 4978860 B2 2012.7.18

ため、本発明では変分ベイズ推定法が用いられる。
【００２３】
第５発明及び第１１発明にあっては、確率モデルを時系列データの入力に応じたスカラ
ー値を出力する判別関数と判別関数の出力値に応じて所定値域の値を出力するロジスティック関数とにより記述する。すなわち、本発明で用いる確率モデルは、複数のカテゴリに分類できるデータの観測過程をモデル化したものであり、このような回帰モデルを用いる
ことにより、観測データが入力された場合、その観測データが複数のカテゴリのうち何れに属するかを判別することができる。
【００２４】
第６発明及び第１２発明にあっては、設定する確率モデルとしてスパースロジスティッ

10

ク回帰モデルを用いるため、確率モデルの解がスパースネスになるという性質が利用され、使用するデータの次元数が大幅に削減される。
【発明の効果】【００２５】
第１発明及び第７発明による場合は、観測データを観測対象者の脳内の複数の電流源における時系列データに写像するフィルタを設定し、設定したフィルタによって写像される
時系列データの入力に対して観測対象者が特定の行動を起こす確率を出力する確率モデルを設定し、新たに取得した観測データを前記フィルタにより電流源における時系列データに変換し、変換した時系列データに基づいて特定の行動を起こす確率を算出するようにしている。したがって、ＥＥＧやＭＥＧなどのセンサ空間における特徴量を利用した従来の20予測手法と異なり、脳内電流源を推定した上で行動予測する本願手法では、外乱などのセンサノイズの影響を最小限に抑えることができ、行動予測を高精度に行うことができる。
【００２６】
第２発明及び第８発明による場合は、予測すべき行動の数だけフィルタを設定し、各フィルタの確率への寄与度を算出するようにしている。したがって、観測データ自体がある
種のカテゴリに分けられる場合に当該モデルを使用することができる。
【００２７】
第３発明及び第９発明による場合は、設定する確率モデルは、学習すべきパラメータを含み、観測対象者が特定の行動を実行した際に観測される観測データに基づいてパラメー
タを学習するようにしている。したがって、事前学習によって得られるパラメータを用い

30て観測対象者が実行する行動を予測することができる。
【００２８】
第４発明及び第１０発明による場合は、パラメータをベイズ推定法を用いて学習する。
パラメータの次元が高次である場合にはモンテカルロ法などの数値解析法は実用的ではな
いが、本発明では変分ベイズ推定法を用いて確率モデルに含まれるパラメータを求めることができる。
【００２９】
第５発明及び第１１発明による場合は、確率モデルを時系列データの入力に応じたスカ
ラー値を出力する判別関数と判別関数の出力値に応じて所定値域の値を出力するロジステ
ィック関数とにより記述する。すなわち、本発明で用いる確率モデルは、複数のカテゴリ40に分類できるデータの観測過程をモデル化したものであり、このような回帰モデルを用いることにより、観測データが入力された場合、その観測データが複数のカテゴリのうち何れに属するかを判別することができる。
【００３０】
第６発明及び第１２発明による場合は、設定する確率モデルとしてスパースロジスティ
ック回帰モデルを用いている。したがって、確率モデルの解がスパースネスになるという性質が利用され、使用するデータの次元数を大幅に削減することができる。その結果、推定精度が劣化せずに済み、高精度に行動予測を行うことができる。
【発明を実施するための最良の形態】【００３１】50

(6)JP 4978860 B2 2012.7.18

以下、本発明をその実施の形態を示す図面に基づいて具体的に説明する。
図１は脳活動を観測する観測手段としてのブレインキャップを示す模式図である。ブレインキャップ１０は、人間の頭部を覆う帽子状のホルダ１０Ａに、第１センサ１１及び第２センサ１２を複数個（例えば、数十個〜数百個）ずつ設けた構成をなしている。これらの第１センサ１１及び第２センサ１２は等ピッチ（例えば、数ミリメートル間隔）で配置されている。
【００３２】
第１センサ１１は、例えば、脳活動の際に生じる電気活動を非侵襲的に計測するための脳波計（ＥＥＧ）センサである。ブレインキャップ１０に配置された第１センサ１１は、脳波計（ＥＥＧ）を構成し、各センサ１１，１１，…が設置された位置にて脳活動に伴う10

脳磁場の経時的変化を計測して出力する。第１センサ１１は時間分解能に優れており、ミ
リ秒単位での計測が可能である。
【００３３】
第２センサ１２は、例えば、ＮＩＲＳの近赤外センサであって、比較的短い波長を持つ
赤外光を出射する発光素子とその赤外光の反射光を受ける受光素子とが一組として構成されており、発光素子からの出射光の脳内での吸収量に基づいて脳血流の状態を非侵襲的に計測する。各センサ１２，１２，…は、自身が設置された各部位における脳血流を計測して出力する。第２センサ１２は、電場又は磁場のように他の領域からの影響を受けないた
め空間分解能に優れており、数ミリメートル〜数十ミリメートル単位での計測が可能である。

20【００３４】
このような第１センサ１１及び第２センサ１２は、小型の構成であっても脳活動を観測することができるため、前述したようなブレインキャップ１０に簡単に取り付けることができ、大型の構成を必要としない。
【００３５】
本実施の形態では、図１に示したブレインキャップ１０を利用して被験者の脳活動を時系列的に観測し、得られた観測データに基づいて被験者の行動を予測する。なお、本明細書における行動とは、外観から観察可能な動作、反応、行為だけでなく、被験者の精神状態も含むものとする。
【００３６】30

図２はブレインキャップ１０の観測データから被験者の行動を予測する行動予測装置の構成を示す模式図である。被験者の頭部に装着されるブレインキャップ１０には、第１セ
ンサ１１及び第２センサ１２が出力する信号を集計する集計装置２０が接続されている。
集計装置２０には、データ処理装置１００と無線通信を行うための通信装置３０が接続されている。本実施の形態では、ブレインキャップ１０の第１センサ１１及び第２センサ１
２が取得したデータを通信装置３０を通じてデータ処理装置１００へ送信し、データ処理装置１００にて被験者の行動を予測する。
【００３７】
データ処理装置１００は、具体的には、パーソナルコンピュータ、ワークステーション
等であり、ブレインキャップ１０に接続された通信装置３０から送信されるデータの解析40を行う。データ処理装置１００は、ＣＰＵ１０１を備えており、このＣＰＵ１０１には、
ＲＯＭ１０３、ＲＡＭ１０４、ＨＤＤ１０５、通信ＩＦ１０６、入力ＩＦ１０７、出力Ｉ
Ｆ１０８等のハードウェアがバス１０２を介して接続されている。ＣＰＵ１０１には、Ｒ
ＯＭ１０３に予め格納された制御プログラムをＲＡＭ１０４上にロードして実行することにより、前述のハードウェア各部を制御し、全体として本発明に係る行動予測装置として機能させる。
【００３８】
ＨＤＤ１０５は、磁気記録媒体を有する記憶装置であり、通信ＩＦ１０６にて受信した
データ、そのデータに基づく解析結果等を記憶する。入力ＩＦ１０７は、データ処理装置
１００の外部にキーボード、マウス等の入力デバイス１１０を接続するためのインタフェ

50

(7)

JP 4978860 B2 2012.7.18

ースであり、解析開始指示、解析結果の出力指示等の操作を受付ける。出力ＩＦ１０８は、液晶ディスプレイ装置、プリンタ装置等の出力デバイス１２０を接続するためのインタ
フェースであり、受信したデータ、解析結果等を出力する。
【００３９】
以下、データ処理装置１００を用いて被験者の行動を予測する手法について具体的に説明する。本実施の形態では、通信ＩＦ１０６を通じて取得する観測データを訓練データとして判別器（データ処理装置１００）の内部パラメータを決定し、新たな観測データ（テストデータ）が入力された場合、決定したパラメータを用いて設定される確率モデルにより、ある特定の行動が起きる確率を算出することによって行動予測を行う。
【００４０】

10

図３はデータ処理装置１００が実行する処理の手順を示すフローチャートである。デー
タ処理装置１００は、まず、通信ＩＦ１０６を通じて訓練データを取得する（ステップＳ１１）。そして、訓練データを用いて階層ベイズ推定法により脳内電流源を特定する（ステップＳ１２）。脳内電流源の特定は、本願発明者らが特願２００５−１０７９５２明細書で開示している手法を用いることができる。すなわち、特願２００５−１０７９５２明細書では、脳内に存する電流源の時空間分布を第１の分解能で推定する第１領域と第１の分解能より低い第２の分解能で推定する第２領域とを設定し、設定した第１及び第２領域での電流源の時空間分布を、対象者の頭蓋外部にて観測した電気現象に基づいて推定するようにしており、高精度に推定したい領域に重点的にパラメータを分配することによってパラメータの増大に伴う推定精度の劣化を防止することができる。

20【００４１】
前述した推定手法では、最適逆フィルタ推定を行うことにより電流源の推定を行っているが、空間逆フィルタがデータ適合型フィルタであることから判別手法には何通りかのバ
リエーションが存在する。本実施の形態では、こういった条件別フィルタを用いた判別問題に対して混合モデルによる重み付き判別手法について説明する。
【００４２】
判別手法としてスパースロジスティック回帰モデルを用いる。ロジスティック回帰モデ
ルは、２カテゴリカルのデータの観測過程をモデル化した確率モデルであり、特徴量ベク
トルｘ∈Ｒd の時に分類ラベルｙ∈｛０，１｝が観測される確率によってモデリングされる。複数のカテゴリを分類する場合はマルチノミナルモデルを用いて以下の議論を拡張す30ることができる。ｉ番目のサンプルを｛ｙi ，ｘi ｝で表し、Ｎサンプル得られたとする。また、Ｎサンプルまとめたラベル及び特徴量をそれぞれｙ＝｛ｙ1，…，ｙN ｝、Ｘ＝
｛ｘ1 ，…，ｘN ｝で表した場合、ロジスティック回帰モデルは次式で表すことができる。
【００４３】【数１】

40【００４４】
つまり、特徴量ｘi がラベル１のカテゴリに属すると判断される確率ｐi は、特徴量ｘiをパラメータθでモデリングされたスカラー値関数ｆ（ｘi ；θ）によってマッピング

し、さらに数式（１）のσ（ｕ）によって値域（０，１）に写像した値によってモデリングされる。数式（１）のσ（ｕ）はロジスティック関数と呼ばれ、ロジスティック回帰と

50

(8)

JP 4978860 B2 2012.7.18

呼ばれる由縁である。スカラー値関数ｆ（ｘ；θ）は判別関数と呼ばれ、ラベルが分からない新しいデータｘが入ってきたときに、ｆ（ｘ；θ）の値と０との大小によって、ラベル１に属するかラベル０に属するかを決定する際に使われる。実用上よく使われるものは線形判別関数であり、このとき、ｆ（ｘ；θ）＝ｘt θで表される。本実施の形態では、
線形判別関数のみを用いる。非線形関数モデルであるカーネル関数を用いた手法も基本的にはパラメータに関して線形であるので、以下同じ議論は可能である。
【００４５】
パラメータθは、最尤推定法によってＮ個の｛ｙi ，ｘi ｝のサンプルから、数式（１）に示されるＰ（ｙ｜Ｘ；θ）のトータルの観測確率が最大となるように推定される。
【００４６】【数２】

10【００４７】
実際の判別問題では、特徴量とラベルが既知の訓練データとから、パラメータθを最尤法により推定値θハットを計算し、ラベルが分かっていないテストデータｘtestに対して
20、【００４８】【数３】【００４９】と判別する。また、ラベルが１であるという確率は、次式で与えられ、判別の確信度を

30

この確率を用いて表すことができる。
【００５０】【数４】【００５１】
スパースロジスティック回帰モデルは、ベイズ推定版のロジスティック回帰モデルであり、数式（１）に加えて、パラメータθの各成分θi の事前分布として平均０、分散αi‑1iの正規分布を仮定し、さらに、逆分散αi も無情報事前分布によってモデル化する。αは対応するパラメータθi がどれだけ推定に重要であるかを決定する変数であり、αiが小さいほど重要度が高くなる。
【００５２】
スパースロジスティック回帰モデルを数式でまとめた場合、以下の式で表すことができる。
【００５３】40

(9)JP 4978860 B2 2012.7.18

【数５】

10
【００５４】
パラメータθ及びハイパーパラメータα＝｛α1 ，…，αd ｝は、以下の式で表される
ベイズ事後確率が最大をとる値として計算される。
【００５５】【数６】

20【００５６】
しかし、数式（６）の右辺の積分は解析的に行うことが困難であり、モンテカルロ法などの数値解析法、又は変分ベイズ法などの近似解法が必要となる。パラメータの次元が高次のときモンテカルロ法は実用的ではないため、本実施の形態では変分ベイズ法に基づいた推定アルゴリズムを用いる。
【００５７】
この推定アルゴリズムをパラメータが収束するまで適用して得られる推定値θハットは30、判別に有効な要素以外の多くが０となるスパースな解をもつ。これは数式（５）で与えられるベイズ事前分布をαi で周辺化した事前分布Ｐ（θi ）＝θi‑1 （ただし、ｉ＝１，…，Ｄ）の形から分かるようにθi →０の極限で無限の確率をとるような事前分布となっているためである。つまり、この事前分布は尤度関数Ｐ（ｙ｜Ｘ；θ）に対して０の値に高い確率を持つような正規化の役割を果たす。
【００５８】
スパースロジスティック回帰モデルの解がスパースネスになるという応用上、非常に重要な情報を与える。ＥＥＧやＭＥＧデータのように膨大な計測データが入力データとなる
場合、そもそもどの特徴量を使ったら良いか分かっていることは稀である。しかも多くの場合、判別精度は適切な特徴量の選択に依存する。線形判別関数ｆ（ｘ；θ）＝ｘt θ＝40

Σθi ｘi におけるスパースなパラメータ推定は、その形をみて分かるように、推定値０に対応する特徴量を使わないことを意味する。すなわち、スパースロジスティック回帰では、パラメータの推定結果として特徴量ベクトルｘの中から自動的に判別に効果的な特徴
量だけを抜き出しているわけである。したがって、初めにラフに大きな次元の初期特徴量ベクトルを入力し、スパースロジステック回帰によって学習すればよいことになる。その結果、構成される判別器は判別に有効な特徴量だけを用いたものとなる。
【００５９】
次に、混合スパースロジスティック回帰モデルについて説明する。混合スパースロジス
ティックモデルは、複数のスパースロジスティック回帰モデルの混合比による重み付きで判別するモデルである。このモデルは、後述する電流源空間での判別の例のように、特徴

50

(10)

JP 4978860 B2 2012.7.18

量自体がある種のカテゴリに分けられ、各モデルの混合比が計算できるときに有用である。
【００６０】
以下では混合分布をＫ個とし、サンプルｉについての混合比をｇi(1)，…，ｇi(K)とした場合、混合スパースロジスティック回帰モデル（確率モデル）を以下の式により設定する（ステップＳ１３）。
【００６１】【数７】
10【００６２】
ここで、Ｐ（ｙ｜Ｘ(k) ；θ(k) ）は入力特徴量の値をＸ(k)とするスパースロジステ
ィック回帰モデル（数式（５）を参照）である。
【００６３】
以下にコンポーネントの数がＫ個の場合のアルゴリズムについて説明する。
20【００６４】(k)

初期パラメータをθi(k)＝１（ただし、ｉ＝１，…，Ｄ、ｋ＝１，…，(k)Ｋ）とする。事前混合比ｇi

＝０，αi
(1)

，…，ｇi

は、条件ごとに空間フィルタを考える場合、データを再構成したときのフィッティングの良さの割合で評価して設定すればよい。
【００６５】
次いで、データ処理装置１００はパラメータの推定を行う（ステップＳ１４）。このとき、事後混合比は次式により与えられる。
【００６６】【数８】30【００６７】
また、推定値θハットは次式により与えられる。
【００６８】【数９】40【００６９】
最大化は各コンポーネントごとに、勾配とヘシアン行列とを用いてニュートン法によって行う。ここで、勾配ｇ(k) 、ヘシアン行列Ｈ(k) 、最大値における事後分散Ｓ(k) は次

50

(11)

JP 4978860 B2 2012.7.18

式で与えられる。
【００７０】【数１０】10【００７１】
次いで、データ処理装置１００は分散の推定を行う（ステップＳ１５）。逆分散αを以下の式により推定する。

20【００７２】【数１１】【００７３】
次いで、データ処理装置１００はパラメータが収束したか否かを判断し（ステップＳ１30６）、パラメータが収束していないと判断した場合（Ｓ１６：ＮＯ）、処理をステップＳ
１４へ戻し、パラメータの値が収束するまで、パラメータの推定（Ｓ１４）及び分散の推定（Ｓ１５）を行う。
【００７４】
パラメータが収束したと判断した場合（Ｓ１６：ＹＥＳ）、行動予測のためにデータ処理装置１００は脳活動の観測データ（テストデータ）を通信ＩＦ１０６を通じて取得し（ステップＳ１７）、前述のステップＳ１１〜Ｓ１６の処理により得られる確率モデルによりトータル確率を算出する（ステップＳ１８）。そして、算出結果に基づいて行動の予測を行う（ステップＳ１９）。
【００７５】

40

次に、ＭＥＧデータへの適用例について説明する。図４は実験手法を説明する説明図である。１７５０ミリ秒間の実験で、プレキュー時間（２５０ミリ秒）、キュー期間（１０００ミリ秒）、刺激期間（５００ミリ秒）からなる。プレキュー期間では、被験者は特に何もせずにスクリーン上の固定点（例えば、黒い小さなサークル）を見続ける。その後、注意を向けるキューとして、黒いサークルがカラー（Ｃ）又は動き（Ｍ）に変化し、被験者は指示された対象に注意を向ける。実際の実験では、１つのセッションでは常に同じキューを表示するようにしているため、キューは被験者が注意を向けるタイミングの指示として使われると理解できる。１０００ミリ秒のキュー期間の後、色（赤または青）のついたランダムドットが外向き又は内向きに５００ミリ秒だけ動く。被験者は指示によって、注意を向けていた属性について反応する。すなわち、色注意なら赤又は青ということにな

50

(12)

JP 4978860 B2 2012.7.18

り、動き注意なら外向き又は内向きということになる。
【００７６】
解析では、瞬きなどのアーティファクトの入っているトライアルの除去、正常に働いていないセンサの除去、１００Ｈｚのデジタルローパスフィルタによる前処理を行った。実際に解析したデータはセンサ数１８６、トライアル数１３６（色注意のトライアル数を６５、動き注意のトライアル数を７１）、被験者が注意を向けているキュー期間は１０００
ミリ秒である。
【００７７】
時刻ｔにＭ個のセンサにおいて計測したＭＥＧデータをｂt ＝（ｂt,1 ，…，ｂt,M ）tとする。センサ空間における判別では、Ｍ個のセンサ全部又は一部を選択し、その時系

10

列を加工して特徴量をつくる。一方、脳皮質における判別では、電流源推定問題を解くことによって得られるセンサ空間から皮質空間への写像（空間フィルタ）をまず適用し、写像した空間における時系列データを用いる。階層ベイズ法では空間フィルタは線形となり、それを表す行列をＦとした場合、ｊt ＝ＦＢt によって得られる脳活動電流源の時系列
データｊt ＝（ｊt,1 ，…，ｊt,L ）t が得られる。このデータから興味ある脳部位を特定し、時系列データの一部を特徴量として用いる。
【００７８】
次に、センサ空間における判別と電流源空間における判別との比較を行う。センサ空間における判別に用いる特徴量として、後頭に位置する２２のチャンネルにおける次数１０の移動平均フィルタで平滑化した後、１００Ｈｚにダウンサンプリングしたデータ点を用20いた。特徴量が［１，−１］の範囲に収まるように訓練データの平均を引き、最大値で割ることによって正規化を行った。
【００７９】
脳皮質上の特徴量として、例えば、前述した手法により同定された６つの興味ある領域（Ｖ４−Ｌｅｆｔ，Ｖ４−Ｒｉｇｈｔ，ＭＴ−Ｌｅｆｔ，ＭＴ−Ｒｉｇｈｔ，ＳＭＧ−Ｌ
ｅｆｔ，Ｐａｒｉｅｔａｌ−Ｒｉｇｈｔ）において、ピーク強度を持つ皮質頂点を選択した。選ばれた皮質点の電流時系列データは、推定によるばらつきを考慮するために事後分
布の平方根で割り、センサ空間における判別と同じく、時間方向には次数１０の移動平均フィルタで平滑化した後１００Ｈｚにダウンサンプリングした。また、特徴量の値が［１，−１］の範囲に収まるように訓練データの平均を引き、絶対値の最大値で割ることによ30って正規化を行った。
【００８０】
階層ベイズ法では空間フィルタは各条件ごとにＦcolor 、Ｆmotionのように計算されるため、電流源に写像した時点ですでにラベルの情報が含まれる。実際の判別では、ラベルの情報は分かっていないので、どちらの空間フィルタを適用するかそれ自体が問題である。そのため、判別問題を解くときは、判別されるべきＭＥＧデータにＦcolor 、Ｆmotion
をそれぞれ適用した結果である電流源それぞれをダウンサンプリングした特徴量を用いる必要がある。そして、それらを用いた判別方法として、混合分布を用いた方法（混合分布法という）、各条件からの電流特徴量を合わせて１つの特徴量として判別する方法（結合法という）が考えられる。

40【００８１】
ここで、空間フィルタを用いた判別問題における注意点を２つ述べる。１つ目は、空間フィルタの構成には訓練データのみから作る点である。クロスバリデーションでは訓練データをさらに、訓練データとテストデータとに分けるが、空間フィルタも分割された訓練データからのみ構成されたフィルタでない場合には公正なバリデーションができないので注意が必要である。空間フィルタの構成に計算量がかかるときには、空間フィルタを作るためのデータを計測するのも１つの手段である。これは、独立成分分析、主成分分析などの空間フィルタを用いる方法すべてに共通するポイントである。２つ目は、階層ベイズ法のように条件別フィルタを計算する場合、空間フィルタ自体にラベルの情報が含まれるので、判別特徴量を作るときはそれぞれのフィルタを適用した結果それぞれを何らかの形で

50

(13)

JP 4978860 B2 2012.7.18

特徴量として用いなければならないという点である。
【００８２】
以下に、センサ空間における判別、電流源空間における判別（混合分布法）、電流源空間における判別（結合法）の結果を示す。図５は判別結果をまとめた図表である。１３６個のトライアル（サンプル）を１０９個の訓練データと、２７個のテストデータとにランダムに分け、訓練データで判別器のパラメータをスパースロジスティック回帰又はその混合モデルで学習し、テストデータを用いて判別精度を評価した。訓練・テストデータの分け方の影響を少なくするために、１０通りの分け方を行い、それぞれについて判別精度を評価した。また、学習後に０でないパラメータ推定値をもつ特徴量の次元を、学習後の特徴量次元として表示している。

10【００８３】
判別精度では電流源空間において混合モデルを用いた方法が平均８９％の判別精度を有し、最も精度が高い。センサ空間における判別では平均８１％であり、電流源空間における判別でも結合法では、センサ空間と同程度の判別精度（８１％）であった。学習後の特徴量次元は、どの手法でも２０前後であり、学習前の２２００次元（すなわち、２００個のセンサ×時間方向のデータ１１個）に比べて、判別に有効でない多くの次元が除去されていることが分かる。このようにスパースロジスティック回帰では、データから自動的に有効なパラメータのみを推定する。
【００８４】
このようにセンサ空間と電流源空間における判別精度の比較を行った結果、電流源空間

20における判別の方法が、信号の生理学的意味から考えて精度が高くなるという予想を支持する結論が得られた。
【００８５】
なお、本実施の形態では、第１センサ１１及び第２センサ１２が出力するデータを用いて脳内の電流源を推定する構成としたが、例えば、第１センサとして脳磁計（ＭＥＧ）センサを用いることもできる。また、第２センサを同時計測するのではなく、事前にｆＭＲＩを用いて計測しても良い。更にｆＭＲＩなどの装置により得られる磁気共鳴画像により電流源の推定を事前に行う構成としてもよい。ｆＭＲＩによる手法では、タスク実行時と安静時との磁気共鳴画像を比較し、磁気共鳴信号の変化を統計分析することにより、脳の活動部位を推定することができる。

30【図面の簡単な説明】【００８６】【図１】脳活動を観測する観測手段としてのブレインキャップを示す模式図である。
【図２】ブレインキャップの観測データから被験者の行動を予測する行動予測装置の構成を示す模式図である。
【図３】データ処理装置が実行する処理の手順を示すフローチャートである。
【図４】実験手法を説明する説明図である。
【図５】判別結果をまとめた図表である。
【符号の説明】【００８７】
１０

ブレインキャップ

２０

集計装置

３０

通信装置

１００

データ処理装置

１０１

ＣＰＵ

１０２

バス

１０３

ＲＯＭ

１０４

ＲＡＭ

１０５

ＨＤＤ

１０６

通信ＩＦ

40

50

(14)
１０７

入力ＩＦ

１０８

出力ＩＦ

１１０

入力デバイス

１２０

出力デバイス【図１】【図２】

JP 4978860 B2 2012.7.18

(15)
【図３】

【図５】

【図４】

JP 4978860 B2 2012.7.18

(16)JP 4978860 B2 2012.7.18

フロントページの続き(51)Int.Cl.

ＦＩ

Ａ６１Ｂ 10/00(2006.01)(72)発明者 吉岡 琢
京都府相楽郡精華町光台二丁目２番地２ 株式会社国際電気通信基礎技術研究所内
審査官 門田 宏
10(56)参考文献 特開２００６−３０５３３４（ＪＰ，Ａ）
特開２００６−２８０８０６（ＪＰ，Ａ）
特表２００２−５２６１８９（ＪＰ，Ａ）
特開２００６−２６３４６０（ＪＰ，Ａ）
川脇大 外５名，fMRI/MEG統合解析による心的な指標運動追跡に関わる皮質電流源の推定，電子情報通信学会技術研究報告，日本，社団法人電子情報通信学会，２００６年１０月 ４日，Vol.
106、No.279，pp.19‑24(58)調査した分野(Int.Cl.，ＤＢ名)
Ａ６１Ｂ

５／０４

Ａ６１Ｂ

５／１４５

Ａ６１Ｂ

− ５／０５

１０／００

ＪＳＴＰｌｕｓ／ＪＭＥＤＰｌｕｓ／ＪＳＴ７５８０（ＪＤｒｅａｍＩＩ）

20