(57)【要約】【課題】知的活動評価のために好適な主指標を選択する。【解決手段】分布推定部１１は、生理指標や観測指標である複数の指標の値と、人間の知的活動を示すラベルとの時系列データの入力を受け、指標または指標の組とラベルとの分布を示す関数を作成する。情報量計算部１２は、分布推定部１１により生成された分布を示す関数から、指標または指標の組とラベルとの情報量を計算する。指標選択部１３は、所望のラベルを最もよく反映し、当該ラベルの推定に最適な主指標として、複数の指標または指標の組の中から、最大の情報量を示す指標または指標の組を選択し、出力する。【選択図】図１
2010022649.tif 000002
【特許請求の範囲】
【請求項１】
人間の生理反応または行動を数値化した複数の種類の指標の中から、人間の知的活動を示す所望のラベルを最もよく反映する主指標を選択する指標選択装置であって、
複数の指標の時系列の値と、当該指標が観測されたときのラベルとの情報から、当該指標又は当該指標の組みと当該ラベルとの分布を得る分布推定部と、
前記分布推定部により得られた前記指標又は前記指標の組みと前記ラベルとの分布から、当該指標または当該指標の組みと、当該ラベルとの組み合わせについての情報量を算出する情報量計算部と、
前記情報量計算部により算出された情報量に応じて前記ラベルの主指標となる前記指標又は前記指標の組みを選択する指標選択部と、
を備えることを特徴とする指標選択装置。
【請求項２】
前記指標または前記指標の組みによって知的活動を評価した場合に期待される誤差を算出する誤差見積部をさらに備え、
前記指標選択部は、前記誤差見積部により算出された誤差に基づいて前記ラベルの主指標となる前記指標又は前記指標の組みを選択する、
ことを特徴とする請求項１に記載の指標選択装置。
【請求項３】
前記情報量計算部によって算出された情報量に基づいて、知的活動の評価の難易度が所定より高い前記指標または前記指標の組みを特定する難易度判定部をさらに備え、
前記指標選択部は、前記難易度判定部により難易度が所定より高いと判断された前記指標または前記指標の組みを除いた前記指標または前記指標の組みを主指標として選択する、
ことを特徴とする請求項１に記載の指標選択装置。
【請求項４】
指標の時系列の値と、当該指標が観測されたときのラベルとの情報を記憶する記憶部と、
前記難易度判定部において、現在の前記指標または前記指標の組み全てについて、知的活動の評価の難易度が所定より高いと判断された場合に新たな前記指標の組みを生成し、前記分布推定部に出力する指標組合せ部とをさらに備え、
前記分布推定部は、前記指標組合せ部から新たな前記指標の組みの入力を受け、前記記憶部内の指標の時系列の値と、当該指標が観測されたときのラベルとの情報を参照して、当該新たな指標の組みと当該ラベルとの分布を算出する、
ことを特徴とする請求項３に記載の指標選択装置。
【請求項５】
人間の生理反応または行動を数値化した複数の種類の指標の中から、人間の知的活動を示す所望のラベルを最もよく反映する主指標を選択する指標選択装置に用いられる指標選択方法であって、
分布推定部が、複数の指標の時系列の値と、当該指標が観測されたときのラベルとの情報から、当該指標又は当該指標の組みと当該ラベルとの分布を得る分布推定過程と、
情報量計算部が、前記分布推定過程において得られた前記指標又は前記指標の組みと前記ラベルとの分布から、当該指標または当該指標の組みと、当該ラベルとの組み合わせについての情報量を算出する情報量計算過程と、
前記情報量計算過程において算出された情報量に応じて前記ラベルの主指標となる前記指標又は前記指標の組みを選択する指標選択過程と、
を有することを特徴とする指標選択方法。
【請求項６】
誤差見積部が、前記指標または前記指標の組みによって知的活動を評価した場合に期待される誤差を算出する誤差見積過程をさらに有し、
前記指標選択過程では、前記誤差見積過程において算出された誤差に基づいて前記ラベルの主指標となる前記指標又は前記指標の組みを選択する、
ことを特徴とする請求項５に記載の指標選択方法。
【請求項７】
難易度判定部が、前記情報量計算過程において算出された情報量に基づいて、知的活動の評価の難易度が所定より高い前記指標または前記指標の組みを特定する難易度判定過程をさらに有し、
前記指標選択過程においては、前記難易度判定過程において難易度が所定より高いと判断された前記指標または前記指標の組みを除いた前記指標または前記指標の組みを主指標として選択する、
ことを特徴とする請求項５に記載の指標選択方法。
【請求項８】
指標組合せ部が、前記難易度判定過程において、現在の前記指標または前記指標の組み全てについて、知的活動の評価の難易度が所定より高いと判断された場合に新たな前記指標の組みを生成し、前記分布推定部に出力する指標組合せ過程と、
前記分布推定部が、前記指標組合せ過程において生成された新たな前記指標の組みの入力を受け、記憶部内の指標の時系列の値と、当該指標が観測されたときのラベルとの情報を参照して、新たな前記指標の組みと当該ラベルとの分布を算出する分布再推定過程とをさらに有し、
前記情報量計算過程においては、前記分布再推定過程において得られた新たな前記指標の組みと前記ラベルとの分布から、新たな前記指標の組みと、当該ラベルとの組み合わせについての情報量を算出し、
前記指標選択過程においては、前記情報量計算過程において算出された情報量に応じて前記ラベルの主指標となる新たな前記指標の組みを選択する、
ことを特徴とする請求項７に記載の指標選択方法。
【請求項９】
請求項１から請求項４のいずれかの項に記載の指標選択装置としてコンピュータを動作させることを特徴とするコンピュータプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、人の知的活動の推定を行うために最適な指標又は指標の組を選択する指標選択装置及び方法ならびにコンピュータプログラムに関する。
【背景技術】
【０００２】
知的活動とは、興味、感情、あるいは印象等の人間の心理的な活動や、認知、記憶あるいは追想等の人間の情報処理活動などを差す。
【０００３】
一方、人間の生理反応は、脳波、心拍変動、血圧、血流、発汗、皮膚電気反応、呼吸、血中酸素飽和度、顔面温度、又は唾液分泌成分等として数値化できる。この数値化された値を総じて生理指標と呼ぶ。
また、人間の行動は、声の調子、非言語表現（溜息等）、視線停留時間、表情変化、又は歩行パターン等として数値化できる。この数値化された値を総じて観測指標と呼ぶ。
このような生理指標、観測指標を算出する測定器は小型で携帯することができるものもある。
【０００４】
従来、人間の知的活動に関するアンケート調査については、十分なテストを行った世界的標準試験法が提案されている。例えば、精神的作業負荷の検査法として、ＮＡＳＡ−ＴＬＸ（Task Load Index）と呼ばれる試験が広く利用されており、このテストによって、知的・知覚的要求、身体的要求、タイムプレッシャー、作業成績、努力、フラストレーションの６つの尺度が得られる（例えば、非特許文献１参照）。
【０００５】
ある信号が得られた時、その分布関数を推定することができる。最も単純には、得られた信号の長さをＮとした時、信号の最小の値から最大の値までを定義域にもつ関数Ｆ（ｘ）を、信号の中でｘよりも小さいものの数を信号の長さＮで除算した値を取るように決めることで得られる。
そこで、ある指標や信号の分布と、生理指標や観測指標の分布とを用いて計算可能な情報量（例えば、非特許文献２参照）を用いることにより、当該指標や信号を最もよく反映する生理指標又は観測指標を知ることができる。
【０００６】
また、シャノンの条件付情報エントロピやフィッシャー情報量により、各指標による知的活動評価の誤り率や二乗誤差の下限値が与えられる（例えば、非特許文献３参照）。
【０００７】
また、指標とラベルの組である標本を機械的に学習器へ入力することにより、指標とラベルの関係を自動分類するルールを自動学習することができる場合がある。このルールを判別ルールと呼ぶ。無関係な指標とラベルの組により標本を構成した場合はこの限りではなく、入力する標本を増やしても判別ルールが定まらない場合がある。
このラベルとして知的活動の状態を示すラベルを用いることにより、知的活動の評価を行う判別ルールを学習することができる。
【０００８】
判別ルールを用いて任意の指標へラベルを付与することができる。このような装置を判別器と呼ぶ。
そして、学習に用いる指標として複数の異なる指標を組み合わせて標本を構成し学習することによって、判別器の精度を向上させることができる場合がある。しかし、全く意味のない指標を組み合わせることで精度が悪化する場合もある。そのような学習器のアルゴリズムとして決定木を自動構築するＣＡＲＴ（classification and regression trees：識別・回帰木）等がある（例えば、非特許文献４参照）。
【非特許文献１】中園  薫、秦  崇洋、安宅  恵理、田中  久弥、長崎  祐二，「ゲーム課題の遂行時における作用負荷の生理指標および心理指標による評価」，信学技報，HIP2007-163(2008-03)，ｐｐ．２５−３０，２００７年
【非特許文献２】Takahiro HADA、Toyonori MUNAKATA，“Filtering for Simple Threshold Systems: Self-Tuning, Mutual Information and Applications” ，IEICE Transactions on Fundamentals of Electronics, Communications and Computer Sciences，Vol. E89-A, No.10，pp. 2566-2574，2006
【非特許文献３】村田昇，“情報理論の基礎”，サイエンス社，日本，２００５年
【非特許文献４】Richard O. Duda、Peter E. Hart、David G. Stork，尾上守夫（監訳），“パターン識別”，第２版，新技術コミュニケ一ションズ，日本，２００１年
【発明の開示】
【発明が解決しようとする課題】
【０００９】
ある知的活動をよく反映する指標（主指標）が分かっている場合、情報量を利用することで、当該指標を最もよく反映する生理指標又は観測指標（副指標）を調べることができ、主指標に代わり副指標をもって知的評価を行うことができる（先願発明：特願２００８−０４６０１４）。
【００１０】
しかしながら、知的活動をよく反映する指標が不明の場合は、主指標を決定することができなかった。
また主指標が分かっていたとしても、生理指標や観測指標のようにノイズや誤りが多く含まれる指標間の学習を行う場合には、ＣＡＲＴのようなアルゴリズムでは過学習してしまう可能性があり、適切な学習を行うためには人間の判断やチューンが必要であるため手間が非常に大きかった。
更に、個人差のある生理指標や観測指標の学習結果がどの程度妥当であるかを判断することが難しかった。そのため、当該人が知的活動をどの程度、生理指標や観測指標に表出しているかが分からず、よって、学習することで知的活動が判別可能であるのかそうでないのかも分からなかった。
【００１１】
本発明は、このような事情を考慮してなされたもので、その目的は、知的活動評価のために好適な指標等を選択するための定量的な尺度を与えることにより、主指標として何を用いるべきかを選択することを可能とする指標選択装置及び方法ならびにコンピュータプログラムを提供することにある。
【課題を解決するための手段】
【００１２】
本発明は、人間の生理反応または行動を数値化した複数の種類の指標の中から、人間の知的活動を示す所望のラベルを最もよく反映する主指標を選択する指標選択装置であって、複数の指標の時系列の値と、当該指標が観測されたときのラベルとの情報から、当該指標又は当該指標の組みと当該ラベルとの分布を得る分布推定部と、前記分布推定部により得られた前記指標又は前記指標の組みと前記ラベルとの分布から、当該指標または当該指標の組みと、当該ラベルとの組み合わせについての情報量を算出する情報量計算部と、前記情報量計算部により算出された情報量に応じて前記ラベルの主指標となる前記指標又は前記指標の組みを選択する指標選択部と、を備えることを特徴とする指標選択装置である。
【００１３】
また、本発明は、上述した指標選択装置において、前記指標または前記指標の組みによって知的活動を評価した場合に期待される誤差を算出する誤差見積部をさらに備え、前記指標選択部は、前記誤差見積部により算出された誤差に基づいて前記ラベルの主指標となる前記指標又は前記指標の組みを選択する、ことを特徴とする。
【００１４】
また、本発明は、上述した指標選択装置において、前記情報量計算部によって算出された情報量に基づいて、知的活動の評価の難易度が所定より高い前記指標または前記指標の組みを特定する難易度判定部をさらに備え、前記指標選択部は、前記難易度判定部により難易度が所定より高いと判断された前記指標または前記指標の組みを除いた前記指標または前記指標の組みを主指標として選択する、ことを特徴とする。
【００１５】
また、本発明は、上述した指標選択装置において、指標の時系列の値と、当該指標が観測されたときのラベルとの情報を記憶する記憶部と、前記難易度判定部において、現在の前記指標または前記指標の組み全てについて、知的活動の評価の難易度が所定より高いと判断された場合に新たな前記指標の組みを生成し、前記分布推定部に出力する指標組合せ部とをさらに備え、前記分布推定部は、前記指標組合せ部から新たな前記指標の組みの入力を受け、前記記憶部内の指標の時系列の値と、当該指標が観測されたときのラベルとの情報を参照して、当該新たな指標の組みと当該ラベルとの分布を算出する、ことを特徴とする。
【００１６】
また、本発明は、人間の生理反応または行動を数値化した複数の種類の指標の中から、人間の知的活動を示す所望のラベルを最もよく反映する主指標を選択する指標選択装置に用いられる指標選択方法であって、分布推定部が、複数の指標の時系列の値と、当該指標が観測されたときのラベルとの情報から、当該指標又は当該指標の組みと当該ラベルとの分布を得る分布推定過程と、情報量計算部が、前記分布推定過程において得られた前記指標又は前記指標の組みと前記ラベルとの分布から、当該指標または当該指標の組みと、当該ラベルとの組み合わせについての情報量を算出する情報量計算過程と、前記情報量計算過程において算出された情報量に応じて前記ラベルの主指標となる前記指標又は前記指標の組みを選択する指標選択過程と、を有することを特徴とする指標選択方法である。
【００１７】
また、本発明は、上述した指標選択方法において、誤差見積部が、前記指標または前記指標の組みによって知的活動を評価した場合に期待される誤差を算出する誤差見積過程をさらに有し、前記指標選択過程では、前記誤差見積過程において算出された誤差に基づいて前記ラベルの主指標となる前記指標又は前記指標の組みを選択する、ことを特徴とする。
【００１８】
また、本発明は、上述した指標選択方法において、難易度判定部が、前記情報量計算過程において算出された情報量に基づいて、知的活動の評価の難易度が所定より高い前記指標または前記指標の組みを特定する難易度判定過程をさらに有し、前記指標選択過程においては、前記難易度判定過程において難易度が所定より高いと判断された前記指標または前記指標の組みを除いた前記指標または前記指標の組みを主指標として選択する、ことを特徴とする。
【００１９】
また、本発明は、上述した指標選択方法において、指標組合せ部が、前記難易度判定過程において、現在の前記指標または前記指標の組み全てについて、知的活動の評価の難易度が所定より高いと判断された場合に新たな前記指標の組みを生成し、前記分布推定部に出力する指標組合せ過程と、前記分布推定部が、前記指標組合せ過程において生成された新たな前記指標の組みの入力を受け、記憶部内の指標の時系列の値と、当該指標が観測されたときのラベルとの情報を参照して、新たな前記指標の組みと当該ラベルとの分布を算出する分布再推定過程とをさらに有し、前記情報量計算過程においては、前記分布再推定過程において得られた新たな前記指標の組みと前記ラベルとの分布から、新たな前記指標の組みと、当該ラベルとの組み合わせについての情報量を算出し、前記指標選択過程においては、前記情報量計算過程において算出された情報量に応じて前記ラベルの主指標となる新たな前記指標の組みを選択する、ことを特徴とする。
【００２０】
また、本発明は、指標選択装置としてコンピュータを動作させることを特徴とするコンピュータプログラムである。
【発明の効果】
【００２１】
本発明によれば、人の知的活動などの所望のラベルの評価を行なうにあたって、定量的な尺度に基づいた最適な指標（例えば、生理指標）を選択することができる。また、本発明によれば、期待される学習誤差を学習の前に見積もることが可能である。加えて、人の知的活動評価の難易度を見積もることを可能とする。
【発明を実施するための最良の形態】
【００２２】
以下、本発明の実施形態を、図面を参照して説明する。
【００２３】
［第１の実施形態］
図１は、本発明の一実施形態による指標選択装置１の構成を示すブロック図である。図１に示すように、本発明の実施形態の指標選択装置１は、分布推定部１１と、情報量計算部１２と、指標選択部１３とから構成される。
【００２４】
分布推定部１１は、指標または指標の組のそれぞれとラベルとの分布関数、又は確率関数、あるいは確率密度関数を推定する。確率変数Ｘがｘ以下の値をとる確率を確率変数Ｘの分布関数という。また、確率関数は、確率変数Ｘに対応した確率を示す関数である。確率密度関数は、分布関数の導関数である。本実施形態では、確率変数として、指標または指標の組みと、１つのラベルとの組み合わせについての確率変数を用いる。
以下、指標または指標の組みを「指標等」と記載する。ここで指標とは、脳波、心拍変動、血圧、血流、発汗、皮膚電気反応、呼吸、血中酸素飽和度、顔面温度、又は唾液分泌成分等として人間の生理反応を数値化した生理指標や、声の調子、非言語表現（溜息等）、視線停留時間、表情変化、又は歩行パターン等として人間の行動を数値化した観測指標である。また、ラベルとは、興味、感情、あるいは印象等の人間の心理的な活動や、認知、記憶あるいは追想等の人間の情報処理活動など、人の知的活動の状態を示すものである。
【００２５】
情報量計算部１２は、分布推定部１１により推定された関数により示される分布に従って、情報量を計算する。情報量とは、ある事象の起こりにくさを示す尺度であり、例えば相互情報量や、シャノンの条件付情報エントロピやフィッシャー情報量等である。相互情報量とは、複数の確率変数についての相互依存の尺度を表す量である。また、条件付情報エントロピとは、ある事象が生じているという条件下における、他の事象の情報量であり、フィッシャー情報量とは、スコアの分布である。
指標選択部１３は、予め与えられた条件から最適な指標等を選択する。
【００２６】
図２は、図１に示す指標選択装置１の処理フローを示す図である。同図において、同一試行中に、複数の指標等の数値とラベルの情報（ＩＮＦ１）が分布推定部１１へ入力される（ステップＳ１１０）。このとき、観測された指標等の数値と、その数値が観測されたときのラベルとが対応可能な情報として入力される。例えば、時系列の各指標等の値と、時系列のラベルが入力される。分布推定部１１は、入力された指標等の値とラベルの情報に基づいて、各指標等とラベルとの結合累積分布関数を計算し、その分布関数をそのまま情報量計算部１２に出力するか、あるいは計算された結合累積分布関数に対して適当に補間を行なった結合累積分布関数を生成し、これを分布関数の情報（ＩＮＦ２）として情報量計算部１２に出力する（ステップＳ１２０）。
【００２７】
ここで、まず、累積分布関数について説明する。例えば、ｘ＝｛４０、４０．５、４１，４１、３９、４０．５｝というデータ列が与えられた場合、ｘ≦３８なるデータは０個、ｘ≦３９なるデータは１個、ｘ≦４０なるデータは２個、ｘ≦４０．５なるデータは４個、ｘ≦４１なるデータは６個であり、これらの数を標本数（ｘの個数）６で除算した値を累積分布関数の値とする。
【００２８】
上記のように、累積分布関数は１つの変数について表したものであるが、本実施形態では、複数の変数についての累積分布関数を使用するため、これを結合累積分布関数と定義する。結合累積分布関数の複数の変数としては、指標等とラベルとが含まれることになる。そこで、まず、分布推定部１１は、指標等の各値が観測された時間、または、所定の周期的な時間における各指標等の値と、対応するラベルとの組み合わせを生成する。そして、例えば、この組み合わせが指標ｙ，ラベルｚである場合、（指標ｙ，ラベルｚ）＝｛（ｙ１，ｚ１）、（ｙ２，ｚ２）、…、（ｙｎ，ｚｎ）｝が生成される。このとき、結合累積分布関数は、（ａ，ｂ）については、ｙ≦ａかつｚ≦ｂを満たす（指標ｙ，ラベルｚ）の個数を、標本数ｎで除算した値をとるような関数となる。また、各指標等の値と、対応するラベルとの組み合わせが、指標ｙと指標ｗの組み，ラベルｚである場合、（指標ｙ，指標ｗ，ラベルｚ）＝｛（ｙ１，ｗ１，ｚ１）、（ｙ２，ｗ２，ｚ２）、…、（ｙｎ，ｗｎ，ｚｎ）｝が生成される。このとき、結合累積分布関数は、（ａ，ｃ，ｂ）については、ｙ≦ａかつｗ≦ｃかつｚ≦ｂを満たす（指標ｙ，指標ｗ，ラベルｚ）の個数を、標本数ｎで除算した値をとるような関数となる。
【００２９】
なお、この処理によって得られる分布関数は階段関数であるが、最も単純な例としては、そのまま階段関数として用いてもよい。あるいは、得られた階段関数に平滑化の処理を施して滑らかな関数として用いてもよい。
【００３０】
情報量計算部１２は、分布推定部１１から出力された複数の指標等についての結合累積分布関数から、目的に応じてシャノンの情報エントロピ、条件付情報エントロピ、相互情報量、又はフィッシャー情報量を算出し、それぞれの指標等とその指標等に対応した当該情報量と（ＩＮＦ３）を指標選択部１３に出力する（ステップＳ１３０）。例えば、情報量計算部１２は、分布関数Ｆ（ｘ）の差分関数であるｆ（ｘ）＝Ｆ（ｘ）−Ｆ（ｘ−Δｘ）を用いて、エントロピ−ｆ（ｘ）×ｌｏｇｆ（ｘ）をｘの全範囲で積分することによって情報エントロピを算出し、情報量として出力する。ただし、実際には、条件付分布関数を用いて計算する。
【００３１】
指標選択部１３は、情報量計算部１２から入力された指標等と情報量の尺度を元に、所望のラベルついて最も情報量の高い指標等を、そのラベルの主指標として選択し、選択した指標等とその情報量（ＩＮＦ４）をディスプレイに表示したり、印刷したり、あるいは、コンピュータ読み書き可能な記録媒体に書き込むなどして出力する（ステップＳ１４０）。
【００３２】
図３は、図１に示した構成において、図２の処理を行うときに用いられる各情報ＩＮＦ１〜４の具体例を示すものである。各情報ＩＮＦ１〜４は、以下を示している。なお、ここでは、指標を複数の生理指標、知的活動ラベルを興味の有無とした場合を例にしている。
【００３３】
ＩＮＦ１：指標選択装置１ヘ入力される指標の値とラベルである。
ＩＮＦ２：分布推定部１１により推定された指標等とラベルの分布である。
ＩＮＦ３：情報量計算部１２により計算された各指標等とラベルの情報量である。
ＩＮＦ４：指標選択部１３により選択された指標等である。
【００３４】
図２のステップＳ１１０において、分布推定部１１は、入力された時系列の指標の値及びラベルから、各指標等及びラベルについてのヒストグラムを作成し、このヒストグラムから結合累積分布関数を作成する。図３（ａ）は、分布推定部１１へ入力された時系列の指標の値とラベル（ＩＮＦ１）である。
【００３５】
指標の時系列の入力データとしては、例えば、被験者に装着した生理センサ等を用いて人間の生理現象を測定することによって得られた時系列のデータを用いることができる。具体例としては、１ｋＨｚ（キロヘルツ）のサンプリングレートでサンプリングされている心電図の電圧値の時系列データなどが挙げられる。
【００３６】
また、ラベルの時系列データとしては、例えば、生理センサを装着した前記被験者に、ボタンなどの入力装置を持たせて外的な刺激等を与え、所定の知的活動の感覚を得たときに、その入力装置を操作するなどして時系列データを取得することができる。具体的には、ある映像を被験者に見せ、怖いと思ったときや、面白いと思ったときにそのラベル（怖い、面白い）に応じたボタンを押し続ける、などである。ここでは、例えば、被験者が興味をもったときに入力装置を操作し続けるようにし、入力装置が操作されているときには興味有のラベル、操作されていないときには興味無のラベルとしている。あるいは、予めその映像が面白いものである、怖いものである、などという情報を入力することでもよい。この場合、ラベル時系列データとしては、実験の最初から最後まで一定してラベルが付与されたものとなる。なお、これに限らず、任意の方法によりラベルの時系列データを入力することが可能である。
【００３７】
図３（ａ）に示す入力データＩＮＦ１の数値（単位：個）は、実際の観測値ではなく、サンプル数を示している。例えば、１ｋＨｚのサンプリングレートで１分間データを収集した場合、指標値は６００００点でサンプリングされるため、６００００個のデータとなる。実際には、データは観測された時点で逐次入力し、処理に必要な数の入力データが蓄積された時点で（例えば、５分間分の入力データが蓄積された時点で）処理を開始するようにすればよい。
なお、心拍間隔（ＲＲＩ）については、データ数が、センサのサンプリングレートによってではなく、利用者の心拍数によって決まるため、数値に「約」をつけている。
【００３８】
分布推定部１１は、図２のステップＳ１２０における処理結果として、図３（ｂ）に示すような、推定された各指標等とラベルの結合累積分布関数（ＩＮＦ２）を情報量計算部１２へ出力する。
図２のステップＳ１３０において、情報量計算部１２では、ＩＮＦ２により示される結合累積分布関数Ｆ（Ｘ，Ｌ）それぞれから、指標等の周辺分布関数Ｐ（Ｘ）とラベルの周辺分布関数Ｐ（Ｌ）を算出し、相互情報量Ｉ（Ｘ；Ｌ）を計算する（Ｘは指標等、Ｌはラベル）。図７（ｃ）は、情報量計算部１２により計算された相互情報量（ＩＮＦ３）であり、指標選択部１３へ入力される。
図２のステップＳ１４０において、指標選択部１３は、ラベルの推定に最適な指標等である主指標として、最大の相互情報量を示す指標等である「血流量」を選択し、出力する。図７（ｄ）は、指標選択部１３により出力される、選択された指標等（ＩＮＦ４）である。
【００３９】
上述するように、本実施形態による指標選択装置１と、先願発明との大きな差異は、主指標が与えられなくても適当な数の入力から作成された分布関数によって、指標等の善し悪しを判定することが可能になることにある。
【００４０】
［第２の実施形態］
図４は、本発明の第２の実施形態による指標選択装置１ａの構成を示すブロック図である。同図に示すように、本実施形態の指標選択装置１ａは、分布推定部１１と、情報量計算部１２と、誤差見積部１４と、指標選択部１３とから構成される。分布推定部１１と、情報量計算部１２と、指標選択部１３は第１の実施形態による指標選択装置１と同様であるため、説明を省略する。誤差見積部１４は、与えられた情報量の種類と値から、誤り率や推定値の分散を計算するものである。
【００４１】
図５は、図４に示す指標選択装置１ａの処理フローを示す図である。同図において、ステップＳ２１０〜Ｓ２３０までの処理は、図２に示す第１の実施形態のステップＳ１１０〜Ｓ１３０と同様の処理を行う。また、ＩＮＦ１〜ＩＮＦ３は第１の実施形態と同様の情報である。
【００４２】
誤差見積部１４は、ステップＳ２３０において情報量計算部１２により計算された情報量に応じて、誤り率や推定値の分散を計算し、指標選択部１３へ出力する（ステップＳ２４０）。例えば、誤差見積部１４は、エントロピの計算式−ｆ（ｘ）×ｌｏｇｆ（ｘ）の逆演算を行なうことによって誤り率を計算し、誤差情報として出力する。相互情報量やシャノンの条件付エントロピの場合、誤差見積部１４は、入力された指標等の値とラベルとから、当該指標等によるラベル推定の誤り率を計算し、指標選択部１３に出力する。また、フィッシャー情報量の場合、誤差見積部１４は、ラベルが何らかの連続値で表現できるならば、当該指標等によるラベルを表現する連続値の推定値の分散を計算し、指標等とそれらに期待される誤差情報とを指標選択部１３に出力する。
【００４３】
指標選択部１３は、誤差見積部１４から入力された指標等と、それらに期待される誤差情報とを元に、所望のラベルついて最も誤差が小さいと期待される指標等を、そのラベルの主指標として選択し、選択した指標等を誤差情報と共に出力する（ステップＳ２５０）。
【００４４】
図６は、図４に示した構成において、図５の処理を行うときに用いられる各情報ＩＮＦ１〜３、５、６の具体例を示すものである。各情報ＩＮＦ１〜３、５、６は、以下を示している。なお、ここでは、指標を複数の生理指標、知的活動ラベルを興味の有無とした場合を例にしている。
【００４５】
ＩＮＦ１：指標選択装置１ａヘ入力される指標の値とラベルである。
ＩＮＦ２：分布推定部１１により推定された指標等とラベルの分布である。
ＩＮＦ３：情報量計算部１２により計算された各指標等とラベルの情報量である。
ＩＮＦ５：誤差見積部１４により計算された、指標等からラベルを推定する際の最低誤り率である。
ＩＮＦ６：指標選択部１３により選択された指標等とその最低誤り率である。
【００４６】
図５のステップＳ２２０において、分布推定部１１は、ステップＳ２１０において入力された時系列の指標の値及びラベルから、各指標等及びラベルについてのヒストグラムを作成し、これを指標等とラベルとの結合確率関数とする。図６（ａ）は、分布推定部１１に入力された時系列の指標の値とラベル（ＩＮＦ１）であり、第１の実施形態におけるＩＮＦ１と同様の方法により収集されたものである。図６（ｂ）は、分布推定部１１により推定され、情報量計算部１２へ入力される結合確率関数（ＩＮＦ２）である。
【００４７】
図５のステップＳ２３０において、情報量計算部１２は、ＩＮＦ２で示される結合確率関数Ｐ（Ｘ，Ｌ）それぞれから指標等の周辺確率関数Ｐ（Ｘ）を算出し、これらから条件付確率関数Ｐ（Ｌ｜Ｘ）を算出する（Ｌはラベル、Ｘは指標等）。図６（ｃ）は、情報量計算部１２により計算された条件付情報エントロピ（ＩＮＦ３）であり、誤差見積部１４へ入力される。
【００４８】
図５のステップＳ２４０において、誤差見積部１４は、条件付情報エントロピから最低誤り率として、最適な推定をした時のビット誤り率を計算する。図６（ｄ）は、誤差見積部により計算された最低誤り率（ＩＮＦ５）であり、指標選択部１３へ入力される。図５のステップＳ２５０において、指標選択部１３は、ラベルの推定に最適な指標等である主指標として、最低誤り率の最も小さい「血流量」を選択し、出力する。図６（ｅ）は、指標選択部１３により出力される、選択された指標等と誤差見積（ＩＮＦ６）である。
【００４９】
上述するように、本実施形態による利点は、第１の実施形態による先願発明との差異に加えて、指標等とラベルとの学習を行う前に、その指標等の理論的に達成可能な誤り率の下限や分散の小ささを、入力から近似的に与えられることにある。
【００５０】
［第３の実施形態］
図７は、本発明の第３の実施形態による指標選択装置１ｂの構成を示すブロック図である。同図に示すように、本実施形態の指標選択装置１ｂは、分布推定部１１と、情報量計算部１２と、難易度判定部１５と、指標選択部１３とから構成される。分布推定部１１と、情報量計算部１２と、指標選択部１３は第１の実施形態による指標選択装置１と同様であるため、説明を省略する。難易度判定部１５は、与えられた情報量の種類と値から、当該人の知的活動評価の難易度を判定するものである。
【００５１】
図８は、図７に示す指標選択装置１ｂの処理フローを示す図である。同図において、ステップＳ３１０〜Ｓ３３０までの処理は、図２に示す第１の実施形態のステップＳ１１０〜Ｓ１３０と同様の処理を行う。また、ＩＮＦ１〜ＩＮＦ３は第１の実施形態と同様の情報である。
【００５２】
難易度判定部１５は、ステップＳ３３０において情報量計算部１２により計算された情報量に応じて、指標等による当該人の知的活動評価の難易度を判定する。例えば、ある指標等について、相互情報量が所定の要求される情報量よりも小さい場合は難易度が高いと判断し、その指標等を情報量計算部１２より入力された情報から排除した結果（ＩＮＦ７）を指標選択部１３へ出力する（ステップ３４０）。要求を満たす指標等がある場合は、その要求を満たす指標等の全てを出力し、要求を満たす指標等がない場合は、要求を満たす指標等がない旨を示す情報（ＩＮＦ７）を指標選択部１３へ出力する。
【００５３】
指標選択部１３は、要求を満たす指標等があった場合、その要求を満たす指標等を全てと、その指標等の情報量（ＩＮＦ８）を出力する（ステップＳ３５０）。当該指標等がなかった場合、その旨を示す情報（ＩＮＦ８）を出力する。
【００５４】
図９は、図７に示した構成において、図８の処理を行うときに用いられる各情報ＩＮＦ１〜３、７、８の具体例を示すものである。各情報ＩＮＦ１〜３、７、８は、以下を示している。なお、ここでは、指標を複数の生理指標、知的活動ラベルを興味の有無とした場合を例にしている。
【００５５】
ＩＮＦ１：指標選択装置１ｂヘ入力される指標の値とラベルである。
ＩＮＦ２：分布推定部１１により推定された指標等とラベルの分布である。
ＩＮＦ３：情報量計算部１２により計算された各指標等とラベルの情報量である。
ＩＮＦ７：難易度判定部１５により判定された難易度である。
ＩＮＦ８：指標選択部１３により選択された指標等である。
【００５６】
図８のステップＳ３２０において、分布推定部１１は、ステップＳ３１０において入力された時系列の指標の値及びラベルから、各指標等及びラベルについてのヒストグラムを作成し、このヒストグラムから結合累積分布関数を作成する。図９（ａ）は、分布推定部１１に入力された時系列の指標の値とラベル（ＩＮＦ１）であり、第１の実施形態におけるＩＮＦ１と同様の方法により収集されたものである。図９（ｂ）は、分布推定部１１により推定され、情報量計算部１２へ入力される結合累積分布関数（ＩＮＦ２）である。
【００５７】
図８のステップＳ３３０において、情報量計算部１２は、ＩＮＦ２で示される結合累積分布関数Ｆ（Ｘ，Ｌ）それぞれから各指標等の周辺分布関数Ｐ（Ｘ）とラベルの周辺分布関数Ｐ（Ｌ）を算出して、相互情報量Ｉ（Ｘ；Ｌ）を計算する（Ｘは指標等、Ｌはラベル）。図９（ｃ）は、情報量計算部１２により計算された相互情報量（ＩＮＦ３）であり、難易度判定部１５へ入力される。
【００５８】
図８のステップＳ３４０において、難易度判定部１５は、各指標等について、入力された相互情報量の値と事前に定義された相互情報量の閾値とを比較し、当該閾値、すなわち、要求される相互情報量を上回るものを「易」、下回るものを「難」と難易度を判定する。図９（ｄ）は、難易度判定部１５により判定された各指標等の難易度（ＩＮＦ７）であり、指標選択部１３へ入力される。
【００５９】
図８のステップＳ３５０において、指標選択部１３は、ラベルの推定に最適な指標等として難易度が「易」である指標等を全て出力するが、難易度が「易」である指標等がない場合は「なし」を出力する。図９（ｅ）が指標選択部１３により出力される指標等の情報（ＩＮＦ８）である。
【００６０】
上述するように、本実施形態による利点は、第１の実施形態による先願発明との差異に加えて、指標等の使用可否を決定する閾値を事前に設定しておくことで、人手を用いずに不要な指標等を全て排除できる点にある。
【００６１】
［第４の実施形態］
図１０は、本発明の第４の実施形態による指標選択装置１ｃの構成を示すブロック図である。図１０に示すように、本実施形態による指標選択装置１ｃは、分布推定部１１と、情報量計算部１２と、難易度判定部１５と、記憶部１６と、指標組合せ部１７と、指標選択部１３とから構成される。分布推定部１１と、情報量計算部１２と、難易度判定部１５と、指標選択部１３とは第３の実施形態による指標選択装置１ｂと同様であるため、説明を省略する。記憶部１６は、入力された情報を記憶し、いつでも当該情報の読み出しが可能な記憶領域である。指標組合せ部１７は、所定の難易度を満たす指標等がない場合に、指標の組みを新たに生成するものである。
【００６２】
図１１は、図１０に示す指標選択装置１ｃの処理フローを示す図である。同図において、ステップＳ４１０〜Ｓ４４０までの処理は、図８に示す第３の実施形態のステップＳ３１０〜Ｓ３４０と同様の処理を行う。また、ＩＮＦ１〜ＩＮＦ３、ＩＮＦ７、ＩＮＦ８は第３の実施形態と同様の情報である。ただし、ステップＳ４２０において、分布推定部１１は、分布の推定を終了した後に、ステップＳ４１０において入力された指標等の値と、ラベルの情報と、推定した分布の情報を記憶部１６へ書き込む。記憶部１６は、書き込まれた情報を蓄積し、情報の読み出し指示を待つ。また、ステップＳ４４０において、難易度判定部１５は、難易度が高いと判断した指標等を、情報量計算部１２より入力された情報から排除した結果（ＩＮＦ７）を指標選択部１３へ出力するとともに、指標組合せ部１７へ出力する。要求を満たす指標等がない場合は、要求を満たす指標等がない旨を示す情報（ＩＮＦ７）を指標選択部１３及び指標組合せ部１７へ出力する。
【００６３】
難易度判定部１５から要求を満たす指標等がない旨を示す情報が指標選択部１３及び指標組合せ部１７に出力された場合（ステップＳ４５０：ＮＯ）、指標選択部１３は、要求を満たす指標等がない旨を示す情報（ＩＮＦ８）を出力する。さらに、指標組合せ部１７は、１つの指標、及び、現在の指標の組みでは、所望のラベルの判定が不可能と判断し、ランダムに指標を選択して２つ以上の指標の新たな組みを生成して分布推定部１１へ出力する（ステップＳ４６０）。
【００６４】
分布推定部１１は、指標組合せ部１７から新たな指標の組みの入力を受けると、指標選択装置１ｃへの入力の有無に関わらず、記憶部１６内の情報を読み出して、ステップＳ４２０と同様に指標等とラベルの分布の関数を推定する。つまり、分布推定部１１は、指標組合せ部１７が生成した新たな指標の組みについて各指標の値域を複数に分割し、分割した値域の各範囲に含まれるサンプル数をカウントして、ステップＳ４２０からの処理を繰り返し、新たな指標の組みについて、第１の実施形態と同様に、分布関数の推定、生成を行ない、情報量計算部１２へ出力する。
【００６５】
そして、ステップＳ４５０において、難易度判定部１５から要求を満たす指標等が指標選択部１３及び指標組合せ部１７に出力された場合（ステップＳ４５０：ＹＥＳ）、指標選択部１３は、難易度判定部１５から入力された、要求を満たす指標等の全てと、その情報量と（ＩＮＦ８）を出力する（ステップＳ４７０）。
【００６６】
図１２は、図１０に示した構成において、図１１の処理を行うときに用いられる各情報各情報ＩＮＦ１〜３、７〜９、１１の具体例を示すものである。各情報ＩＮＦ１〜３、７〜９、１１は、以下を示している。ここでは、指標を複数の生理指標、知的活動ラベルを興味の有無とした場合を例にしている。
【００６７】
ＩＮＦ１：指標選択装置１ａヘ入力される指標の値とラベルである。
ＩＮＦ２：分布推定部１１により推定された指標等とラベルの分布である。
ＩＮＦ３：情報量計算部１２により計算された各指標等とラベルの情報量である。
ＩＮＦ７：難易度判定部１５により判定された難易度である。
ＩＮＦ８：指標選択部１３により選択された指標等である。
ＩＮＦ９：記憶部１６の蓄積情報である。
ＩＮＦ１１：分布推定部１１により生成された指標の組みの情報である。
【００６８】
図１１のステップＳ４２０において、分布推定部１１は、ステップＳ４１０において入力された時系列の指標の値及びラベルから、各指標等及びラベルについてのヒストグラムを作成し、このヒストグラムから結合累積分布関数を作成する。図１２（ａ）は、分布推定部１１に入力された時系列の指標とラベルの情報（ＩＮＦ１）であり、第１の実施形態におけるＩＮＦ１と同様の方法により収集されたものである。図１２（ｂ）は、分布推定部１１により推定され、情報量計算部１２へ入力される結合累積分布関数の情報（ＩＮＦ２）である。
【００６９】
図１２（ｃ）に示すように、記憶部１６へは、分布推定部１１へ入力された時系列の指標とラベルの情報が書き込まれる。記憶部１６は、情報の読み出しの指示を受けるまで入力された情報を保存し、読み出しの指示を受けると、保持していた情報を出力する。
【００７０】
図１１のステップＳ４３０において、情報量計算部１２は、ＩＮＦ２で示される結合累積分布関数Ｆ（Ｘ，Ｌ）それぞれから指標等の周辺分布関数Ｐ（Ｘ）とラベルの周辺分布関数Ｐ（Ｌ）を算出して、相互情報量Ｉ（Ｘ；Ｌ）を計算する（Ｘは指標等、Ｌはラベル）。図１２（ｄ）は、情報量計算部１２により計算された相互情報量（ＩＮＦ３）であり、難易度判定部１５へ入力される。
【００７１】
図１１のステップＳ４４０において、難易度判定部１５は、各指標等について、入力された相互情報量の値と事前に定義された相互情報量の閾値とを比較し、当該閾値、すなわち、要求される相互情報量を上回るものを「易」、下回るものを「難」と難易度を判定する。図１１（ｅ）は、難易度判定部１５により判定された各指標等の難易度（ＩＮＦ７）であり、指標組合せ部１７及び指標選択部１３へ入力される。
【００７２】
指標選択部１３は、図１１のステップＳ４７０において、ラベルの推定に最適な指標等として難易度が「易」である指標等全てを出力し、全ての指標等の難易度が「難」であり、最適な指標等がない場合は、図１１のステップＳ４６０において「なし」を示す情報を出力する。図１２（ｆ）は、最適な指標等がない場合に出力される指標等の情報（ＩＮＦ８）である。
【００７３】
指標組合せ部１７は、全ての指標等の難易度が「難」である場合、図１１のステップＳ４６０において、指標等を組み合わせた新たな組みを生成し、指標選択装置１へ出力する。ただし、既に生成済みの指標の組みは作成せず、また組み合わせはランダムに生成する。図１２（ｇ）は、指標組合せ部１７により生成され、出力される指標の組みを示す情報（ＩＮＦ１１）である。
【００７４】
分布推定部１１は、指標組合せ部１７から新たな指標の組みの情報（ＩＮＦ１１）を受けとった場合、図１１のステップＳ４２０に戻り、新しい指標の組みそれぞれとラベルの結合累積分布関数を作成するため、記憶部１６より蓄積情報（ＩＮＦ９）を読み出し、結合累積分布関数を再推定し、情報量計算部１２へ入力する。図１２（ｈ）が再推定されて出力される結合累積分布関数（ＩＮＦ２）である。
【００７５】
情報量計算部１２は、１回目の入力を受けた時と全く同様に動作し、相互情報量を計算する。図１２（ｉ）は、情報量計算部１２により計算された相互情報量を示す情報（ＩＮＦ３）であり、難易度判定部１５へ入力される。難易度判定部１５は、１回目の入力を受けた時と全く同様に動作し、難易度を判定する。図１２（ｊ）は、難易度判定部１５により判定された難易度（ＩＮＦ７）であり、指標選択部１３及び指標組合せ部１７へ入力される。指標選択部１３は、１回目の入力を受けた時と全く同様に動作する。ここでは、難易度が「易」と判定される指標の組があるため、図１２（ｋ）のような指標の組（ＩＮＦ８）が出力される。また、指標組合せ部１７は、１回目の入力を受けた時と全く同様に動作するが、ここでは、難易度が「易」と判定される指標の組があるため、何もせずに停止する。
【００７６】
本実施形態の長所は、上述した第３の実施形態の利点に加えて、単純な指標の組だけでは知的活動評価が推定できない場合であっても、自律的に当該評価の推定が可能な指標の組みを探索し、要求を満たす組が存在すれば、適切な結果を出力する点にある。
【００７７】
以上説明した本発明の実施形態によれば、知的活動評価のための主指標が不明の状況でも、知的活動評価に好適な指標等を示すことができ、更に学習前にその指標の誤差を見積もることができ、加えて当該人の知的活動評価の難易度を見積もることができる。
【００７８】
なお、上述の指標選択装置１、１ａ、１ｂ、１ｃは、内部にコンピュータシステムを有している。そして、上述した分布推定部１１、情報量計算部１２、指標選択部１３、誤差見積部、難易度判定部１５、指標組合せ部１７の動作の過程は、プログラムの形式でコンピュータ読み取り可能な記録媒体に記憶されており、このプログラムをコンピュータシステムが読み出して実行することによって、上記処理が行われる。ここでいうコンピュータシステムとは、ＣＰＵ及び各種メモリやＯＳ、周辺機器等のハードウェアを含むものである。
【００７９】
また、「コンピュータシステム」は、ＷＷＷシステムを利用している場合であれば、ホームページ提供環境（あるいは表示環境）も含むものとする。
また、「コンピュータ読み取り可能な記録媒体」とは、フレキシブルディスク、光磁気ディスク、ＲＯＭ、ＣＤ−ＲＯＭ等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置のことをいう。さらに「コンピュータ読み取り可能な記録媒体」とは、インターネット等のネットワークや電話回線等の通信回線を介してプログラムを送信する場合の通信線のように、短時間の間、動的にプログラムを保持するもの、その場合のサーバやクライアントとなるコンピュータシステム内部の揮発性メモリのように、一定時間プログラムを保持しているものも含むものとする。また上記プログラムは、前述した機能の一部を実現するためのものであっても良く、さらに前述した機能をコンピュータシステムにすでに記録されているプログラムとの組み合わせで実現できるものであっても良い。
【図面の簡単な説明】
【００８０】
【図１】本発明の第１の実施形態による指標選択装置の構成を示すブロック図である。
【図２】同実施形態における指標選択装置の処理フローを示す図である。
【図３】同実施形態における機能ブロック間のインタフェース例を示す図である。
【図４】本発明の第２の実施形態による指標選択装置の構成を示すブロック図である。
【図５】同実施形態における指標選択装置の処理フローを示す図である。
【図６】同実施形態における機能ブロック間のインタフェース例を示す図である。
【図７】本発明の第３の実施形態による指標選択装置の構成を示すブロック図である。
【図８】同実施形態における指標選択装置の処理フローを示す図である。
【図９】同実施形態における機能ブロック間のインタフェース例を示す図である。
【図１０】本発明の第４の実施形態による指標選択装置の構成を示すブロック図である。
【図１１】同実施形態における指標選択装置の処理フローを示す図である。
【図１２】同実施形態における機能ブロック間のインタフェース例を示す図である。
【符号の説明】
【００８１】
１、１ａ、１ｂ、１ｃ…指標選択装置
１１…分布推定部
１２…情報量計算部
１３…指標選択部
１４…誤差見積部
１５…難易度判定部
１６…記憶部
１７…指標組合せ部
【図１】
2010022649.tif 000003
【図２】
2010022649.tif 000004
【図３】
2010022649.tif 000005
【図４】
2010022649.tif 000006
【図５】
2010022649.tif 000007
【図６】
2010022649.tif 000008
【図７】
2010022649.tif 000009
【図８】
2010022649.tif 000010
【図９】
2010022649.tif 000011
【図１０】
2010022649.tif 000012
【図１１】
2010022649.tif 000013
【図１２】
2010022649.tif 000014