ホーム > 特許ランキング > 国立大学法人大阪大学
知財求人 - 知財ポータルサイト「IP Force」
▶ 国立大学法人大阪大学の特許一覧 ▶ ＰａＭｅＬａ株式会社の特許一覧
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2023-01-04
(45)【発行日】2023-01-13
(54)【発明の名称】快不快の判別
(51)【国際特許分類】
A61B 5/377 20210101AFI20230105BHJP
A61B 5/372 20210101ALI20230105BHJP
A61B 5/16 20060101ALI20230105BHJP
【ＦＩ】
A61B5/377
A61B5/372
A61B5/16 110
【請求項の数】
15
(21)【出願番号】P 2020018127
(22)【出願日】2020-02-05
(62)【分割の表示】P 2018560695の分割
【原出願日】2018-07-27
(65)【公開番号】P2020073084
(43)【公開日】2020-05-14
【審査請求日】2021-07-26
(31)【優先権主張番号】P 2017146553
(32)【優先日】2017-07-28
(33)【優先権主張国・地域又は機関】JP
【国等の委託研究の成果に係る記載事項】（出願人による申告）革新的研究開発推進プログラム（ＩｍＰＡＣＴ）「対話健康支援ロボティクス」に係る委託業務、産業技術力強化法第１７条の適用を受ける特許出願
(73)【特許権者】
【識別番号】504176911
【氏名又は名称】国立大学法人大阪大学
(73)【特許権者】
【識別番号】517249071
【氏名又は名称】ＰａＭｅＬａ株式会社
(74)【代理人】
【識別番号】100131705
【弁理士】
【氏名又は名称】新山 雄一
(74)【代理人】
【識別番号】100190621
【弁理士】
【氏名又は名称】崎間 伸洋
(72)【発明者】
【氏名】中江 文
(72)【発明者】
【氏名】曽雌 崇弘
【審査官】門田 宏
(56)【参考文献】
【文献】国際公開第２０１６／１３６３６１（ＷＯ，Ａ１）
【文献】特開２０１５－０９５２１２（ＪＰ，Ａ）
【文献】国際公開第２０１１／１５５１９６（ＷＯ，Ａ１）
【文献】国際公開第２０１２／０７０２１２（ＷＯ，Ａ１）
【文献】国際公開第２０１８／０３８１２１（ＷＯ，Ａ１）
【文献】国際公開第２０１４／１５１８７４（ＷＯ，Ａ１）
【文献】特開２０１３－２５５７４２（ＪＰ，Ａ）
【文献】米田徹，痛みを誘発する電気的な新奇刺激に対する事象関連電位の測定，電子情報通信学会技術研究報告，Vol.98 No.400，日本，1998年11月17日，Page 7-12
【文献】青樹秀樹，痛覚刺激時における脳波ゆらぎについて，電子情報通信学会技術研究報告，Vol.95 No.356，日本，1995年11月16日，Page 71-78
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ａ６１Ｂ ５／３６９ － ５／３８６
Ａ６１Ｂ ５／１６
Ａ６１Ｂ １０／００
ＪＳＴＰｌｕｓ／ＪＭＥＤＰｌｕｓ／ＪＳＴ７５８０（ＪＤｒｅａｍＩＩＩ）
(57)【特許請求の範囲】
【請求項1】
対象の快不快を判定する手段をコンピュータにより生成する方法であって、前記コンピュータは、演算手段を含み、
ａ）被験者に、同じ環境下で異なる刺激を付与することで、前記演算手段が、センサーにより前記被験者について計測された各々の脳波データまたはその分析データを入手する工程と、
ｂ）前記演算手段が、ａ）の工程で得られた前記脳波データまたはその分析データの相違と、前記刺激に関する前記被験者の反応とを関連付けを実行する工程と、
ｃ）前記関連付けに基づいて、前記演算手段が、前記被験者の快不快の判定のためのモデル作成処理を実行して、前記対象の快不快を判定するための快不快判定手段を生成する工程と
を備え、
前記快不快判定手段を生成する工程は、
ｉ）前記脳波データに対して、瞬目および筋電位によるノイズ除去と、前記刺激の付与前の所定期間のデータを用いた標準化を実施する前処理と、
ｉｉ）前記前処理後のデータについて、第１の刺激提示後の所定期間の快状態の第１の振幅データと第２の刺激提示後の所定期間の不快状態の第２の振幅データとを用いて、ロジスティック回帰モデルにより前記快不快判定手段を生成する処理とを含む、方法。
【請求項2】
対象の快不快をコンピュータにより判定する方法であって、前記コンピュータは、演算手段を含み、
ｄ）被験者に同じ環境下で異なる刺激を付与する試験により取得された、モデル作成用の脳波データまたはその分析データの相違と前記刺激に関する前記被験者の反応との関連付けに基づいて、前記被験者の快不快の判定のモデル作成処理によって作成された、前記対象の快不快を判定するための快不快判定手段を提供する工程と、
ｅ）前記演算手段が、センサーにより計測された前記対象からの脳波データまたはその分析データを入手し、前記快不快判定手段に適用して、前記対象の快不快を判定する工程と
を備え、
前記快不快判定手段は、
ｉ）前記脳波データに対して、瞬目および筋電位によるノイズ除去と、前記刺激の付与前の所定期間のデータを用いた標準化を実施する前処理を実行するものであり、
ｉｉ）前記前処理後のデータについて、刺激提示後の所定期間の快状態の第１の振幅データと不快状態の第２の振幅データとを用いて、ロジスティック回帰モデルにより作成される、ことを特徴とする、方法。
【請求項3】
前記前処理は、
前記センサーの複数の電極からの脳波信号に対する主成分分析による瞬目成分の除去および所定周波数についての高域カットフィルターによる筋電位ノイズ除去の処理と、
前記刺激の付与前の所定期間に計測される脳波データを用いたｚ値化による標準化処理を含む、請求項１または２記載の方法。
【請求項4】
前記演算手段が、前記モデル作成処理の前に、前記複数の被験者についての前記第１の振幅データおよび前記第２の振幅データの統計的な分布に基づいて、自己複製処理を実行し、前記モデル作成処理のためのデータの数を増加させる工程をさらに含む、請求項３に記載の方法。
【請求項5】
前記演算手段は、前記自己複製処理を、前記統計的な分布に基づく正規乱数を発生させることによって実行する、請求項４に記載の方法。
【請求項6】
前記演算手段は、前記自己複製処理を、前記統計的な分布に基づくピアソンシステム乱数を発生させることによって実行する、請求項４に記載の方法。
【請求項7】
前記快不快判定手段の生成におけるモデル作成処理は、Ｌ１正則化およびベイズ最適化によるハイパーパラメータの決定処理を含む、請求項５または６に記載の方法。
【請求項8】
前記快不快判定手段の生成におけるモデル作成処理は、Ｌ１正則化およびグリッドサーチによるハイパーパラメータの決定処理を含む、請求項５または６に記載の方法。
【請求項9】
対象の快不快を判定する手段を生成する装置であって、
Ａ）被験者に、同じ環境下で異なる刺激を付与することで、センサーにより前記被験者について計測された各々の脳波データまたはその分析データを入手するデータ入手部と、
Ｂ）前記データ入手部で得られた前記脳波データまたはその分析データの相違と、前記刺激に関する前記被験者の反応とを関連付ける処理部と、
Ｃ）前記関連付けに基づいて、前記被験者の快不快の判定のためのモデル作成処理を実行して、前記対象の快不快を判定するための快不快判定手段を生成する判定手段生成部と
を備え、
前記判定手段生成部は、
ｉ）前記脳波データに対して、前処理として、瞬目および筋電位によるノイズ除去と、前記刺激の付与前の所定期間のデータを用いた標準化を実施する前処理部と、
ｉｉ）前記前処理後のデータについて、第１の刺激提示後の所定期間の快状態の第１の振幅データと第２の刺激提示後の所定期間の不快状態の第２の振幅データとを用いて、ロジスティック回帰モデルにより前記快不快判定手段を生成するモデル作成処理部とを含む、装置。
【請求項10】
対象の快不快を判定する装置であって、
Ｄ）被験者に同じ環境下で異なる刺激を付与する試験により取得された、モデル作成用の脳波データまたはその分析データの相違と前記刺激に関する前記被験者の反応との関連付けに基づいて、前記被験者の快不快の判定のモデル作成処理によって作成された、前記対象の快不快を判定するための快不快判定手段を提供する判定手段提供部と、
Ｅ）前記対象から試験用の脳波データまたはその分析データを入手し、前記快不快判定手段に適用し、前記対象の快不快を判定する判定部と
を備え、
前記快不快判定手段は、
前記脳波データに対して、前処理として、瞬目および筋電位によるノイズ除去と、前記刺激の付与前の所定期間のデータを用いた標準化を実施する前処理部を含み、
前記前処理後のデータについて、刺激提示後の所定期間の快状態の第１の振幅データと不快状態の第２の振幅データとを用いて作成された、ロジスティック回帰モデルである、装置。
【請求項11】
前記前処理は、
前記センサーの複数の電極からの脳波信号に対する主成分分析による瞬目成分の除去および所定周波数についての高域カットフィルターによる筋電位ノイズ除去の処理と、
前記刺激の付与前の所定期間に計測される脳波データを用いたｚ値化による標準化処理を含む、請求項９または１０記載の方法。
【請求項12】
対象の快不快を判定する手段を生成する方法をコンピュータに実行させるためのプログラムであって、前記コンピュータは、演算手段を含み、前記方法は、
ａ）被験者に、同じ環境下で異なる刺激を付与することで、前記演算手段が、センサーにより前記被験者について計測された各々の脳波データまたはその分析データを入手するステップと、
ｂ）前記演算手段が、ａ）の工程で得られた前記脳波データまたはその分析データの相違と、前記刺激に関する前記被験者の反応とを関連付けを実行するステップと、
ｃ）前記関連付けに基づいて、前記演算手段が、前記被験者の快不快の判定のためのモデル作成処理を実行して、前記対象の快不快を判定するための快不快判定手段を生成するステップと
を備え、
前記快不快判定手段を生成するステップは、
ｉ）前記脳波データに対して、瞬目および筋電位によるノイズ除去と、前記刺激の付与前の所定期間のデータを用いた標準化を実施する前処理と、
ｉｉ）前記前処理後のデータについて、第１の刺激提示後の所定期間の快状態の第１の振幅データと第２の刺激提示後の所定期間の不快状態の第２の振幅データとを用いて、ロジスティック回帰モデルにより前記快不快判定手段を生成する処理とを含む、プログラム。
【請求項13】
対象の快不快を判定する方法をコンピュータに実行させるためのプログラムであって、前記コンピュータは、演算手段を含み、前記方法は、
ｄ）被験者に同じ環境下で異なる刺激を付与する試験により取得された、モデル作成用の脳波データまたはその分析データの相違と前記刺激に関する前記被験者の反応との関連付けに基づいて、前記被験者の快不快の判定のモデル作成処理によって作成された、前記対象の快不快を判定するための快不快判定手段を提供するステップと、
ｅ）前記演算手段が、センサーにより計測された前記対象からの脳波データまたはその分析データを入手し、前記快不快判定手段に適用して、前記対象の快不快を判定するステップと
を備え、
前記快不快判定手段は、
ｉ）前記脳波データに対して、瞬目および筋電位によるノイズ除去と、前記刺激の付与前の所定期間のデータを用いた標準化を実施する前処理を実行し、
ｉｉ）前記前処理後のデータについて、刺激提示後の所定期間の快状態の第１の振幅データと不快状態の第２の振幅データとを用いて、ロジスティック回帰モデルを作成する処理により生成されたものである、プログラム。
【請求項14】
対象の快不快を判定する手段を生成する方法をコンピュータに実行させるための請求項１２に記載のプログラムを記録した記録媒体。
【請求項15】
対象の快不快を判定する方法をコンピュータに実行させるための請求項１３に記載のプログラムを記録した記録媒体。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、脳波を用いて快不快を判別する技術に関する。より詳細には、本発明は、同一の環境における異なる刺激での脳波データまたはその分析データに基づいて、快不快を判別することに関し、例えば、痛みがある場合にそれが心地よい痛みなのか不快な痛みなのかを判別することができる。
【背景技術】
【0002】
種々の感覚は、一方向のベクトルで表されることが多く、例えば、疼痛をみると、痛いかまたは痛くないかで区別することが多い。しかし、痛い感覚でも、それが鍼灸治療におけるもののように心地よい（快い）痛みなのか、それとも、不快な痛みなのか、簡単に区別できるものではない。
【発明の概要】
【課題を解決するための手段】
【0003】
本発明者らは鋭意研究をした結果、種々の感覚について、それが心地よい（快い）感覚なのか、またはそうではない不快な感覚なのかを判別することができる判別技術を見出した。
【0004】
本発明は、例えば、以下を提供する。
（項目１）この発明の１つの局面に従うと、対象の快不快を判定する手段をコンピュータにより生成する方法であって、前記コンピュータは、演算手段を含み、
ａ）被験者に、同じ環境下で異なる刺激を付与することで、前記演算手段が、センサーにより前記被験者について計測された各々の脳波データまたはその分析データを入手する工程と、
ｂ）前記演算手段が、ａ）の工程で得られた前記脳波データまたはその分析データの相違と、前記刺激に関する前記被験者の反応とを関連付けを実行する工程と、
ｃ）前記関連付けに基づいて、前記演算手段が、前記被験者の快不快の判定のためのモデル作成処理を実行して、前記対象の快不快を判定するための快不快判定手段を生成する工程と
を備え、
前記快不快判定手段を生成する工程は、
ｉ）前記脳波データに対して、瞬目および筋電位によるノイズ除去と、前記刺激の付与前の所定期間のデータを用いた標準化を実施する前処理と、
ｉｉ）前記前処理後のデータについて、第１の刺激提示後の所定期間の快状態の第１の振幅データと第２の刺激提示後の所定期間の不快状態の第２の振幅データとを用いて、ロジスティック回帰モデルにより前記快不快判定手段を生成する処理とを含む、方法。
（項目２）この発明の他の局面に従うと、対象の快不快をコンピュータにより判定する方法であって、前記コンピュータは、演算手段を含み、
ｄ）被験者に同じ環境下で異なる刺激を付与する試験により取得された、モデル作成用の脳波データまたはその分析データの相違と前記刺激に関する前記被験者の反応との関連付けに基づいて、前記被験者の快不快の判定のモデル作成処理によって作成された、前記対象の快不快を判定するための快不快判定手段を提供する工程と、
ｅ）前記演算手段が、センサーにより計測された前記対象からの脳波データまたはその分析データを入手し、前記快不快判定手段に適用して、前記対象の快不快を判定する工程と
を備え、
前記快不快判定手段は、
ｉ）前記脳波データに対して、瞬目および筋電位によるノイズ除去と、前記刺激の付与前の所定期間のデータを用いた標準化を実施する前処理を実行するものであり、
ｉｉ）前記前処理後のデータについて、刺激提示後の所定期間の快状態の第１の振幅データと不快状態の第２の振幅データとを用いて、ロジスティック回帰モデルにより作成される、ことを特徴とする、方法。
（項目３）好ましくは、前記前処理は、
前記センサーの複数の電極からの脳波信号に対する主成分分析による瞬目成分の除去および所定周波数についての高域カットフィルターによる筋電位ノイズ除去の処理と、
前記刺激の付与前の所定期間に計測される脳波データを用いたｚ値化による標準化処理を含む。
（項目４）好ましくは、前記演算手段が、前記モデル作成処理の前に、複数の被験者についての前記第１の振幅データおよび前記第２の振幅データの統計的な分布に基づいて、自己複製処理を実行し、前記モデル作成処理のためのデータの数を増加させる工程をさらに含む。
（項目５）好ましくは、前記演算手段は、前記自己複製処理を、前記統計的な分布に基づく正規乱数を発生させることによって実行する。
（項目６）好ましくは、前記演算手段は、前記自己複製処理を、前記統計的な分布に基づくピアソンシステム乱数を発生させることによって実行する。
（項目７）好ましくは、前記快不快判定手段の生成におけるモデル作成処理は、Ｌ１正則化およびベイズ最適化によるハイパーパラメータの決定処理を含む。
（項目８）好ましくは、前記快不快判定手段の生成におけるモデル作成処理は、Ｌ１正則化およびグリッドサーチによるハイパーパラメータの決定処理を含む。
（項目９）この発明のさらに他の局面に従うと、対象の快不快を判定する手段を生成する装置であって、
Ａ）被験者に、同じ環境下で異なる刺激を付与することで、センサーにより前記被験者について計測された各々の脳波データまたはその分析データを入手するデータ入手部と、
Ｂ）前記データ入手部で得られた前記脳波データまたはその分析データの相違と、前記刺激に関する前記被験者の反応とを関連付ける処理部と、
Ｃ）前記関連付けに基づいて、前記被験者の快不快の判定のためのモデル作成処理を実行して、前記対象の快不快を判定するための快不快判定手段を生成する判定手段生成部と
を備え、
前記判定手段生成部は、
ｉ）前記脳波データに対して、前処理として、瞬目および筋電位によるノイズ除去と、前記刺激の付与前の所定期間のデータを用いた標準化を実施する前処理部と、
ｉｉ）前記前処理後のデータについて、第１の刺激提示後の所定期間の快状態の第１の振幅データと第２の刺激提示後の所定期間の不快状態の第２の振幅データとを用いて、ロジスティック回帰モデルにより前記快不快判定手段を生成するモデル作成処理部とを含む、装置。
（項目１０）この発明のさらに他の局面に従うと、対象の快不快を判定する装置であって、
Ｄ）被験者に同じ環境下で異なる刺激を付与する試験により取得された、モデル作成用の脳波データまたはその分析データの相違と前記刺激に関する前記被験者の反応との関連付けに基づいて、前記被験者の快不快の判定のモデル作成処理によって作成された、前記対象の快不快を判定するための快不快判定手段を提供する判定手段提供部と、
Ｅ）前記対象から試験用の脳波データまたはその分析データを入手し、前記快不快判定手段に適用し、前記対象の快不快を判定する判定部と
を備え、
前記快不快判定手段は、
前記脳波データに対して、前処理として、瞬目および筋電位によるノイズ除去と、前記刺激の付与前の所定期間のデータを用いた標準化を実施する前処理部を含み、
前記前処理後のデータについて、刺激提示後の所定期間の快状態の第１の振幅データと不快状態の第２の振幅データとを用いて作成された、ロジスティック回帰モデルである、装置。
（項目１１）好ましくは、前記前処理は、
前記センサーの複数の電極からの脳波信号に対する主成分分析による瞬目成分の除去および所定周波数についての高域カットフィルターによる筋電位ノイズ除去の処理と、
前記刺激の付与前の所定期間に計測される脳波データを用いたｚ値化による標準化処理を含む。
（項目１２）この発明のさらに他の局面に従うと、対象の快不快を判定する手段を生成する方法をコンピュータに実行させるためのプログラムであって、前記コンピュータは、演算手段を含み、前記方法は、
ａ）被験者に、同じ環境下で異なる刺激を付与することで、前記演算手段が、センサーにより前記被験者について計測された各々の脳波データまたはその分析データを入手するステップと、
ｂ）前記演算手段が、ａ）の工程で得られた前記脳波データまたはその分析データの相違と、前記刺激に関する前記被験者の反応とを関連付けを実行するステップと、
ｃ）前記関連付けに基づいて、前記演算手段が、前記被験者の快不快の判定のためのモデル作成処理を実行して、前記対象の快不快を判定するための快不快判定手段を生成するステップと
を備え、
前記快不快判定手段を生成するステップは、
ｉ）前記脳波データに対して、瞬目および筋電位によるノイズ除去と、前記刺激の付与前の所定期間のデータを用いた標準化を実施する前処理と、
ｉｉ）前記前処理後のデータについて、第１の刺激提示後の所定期間の快状態の第１の振幅データと第２の刺激提示後の所定期間の不快状態の第２の振幅データとを用いて、ロジスティック回帰モデルにより前記快不快判定手段を生成する処理とを含む、プログラム。
（項目１３）この発明のさらに他の局面に従うと、対象の快不快を判定する方法をコンピュータに実行させるためのプログラムであって、前記コンピュータは、演算手段を含み、前記方法は、
ｄ）被験者に同じ環境下で異なる刺激を付与する試験により取得された、モデル作成用の脳波データまたはその分析データの相違と前記刺激に関する前記被験者の反応との関連付けに基づいて、前記被験者の快不快の判定のモデル作成処理によって作成された、前記対象の快不快を判定するための快不快判定手段を提供するステップと、
ｅ）前記演算手段が、センサーにより計測された前記対象からの脳波データまたはその分析データを入手し、前記快不快判定手段に適用して、前記対象の快不快を判定するステップと
を備え、
前記快不快判定手段は、
ｉ）前記脳波データに対して、瞬目および筋電位によるノイズ除去と、前記刺激の付与前の所定期間のデータを用いた標準化を実施する前処理を実行し、
ｉｉ）前記前処理後のデータについて、刺激提示後の所定期間の快状態の第１の振幅データと不快状態の第２の振幅データとを用いて、ロジスティック回帰モデルを作成する処理により生成されたものである、プログラム。
（項目１４）好ましくは、対象の快不快を判定する手段を生成する方法をコンピュータに実行させるための請求項１２に記載のプログラムを記録した記録媒体。
（項目１５） 対象の快不快を判定する方法をコンピュータに実行させるための請求項１３に記載のプログラムを記録した記録媒体。
（項目１）この発明の１つの局面に従うと、対象の快不快を判定する手段をコンピュータにより生成する方法であって、前記コンピュータは、演算手段を含み、
ａ）被験者に、同じ環境下で異なる刺激を付与することで、前記演算手段が、センサーにより前記被験者について計測された各々の脳波データまたはその分析データを入手する工程と、
ｂ）前記演算手段が、ａ）の工程で得られた前記脳波データまたはその分析データの相違と、前記刺激に関する前記被験者の反応とを関連付けを実行する工程と、
ｃ）前記関連付けに基づいて、前記演算手段が、前記被験者の快不快の判定のためのモデル作成処理を実行して、前記対象の快不快を判定するための快不快判定手段を生成する工程と
を備え、
前記快不快判定手段を生成する工程は、
ｉ）前記脳波データに対して、瞬目および筋電位によるノイズ除去と、前記刺激の付与前の所定期間のデータを用いた標準化を実施する前処理と、
ｉｉ）前記前処理後のデータについて、第１の刺激提示後の所定期間の快状態の第１の振幅データと第２の刺激提示後の所定期間の不快状態の第２の振幅データとを用いて、ロジスティック回帰モデルにより前記快不快判定手段を生成する処理とを含む、方法。
（項目２）この発明の他の局面に従うと、対象の快不快をコンピュータにより判定する方法であって、前記コンピュータは、演算手段を含み、
ｄ）被験者に同じ環境下で異なる刺激を付与する試験により取得された、モデル作成用の脳波データまたはその分析データの相違と前記刺激に関する前記被験者の反応との関連付けに基づいて、前記被験者の快不快の判定のモデル作成処理によって作成された、前記対象の快不快を判定するための快不快判定手段を提供する工程と、
ｅ）前記演算手段が、センサーにより計測された前記対象からの脳波データまたはその分析データを入手し、前記快不快判定手段に適用して、前記対象の快不快を判定する工程と
を備え、
前記快不快判定手段は、
ｉ）前記脳波データに対して、瞬目および筋電位によるノイズ除去と、前記刺激の付与前の所定期間のデータを用いた標準化を実施する前処理を実行するものであり、
ｉｉ）前記前処理後のデータについて、刺激提示後の所定期間の快状態の第１の振幅データと不快状態の第２の振幅データとを用いて、ロジスティック回帰モデルにより作成される、ことを特徴とする、方法。
（項目３）好ましくは、前記前処理は、
前記センサーの複数の電極からの脳波信号に対する主成分分析による瞬目成分の除去および所定周波数についての高域カットフィルターによる筋電位ノイズ除去の処理と、
前記刺激の付与前の所定期間に計測される脳波データを用いたｚ値化による標準化処理を含む。
（項目４）好ましくは、前記演算手段が、前記モデル作成処理の前に、複数の被験者についての前記第１の振幅データおよび前記第２の振幅データの統計的な分布に基づいて、自己複製処理を実行し、前記モデル作成処理のためのデータの数を増加させる工程をさらに含む。
（項目５）好ましくは、前記演算手段は、前記自己複製処理を、前記統計的な分布に基づく正規乱数を発生させることによって実行する。
（項目６）好ましくは、前記演算手段は、前記自己複製処理を、前記統計的な分布に基づくピアソンシステム乱数を発生させることによって実行する。
（項目７）好ましくは、前記快不快判定手段の生成におけるモデル作成処理は、Ｌ１正則化およびベイズ最適化によるハイパーパラメータの決定処理を含む。
（項目８）好ましくは、前記快不快判定手段の生成におけるモデル作成処理は、Ｌ１正則化およびグリッドサーチによるハイパーパラメータの決定処理を含む。
（項目９）この発明のさらに他の局面に従うと、対象の快不快を判定する手段を生成する装置であって、
Ａ）被験者に、同じ環境下で異なる刺激を付与することで、センサーにより前記被験者について計測された各々の脳波データまたはその分析データを入手するデータ入手部と、
Ｂ）前記データ入手部で得られた前記脳波データまたはその分析データの相違と、前記刺激に関する前記被験者の反応とを関連付ける処理部と、
Ｃ）前記関連付けに基づいて、前記被験者の快不快の判定のためのモデル作成処理を実行して、前記対象の快不快を判定するための快不快判定手段を生成する判定手段生成部と
を備え、
前記判定手段生成部は、
ｉ）前記脳波データに対して、前処理として、瞬目および筋電位によるノイズ除去と、前記刺激の付与前の所定期間のデータを用いた標準化を実施する前処理部と、
ｉｉ）前記前処理後のデータについて、第１の刺激提示後の所定期間の快状態の第１の振幅データと第２の刺激提示後の所定期間の不快状態の第２の振幅データとを用いて、ロジスティック回帰モデルにより前記快不快判定手段を生成するモデル作成処理部とを含む、装置。
（項目１０）この発明のさらに他の局面に従うと、対象の快不快を判定する装置であって、
Ｄ）被験者に同じ環境下で異なる刺激を付与する試験により取得された、モデル作成用の脳波データまたはその分析データの相違と前記刺激に関する前記被験者の反応との関連付けに基づいて、前記被験者の快不快の判定のモデル作成処理によって作成された、前記対象の快不快を判定するための快不快判定手段を提供する判定手段提供部と、
Ｅ）前記対象から試験用の脳波データまたはその分析データを入手し、前記快不快判定手段に適用し、前記対象の快不快を判定する判定部と
を備え、
前記快不快判定手段は、
前記脳波データに対して、前処理として、瞬目および筋電位によるノイズ除去と、前記刺激の付与前の所定期間のデータを用いた標準化を実施する前処理部を含み、
前記前処理後のデータについて、刺激提示後の所定期間の快状態の第１の振幅データと不快状態の第２の振幅データとを用いて作成された、ロジスティック回帰モデルである、装置。
（項目１１）好ましくは、前記前処理は、
前記センサーの複数の電極からの脳波信号に対する主成分分析による瞬目成分の除去および所定周波数についての高域カットフィルターによる筋電位ノイズ除去の処理と、
前記刺激の付与前の所定期間に計測される脳波データを用いたｚ値化による標準化処理を含む。
（項目１２）この発明のさらに他の局面に従うと、対象の快不快を判定する手段を生成する方法をコンピュータに実行させるためのプログラムであって、前記コンピュータは、演算手段を含み、前記方法は、
ａ）被験者に、同じ環境下で異なる刺激を付与することで、前記演算手段が、センサーにより前記被験者について計測された各々の脳波データまたはその分析データを入手するステップと、
ｂ）前記演算手段が、ａ）の工程で得られた前記脳波データまたはその分析データの相違と、前記刺激に関する前記被験者の反応とを関連付けを実行するステップと、
ｃ）前記関連付けに基づいて、前記演算手段が、前記被験者の快不快の判定のためのモデル作成処理を実行して、前記対象の快不快を判定するための快不快判定手段を生成するステップと
を備え、
前記快不快判定手段を生成するステップは、
ｉ）前記脳波データに対して、瞬目および筋電位によるノイズ除去と、前記刺激の付与前の所定期間のデータを用いた標準化を実施する前処理と、
ｉｉ）前記前処理後のデータについて、第１の刺激提示後の所定期間の快状態の第１の振幅データと第２の刺激提示後の所定期間の不快状態の第２の振幅データとを用いて、ロジスティック回帰モデルにより前記快不快判定手段を生成する処理とを含む、プログラム。
（項目１３）この発明のさらに他の局面に従うと、対象の快不快を判定する方法をコンピュータに実行させるためのプログラムであって、前記コンピュータは、演算手段を含み、前記方法は、
ｄ）被験者に同じ環境下で異なる刺激を付与する試験により取得された、モデル作成用の脳波データまたはその分析データの相違と前記刺激に関する前記被験者の反応との関連付けに基づいて、前記被験者の快不快の判定のモデル作成処理によって作成された、前記対象の快不快を判定するための快不快判定手段を提供するステップと、
ｅ）前記演算手段が、センサーにより計測された前記対象からの脳波データまたはその分析データを入手し、前記快不快判定手段に適用して、前記対象の快不快を判定するステップと
を備え、
前記快不快判定手段は、
ｉ）前記脳波データに対して、瞬目および筋電位によるノイズ除去と、前記刺激の付与前の所定期間のデータを用いた標準化を実施する前処理を実行し、
ｉｉ）前記前処理後のデータについて、刺激提示後の所定期間の快状態の第１の振幅データと不快状態の第２の振幅データとを用いて、ロジスティック回帰モデルを作成する処理により生成されたものである、プログラム。
（項目１４）好ましくは、対象の快不快を判定する手段を生成する方法をコンピュータに実行させるための請求項１２に記載のプログラムを記録した記録媒体。
（項目１５） 対象の快不快を判定する方法をコンピュータに実行させるための請求項１３に記載のプログラムを記録した記録媒体。
【0005】
本発明において、上記１または複数の特徴は、明示された組み合わせに加え、さらに組み合わせて提供されうることが意図される。本発明のさらなる実施形態および利点は、必要に応じて以下の詳細な説明を読んで理解すれば、当業者に認識される。
【発明の効果】
【0006】
本発明は、快不快を判別することができる。本発明はまた疼痛についてみると、心地よい痛みと、不快な痛みとを判別することができ、よりきめ細やかで主観に合致した治療や手術を行うことができ、医療関連産業において有用性がある。
【図面の簡単な説明】
【0007】
【図1】図１は、異なる刺激タイプによる疼痛不快度の違いを判別するための実験パラダイムを示す。不快疼痛刺激として低温刺激（レベル３）の—１０℃を使い、対応する不快でない疼痛刺激として、強度レベルを個人ごとに統制した電気刺激（レベル３）を用いた。各レベルは３刺激含み、１５秒間続いた。
【図2】図２は、不快疼痛の判別分析のプロセスを示す。判別分析には、サポートベクターマシーンを用いて行った。より詳細には、ＳＶＭ－ＲＦＥ（Ｓｕｐｐｏｒｔ ｖｅｃｔｏｒ ｍａｃｈｉｎｅ ｒｅｃｕｒｓｉｖｅ ｆｅａｔｕｒｅ ｅｌｉｍｉｎａｔｉｏｎ）を用いた。最初に、２４個の脳波特徴量のランキングを行い、上位から１個ずつ特徴量を増やしながら、１個抜き交差検証により、判別精度を調べた。最後に、最も判別精度が高い特徴量群を用いて、ランダム化テストを行い、偶然レベルの判別精度を調べた。
【図5】図５は、上位５個の特徴量を用いたＳＶＭの判別分析によるランダム化テストの結果を示す。判別ラベルを１０００回ランダム化し、その度に、1個抜き交差検証により判別性を参集した。実際の「５６．０９８％」の判別精度は、ランダム分布の上位５％以内を占め、有意な判別精度を実現していることを示す。
【図6】図６は、同一高温刺激の痛み文脈効果（主観評価）を示す。同じ温度刺激（４０℃）にかかわらず、痛みの文脈（環境）が変わると不快度（主観評価）が変化することが示された。すなわち、３８℃の「痛くない」刺激と一緒に出た場合、４８℃の「痛い」刺激と一緒に出た場合よりも、４０℃での刺激で不快度が高くなることが示された。
【図7】図７は、同一高温刺激の痛み文脈効果（脳波振幅）を示す。同じ温度刺激（４０℃）にかかわらず、痛みの文脈(環境）が変わると脳波振幅が変化することが示された。４８℃の痛い刺激と一緒に出た場合よりも、脳波振幅が小さくなる。
【図8】図８は、不快疼痛判別器（シグモイド関数）を示す。変曲点の閾値を用いて、不快度の異なる痛みを分ける。＞１．０５５５は、不快度が低いと判断し、≦１．０５５５は不快度が高いと判断する。全体の判別精度は、６４％であった。内訳をみると、不快度低の判別精度は、６８％であり、不快度高の判別精度は６０％であった。
【図9】図９は、不快疼痛判別器（シグモイド関数）の応用例を示す。前半は、痛みモニタリング初期の不快疼痛レベル強状態を保ち、後半は、痛みモニタリング中・後期の不快疼痛レベル弱状態を示す。１５秒以上の脳活動平均値群（Ａ、Ｂ）を取得し統計検定による特徴量変化の比較を行い、不快度緩和の客観的評価を行う。同じ物理量の刺激を異なる痛み環境で提示すると、「文脈依存型リファレンス検査手法」により、より厳密な不快疼痛評価を実現する。
【図10】図１０は、心理ストレスパラダイムの実験である。認知課題として、カラーストループ課題（認知的競合課題）を付与し、フォントカラー課題：「フォントの色」を答えることを３ブロック（各５０試行）行う。図５の上に示すのは一致条件（競合無し、フォントカラーと文字情報が一致する）、不一致条件（競合あり、難しい、フォントカラーと文字情報が不一致である）の実験を行うスキームが示される。
【図11】図１１は、心理ストレスパラダイムにおいて、行動データ（回答に要した反応時間）を示す。左が心理ストレスなし（第三者による課題中の監視がない）、右が心理ストレスあり（第三者による監視がある）である。難しい不一致条件の方が回答時間が長くなっている。心理的ストレスあり、なしの反応時間にはどの条件間でも有意差は見られなかった。したがって、表立って、ストレスありなしの影響は見られないと結論付けられた。
【図12】図１２は、脳波データ（事象関連電位：Ｆ３の左前頭）の結果を示す。実線は不一致条件で、破線は一致条件（比較基準条件）である。上（Ａ)は快条件（心理ストレスなし）を示し、下（Ｂ)は不快条件（心理ストレスあり）を示す。快条件（心理ストレスなし）では、不一致条件の波形が陽性方向にシフトしていた。不快条件（心理ストレスあり）では、不一致条件の波形が陰性方向にシフトし続けた。例として、陰性電位の３００～８００ミリ秒までの占有率を計算し特徴量とした。占有率（％）＝陰性電位時間／全体時間×１００で示す。
【図13】図１３は、判別分析（シグモイド関数）である。図１２の結果をシグモイドフィッティングした結果である。潜在心理不快判別器としてその閾値を求めたところ、閾値＝４３．５であった。これ以下だとストレスなし、これを超えるとストレスありと判断する。ストレスなしの判別精度：７３．１％、ストレスありの判別精度：６１．５％であった。全体：６７．６％であった。
【図21】図２１は、個人に合わせた快不快判別モデルの作成について、サンプル増幅による自己複製型特徴量作成を実施した例である、被検者1名の熱痛み主観評価の例を示す。本実施例で用いた被検者(ID185)は、熱刺激が40℃～46℃までは、ほとんど痛みを感じなかったが、レベル5の48℃で急に痛みの不快度(耐えられないほど不快)が上昇し、レベル6において、不快度はスケールの天井に達した。したがって、判別モデル作成は、レベル１とレベル６の最も遠い条件を用いて行った。
【図22】図２２は、実データ（熱刺激レベル１と６）のサンプル分布特性を示す。左は熱刺激レベル1(40℃)のサンプル分布(n = 30003)を示す。右は、熱刺激レベル6(50℃)のサンプル分布(n= 30003)を示す。熱刺激レベル6の分布特性が、より広い範囲で右方向に広がっていることが分かる。各グラフは、y軸がサンプル数、ｘ軸が、平均電位（絶対値、ｚ値化）を示す。
【図23】図２３は、自己複製型特徴量のサンプル分布の例（分布1；正規乱数増幅）を示す。レベル１と6の実測サンプルの分布特性をもつ正規乱数(n=10000)を発生し、「自己複製型特徴量」を作成した。グラフは、Fp1電極における熱刺激レベル1(40℃)のサンプル分布(n = 10003)は、平均:-0.27、SD:0.54、ｙ軸がサンプル数、ｘ軸が、平均電位（絶対値、ｚ値化）を示す。Fp1電極における熱刺激レベル6(50℃)のサンプル分布(n =10000)は、平均:1.10、SD: 1.74である。オーバーラップするエリアが狭いほど、判別には有利である。
【図24】図２４は、個人判別モデルの汎化例（汎化1；正規乱数増幅)を示す。個人（ID185）判別モデル用特徴量係数の入手を行い、判別モデルの決定し、異なる被検者（n= 169）のレベル１とレベル6の痛み判別推定を行った。推定結果は、全体の判別精度:70.4% (SD21.8)、判別精度70以上: 78名（46.2％）、判別精度50以上:110名（65.1％）であった。
【図25】図２５は、自己複製型特徴量のサンプル分布の例（サンプル分布2；ピアソンシステム乱数増幅)を示す。縦軸はサンプル数、横軸は平均電位（絶対値、ｚ値化）である。レベル１と6の実測30,000サンプルの分布特性を持つピアソンシステム乱数(n=10000)を発生し、「自己複製型特徴量」を作成した。Fp1電極における熱刺激レベル1(40℃)のサンプル分布(n = 10,000)であった。また、平均:-0.27、SD:0.54であり、歪度: 1.10、尖度：2.62であった。Fp1電極における熱刺激レベル6(50℃)のサンプル分布(n= 10,000)であり、平均:1.10、SD: 1.74、歪度: 1.77、尖度:6.35であった。オーバーラップするエリアが狭いほど、判別には有利である。
【図26】図２６は、個人判別モデルの汎化例（汎化2；ピアソンシステム乱数増幅）を示す。個人（ID185）判別モデル用特徴量係数を入手し、判別モデルを決定し、異なる被検者（n= 169）のレベル1とレベル6の痛み判別推定を行う。推定結果は、全体の判別精度:71.7% (SD21.0)、判別精度70以上: 78名（46.2％）、判別精度50以上:114名（67.5％）であった。
【図27】図２７は、自己複製型特徴量のサンプル分布の例（サンプル分布3；ピアソンシステム乱数増幅)を示す。少ない観測サンプル（各レベル１０サンプル）の分布特性から、ピアソンシステム乱数(n=10000)を発生し、「自己複製型特徴量」を作成した。Fp1電極における熱刺激レベル1(40℃)のサンプル分布(n = 10,000)は、平均: -0.27、 SD:0.7、歪度:1.90、尖度：6.22であった。Fp1電極における熱刺激レベル6(50℃)のサンプル分布(n = 10,000)は、平均:1.21、SD: 3.15、歪度:3.44、尖度:15.97であった。オーバーラップするエリアが狭いほど、判別には有利である。
【発明を実施するための形態】
【0008】
以下、本発明を説明する。本明細書の全体にわたり、単数形の表現は、特に言及しない限り、その複数形の概念をも含むことが理解されるべきである。従って、単数形の冠詞（例えば、英語の場合は「ａ」、「ａｎ」、「ｔｈｅ」など）は、特に言及しない限り、その複数形の概念をも含むことが理解されるべきである。また、本明細書において使用される用語は、特に言及しない限り、当該分野で通常用いられる意味で用いられることが理解されるべきである。したがって、他に定義されない限り、本明細書中で使用されるすべての専門用語および科学技術用語は、本発明の属する分野の当業者によって一般的に理解されるのと同じ意味を有する。矛盾する場合、本明細書（定義を含めて）が優先する。
【0009】
（定義）
最初に本発明において使用される用語および一般的な技術を説明する。
最初に本発明において使用される用語および一般的な技術を説明する。
【0010】
本明細書において、「対象」（英文ではｏｂｊｅｃｔ）とは、患者（ｐａｔｉｅｎｔ）または被検者もしくは被験者（ｓｕｂｊｅｃｔ）と同義に用いられ、疼痛測定および脳波測定などの本開示の技術が対象とする任意の生体または動物をいう。対象としては、好ましくは、ヒトであるがこれに限定されない。本明細書において、疼痛の推定を行う場合、「推定対象」とすることがあるが、これは対象などと同じ意味である。
【0011】
本明細書において「脳波」は当該分野で通常用いられるのと同義であり、頭皮上に１対の電極を置いて脳の神経活動にともなう電位差によって発生する電流をいう。脳波には、電流の時間的変化を導出記録した脳電図（ｅｌｅｃｔｒｏｅｎｃｅｐｈａｌｏｇｒａｍ，ＥＥＧ）を包含する。安静時には振幅約５０μＶ、周波数１０Ｈｚ前後の波が主成分をなすとされる。これをα波という。精神活動時にはα波は抑制され、振幅の小さい１７～３０Ｈｚの速波が現われるとされ、これはβ波という。浅い睡眠の時期にはα波はしだいに減少して４～８Ｈｚのθ波が現われるとされる。深い睡眠中は１～４Ｈｚのδ波が現われるとされる。これらの脳波は特定の振幅および周波数で記述することができ、本発明では、振幅の解析が重要でありうる。
【0012】
本明細書において「脳波データ」には、脳波に関する任意のデータであり（「脳活動量」、「脳特徴量」等ともいう）、振幅データ（ＥＥＧ振幅、周波数特性など）が含まれる。これらの脳波データを分析した「分析データ」は、脳波データと同様に用いることができることから、本明細書では、「脳波データまたはその分析データ」とまとめて呼ぶことがある。分析データとしては、例えば、脳波データの平均振幅やピーク振幅（例えば、Ｆｚ、Ｃｚ、Ｃ３、Ｃ４）、周波数パワー（例えば、Ｆｚ（δ）、Ｆｚ（θ）、Ｆｚ（α）、Ｆｚ（β）、Ｆｚ（γ）、Ｃｚ（δ）、Ｃｚ（θ）、Ｃｚ（α）、Ｃｚ（β）、Ｃｚ（γ）、Ｃ３（δ）、Ｃ３（θ）、Ｃ３（α）、Ｃ３（β）、Ｃ３（γ）、Ｃ４（δ）、Ｃ４（θ）、Ｃ４（α）、Ｃ４（β）、Ｃ４（γ）など）等を挙げることができる。もちろん、脳波データまたはその解析データとして通常使用される他のデータを排除するものではない。
【0013】
本明細書において「振幅データ」とは、「脳波データ」の一種であり、脳波の振幅のデータをいう。単に「振幅」ということもあり、「ＥＥＧ振幅」ともいう。このような振幅データは、脳活動の指標であることから、「脳活動データ」「脳活動量」などと称されることもある。振幅データは、脳波の電気信号を測定することによって得ることができ、電位（μＶ等で表示され得る）で表示される。振幅データとしては、平均振幅を使用することができるがこれに限定されない。
【0014】
本明細書において、「周波数パワー」とは、波形の周波数成分をエネルギーとして表したものであり、パワースペクトルともいう。周波数パワーは、高速フーリエ変換（ＦＦＴ）（離散フーリエ変換（ＤＦＴ）を計算機上で高速に計算するアルゴリズム）を利用して、時間領域のノイズの含まれる信号内に埋もれた信号の周波数成分を抽出し、これを計算することで算出することができる。信号のＦＦＴは、例えば、ＭＡＴＬＡＢの関数ｐｅｒｉｏｄｇｒａｍを用いて、その出力を正規化してパワースペクトル密度（ＰＳＤ）、またはパワーの測定元となるパワースペクトルを算出することができる。ＰＳＤは、時間信号のパワーが周波数についてどのように分布しているかを示すもので、単位はワット／Ｈｚである。ＰＳＤの各点を、その点が定義された周波数範囲にわたって（つまり、ＰＳＤの分解能帯域幅にわたって）積分して、パワースペクトルを計算する。パワースペクトルの単位はワットである。パワーの値は、周波数範囲にわたって積分せずに、直接パワースペクトルから読み取ることができる。ＰＳＤもパワースペクトルも実数なので、位相情報は何も含まれていないことになる。このように周波数パワーの計算はＭＡＴＬＡＢの標準的な機能で算出することができる。
【0015】
本明細書において「快」「不快」とは、行動を理解するための最も基本的な心的属性の1つであり、快をもたらす刺激には接近するが、不快をもたらす刺激からは遠ざかろうとする属性である。動物は快をもたらす刺激を獲得しようと接近するが、不快をもたらす刺激からは回避したり、不快な状態を維持する刺激からは逃避、もしくは、不快な状態を解消するような刺激を得ようと行動する。これらの接近・回避・逃避行動は環境に適応し、生存確率を高めるための基本的な行動の原理である。
【0016】
このように「快」と「不快」とは、対をなす概念である。
【0017】
本明細書において「ストレス」とは、外部からのさまざまな刺激（ストレッサー）によって心や体に負担がかかることで、心身に歪みが生じることをストレスという。ストレスは、不眠やうつ、胃痛や頭痛、さらには胃・十二指腸潰瘍など、心と体に不調を引き起こす。ストレスには、快いストレスである快ストレスと、不快なストレスである不快ストレスとがある。本明細書では、ストレスは快不快と同義で用いられる。ストレスは「身体的ストレス」と「精神的ストレス」に分けることができ、「身体的ストレス」は外からの刺激によって起きる「外的ストレス」と、特定の環境に置かれた場合に自分自身の内面から起きてくる「内的ストレス」に分けることができ、「精神的ストレス」は社会生活の中から生まれてくる「社会的ストレス」と、自分自身の心理的な要素から生まれてくる「心理的ストレス」に分けることができる。本明細書でいう心理学的パラダイムは、対象が特定の環境に置かれたとき、自身が持つ性格や心理状態に応じて、心理的ストレスを顕在的、もしくは潜在的に感じるように設定された検査方法である。
【0018】
本明細書において「痛み」および「疼痛」は同義であり、身体部分に傷害・炎症など一般に強い侵害のあるとき、これを刺激として生ずる感覚をいう。ヒトでは、強い不快感情を伴う感覚として一般感覚にも含まれる。加えて、皮膚痛覚などはある程度は外部受容の性格も備え、他の皮膚感覚や味覚と協同して、外物の硬さ・鋭さ・熱さ（熱痛）・冷たさ（冷痛）・辛さなどの質の判断に役立つとされる。ヒトの痛覚は皮膚・粘膜以外に身体のほとんどあらゆる部分（例えば、胸膜、腹膜、内臓（内臓痛覚、脳を除く）、歯、眼および耳など）に起こり得、いずれも脳において脳波またはその変動として感知され得る。この他、内臓痛に代表される内部痛覚もまた、痛覚に包含される）。内臓痛に対して上述した痛覚は体性痛という。体性痛および内臓痛に加えて、実際に障害されている部位と異なる部位の表面が痛くなるような現象である「関連痛」という痛覚も報告されており、本発明は、これらの多様な疼痛タイプを快不快という観点で分類することができる。
【0019】
痛覚には、感受性（痛閾）に個人差があり、痛刺激の起こり方や受容器部位の相違により、質的相違があり、鈍痛や鋭利痛などの分類があるが、本開示ではいずれの種類の痛覚でも測定、推定および分類することができる。また、速い痛覚（Ａ痛覚）および遅い痛覚（Ｂ痛覚）、（速い）局所的痛みおよび（遅い）瀰漫性痛みにも対応可能である。本発明は、痛覚異常過敏などの痛覚の異常症などにも対応し得る。痛みを伝える末梢神経には「Ａδ繊維」と「Ｃ繊維」の２つの神経繊維が知られており、例えば手をたたくと、始めの痛みはＡδ繊維の伝導により、局在が明確な鋭い痛み（一次痛；鋭利痛）が伝わる。その後、Ｃ繊維の伝導により、局在が不明確なじんじんとした痛み（二次痛；鈍痛）を感じるとされている。痛みは４～６週間以内持続する「急性疼痛」、と４～６週間以上持続する「慢性疼痛」に分類される。痛みは、脈拍や体温、血圧、呼吸と並ぶ重要なバイタルサインであるが、客観的データとして表示することは難しい。代表的な痛みスケールＶＡＳ（Vｉｓｕａｌ Aｎａｌｏｇｕｅ Sｃａｌｅ）やｆａｃｅｓ ｐａｉｎ ｒａｔｉｎｇ ｓｃａｌｅは主観的な評価法であり、患者間の痛みを比較することはできない。他方で、本発明者は、痛みの客観的評価のための指標として、末梢循環系の影響を受けにくい脳波に着目し、その痛み刺激に対する振幅／潜時の変化を観測すれば、痛みの種類（快不快）をも分類することが導かれた。瞬間刺激も持続刺激もこの分類が可能である。
【0020】
本発明では、強度自体よりも「治療が必要な」痛み（不快な痛み）かどうかということが区別できることが重要な点の一つである。したがって「治療」という概念を軸に「痛み」の類別化を明確にできることも重要である。
【0021】
本明細書において「主観的疼痛感覚レベル」とは、対象が有する疼痛感覚のレベルをいい、コンピュータ化された可視化アナログスケール（ＣＯＶＡＳ）等の慣用技術または他の公知技術、例えば、Ｓｕｐｐｏｒｔ Ｔｅａｍ Ａｓｓｅｓｍｅｎｔ Ｓｃｈｅｄｕｌｅ（ＳＴＡＳ－Ｊ）、Ｎｕｍｅｒｉｃａｌ Ｒａｔｉｎｇ Ｓｃａｌｅ（ＮＲＳ）、Ｆａｃｅｓ Ｐａｉｎ Ｓｃａｌｅ（ＦＰＳ）、Ａｂｂｅｙ ｐａｉｎ ｓｃａｌｅ（Ａｂｂｅｙ）、Ｃｈｅｃｋｉｎｌｉｓｔ ｏｆ Ｎｏｎｖｅｒｂａｌ Ｐａｉｎ Ｉｎｄｉｃａｔｏｒｓ（ＣＮＰＩ）、Ｎｏｎ－ｃｏｍｍｕｎｉｃａｔｉｖｅ Ｐａｔｉｅｎｔ‘ｓ Ｐａｉｎ Ａｓｓｅｓｓｍｅｎｔ Ｉｎｓｔｒｕｍｅｎｔ（ＮＯＰＰＡＩＮ）、Ｄｏｌｏｐｌｕｓ２などで表現することができる。
【0022】
本明細書において「条件」は、ある物事が成り立つ、または起こるもととなる事柄のうち、その原因ではないが、それを制約するものをいう。本明細書では、条件は「刺激」および「環境」を含む。
【0023】
本明細書において「刺激」とは、対象に対して何らかの反応を生じさせるものをいい、対象が生物の場合、生物やそのある部分の生理学的活性に，一時的な変化をもたらす要因をいう。理論に束縛されることは望まないが、例示的に説明すると、「環境」がその対象を空間的に包含し、何ら直接の焦点化した作用（例えば、身体部位への物理的刺激）を及ぼさないのに対して、「刺激」はその対象に対してそのような作用を及ぼす。
【0024】
「刺激」の具体的な例として例示される痛覚に関した事象の場合は痛覚を生じさせ得る任意の刺激を含む。例えば、電気刺激、冷刺激、熱刺激、物理的刺激、化学的刺激などが含まれる。本発明において、刺激はどのような刺激であってもよい。刺激の評価は、例えば、コンピュータ化された可視化アナログスケール（ＣＯＶＡＳ）等の慣用技術または他の公知技術、例えば、Ｓｕｐｐｏｒｔ Ｔｅａｍ Ａｓｓｅｓｓｍｅｎｔ Ｓｃｈｅｄｕｌｅ（ＳＴＡＳ－Ｊ）、Ｎｕｍｅｒｉｃａｌ Ｒａｉｎｇ Ｓｃａｌｅ（ＮＲＳ）、Ｆａｃｅｓ Ｐａｉｎ Ｓｃａｌｅ（ＦＰＳ）、Ａｂｂｅｙ ｐａｉｎ ｓｃａｌｅ（Ａｂｂｅｙ）、Ｃｈｅｃｋｌｉｓｔ ｏｆ Ｎｏｎｖｅｒｂａｌ Ｐａｉｎ Ｉｎｄｉｃａｔｏｒｓ（ＣＮＰＩ）、Ｎｏｎ－ｃｏｍｍｕｎｉｃａｔｉｖｅ Ｐａｔｉｅｎｔ‘ｓ Ｐａｉｎ Ａｓｓｅｓｓｍｅｎｔ Ｉｎｓｔｒｕｍｅｎｔ（ＮＯＰＰＡＩＮ）、Ｄｏｌｏｐｌｕｓ２などを用いて、主観的疼痛感覚レベルと対応させることができる。刺激の強度として採用され得る値としては、例えば、侵害受容閾値（ｎｏｃｉｃｅｐｔｉｖｅ ｔｈｒｅｓｈｏｌｄ；侵害受容線維に神経インパルスを生じさせる閾値）、痛覚判別閾値（ｐａｉｎ ｄｅｔｅｃｔｉｏｎ ｔｈｒｅｓｈｏｌｄ；ヒトが痛みとして知覚できる侵害刺激の強度）および痛覚許容限度閾値（ｐａｉｎ ｔｏｌｅｒａｎｃｅ ｔｈｒｅｓｈｏｌｄ；ヒトが実験的に許容できる侵害刺激の中でもっとも強い刺激強度）等を挙げることができる。
【0025】
心理学的な条件の場合、刺激としては、例えば、五感（視覚、聴覚、味覚、触覚、嗅覚）などで感じることができ、脳において情報処理がなされる任意の要因や、社会的ストレスなどの精神的に感知することができる任意の要因などを挙げることができる。
【0026】
本明細書において「文脈」（ｃｏｎｔｅｘｔ）および「環境」（ｅｎｖｉｒｏｎｍｅｎｔ）は同じ意味で用いられ、条件の一つであり、対象に対して何ら直接の焦点化した作用を及ぼさないが、その対象を取り囲む任意の条件をいい、刺激を対象に提示する際の、周囲の状況を言う。例えば、対象がおかれる空間に言及する場合は、温度、湿度、密度、明るさ、広さ、傾き、高さなどの空間特性を規定する多様なパラメータにより決まる３次元空間である。また、身体レベルでは、衣服や装飾品も身体空間を規定するパラメータであると理解できる。車いすを使う人の場合に見られるように、身体空間は外的環境に拡張する場合もある。本発明では、最初に、異なる刺激が不快度の違いを生じる基本的な状況に加え、対象がおかれる環境の違いや対象が受ける刺激の提示文脈の違いによって快不快が異なって感知されること、およびそれを信号解析結果として提示することができることが判明した。
【0027】
本明細書において快不快の「分類」は、種々の観点（例えば、痛覚の他、心理的なもの）で行うことができる。代表的な例としては、対象の痛みが「快」なのか「不快」なのかを分類することを含むが、これは、「我慢できる」痛みと「我慢できない」痛みとに分類する手法も想定され得る。
【0028】
本明細書において「陰性効果」とは、比較の基準となる条件の電位活動に比べて、関心対象となる検査条件の電位活動がマイナス方向にシフトすることを意味する。本明細書においては、基準となるストレスがない状態に比較して、顕在的、もしくは潜在的ストレスがある場合に、電位がマイナスの向きにシフトすることをいう。また、逆に、ストレスがない状態が検査の関心である場合は、陽性電位の増大が指標となり、陰性効果と陽性効果は、検査の基準条件の設定により言及の仕方が変わる。
【0029】
本明細書において「陰性占有率」または「占有率」とは、陰性電位の総時間（もしくは陰性電位データポイント数）／全体時間（もしくは、総データポイント数）×１００で算出する。占有率は、心理ストレスの指標であることが見いだされた。ここで時間としては、刺激後３００～８００ミリ秒を代表的に用いることができるが、これに限定されず、刺激後１５０ミリ秒より後の範囲、刺激後３００ミリ秒より後の範囲、刺激後３００ミリ秒～８００ミリ秒の範囲、刺激後３００ミリ秒～６００ミリ秒の範囲などを指すことができる。占有率は、電位効果の持続時間と類似するが、途中、効果が消失し、すぐに再開するような場合、効果の時間的特性を表現するのに適している。また、上記の３００～６００ミリ秒は、意識的判断などが収束する場合、陽性電位の効果が見られるが、認知作業の負荷が高い場合や抑制機能が続いている場合には、陰性効果が持続して見られ、心理ストレスの指標として用いるのに客観的な根拠がある（Ｋｉｎｇ ＪＷ， Ｋｕｔａｓ Ｍ．”Ｗｈｏ ｄｉｄ ｗｈａｔ ａｎｄ ｗｈｅｎ？ Ｕｓｉｎｇ ｗｏｒｄ－ ａｎｄ ｃｌａｕｓｅ－ｌｅｖｅｌ ＥＲＰｓ ｔｏ ｍｏｎｉｔｏｒ ｗｏｒｋｉｎｇ ｍｅｍｏｒｙ ｕｓａｇｅ ｉｎ ｒｅａｄｉｎｇ．”Ｊ． Ｃｏｇｎ． Ｎｅｕｒｏｓｃｉ． ７（３）： ３７６－３９５,１９９５）。
【0030】
本明細書において「自己複製」とは、あるサンプルやデータについて、そのデータのサンプル数をもとのデータに基づいて複製することをいう。自己複製の結果サンプル数が大幅に増加することになる。自己複製は、例えば、サンプルの分布特性（平均、ＳＤ）を用いて正規乱数やピアソンシステム乱数増幅などの乱数を発生させることによって、実現することができ、例えば、１００００サンプルの「自己複製型特徴量」を生成することができる。複製量は、１０サンプル、５０サンプル、１００サンプル、５００サンプル、１０００サンプル、５０００サンプル、１００００サンプル、２００００サンプル、５００００サンプル、１０００００サンプルなど任意の倍数が使用され得る。したがって、自己複製は「サンプル増幅法」とも呼ばれ、「サンプル増幅法」は、サンプルが少ない場合でも、その分布特性を用いて、大幅にサンプル数を増加させる技法である。例えば、実施例で用いられた事象関連電位は、通常、20、30、４０、５０、もしくは100回などの比較的少ない刺激提示数を用いるが、加算平均すれば１電極につき1脳波データにとどまる。しかしながら、本発明のサンプル増幅法を用いれば、20回の刺激提示で得た1電極の２０個の脳波データを、その分布特性（平均、散布度、尖度、歪度）に基づき、原理上無限にサンプル増幅させることができる。これに、機械学習（例えば、本実施例におけるLASSO正則化やSVM）により判別モデルのハイパーパラメタ（ラムダやCost、γ）を決定し、特徴量（例えば、振幅5個）の重みづけ係数とモデル切片を決定することができる。
【0031】
（好ましい実施形態）
以下に本発明の好ましい実施形態を説明する。以下に提供される実施形態は、本発明のよりよい理解のために提供されるものであり、本発明の範囲は以下の記載に限定されるべきでないことが理解される。従って、当業者は、本明細書中の記載を参照して、本発明の範囲内で適宜改変を行うことができることは明らかである。また、本発明の以下の実施形態は単独でも使用されあるいはそれらを組み合わせて使用することができることが理解される。
以下に本発明の好ましい実施形態を説明する。以下に提供される実施形態は、本発明のよりよい理解のために提供されるものであり、本発明の範囲は以下の記載に限定されるべきでないことが理解される。従って、当業者は、本明細書中の記載を参照して、本発明の範囲内で適宜改変を行うことができることは明らかである。また、本発明の以下の実施形態は単独でも使用されあるいはそれらを組み合わせて使用することができることが理解される。
【0032】
また、以下で説明する実施の形態は、いずれも包括的または具体的な例を示すものである。以下の実施の形態で示される数値、形状、材料、構成要素、構成要素の配置位置及び接続形態、ステップ、ステップの順序などは、一例であり、請求の範囲を限定する主旨ではない。また、以下の実施の形態における構成要素のうち、最上位概念を示す独立請求項に記載されていない構成要素については、任意の構成要素として説明される。
【0033】
（快不快分類法）
1つの局面において、本発明は、対象のストレスまたは快不快を判定する手段または判定値（機械学習などで判定器や分類器等を作成するときは、ｋｅｒｎｅｌと呼ばれる関数を用いた数式モデルである。ここに分類ラベルと特徴量を投入し、パラメータを決定する。このモデルが生み出す推定値を閾値、すなわち分類値でカテゴリ化（「１」、「０」など）し、実際の分類ラベルと照合して、判別精度を算出する。従って、手段、もしくは数式モデルと分類値は、プロセスと結果の関係にあるといえる。）を生成する方法であって、ａ）対象を少なくとも２つの環境下で、同一、もしくは強度などが同程度の刺激を付与し、各々の脳波データまたはその分析データを入手する工程と、ｂ）該少なくとも２つの環境下で得られた該脳波データまたはその分析データの相違と、該刺激に関する該対象の反応とを関連付ける工程と、ｃ）該関連付けに基づいて、対象のストレスまたは快不快を判定する快不快判定手段を生成する工程とを包含する、方法を提供する。
1つの局面において、本発明は、対象のストレスまたは快不快を判定する手段または判定値（機械学習などで判定器や分類器等を作成するときは、ｋｅｒｎｅｌと呼ばれる関数を用いた数式モデルである。ここに分類ラベルと特徴量を投入し、パラメータを決定する。このモデルが生み出す推定値を閾値、すなわち分類値でカテゴリ化（「１」、「０」など）し、実際の分類ラベルと照合して、判別精度を算出する。従って、手段、もしくは数式モデルと分類値は、プロセスと結果の関係にあるといえる。）を生成する方法であって、ａ）対象を少なくとも２つの環境下で、同一、もしくは強度などが同程度の刺激を付与し、各々の脳波データまたはその分析データを入手する工程と、ｂ）該少なくとも２つの環境下で得られた該脳波データまたはその分析データの相違と、該刺激に関する該対象の反応とを関連付ける工程と、ｃ）該関連付けに基づいて、対象のストレスまたは快不快を判定する快不快判定手段を生成する工程とを包含する、方法を提供する。
【0034】
別の局面において、本発明はまた、「判定値（手段）」生成～「実際の判定」までのすべての段階を含む方法を提供する。すなわち、本発明は、対象のストレスまたは快不快を判定する方法であって、ａ）対象を少なくとも２つの環境下で、同一、もしくは強度などが同程度の刺激を付与し、各々のモデル用の脳波データまたはその分析データを入手する工程と、ｂ）該少なくとも２つの環境下で得られた該モデル用の脳波データまたはその分析データの相違と、該刺激に関する該対象の反応とを関連付ける工程と、ｃ）該関連付けに基づいて、対象のストレスまたは快不快を判定する快不快判定手段を生成する工程とｄ）該対象から試験用の脳波データまたはその分析データを入手し、該快不快判定手段に適用し、該対象の快不快を判定する工程とを包含する、方法を提供する。
【0035】
別の局面では、本発明は、「判定値」はすでに出ている状態で、医療機器に実装する段階を含む方法を提供する。すなわち、本発明は、対象のストレスまたは快不快を判定する方法であって、ｃ）少なくとも２つの環境下での試験に基づく関連付けに基づいた、対象のストレスまたは快不快を判定する快不快判定手段を提供する工程とｄ）該対象から試験用の脳波データまたはその分析データを入手し、該快不快判定手段に適用し、該対象の快不快を判定する工程とを包含する、方法を提供する。
【0036】
別の局面では、本発明は、対象のストレスまたは快不快を判定する手段または値を生成する装置であって、Ａ）対象を少なくとも２つの環境下で、同一、もしくは強度などが同程度の刺激を付与し、各々の脳波データまたはその分析データを入手するデータ入手部と、Ｂ）該少なくとも２つの環境下で得られた該脳波データまたはその分析データの相違と、該刺激に関する該対象の反応とを関連付ける処理部と、Ｃ）該関連付けに基づいて、対象のストレスまたは快不快を判定する快不快判定手段または値を生成する判定手段生成部とを包含する、装置を提供する。本発明のＡ）、Ｂ）およびＣ）は、それぞれ工程ａ）、ｂ）、およびｃ）を実現するように構成される。
【0037】
別の局面では、本発明は、対象のストレスまたは快不快を判定する装置であって、Ａ）対象を少なくとも２つの環境下で同一、もしくは強度などが同程度の刺激を付与し、各々のモデル用の脳波データまたはその分析データを入手するデータ入手部と、Ｂ）該少なくとも２つの環境下で得られた該モデル用の脳波データまたはその分析データの相違と、該刺激に関する該対象の反応とを関連付ける処理部と、Ｃ）該関連付けに基づいて、対象のストレスまたは快不快を判定する快不快判定手段または値を生成する判定手段提供部とＤ）該対象から試験用の脳波データまたはその分析データを入手し、該快不快判定手段または値に適用し、該対象の快不快を判定する判定部とを包含する、装置を提供する。本発明のＡ）、Ｂ）、Ｃ）、およびＤ）は、それぞれ工程ａ）、ｂ）、ｃ）、およびｄ）を実現するように構成される。
【0038】
別の局面では、本発明は、対象のストレスまたは快不快を判定する装置であって、Ｃ）少なくとも２つの環境下での試験に基づく相関付けに基づいた、対象のストレスまたは快不快を判定する快不快判定手段または値を提供する判定手段提供部とＤ）該対象から試験用の脳波データまたはその分析データを入手し、該快不快判定手段に適用し、該対象の快不快を判定する判定部とを包含する、装置を提供する。本発明のＣ）、およびＤ）は、それぞれ工程ｃ）、およびｄ）を実現するように構成される。
【0039】
Ａ）に含まれる刺激機能は、刺激強度を複数種類提供することができる手段や機能を有していてもよい。そして、そのような刺激を、対象に対して付与することができるように構成される。Ａ）に含まれる脳波データ（例えば、振幅データ）の取得機能は、対象の脳波データを取得するように構成される。Ａ）はこれら２つの機能を別々の部分として実現してもよく、一体の部分として実現してもよい。また、これら以外の機能を持っていてもよい。
【0040】
Ｂ）では、シグモイド曲線のフィッティングなどの任意の関数フィッティングおよび２つ以上の文脈（環境）の相違と測定結果との関連付けの解析を行う機能を有する。
【0041】
Ｃ）は快不快判定手段または値の生成を行う機能を有し得る。
【0042】
Ｄ）では、測定用の脳波等を入手し、これをＣ）で得られた快不快判定手段または値に当てはめて快不快を判定する機能を有し得る。
【0043】
Ａ）、Ｂ）、Ｃ）およびＤ）の少なくとも２つの機能は、別々の装置、デバイス、ＣＰＵまたは端末などで実現されてもよく、１つの部分として実現されてもよい。通常は、１つのＣＰＵまたは計算装置において、これらの計算を実現するプログラムが組み込まれたか組み込まれ得る構成となっている。
【0044】
以下各工程を説明する。
【0045】
以下スキーム図を用いて快不快手法を説明する（図１５）。
【0046】
ａ）対象を少なくとも２つの環境下で、同一、もしくは強度などが同程度の刺激を付与し、各々のモデル用の脳波データまたはその分析データを入手する工程（Ｓ１００～Ｓ２００）について：
この工程では、対象に対して少なくとも２つの環境（好ましくは異なる環境）下で、同一（もしくは同一種類）、または同程度の刺激（例えば、痛み刺激）を複数回付与し、各々のモデル用の脳波データまたはその分析データを測定または入手する。または複数のレベルの刺激強度で該推定対象を刺激する工程（Ｓ１００）では、推定対象が複数のレベル（強さまたは大きさ）の複数の刺激（例えば、冷温刺激、電気刺激など）から選択して刺激され、該刺激強度に対応する該推定対象の脳波データ（脳活動データ、脳活動量などともいう（Ｓ１００）。例えば、振幅データ（ＥＥＧ振幅）、周波数特性などを含む）を取得する（Ｓ２００）。このような脳波データは、当該分野で周知の任意の手法を用いて取得することができる。脳波データは、脳波の電気信号を測定することによって得ることができ、振幅データなどとして電位（μＶ等で表示され得る）で表示される。周波数特性はパワースペクトル密度などで表示される。
この工程では、対象に対して少なくとも２つの環境（好ましくは異なる環境）下で、同一（もしくは同一種類）、または同程度の刺激（例えば、痛み刺激）を複数回付与し、各々のモデル用の脳波データまたはその分析データを測定または入手する。または複数のレベルの刺激強度で該推定対象を刺激する工程（Ｓ１００）では、推定対象が複数のレベル（強さまたは大きさ）の複数の刺激（例えば、冷温刺激、電気刺激など）から選択して刺激され、該刺激強度に対応する該推定対象の脳波データ（脳活動データ、脳活動量などともいう（Ｓ１００）。例えば、振幅データ（ＥＥＧ振幅）、周波数特性などを含む）を取得する（Ｓ２００）。このような脳波データは、当該分野で周知の任意の手法を用いて取得することができる。脳波データは、脳波の電気信号を測定することによって得ることができ、振幅データなどとして電位（μＶ等で表示され得る）で表示される。周波数特性はパワースペクトル密度などで表示される。
【0047】
好ましい実施形態では、本発明を実施するために、脳波データは、１）できるだけ少ない電極（２つ程度）で、２）毛髪のある頭皮を極力避け、３）寝ている状態でも記録できるような簡便な方法で実施することが好ましいが、必要に応じて電極の数は増加させてもよい（例えば、３つ、４つ、５つなどでもよい）。
【0048】
あるいは、モデル用の脳波データまたはその分析データ（例えば、本実施例における事象関連電位や、誘発電位）を、自己複製させて、該脳波データまたはその分析データの数を増加させてもよい。自己複製は、分布特性に基づいて実施されうる。自己複製はまた、正規乱数、ピアソンシステム乱数増幅などの乱数を発生させることによって実施され得る。あるいは、自己複製は、サンプルの分布特性（平均、標準偏差（ＳＤ））を用いて正規乱数を発生させ、適切な数のサンプル（例えば、１万、２万等）の「自己複製型特徴量」を増幅作成することができる。「サンプル増幅法」は、サンプルが少ない場合でも、その分布特性を用いて、大幅にサンプル数を増加させる技法であるといえる。
【0049】
ｂ）該少なくとも２つの環境下で得られた該モデル用の脳波データまたはその分析データの相違と、該刺激に関する該対象の反応とを関連付ける工程（Ｓ３００）について：
ａ）のようにして得られた脳波データまたはその分析データは、適宜の手法に基づいて、相違点を、刺激タイプや刺激呈示環境などの、刺激および環境に関するパラメータを含む条件パラメータ（例えば、痛み刺激不快度など）と関連付ける。
ａ）のようにして得られた脳波データまたはその分析データは、適宜の手法に基づいて、相違点を、刺激タイプや刺激呈示環境などの、刺激および環境に関するパラメータを含む条件パラメータ（例えば、痛み刺激不快度など）と関連付ける。
【0050】
Ｓ３００は、条件（あるいは刺激）の違い、環境（あるいは文脈）の違いに基づいて、快不快の違いを設定し、それに関連する特徴量を見つけ出す工程である。
【0051】
Ｓ４００は、その特徴量を用いて、条件の違いにラベルをつけて、シグモイドフィッティングや機械学習などを行う工程である。より詳細には以下のとおりである。
【0052】
Ｓ３００における、不快度に関わる手続きや環境の違いを脳波特徴量に関連付ける手続きは、例えば、以下のように行われる。日常環境で、不快度の異なる類似した環境ＡとＢがある場合、環境構成に関わる要因をできる限り枚挙する。その要因を環境ＡとＢの間で統計的に比較し、有意な違いがあるかどうかを検討する。もし有意な違いが見られたら、環境ＡとＢの不快度の違いに関わる有力候補とする。脳波特徴量は、振幅、潜時、効果の持続時間、分布、周波数パワーなど、時間的、空間的、もしくは両者の相互作用から構成される複雑な特性を持つ。したがって、環境ＡとＢの違いに関連付けられる特徴量が直観的に見つかる保証はない。したがって、時間的、空間的に分析して部分に分解し、環境ＡとＢに関わる特徴量を統計的に比較したり（ｔ検定や分散分析）、連続的な関係（相関関係や回帰）を調べたりすることで、刺激や環境条件と特徴量の関係を特定する。
【0053】
この工程では、工程ａで得られた脳波データを、フィルター処理、眼球運動補正、アーチファクト除去などの基本的な信号処理を行った後、条件パラメータに関連付けて該当部分の信号を抽出し、脳波特徴量を作成する（Ｓ３００）。これには、平均値（算術平均、幾何平均）、他の代表値（中央値、最頻値）、エントロピー、周波数パワー、ウェーブレット、もしくは平均、ならびに単一試行の事象関連電位の成分などが含まれる。
【0054】
ｃ）該関連付けに基づいて、対象のストレスまたは快不快を判定する快不快判定手段を生成する工程について（Ｓ４００）：
ｂ）で算出された環境および刺激に関する条件パラメータに関連づいた脳波特徴量を用い、フィッティングで得られたモデル曲線において、閾値や判定指数を設定する工程である。閾値電位や（陰性）占有率などの数値で閾値を設定して、これを判定指数とすることができる。
ｂ）で算出された環境および刺激に関する条件パラメータに関連づいた脳波特徴量を用い、フィッティングで得られたモデル曲線において、閾値や判定指数を設定する工程である。閾値電位や（陰性）占有率などの数値で閾値を設定して、これを判定指数とすることができる。
【0055】
Ｓ４００では、Ｓ３００における関連付けで特定された特徴量を用いて、既存、もしくは未知の刺激や環境を判別・推定するモデルを作成する。例えば、上記の特徴量の統計検定では、サンプル数が増えれば、環境ＡとＢの違いが小さかったり、なかったり、逆転しているようなサンプル数が増えても、「条件間に有意差がある」と判断される可能性が上がる。しかしながら、どの程度、該特徴量がサンプルを判別できるかという観点は含まれていない。しかしながら、痛みや心理的ストレスの判別では、個人が痛みやストレスを感じているかをできる限り精度高く区別できることが重要であり、統計的有意差の検出とは異なる有効性が要求され、また、目的とする。
【0056】
脳波特徴量を用いて、条件パラメータに合わせて２、３、もしくは３分類以上を行うために判別／推定モデルを作成するが、一つの方法としてはプロット図を作成し、シグモイド関数パターンなどの適宜なフィッティング関数に当てはめ（フィッティングさせ）る。当てはめは、当該分野で公知の任意の手法を用いて行うことができる。このような具体的なフィッティング関数としては、ボルツマン関数、二重ボルツマン関数、Ｈｉｌｌ関数、ロジスティック用量応答、シグモイドリチャード関数、シグモイドワイブル関数などを挙げることができるが、これらに限定されるものではない。これらのうち、標準ロジスティック関数は、シグモイド関数と呼ばれ、標準関数、もしくは変形が一般的であり、好ましい。
【0057】
必要に応じて該シグモイド関数パターンなどの適宜の関数パターンへの当てはめの回帰係数が所定以上の場合、該シグモイド曲線等に基づいて、快不快の判定のための閾値を設定することができる。ここでは、シグモイド曲線の場合変曲点（変調点）に基づいて生成することができるがこれに限定されない。必要に応じて、快不快のレベルの分類が最大限になるように疼痛分類値を較正（キャリブレーション）してもよい。閾値は、快不快レベルの算出または分類へと応用することができ、治療の効果の判定に利用することができる。
【0058】
したがって、１つの具体的な実施形態では、前記関連付けは、条件、例えば前記環境および前記刺激に基づいて、快不快の相違を設定し、該相違に関する特徴量を見出すことを含み、および前記快不快判定手段の生成は、該特徴量を用いて、該刺激の相違に識別するラベルを付すことを含む。特定の実施形態では、前記快不快判定手段の生成は、シグモイドフィッティングまたは機械学習によって達成される。
【0059】
１つの実施形態では、前記ストレスまたは快不快は、１）前記対象に与えられる刺激および前記環境の両方が異なるか、２）前記対象に与えられる刺激が異なるが前記環境が同じであるか、または３）前記対象に与えられる刺激は同じであるが前記環境が異なる。より具体的に説明すると、１）は、例えば、提示刺激のタイプも刺激が提示される空間も異なり、その違いが不快度と関連する特性に基づくものであり、２）は、例えば、提示刺激タイプや刺激の強度が異なるが、提示される文脈（環境）が異なり、それが不快度の違いに関連するものであり、３）は、例えば、刺激の提示方法や課題内容は全く同じであるが、対象が配置される空間特性が異なり、それが不快度の違いに関連するものである。
【0060】
実際の医療装置では、ａ）～ｃ）を行うように構成されてもよいが、予めこの判定器または判定値を設定しておいてもよく、その場合、ｃ）工程は、ｃ’）少なくとも２つの環境下での試験に基づく関連付けに基づいた、対象のストレスまたは快不快を判定する快不快判定手段を提供する工程であってもよい。
【0061】
同一の被験者が対象者になる場合、以前の快不快判定手段（値など）を用いて、判定器や判定値を継承、もしくは更新するような工程を含んでいてもよい。
【0062】
ｄ）該対象から試験用の、未知の状態に関わる脳波データまたはその分析データを入手し（Ｓ４５０）、該快不快判定手段に適用し、該対象の快不快を判定する工程について（Ｓ５００）：
この工程では、判定器や閾値をもとに、実際に測定した対象の未知の状態に関わる測定値、例えば、脳波データまたは分析データから、その判定器や閾値に対応する数値を算出し、これを判定器や閾値と比較して快不快の有無またはそのレベルを決定する。
この工程では、判定器や閾値をもとに、実際に測定した対象の未知の状態に関わる測定値、例えば、脳波データまたは分析データから、その判定器や閾値に対応する数値を算出し、これを判定器や閾値と比較して快不快の有無またはそのレベルを決定する。
【0063】
該対象の脳波データ（例えば振幅データ）を得る工程である（Ｓ４５０）。この工程は、何らかの刺激を受けているあるいは処置されているかどうかにかかわらず、測定を意図する対象から、対象の状態が未知の脳波データを得る工程であり、脳波データを取得できる手法であれば、どのような手法であってもよい。ａ）工程で使用される脳波データ取得と同じ手法を用いることができ、通常は同一の手法を用いる。そして、Ｓ５００にあるように、快不快判定器や判定値に適用し、該対象の快不快を判定する。所定の快不快判定手段または値は、対象に関して判別／推定するレベルと紐付けられて「不快度（快不快）判定器」、もしくは「不快度（快不快）判定予測器」と呼ばれる。閾値よりも不快側の数値では不快と判定または予測し、快側での数値では快と判定または予測する。
【0064】
１つの実施形態では、脳波データ（例えば、振幅データ）の前記快不快判定手段または値への当てはめは、平均値で行うことができる。このような平均値は、１５秒～２００秒、さらに数時間にわたるデータ記録の場合、２００秒を超える（例えば、３００秒、５００秒、６００秒、９００秒、１２００秒など）平均値であってもよい。短時間（例えば、１秒以下）の外的、もしくは内的出来事に同期する事象関連脳活動の場合は、出来事が発生した後、１秒以内の平均値や代表値や時系列変化データなどであってもよい。また、判別器や判定値において用いられるデータは、標準化や正規化したものが望ましい。また、異なる日時のデータを比較する場合は、共通のベースライン環境や刺激タイプを設定し、それに対する変動や効果を算出して判別器や判定値に関連付けることが望ましい。１つの実施形態では、前記関連付けは、前記のシグモイド関数を判別器として用いることで実現され得る。
【0065】
別の実施形態では、判別／推定する未知のストレスまたは快不快は、前記対象が痛みを感じているときのものである。
【0066】
さらに別の実施形態では、「痛み」の物理的強度と、「ストレス（快不快）」が別パラメータで表示できる。この実施形態では、二次元パラメータ化することができる。詳細に説明すると、前記快不快判定手段は、前記対象の疼痛の強度と、前記対象のストレスまたは快不快の物理的強度によらないレベルとを識別する。ここでは、痛みを感じる痛覚や色を感じる色覚としての物理量とは別のストレスまたは快不快のレベルを示すパラメータが存在することが見いだされた。物理的なパラメータ強度以外の指数として、刺激を同じ、もしくは出来る限り統制するが、ストレスと不快度に影響を与える環境パラメータを変えることにより、このような快不快またはストレス指数を算出することができる。これらの指数は、例えば、心理パラダイムにおける陰性電位活動の特定の時間範囲における占有率などを挙げることができるが、この前提として、陰性脳波活動が認知的負荷（例．作業記憶負荷）や抑制機能の増大により振幅が増大したり、効果時間が持続するような性質に基づく。
【0067】
１つの実施形態では、前記脳波データまたはその分析データは、データ記録位置として国際１０－２０基準、もしくはその拡張基準に準拠したＦ３、Ｆ４、Ｃ３、Ｃ４、Ｐ３、Ｐ４などの前頭部から頭頂部、さらには後頭部にかけた頭皮上位置を電極位置として含み、あるいは、特定の均一な距離（例．２．５ｃｍなど）でおかれた位置を網羅してもよく、記録、ならびに解析の時間帯は、短時間の事象関連電位活動の場合は、０－１００、１００－２００、２００－３００、３００－４００、４００－５００、５００－６００、６００－７００、７００－８００ミリ秒（ｍｓ）、あるいは、より細かい時間区分（１０ミリ秒など）や、より長い時間帯（時に、数秒間に及ぶ）であってもよい。前記脳波データまたはその分析データは、これらの組合せから選択される、少なくとも１つの脳波特徴量を含む。
【0068】
さらに別の実施形態では、前記脳波特徴量は、Ｆｐ１、Ｆｐ２、Ｆｐｚ、Ｆ３、Ｆ４、Ｆｚ、Ｃ３、Ｃ４、Ｃｚ、Ｐ３、Ｐ４およびＰｚからなる群より選択される少なくとも１つを含み、例えば、平均振幅Ｆｚ、Ｃ３およびＣ４；ならびに周波数Ｆｚ（δ）、Ｆｚ（β）、Ｃｚ（δ）、Ｃ３（θ）およびＣ４（β）を含む。Ｃｚ（振幅），Ｃ３（α），Ｃｚ（β），Ｆｚ（δ），Ｃｚ（γ）を含むことが好ましいがこれに限定されない。
【0069】
図１６には、本発明の装置の模式図が記載される。このうちこの実施形態は快不快判定器（手段）の生成の場合で、１０００～３０００が関与する。刺激呈示部１０００はＡ）に該当し、これは、刺激が提示される環境や刺激タイプに関する情報が脳波データ取得部２０００、および快不快判定値生成部３０００に連絡される。脳波データ取得部２０００では、刺激呈示部から対象（１５００）に発出された刺激に同期した脳波データを得られるように、対象（１５００）に連結されたか、またはされ得る脳波計を備えるか脳波計に連結されるように構成される（２５００）。
【0070】
図１７は、１つの実施形態の快不快判定システム５１００の機能構成を示すブロック図である（なお、この構成図のいくつかは任意の構成部分であり省略され得ることに留意されたい）。システム５１００は、脳波測定部５２００を備え、この脳波測定部は脳波記録センサー５２５０および必要に応じて脳波増幅部５２７０を内部に備えるか、外部で接続し、快不快判定装置５３００にて快不快の信号処理および判別／推定が行われる。快不快判定装置５３００では、脳波信号処理部５４００で脳波信号を処理し、（必要な場合は、脳波特徴量抽出部５５００で脳波特徴量を抽出、ときにサンプル増幅し）快不快判定部５６００で快不快、または不快度を判別／推定し、（必要に応じて）判別レベル可視化部５８００で不快度を可視化する。また、内部、もしくは外部に、刺激装置部５９００を備え、この刺激装置部５９００は、対象の快不快および不快判別器作成と実際の未知の不快度レベルの判別のために刺激情報（刺激タイプや環境情報などを）を送る。刺激装置部５９００は、刺激呈示部５９２０の他必要に応じて刺激情報可視化部５９６０を備え、刺激や環境に係る画像や数字などの情報を表示してもよい。また、快不快判定システムは、判別器や判定値を生成する生成部５７００を装置５３００の外部、もしくは内部に備えていてもよい。
【0071】
このように快不快判定システム５１００は、脳波測定部５２００と、快不快判定装置５３００とを備え、必要に応じて刺激装置部５９００を備える。快不快判定装置５３００は、例えば、プロセッサ及びメモリを備えるコンピュータによって実現される。この場合、快不快判定装置５３００は、メモリに格納されたプログラムがプロセッサによって実行されたときに、プロセッサを、必要に応じて脳波増幅部５２７０、脳波信号処理部５４００、(必要に応じて）快不快判定部５６００、(必要に応じて）判別レベル可視化部５８００などとして機能させる。必要に応じて刺激や環境情報の可視化もさせる。また、本発明のシステム５１００または快不快判定装置５３００は、例えば、専用電子回路によって実現されてもよい。専用電子回路は、１つの集積回路であってもよいし、複数の電子回路であってもよい。脳波データ取得部および快不快判定値生成部は、この快不快判定装置と同様の構成をとっていてもよい。
【0072】
脳波測定部５２００は、脳波計（脳波記録センサー５２５０）を介して推定対象から複数回の脳波測定を行うことにより複数の脳波データを取得する。推定対象とは、刺激や環境によって脳波に変化が生じる生体であり、ヒトに限定される必要はない。
【0073】
快不快判定部５６００は、判定値を用いて不快度を判別／推定する。判別器や判定値が外部、もしくは内部であらかじめ生成されない場合は、生成も行う。判別器や判定値を生成する部門は、快不快判定値生成部５７００として装置５３００の外部、もしくは内部に備えることができる。不快度判定値は、複数の脳波データの振幅から、不快度を推定または分類するためのものである。つまり、快不快判定部５６００、もしくは快不快判定値生成部５７００は、脳波データから対象の不快度を推定または分類するための判定値を生成することができる。
【0074】
脳波記録センサー５２５０は、推定対象の脳内で発生する電気活動を頭皮上の電極で計測する。そして、脳波記録センサー５２５０は、計測結果である脳波データを出力する。脳波データは必要に応じて増幅することができる。
【0075】
図１６に基づいてさらに説明する。判定部まで備える局面を説明する。図１６では、快不快判定部４０００の他、脳波データ取得部２０００が参照される。点線は判別モデル作成の手順を示し、実線は、実際の疼痛レベルの判別／推定を行う手順を示す。この場合、（快不快判定値生成）の節で説明したように、脳波データは、対象１５００から脳波計を介して取得することができる。すなわち、脳波データ取得部２０００は、対象１５００に連結され得るように構成され、脳波データ取得部２０００では、対象（１５００）から得られる脳波データを取得するように、対象（１５００）に連結されたかされ得る脳波計を備えるか脳波計に連結されるように構成される（２５００）。快不快判定部４０００では、予め快不快判定値が格納されているか、別途生成したデータを受容するように構成され、必要に応じて参照することができるように構成される。そのような連結の構成は、有線であっても無線であってもよい。あらかじめ格納される快不快判定値は、快不快判定値生成部３０００において、例えば、特徴量の判別器（シグモイド関数フィッティングなど）に基づき生成される。
【0076】
図１７は、１つの実施形態の快不快判定システム５１００の機能構成を示すブロック図である。システム５１００は、脳波測定部５２００を備え、この脳波測定部は脳波記録センサー５２５０および必要に応じて脳波増幅部５２７０を内部に備えるか、外部で接続し、快不快判定装置５３００にて痛みの信号処理および判別／推定が行われる。快不快判定装置５３００では、脳波信号処理部５４００で脳波信号を処理し、(必要に応じて）快不快判定部５６００で痛みを推定／判別し、(必要に応じて）判別レベル可視化部５８００で痛みを可視化する。また、内部、もしくは外部に、刺激装置部５９００を備え、この刺激装置部５９００は、対象の不快度判別器作成のために寄与する。判定値は、快不快判定値生成部５７００であらかじめ作成されてもよい。
【0077】
このように快不快判定システム５１００は、脳波測定部５２００と、快不快判定装置５３００とを備える。快不快判定装置５３００は、例えば、プロセッサ及びメモリを備えるコンピュータによって実現される。この場合、快不快判定装置５３００は、メモリに格納されたプログラムがプロセッサによって実行されたときに、プロセッサを必要に応じて脳波増幅部５２７０、脳波信号処理部５４００、(必要に応じて）快不快判定部５６００、(必要に応じて）判別レベル可視化部５８００などとして機能させる。必要に応じてリファレンス刺激発生および可視化もさせることができる。また、本発明のシステム５１００または装置５３００は、例えば、専用電子回路によって実現されてもよい。専用電子回路は、１つの集積回路であってもよいし、複数の電子回路であってもよい。脳波データ測定部および快不快判定値生成部３０００（図１６参照）は、この疼痛推定装置と同様の構成をとっていてもよいし、外部に構成されてもよい。
【0078】
脳波測定部５２００は、脳波計（脳波記録センサー５２５０）を介して推定対象から複数回の脳波測定を行うことにより複数の脳波データを取得する。推定対象とは、痛みによって脳波に変化が生じる生体であり、ヒトに限定される必要はない。
【0079】
快不快判定部５６００は、快不快判定値生成部３０００（図１６参照）により作成された疼痛分類値に基づいて、複数の脳波データの振幅から、痛みの大きさを推定または分類する。つまり、快不快判定部５６００は、判定値に基づいて、脳波データから対象の疼痛を推定または分類する。
【0080】
脳波記録センサー５２５０は、推定対象の脳内で発生する電気活動を頭皮上の電極で計測する。そして脳波記録センサー５２５０は、計測結果である脳波データを出力する。脳波データは必要に応じて増幅することができる。
【0081】
次に、以上のように構成された装置の処理または方法について説明する。図１５は、一連の処理を示すフローチャートである。この局面では、Ｓ４００～Ｓ６００までが関与しうる。Ｓ４００で、刺激や環境情報（条件パラメータ）を用いて不快度判定値（不快度判定器／不快度判定器ともいう）が生成された後の工程である。あるいは、別途快不快度判定値が入手可能な場合（以前に取得して格納してある場合など）である。
【0082】
この不快度判定値は、作成後、快不快判定部４０００（図１６参照）に予め格納されていてもよく、また、値データを受容できるように快不快判定部４０００が構成されていてもよい。あるいは快不快判定値生成部３０００が付属する場合は、その生成部に格納されてもよく、別途記録媒体が配置されていてもよい。通信でこの値を受容することもできる。
【0083】
次に、対象から脳波データを取得する（Ｓ４５０）（図１５参照）。この脳波データの取得はＳ２００で説明したものと同様の技術を用いることができ、同様の実施形態を採用することができるが、Ｓ２００と同じ装置またはデバイスを用いる必要はなく、異なっていても同じであってもよい。
【0084】
次に、Ｓ４５０で得られた脳波データ（例えば、振幅データ）を、不快度判定値に当てはめて、その脳波データに対応する不快度の判別／推定を行う（Ｓ５００）（図１５参照）。このような不快度判定は、予め所定の値が出た場合に特定の文言（好き、嫌い等）を表示するあるいは発声するように構成されてもよく、実際の値と不快度判定値とを並列して表示させユーザー（臨床医）が検討できるようにしてもよい。
【0085】
図１８は、図１７のブロック図に、判別器生成（例えば、シグモイド関数フィッティング）の過程を含めて拡張し、作業内容を記入した快不快判定システム５１００の図である。システム５１００は、脳波測定部５２００を備え、脳波計５２２０に連結され、収集した脳波データから特徴量抽出部５５００で平均値などの脳波特徴量を必要に応じて得たり、サンプル増幅する。あらかじめ不快度判定値が生成される場合は、快不快判定装置５３００の外部、もしくは内部にある快不快判定値生成部３０００で、判別器、例えば、シグモイド関数フィッティングにより、不快度判定値が生成され、快不快判定部５６００に送られ、格納される。実際の未知の刺激タイプ、もしくは環境の不快度判定では、刺激呈示部５９２０の刺激呈示、もしくは表示と同期した脳波データが脳波計５２２０から測定部５２００に送られた後、特徴量抽出部５５００で脳波特徴量が作られ、快不快判定部５６００に送られ、不快度判定値を用いて、未知の刺激や環境の不快度を判別／推定する。(必要に応じて）判別レベル可視化部５８００で不快度を可視化する。このような一連の過程は、プロセッサ及びメモリを備えたコンピュータや携帯端末で実現されたり、専用電子回路によって実現されてもよい。専用電子回路は、１つの集積回路であってもよいし、複数の電子回路であってもよい。また、ソフトウェアで実現され、必要なハードウェアを制御することで実現しても良い。
【0086】
（心理ストレス特徴量）
１つの局面において、本発明は、潜時における波形について、基準波形と比較した場合の陰性レベルによって快不快を判定する方法を提供する。
１つの局面において、本発明は、潜時における波形について、基準波形と比較した場合の陰性レベルによって快不快を判定する方法を提供する。
【0087】
１つの実施形態では、前記陰性レベルは、対象が処理する刺激や環境に依存するが、刺激後１５０ミリ秒より後の範囲の波形を基準とし、好ましくは、刺激後３００ミリ秒より後の範囲の波形を基準とし、刺激後３００ミリ秒～８００ミリ秒の範囲の波形を基準とし、あるいは、刺激後３００ミリ秒～６００ミリ秒の範囲の波形を基準とし、好ましくは、刺激後３００ミリ秒～８００ミリ秒の範囲の（陰性）占有率を基準とすることが有利であり得る。また、占有率のほかに、追加または代替として、特定の時間範囲の平均振幅、振幅の積分値、ピーク振幅、ピーク潜時、効果の持続時間、周波数パワーを含む特徴量も使うことができる。
【0088】
この陰性レベルで規定された前記快不快判定手段を用いて、心理ストレスを判定することができる。
【0089】
（他の実施の形態）
以上、本発明の１つまたは複数の態様に係る痛み推定装置について、実施の形態に基づいて説明したが、本発明は、この実施の形態に限定されるものではない。本発明の趣旨を逸脱しない限り、当業者が思いつく各種変形を本実施の形態に施したものや、異なる実施の形態における構成要素を組み合わせて構築される形態も、本発明の１つまたは複数の態様の範囲内に含まれてもよい。
以上、本発明の１つまたは複数の態様に係る痛み推定装置について、実施の形態に基づいて説明したが、本発明は、この実施の形態に限定されるものではない。本発明の趣旨を逸脱しない限り、当業者が思いつく各種変形を本実施の形態に施したものや、異なる実施の形態における構成要素を組み合わせて構築される形態も、本発明の１つまたは複数の態様の範囲内に含まれてもよい。
【0090】
例えば、上記各実施の形態において、脳波データの振幅値として、ピークピーク値が用いられうるが、これに限定されない。例えば、振幅値として、単なるピーク値が用いられてもよい。
【0091】
なお、上記実施形態では、脳波振幅の上限値Ａｍａｘに対応する不快度の大きさの値Ｐｍａｘが１となり、脳波振幅の上限値Ａｍｉｎに対応する痛みの大きさの値Ｐｍｉｎが０となるように、不快度の大きさの値の範囲を設定したが、これに限定されない。例えば、０～１００で痛みの大きさが表されてもよい。この場合、快不快判定部５６００は、不快度の大きさの値Ｐｘを以下の式で推定すればよい。
【0092】
Ｐｘ＝Ｐｍａｘ×（Ａｘ－Ａｍｉｎ）／（Ａｍａｘ－Ａｍｉｎ）
また、上記では、複数の脳波データを分析することにより快不快判定値の生成例として、曲線あてはめ（ｃｕｒｖｅ ｆｉｔｔｉｎｇ）を説明したが、これに限定されない。また、脳波振幅の上限値は、予め定められた値が用いられてもよい。例えば、予め定められた値は、例えば５０μＶ～１００μＶであり、実験的又は経験的に定められればよい。このように通常の解析では、アーチファクト除去方法として、プラスマイナス５０μＶから１００μＶぐらいのデータを排除するが、このようなアーチファクト除去は必要に応じて本発明でも実施してもよい。
また、上記では、複数の脳波データを分析することにより快不快判定値の生成例として、曲線あてはめ（ｃｕｒｖｅ ｆｉｔｔｉｎｇ）を説明したが、これに限定されない。また、脳波振幅の上限値は、予め定められた値が用いられてもよい。例えば、予め定められた値は、例えば５０μＶ～１００μＶであり、実験的又は経験的に定められればよい。このように通常の解析では、アーチファクト除去方法として、プラスマイナス５０μＶから１００μＶぐらいのデータを排除するが、このようなアーチファクト除去は必要に応じて本発明でも実施してもよい。
【0093】
また、刺激呈示部５９２０（図１７参照）が推定対象５０９９に与える刺激は、刺激タイプや呈示環境に応じて推定対象５０９９が感じる不快度の大きさが変わるのであれば、どのような種類の刺激が与えられてもよいが、刺激の物理的強度から分離された不快度を判定する場合は、特に、同じ刺激や極力類似した刺激が異なる環境や条件で呈示されることが望ましい。これを「文脈依存型不快度検出方法」もしくは「文脈依存型リファレンス検査手法」と呼ぶ。例えば、本発明の実施例で示すように、１）同じ高温刺激（例えば、４０℃）を、ことなる刺激呈示文脈（３８℃文脈、もしくは４８℃文脈）で呈示し、異なる不快度を検出する手法や、２）同じ認知処理を異なる状況（例えば、監視あり、なし）で行い、潜在的な不快度を検出する方法などがある。
【0094】
また、上記各実施の形態における快不快判定装置が備える構成要素の一部または全部は、１個のシステムＬＳＩ（Ｌａｒｇｅ Ｓｃａｌｅ Ｉｎｔｅｇｒａｔｉｏｎ：大規模集積回路）から構成されているとしてもよい。例えば、図１７に示すように、快不快判定装置５３００は、必要に応じて測定部５２００と必要に応じて刺激呈示部５９２０とを有するシステムＬＳＩから構成されてもよい。
【0095】
システムＬＳＩは、複数の構成部を１個のチップ上に集積して製造された超多機能ＬＳＩであり、具体的には、マイクロプロセッサ、ＲＯＭ（Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）、ＲＡＭ（Ｒａｎｄｏｍ Ａｃｃｅｓｓ Ｍｅｍｏｒｙ）などを含んで構成されるコンピュータシステムである。ＲＯＭには、コンピュータプログラムが記憶されている。前記マイクロプロセッサが、コンピュータプログラムに従って動作することにより、システムＬＳＩは、その機能を達成する。
【0096】
なお、ここでは、システムＬＳＩとしたが、集積度の違いにより、ＩＣ、ＬＳＩ、スーパーＬＳＩ、ウルトラＬＳＩと呼称されることもある。また、集積回路化の手法はＬＳＩに限るものではなく、専用回路または汎用プロセッサで実現してもよい。ＬＳＩ製造後に、プログラムすることが可能なＦＰＧＡ（Ｆｉｅｌｄ Ｐｒｏｇｒａｍｍａｂｌｅ Ｇａｔｅ Ａｒｒａｙ）、あるいはＬＳＩ内部の回路セルの接続や設定を再構成可能なリコンフィギュラブル・プロセッサを利用してもよい。
【0097】
さらには、半導体技術の進歩または派生する別技術によりＬＳＩに置き換わる集積回路化の技術が登場すれば、当然、その技術を用いて機能ブロックの集積化を行ってもよい。バイオ技術の適用等が可能性としてありえる。
【0098】
また、本発明の一態様は、このような快不快判定値生成、快不快判定装置だけではなく、痛み推定装置に含まれる特徴的な構成部をステップとする疼痛分類値生成、痛み判別・分類方法であってもよい。また、本発明の一態様は、快不快判定値生成、快不快判定方法に含まれる特徴的な各ステップをコンピュータに実行させるコンピュータプログラムであってもよい。また、本発明の一態様は、そのようなコンピュータプログラムが記録された、コンピュータ読み取り可能な非一時的な記録媒体であってもよい。
【0099】
なお、上記各実施の形態において、各構成要素は、専用のハードウェアで構成されるか、各構成要素に適したソフトウェアプログラムを実行することによって実現されてもよい。各構成要素は、ＣＰＵまたはプロセッサなどのプログラム実行部が、ハードディスクまたは半導体メモリなどの記録媒体に記録されたソフトウェアプログラムを読み出して実行することによって実現されてもよい。ここで、上記各実施の形態の痛み推定装置などを実現するソフトウェアは、本明細書において上述したプログラムであり得る。
【0100】
本明細書において「または」は、文章中に列挙されている事項の「少なくとも１つ以上」を採用できるときに使用される。「もしくは」も同様である。本明細書において「２つの値の範囲内」と明記した場合、その範囲には２つの値自体も含む。
【0101】
本明細書において引用された、科学文献、特許、特許出願などの参考文献は、その全体が、各々具体的に記載されたのと同じ程度に本明細書において参考として援用される。
【0102】
以上、本発明を、理解の容易のために好ましい実施形態を示して説明してきた。以下に、実施例に基づいて本発明を説明するが、上述の説明および以下の実施例は、例示の目的のみに提供され、本発明を限定する目的で提供したのではない。従って、本発明の範囲は、本明細書に具体的に記載された実施形態にも実施例にも限定されず、特許請求の範囲によってのみ限定される。
【実施例】
【0103】
以下に実施例を記載する。以下の実施例で用いる対象の取り扱いは、必要な場合、大阪大学において規定される基準を遵守し、臨床研究が関係する場合はヘルシンキ宣言およびＩＣＨ－ＧＣＰに準拠して行った。
【0104】
（実施例１：異なる不快度を持つ、同程度の疼痛強度を持つ刺激タイプの判別）
本実施例では、不快度が異なる２タイプの痛み刺激の判別を機械学習（ＳＶＭ－ＲＦＥ）により行った。２タイプの刺激は、不快でなく、マッサージで使われるような電気刺激と、不快度が高い低温刺激を用いた。両者は刺激強度が主観報告により統制された。
本実施例では、不快度が異なる２タイプの痛み刺激の判別を機械学習（ＳＶＭ－ＲＦＥ）により行った。２タイプの刺激は、不快でなく、マッサージで使われるような電気刺激と、不快度が高い低温刺激を用いた。両者は刺激強度が主観報告により統制された。
【0105】
（方法）
（参加者）
２０代から７０代の４１名の健常な成人被験者の同じグループが、低温刺激パラダイムと電気刺激パラダイムの実験に参加した。参加者らは、実験前にインフォームドコンセントに同意した。全参加者は、神経性および／もしくは精神性の病気、または臨床薬物治療条件下における急性および／もしくは慢性疼痛を経験していないことを自己報告した。本実施例を、大阪大学医学部付属病院の倫理委員会ならびにヘルシンキ宣言による承認下で実施した。
（参加者）
２０代から７０代の４１名の健常な成人被験者の同じグループが、低温刺激パラダイムと電気刺激パラダイムの実験に参加した。参加者らは、実験前にインフォームドコンセントに同意した。全参加者は、神経性および／もしくは精神性の病気、または臨床薬物治療条件下における急性および／もしくは慢性疼痛を経験していないことを自己報告した。本実施例を、大阪大学医学部付属病院の倫理委員会ならびにヘルシンキ宣言による承認下で実施した。
【0106】
（実験刺激、ならびに手続）
図１に実験方法の概要を示す。温度刺激システム(Ｐａｔｈｗａｙ：Ｍｅｄｏｃ Ｃｏ．，Ｌｔｄ．，ＲａｍａｔＹｉｓｈａｉ， Ｉｓｒａｅｌ)を使用して、低温刺激を参加者の右前腕に与えた。低温刺激は、３レベルの温度強度（１０℃、０℃、－１０℃）を包含していた。各温度レベルは３回の刺激からなり、刺激間間隔（ＩＳＩ）は２０秒であった。各刺激は、５秒間持続するプラトーを有し、基準温度（３５℃）からの上昇および下降待機時間は、約１０秒であった。各レベルにおける３回の刺激の後、ブロック間間隔は１００秒間持続した。参加者は、コンピュータ化された可視化アナログスケール（ＣＯＶＡＳ）上で、０から１００の範囲（０：「無痛」；１００：「耐えられない疼痛」）で疼痛強度を連続的に評価した。ＣＯＶＡＳデータを刺激強度の変更と同時に記録した。電気刺激パラダイムでは、最初に、個人ごとに、低温刺激１０℃、０℃、－１０℃に対応する電気刺激強度を知覚・痛覚定量分析装置（ＰＡＩＮＶＩＳＩＯＮ ＣＯ．， Ｌｔｄ．，Ｏｓａｋａ， Ｊａｐａｎ）を用いて特定した。その後、電気刺激の弱・中・強の３レベル、各３刺激を低温刺激と同じ呈示方法で呈示した。参加者は、ＣＯＶＡＳを用いて痛みレベルの主観評価を刺激呈示と並行して行った。
図１に実験方法の概要を示す。温度刺激システム(Ｐａｔｈｗａｙ：Ｍｅｄｏｃ Ｃｏ．，Ｌｔｄ．，ＲａｍａｔＹｉｓｈａｉ， Ｉｓｒａｅｌ)を使用して、低温刺激を参加者の右前腕に与えた。低温刺激は、３レベルの温度強度（１０℃、０℃、－１０℃）を包含していた。各温度レベルは３回の刺激からなり、刺激間間隔（ＩＳＩ）は２０秒であった。各刺激は、５秒間持続するプラトーを有し、基準温度（３５℃）からの上昇および下降待機時間は、約１０秒であった。各レベルにおける３回の刺激の後、ブロック間間隔は１００秒間持続した。参加者は、コンピュータ化された可視化アナログスケール（ＣＯＶＡＳ）上で、０から１００の範囲（０：「無痛」；１００：「耐えられない疼痛」）で疼痛強度を連続的に評価した。ＣＯＶＡＳデータを刺激強度の変更と同時に記録した。電気刺激パラダイムでは、最初に、個人ごとに、低温刺激１０℃、０℃、－１０℃に対応する電気刺激強度を知覚・痛覚定量分析装置（ＰＡＩＮＶＩＳＩＯＮ ＣＯ．， Ｌｔｄ．，Ｏｓａｋａ， Ｊａｐａｎ）を用いて特定した。その後、電気刺激の弱・中・強の３レベル、各３刺激を低温刺激と同じ呈示方法で呈示した。参加者は、ＣＯＶＡＳを用いて痛みレベルの主観評価を刺激呈示と並行して行った。
【0107】
（ＥＥＧ分析）
（振幅の特徴量抽出）
低温、および電気刺激条件における連続的なＥＥＧデータは、眼球運動ノイズ（ＥＯＧ）を減衰するため、以下の回帰フィルターを適用した：
［数１］
生ＥＥＧ＝β×ＥＯＧ＋Ｃ
推定ＥＥＧ＝生ＥＥＧ－β×ＥＯＧ
β：回帰係数
Ｃ：切片
推定ＥＥＧ：推定されたＥＥＧ
Ｆｐ１は、左眼に最も近接しており、眼球運動に強く影響を受けるため、Ｆｐ１データをＥＯＧデータとして使用した。ＥＯＧの補正後、刺激呈示前５秒から刺激呈示後１５秒のエポック波形を、各レベルの刺激ごとに切り出した。ベースライン補正を刺激呈示前の平均電位を用いて行った後、±１００μＶでアーチファクト除去を行った。電位を絶対化した後、最大振幅で標準化を行い、刺激呈示後１５秒間の平均振幅をレベルごとに求め、振幅の特徴量とした（４個：Ｆｚ，Ｃｚ，Ｃ３，Ｃ４）。快不快の特徴量として、低温刺激、電気刺激ともに、確実に「痛い」という主観報告をされた強度レベル３だけを用いた。
（振幅の特徴量抽出）
低温、および電気刺激条件における連続的なＥＥＧデータは、眼球運動ノイズ（ＥＯＧ）を減衰するため、以下の回帰フィルターを適用した：
［数１］
生ＥＥＧ＝β×ＥＯＧ＋Ｃ
推定ＥＥＧ＝生ＥＥＧ－β×ＥＯＧ
β：回帰係数
Ｃ：切片
推定ＥＥＧ：推定されたＥＥＧ
Ｆｐ１は、左眼に最も近接しており、眼球運動に強く影響を受けるため、Ｆｐ１データをＥＯＧデータとして使用した。ＥＯＧの補正後、刺激呈示前５秒から刺激呈示後１５秒のエポック波形を、各レベルの刺激ごとに切り出した。ベースライン補正を刺激呈示前の平均電位を用いて行った後、±１００μＶでアーチファクト除去を行った。電位を絶対化した後、最大振幅で標準化を行い、刺激呈示後１５秒間の平均振幅をレベルごとに求め、振幅の特徴量とした（４個：Ｆｚ，Ｃｚ，Ｃ３，Ｃ４）。快不快の特徴量として、低温刺激、電気刺激ともに、確実に「痛い」という主観報告をされた強度レベル３だけを用いた。
【0108】
（周波数パワーの特徴量抽出）
周波数解析では、ＥＯＧ補正処理を全体ＥＥＧデータに行った後、各レベルの刺激ごとに、刺激開始から刺激呈示後１５秒までの脳波データを切り出した。フーリエ変換を行った後、周波数パワー（実数部のｌｏｇ１０変換データ）を算出した。δ（１－３Ｈｚ）、θ（４－７Ｈｚ）、α（８－１３Ｈｚ）、β（１４－３０Ｈｚ）、γ（３１－１００Ｈｚ）ごとに、パワーの平均値をレベルごとに算出し、全レベルのデータを用いて個人ごとに最大値で標準化し、周波数特徴量とした（２０個：４電極×５帯域）。快不快の特徴量として、低温と電気刺激ともに、確実に「痛い」という主観報告がなされたレベル３だけを用いた。
周波数解析では、ＥＯＧ補正処理を全体ＥＥＧデータに行った後、各レベルの刺激ごとに、刺激開始から刺激呈示後１５秒までの脳波データを切り出した。フーリエ変換を行った後、周波数パワー（実数部のｌｏｇ１０変換データ）を算出した。δ（１－３Ｈｚ）、θ（４－７Ｈｚ）、α（８－１３Ｈｚ）、β（１４－３０Ｈｚ）、γ（３１－１００Ｈｚ）ごとに、パワーの平均値をレベルごとに算出し、全レベルのデータを用いて個人ごとに最大値で標準化し、周波数特徴量とした（２０個：４電極×５帯域）。快不快の特徴量として、低温と電気刺激ともに、確実に「痛い」という主観報告がなされたレベル３だけを用いた。
【0109】
（快不快判別分析）
参加者４１名の不快刺激（低温刺激）と快刺激（電気刺激）の判別を、Ｓｕｐｐｏｒｔ Ｖｅｃｔｏｒ Ｍａｃｈｉｎｅ Ｒｅｃｕｒｓｉｖｅ Ｆｅａｔｕｒｅ Ｅｌｉｍｉｎａｔｉｏｎ（ＳＶＭ－ＲＦＥ）を用いて行った（Ｇｕｙｏｎ Ｉ， Ｗｅｓｔｏｎ Ｊ， Ｂａｒｎｈｉｌｌ Ｓ， ＆ Ｖａｐｎｉｋ Ｖ． Ｇｅｎｅ ｓｅｌｅｃｔｉｏｎ ｆｏｒ ｃａｎｃｅｒ ｃｌａｓｓｉｆｉｃａｔｉｏｎ ｕｓｉｎｇ Ｓｕｐｐｏｒｔ Ｖｅｃｔｏｒ Ｍａｃｈｉｎｅ． Ｍａｃｈｉｎｅ Ｌｅａｒｎｉｎｇ ４６， ３８９－４２２， ２００２）（図２）。図２に示すように、特徴量は「２４個」（個人内で標準化済み）利用し、不快レベルは、「２レベル（不快・不快でない）」に設定した。特徴のランキングとして、２４個の脳波特徴量のランキングを行った。判別アルゴリズムの決定は、サポートベクターマシーン（ＳＶＭ）を用いてデータを「１個抜き交差検証」を行い、最も判別精度が高い特徴の組み合わせを見つけることで行った。ランダム化テスト（１０００回）は、最も判別精度が高い特徴量数を用いて、不快度のラベルをランダム化し、偶然レベルの判別精度を算出することで行った（上位９５％以上を占める場合、有意な判別精度とした）。データ解析には統計ソフトウェアパッケージのＲ，ならびにＳＶＭ－ＲＦＥ のＲ－ｃｏｄｅ（ｈｔｔｐ：／／ｇｉｔｈｕｂ．ｃｏｍ／ｊｏｈｎｃｏｌｂｙ／ＳＶＭ－ＲＦＥ）を用いた。図３に示したとおり、ＳＶＭ－ＲＦＥのおおまかな流れは、１）訓練データを用いて判別器を訓練する、２）特徴量をランキングする、３）貢献度が最も低い特徴量を除く、というプロセスを、最後の特徴量が無くなるまで繰り返す過程を含み、詳細な流れは、図３Ｂのとおりであり、具体的には、以下に示すとおりである。
用語
ａ）訓練データ：Ｓａｍｐｌｅ＝［ｘ１，ｘ２，．．．，ｘｋ，．．．，ｘ］
ｂ）分類ラベル：Ｃｌａｓｓ＝［１，０，１，０，．．．，ｙｋ，．．．，ｙ］
ｃ）ＲＦＥにより生き残った特徴量の集合 ：Ｓｆｅａｔｕｒｅ＝［１，２，．．．，ｎ］
ｄ）特徴量のランキング：Ｒａｎｋ＝［］
手順
・訓練データを生き残った特徴量に限定する：Ｘ＝Ｓａｍｐｌｅ（：，Ｓｆｅａｔｕｒｅ）
・判別器を訓練する（ＳＶＭ）：Ｃｌａｓｓｉｆｉｅｒ＝ＳＶＭ（Ｘ，Ｃｌａｓｓ）
・特徴量の重み付け係数を計算する：Ｗｅｉｇｈｔ＝Σαｋｙｋｘｋ
・ランキング基準をすべての特徴量において計算する：Ｃｒｉｔｅｒｉｏｎｉ＝（Ｗｅｉｇｈｔｉ）２，ｆｏｒ ａｌｌｉ
・もっとも小さいランキング基準の特徴量を見つける：Ｆ＝ａｒｇｍａｉｎ（Ｃｒｉｔｅｒｉｏｎ）
・生き残った特徴量のランキングリストを刷新する：Ｒａｎｋ＝［Ｓｆｅａｔｕｒｅ（Ｆ），Ｒａｎｋ］
・もっとも小さいランキング基準の特徴量を除く：Ｓｆｅａｔｕｒｅ＝Ｓｆｅａｔｕｒｅ（１：ｆ－１，ｆ＋１：ｌｅｎｇｔｈ（Ｓｆｅａｔｕｒｅ））
・特徴量ランキングリストを出力 ：Ｆｅａｔｕｒｅ ｒａｎｋｅｄ ｌｉｓｔ ｒ。
参加者４１名の不快刺激（低温刺激）と快刺激（電気刺激）の判別を、Ｓｕｐｐｏｒｔ Ｖｅｃｔｏｒ Ｍａｃｈｉｎｅ Ｒｅｃｕｒｓｉｖｅ Ｆｅａｔｕｒｅ Ｅｌｉｍｉｎａｔｉｏｎ（ＳＶＭ－ＲＦＥ）を用いて行った（Ｇｕｙｏｎ Ｉ， Ｗｅｓｔｏｎ Ｊ， Ｂａｒｎｈｉｌｌ Ｓ， ＆ Ｖａｐｎｉｋ Ｖ． Ｇｅｎｅ ｓｅｌｅｃｔｉｏｎ ｆｏｒ ｃａｎｃｅｒ ｃｌａｓｓｉｆｉｃａｔｉｏｎ ｕｓｉｎｇ Ｓｕｐｐｏｒｔ Ｖｅｃｔｏｒ Ｍａｃｈｉｎｅ． Ｍａｃｈｉｎｅ Ｌｅａｒｎｉｎｇ ４６， ３８９－４２２， ２００２）（図２）。図２に示すように、特徴量は「２４個」（個人内で標準化済み）利用し、不快レベルは、「２レベル（不快・不快でない）」に設定した。特徴のランキングとして、２４個の脳波特徴量のランキングを行った。判別アルゴリズムの決定は、サポートベクターマシーン（ＳＶＭ）を用いてデータを「１個抜き交差検証」を行い、最も判別精度が高い特徴の組み合わせを見つけることで行った。ランダム化テスト（１０００回）は、最も判別精度が高い特徴量数を用いて、不快度のラベルをランダム化し、偶然レベルの判別精度を算出することで行った（上位９５％以上を占める場合、有意な判別精度とした）。データ解析には統計ソフトウェアパッケージのＲ，ならびにＳＶＭ－ＲＦＥ のＲ－ｃｏｄｅ（ｈｔｔｐ：／／ｇｉｔｈｕｂ．ｃｏｍ／ｊｏｈｎｃｏｌｂｙ／ＳＶＭ－ＲＦＥ）を用いた。図３に示したとおり、ＳＶＭ－ＲＦＥのおおまかな流れは、１）訓練データを用いて判別器を訓練する、２）特徴量をランキングする、３）貢献度が最も低い特徴量を除く、というプロセスを、最後の特徴量が無くなるまで繰り返す過程を含み、詳細な流れは、図３Ｂのとおりであり、具体的には、以下に示すとおりである。
用語
ａ）訓練データ：Ｓａｍｐｌｅ＝［ｘ１，ｘ２，．．．，ｘｋ，．．．，ｘ］
ｂ）分類ラベル：Ｃｌａｓｓ＝［１，０，１，０，．．．，ｙｋ，．．．，ｙ］
ｃ）ＲＦＥにより生き残った特徴量の集合 ：Ｓｆｅａｔｕｒｅ＝［１，２，．．．，ｎ］
ｄ）特徴量のランキング：Ｒａｎｋ＝［］
手順
・訓練データを生き残った特徴量に限定する：Ｘ＝Ｓａｍｐｌｅ（：，Ｓｆｅａｔｕｒｅ）
・判別器を訓練する（ＳＶＭ）：Ｃｌａｓｓｉｆｉｅｒ＝ＳＶＭ（Ｘ，Ｃｌａｓｓ）
・特徴量の重み付け係数を計算する：Ｗｅｉｇｈｔ＝Σαｋｙｋｘｋ
・ランキング基準をすべての特徴量において計算する：Ｃｒｉｔｅｒｉｏｎｉ＝（Ｗｅｉｇｈｔｉ）２，ｆｏｒ ａｌｌｉ
・もっとも小さいランキング基準の特徴量を見つける：Ｆ＝ａｒｇｍａｉｎ（Ｃｒｉｔｅｒｉｏｎ）
・生き残った特徴量のランキングリストを刷新する：Ｒａｎｋ＝［Ｓｆｅａｔｕｒｅ（Ｆ），Ｒａｎｋ］
・もっとも小さいランキング基準の特徴量を除く：Ｓｆｅａｔｕｒｅ＝Ｓｆｅａｔｕｒｅ（１：ｆ－１，ｆ＋１：ｌｅｎｇｔｈ（Ｓｆｅａｔｕｒｅ））
・特徴量ランキングリストを出力 ：Ｆｅａｔｕｒｅ ｒａｎｋｅｄ ｌｉｓｔ ｒ。
【0110】
本実施例においては、快４１ラベル、不快４１ラベルの８２サンプルと２４特徴量を用いて判別分析を行った。放射基底（ガウス）関数をｋｅｒｎｅｌとして用い、データは正規化した。
［数２］
放射基底関数：
Ｇ（ｘ１,ｘ２）＝ｅｘｐ（－|ｘ１－ｘ２|２）
Ｇ＝ガウス関数
ｘ＝データポイント
ｅｘｐ＝指数関数。
［数２］
放射基底関数：
Ｇ（ｘ１,ｘ２）＝ｅｘｐ（－|ｘ１－ｘ２|２）
Ｇ＝ガウス関数
ｘ＝データポイント
ｅｘｐ＝指数関数。
【0111】
（結果と考察）
図４に２４特徴量のランキングリストと上位特徴量から順番に1個ずつ増やしていった
場合の８２サンプルの快不快判別精度を示す。頭皮上中央部付近の電極における振幅と周波数の特徴量が上位１０特徴量のうち８割を占めており、上位から、Ｃｚ（振幅）、Ｃ３（α）、Ｃｚ（β）、Ｆｚ（δ）、Ｃｚ（γ）、Ｃ４（β）、Ｃ３（δ）、Ｃｚ（α）、Ｆｚ（α）、Ｃ３（θ）であった。判別精度がもっとも高い判別器は、上位５特徴のＣｚ（振幅），Ｃ３（α），Ｃｚ（β），Ｆｚ（δ），Ｃｚ（γ）を用いたもので、「５６．０９８％」の判別精度であった。この判別精度が、偶然レベルではなく、統計的に有意であることを確認するために、快不快の判別ラベルを１０００回ランダム化し、上位５特徴を使って訓練を行った。図５に、判別精度のランダム化分布を示す。「５６．０９８％」の実際の判別精度は、ランダム分布の上位９５％に含まれ、統計的に有意であることが示された。
図４に２４特徴量のランキングリストと上位特徴量から順番に1個ずつ増やしていった
場合の８２サンプルの快不快判別精度を示す。頭皮上中央部付近の電極における振幅と周波数の特徴量が上位１０特徴量のうち８割を占めており、上位から、Ｃｚ（振幅）、Ｃ３（α）、Ｃｚ（β）、Ｆｚ（δ）、Ｃｚ（γ）、Ｃ４（β）、Ｃ３（δ）、Ｃｚ（α）、Ｆｚ（α）、Ｃ３（θ）であった。判別精度がもっとも高い判別器は、上位５特徴のＣｚ（振幅），Ｃ３（α），Ｃｚ（β），Ｆｚ（δ），Ｃｚ（γ）を用いたもので、「５６．０９８％」の判別精度であった。この判別精度が、偶然レベルではなく、統計的に有意であることを確認するために、快不快の判別ラベルを１０００回ランダム化し、上位５特徴を使って訓練を行った。図５に、判別精度のランダム化分布を示す。「５６．０９８％」の実際の判別精度は、ランダム分布の上位９５％に含まれ、統計的に有意であることが示された。
【0112】
以上の結果から、脳波特徴量を用いて判別器を作成し、快不快の異なる刺激タイプを判別することには意義があり、さらに、貢献度が高い特徴量を見つけて投入すれば、判別器の精度がさらに向上する余地があることを示唆している。
【0113】
（実施例２：異なる不快度を持つ、同じ高温刺激の判別）
本実施例では、高温疼痛刺激での快不快を分析した。特に、刺激呈示の文脈を変えることで、同じ高温刺激の快不快評価を変動させ、関連する脳波特徴量とシグモイド関数を用いて、快不快の判別を行った。これにより、従来の単なるデータの関連付けでは難しかった痛みでの、刺激強度が統制された快不快の判定が可能であることを実証する。
本実施例では、高温疼痛刺激での快不快を分析した。特に、刺激呈示の文脈を変えることで、同じ高温刺激の快不快評価を変動させ、関連する脳波特徴量とシグモイド関数を用いて、快不快の判別を行った。これにより、従来の単なるデータの関連付けでは難しかった痛みでの、刺激強度が統制された快不快の判定が可能であることを実証する。
【0114】
（方法）
（参加者）
２０代から３０代の２５名の健常な成人被験者が、高温刺激を用いた疼痛実験に、異なる日に２回参加した。参加者らは、実験前にインフォームドコンセントに同意した。全参加者は、神経性および／もしくは精神性の病気、または臨床薬物治療条件下における急性および／もしくは慢性疼痛を経験していないことを自己報告した。本実施例を、大阪大学医学部付属病院の倫理委員会ならびにヘルシンキ宣言による承認下で実施した。
（参加者）
２０代から３０代の２５名の健常な成人被験者が、高温刺激を用いた疼痛実験に、異なる日に２回参加した。参加者らは、実験前にインフォームドコンセントに同意した。全参加者は、神経性および／もしくは精神性の病気、または臨床薬物治療条件下における急性および／もしくは慢性疼痛を経験していないことを自己報告した。本実施例を、大阪大学医学部付属病院の倫理委員会ならびにヘルシンキ宣言による承認下で実施した。
【0115】
（手順）
温度刺激システム（Ｐａｔｈｗａｙ：Ｍｅｄｏｃ Ｃｏ．，Ｌｔｄ．，ＲａｍａｔＹｉｓｈａｉ， Ｉｓｒａｅｌ）を使用して、高温刺激を参加者の右前腕に与えた。刺激提示条件は、２つ用意した。条件１は、疼痛弱文脈条件であり、３６℃（ベース温度）、４０℃、３８℃を含んでいた。３６℃のベース区間が最初に１分あり、４０℃、３８℃の刺激区間がランダムに３回ずつ発生した。４０℃と３８℃の刺激区間は、１５秒の刺激（立ち上がり、立下り５秒）を５回連続呈示した。疼痛強文脈条件は、３６℃（ベース温度）、４０℃、４８℃を含んでいた。３８℃刺激ブロックを４８℃ブロックに置き換えることで、刺激提示文脈を変化させた。疼痛弱文脈条件と同様に、３６℃のベース区間が最初に１分あり、４０℃、４８℃の刺激区間がランダムに３回ずつ発生した。４０℃と４８℃の刺激区間は、１５秒の刺激（立ち上がり、立下り５秒）を５回連続呈示した。参加者は、コンピュータ化された可視化アナログスケール（ＣＯＶＡＳ）上で、０から１００の範囲（０：「無痛」；１００：「耐えられない疼痛」）で疼痛強度を連続的に評価した。ＣＯＶＡＳデータを刺激強度の変更と同時に記録した。
温度刺激システム（Ｐａｔｈｗａｙ：Ｍｅｄｏｃ Ｃｏ．，Ｌｔｄ．，ＲａｍａｔＹｉｓｈａｉ， Ｉｓｒａｅｌ）を使用して、高温刺激を参加者の右前腕に与えた。刺激提示条件は、２つ用意した。条件１は、疼痛弱文脈条件であり、３６℃（ベース温度）、４０℃、３８℃を含んでいた。３６℃のベース区間が最初に１分あり、４０℃、３８℃の刺激区間がランダムに３回ずつ発生した。４０℃と３８℃の刺激区間は、１５秒の刺激（立ち上がり、立下り５秒）を５回連続呈示した。疼痛強文脈条件は、３６℃（ベース温度）、４０℃、４８℃を含んでいた。３８℃刺激ブロックを４８℃ブロックに置き換えることで、刺激提示文脈を変化させた。疼痛弱文脈条件と同様に、３６℃のベース区間が最初に１分あり、４０℃、４８℃の刺激区間がランダムに３回ずつ発生した。４０℃と４８℃の刺激区間は、１５秒の刺激（立ち上がり、立下り５秒）を５回連続呈示した。参加者は、コンピュータ化された可視化アナログスケール（ＣＯＶＡＳ）上で、０から１００の範囲（０：「無痛」；１００：「耐えられない疼痛」）で疼痛強度を連続的に評価した。ＣＯＶＡＳデータを刺激強度の変更と同時に記録した。
【0116】
（ＥＥＧデータ記録）
市販されているＢｉｏ－Ａｍｐｌｉｆｉｅｒ（ＥＥＧ１２００; Ｎｉｈｏｎ Ｋｏｄｅｎ）を使用して、７つの頭皮Ａｇ／ＡｇＣｌ電極（Ｆｐ１、Ｆｐ２、Ｆ３、Ｆ４、Ｃ３、Ｃ４、Ｐｚ）からＥＥＧを記録した。最前部の電極であるＦｐ１とＦｐ２を眼球電位（ＥＯＧ）補正のために使用した。脳波導出用の基準電極を両方の耳たぶに装着し、アース電極を額の中央部に置いた。サンプリングレートは１，０００Ｈｚで、０．３～１２０Ｈｚの範囲の帯域で増幅した。全電極のインピーダンスは１５ｋΩ未満であった。
市販されているＢｉｏ－Ａｍｐｌｉｆｉｅｒ（ＥＥＧ１２００; Ｎｉｈｏｎ Ｋｏｄｅｎ）を使用して、７つの頭皮Ａｇ／ＡｇＣｌ電極（Ｆｐ１、Ｆｐ２、Ｆ３、Ｆ４、Ｃ３、Ｃ４、Ｐｚ）からＥＥＧを記録した。最前部の電極であるＦｐ１とＦｐ２を眼球電位（ＥＯＧ）補正のために使用した。脳波導出用の基準電極を両方の耳たぶに装着し、アース電極を額の中央部に置いた。サンプリングレートは１，０００Ｈｚで、０．３～１２０Ｈｚの範囲の帯域で増幅した。全電極のインピーダンスは１５ｋΩ未満であった。
【0117】
（ＥＥＧ分析）
連続的なＥＥＧデータを、ＥＯＧノイズを減衰するため、以下の回帰フィルターを生のＥＥＧデータに適用した：
［数３］
生ＥＥＧ＝β×（Ｆｐ１＋Ｆｐ２）＋Ｃ
推定ＥＥＧ＝生ＥＥＧ－β×（Ｆｐ１＋Ｆｐ２）
β：回帰係数
Ｃ：切片
推定ＥＥＧ：推定されたＥＥＧ
Ｆｐ１とＦｐ２データは、垂直眼球運動、もしくは瞬目活動の電位を増幅するために加算して用いた。ＶＥＯＧ補正後、ハムノイズ（６０Ｈｚ）を省くために、ノッチフィルターを全ＥＥＧデータにかけた。刺激開始から１５秒間の脳波データを３６℃（疼痛弱、強文脈：１エポック）、４０℃（疼痛弱、強文脈：５エポック×２ブロック）、３８℃（疼痛弱文脈：５エポック×２ブロック）、４８℃（疼痛強文脈：５エポック×２ブロック）の各刺激に対して抽出した。振幅を絶対値化した後、１００μＶを超える電位をアーチファクト除去し、最大振幅を用いて標準化した。刺激ごとに時間方向で振幅（絶対値化、ならびに標準化済）の平均値を求め、さらに、各温度刺激で全体平均振幅を求めた。
連続的なＥＥＧデータを、ＥＯＧノイズを減衰するため、以下の回帰フィルターを生のＥＥＧデータに適用した：
［数３］
生ＥＥＧ＝β×（Ｆｐ１＋Ｆｐ２）＋Ｃ
推定ＥＥＧ＝生ＥＥＧ－β×（Ｆｐ１＋Ｆｐ２）
β：回帰係数
Ｃ：切片
推定ＥＥＧ：推定されたＥＥＧ
Ｆｐ１とＦｐ２データは、垂直眼球運動、もしくは瞬目活動の電位を増幅するために加算して用いた。ＶＥＯＧ補正後、ハムノイズ（６０Ｈｚ）を省くために、ノッチフィルターを全ＥＥＧデータにかけた。刺激開始から１５秒間の脳波データを３６℃（疼痛弱、強文脈：１エポック）、４０℃（疼痛弱、強文脈：５エポック×２ブロック）、３８℃（疼痛弱文脈：５エポック×２ブロック）、４８℃（疼痛強文脈：５エポック×２ブロック）の各刺激に対して抽出した。振幅を絶対値化した後、１００μＶを超える電位をアーチファクト除去し、最大振幅を用いて標準化した。刺激ごとに時間方向で振幅（絶対値化、ならびに標準化済）の平均値を求め、さらに、各温度刺激で全体平均振幅を求めた。
【0118】
（主観評価の分析）
脳波データと同様に、３６℃（疼痛弱、強文脈：１エポック）、４０℃（疼痛弱、強文脈：５エポック×２ブロック）、３８℃（疼痛弱文脈：５エポック×２ブロック）、４８℃（疼痛強文脈：５エポック×２ブロック）の各刺激に対して抽出した。抽出した時間幅は、主観評価は痛み刺激提示後遅延して変化するので、刺激提示後１０秒から次の刺激提示までとし、その区間の最大値を評価得点とし、各温度条件で全体平均得点を算出した。
脳波データと同様に、３６℃（疼痛弱、強文脈：１エポック）、４０℃（疼痛弱、強文脈：５エポック×２ブロック）、３８℃（疼痛弱文脈：５エポック×２ブロック）、４８℃（疼痛強文脈：５エポック×２ブロック）の各刺激に対して抽出した。抽出した時間幅は、主観評価は痛み刺激提示後遅延して変化するので、刺激提示後１０秒から次の刺激提示までとし、その区間の最大値を評価得点とし、各温度条件で全体平均得点を算出した。
【0119】
（統計解析）
統計検定を行う前に、疼痛弱文脈と疼痛強文脈の比較ベースラインをそろえる為に、３６℃の脳波振幅、ならびに主観評価スコアをもとに他の温度条件の数値をベースライン補正した。主観評価の文脈効果を調べる為に、温度強度が同じ、すなわち４０℃条件の評価スコアを弱と強文脈間でｔ検定を用いて比較した。同様に、脳波振幅はＦ３とＦ４の平均値をもとめ、４０℃条件の平均振幅を弱と強文脈間でｔ検定を用いて比較した。
統計検定を行う前に、疼痛弱文脈と疼痛強文脈の比較ベースラインをそろえる為に、３６℃の脳波振幅、ならびに主観評価スコアをもとに他の温度条件の数値をベースライン補正した。主観評価の文脈効果を調べる為に、温度強度が同じ、すなわち４０℃条件の評価スコアを弱と強文脈間でｔ検定を用いて比較した。同様に、脳波振幅はＦ３とＦ４の平均値をもとめ、４０℃条件の平均振幅を弱と強文脈間でｔ検定を用いて比較した。
【0120】
（判別モデル）
脳波振幅による不快度判別器をシグモイド関数フィッティングにより作成したあと、変調点の数値を算出し、不快度判定値とした。その判定値に基づき、異なる文脈における４０℃刺激の判別を行った。
脳波振幅による不快度判別器をシグモイド関数フィッティングにより作成したあと、変調点の数値を算出し、不快度判定値とした。その判定値に基づき、異なる文脈における４０℃刺激の判別を行った。
【0121】
（結果）
図６に主観評価の結果を示す。同じ４０℃刺激でも、３８℃の疼痛弱文脈と４８℃の疼痛強文脈では、３６℃を共通のベースラインとした場合、主観評価の程度が有意に異なっていた。すなわち、４０℃刺激は、３８℃の疼痛弱刺激と同じ提示文脈で呈示された場合、４８℃の疼痛強文脈に比べてより不快度が高いと感じられることが観察された。
図６に主観評価の結果を示す。同じ４０℃刺激でも、３８℃の疼痛弱文脈と４８℃の疼痛強文脈では、３６℃を共通のベースラインとした場合、主観評価の程度が有意に異なっていた。すなわち、４０℃刺激は、３８℃の疼痛弱刺激と同じ提示文脈で呈示された場合、４８℃の疼痛強文脈に比べてより不快度が高いと感じられることが観察された。
【0122】
また、前頭電極の４０℃刺激に係る脳波振幅も、弱と強文脈では有意差があり、強文脈の方が高い振幅を示していた（図７）。４８℃条件の高い脳活動性が、文脈効果により４０℃条件における高い活動性を惹起し、その高い活動性が評価基準に影響を与え、主観評価を相対的に下げたことが示唆される。
【0123】
図８に、シグモイド関数を用いた、不快刺激２５サンプルと快刺激２５サンプルの判別分析の結果を示す。フィッティングにより、以下のような判別器（シグモイド関数）が得られた。
［数４］
ｙ＝０．９２＋０．２７７４／（１＋１０（１９．６７－ｘ）×８８．１２）
変調点の閾値は「１．０５５５」となり、特徴量が閾値より高い場合が「不快度が低い」と判別され、閾値より低い場合が「不快度が高い」と判別される。この判別基準に基づき、全５０サンプルを判別したところ、全体のエラー率は「３６％」であり、判別精度は６０～７０％の範囲に収まる結果となっていた。シンプルな二値判別関数により、物理量が同じ刺激に関する不快度の違いが、偶然レベルを超えた判別結果を示したことは本発明の特筆すべき結果であり、同じ物理量の刺激を異なる呈示文脈で提示する、文脈依存型リファレンス検査手法により、より厳密な不快疼痛評価を実現可能となる。このように、不快度評価は、脳波振幅からは読み取れない実際の被験者の不快度を精確に判断することができることが示された。
［数４］
ｙ＝０．９２＋０．２７７４／（１＋１０（１９．６７－ｘ）×８８．１２）
変調点の閾値は「１．０５５５」となり、特徴量が閾値より高い場合が「不快度が低い」と判別され、閾値より低い場合が「不快度が高い」と判別される。この判別基準に基づき、全５０サンプルを判別したところ、全体のエラー率は「３６％」であり、判別精度は６０～７０％の範囲に収まる結果となっていた。シンプルな二値判別関数により、物理量が同じ刺激に関する不快度の違いが、偶然レベルを超えた判別結果を示したことは本発明の特筆すべき結果であり、同じ物理量の刺激を異なる呈示文脈で提示する、文脈依存型リファレンス検査手法により、より厳密な不快疼痛評価を実現可能となる。このように、不快度評価は、脳波振幅からは読み取れない実際の被験者の不快度を精確に判断することができることが示された。
【0124】
図９に、本実施例の応用例を示す。臨床現場において、持続する不快疼痛をモニタリングし鎮静化を行う場合、図に示したように、不快度強状態から不快度弱状態に漸次移行することが予測される。本発明における不快度判別器は、この不快度の推移を判別／推定するために、以下のようなプロセスを用いて行う。
１．侵襲性が低い刺激（電気や熱刺激）を一定の強度、間隔で患者に背景的に提示する。２．一定の時間幅（例えば、１５秒）の平均振幅（絶対値化、最初の刺激の振幅などを用いて標準化を行う）を時間変化に応じて算出し続け、トレンド表示する。３．複数の平均振幅を時間的な前後群（両群１０個など）に分け、t検定などで比較し、有意な変化パターン時点を見つける。
１．侵襲性が低い刺激（電気や熱刺激）を一定の強度、間隔で患者に背景的に提示する。２．一定の時間幅（例えば、１５秒）の平均振幅（絶対値化、最初の刺激の振幅などを用いて標準化を行う）を時間変化に応じて算出し続け、トレンド表示する。３．複数の平均振幅を時間的な前後群（両群１０個など）に分け、t検定などで比較し、有意な変化パターン時点を見つける。
【0125】
このような未知の不快疼痛レベルの推移の客観的評価方法は、刺激提示文脈を利用した手法であることから、「文脈依存型リファレンス検査手法」と呼びうる。
【0126】
（実施例３：異なる文脈における認知課題試行中の潜在的ストレスの判別）
本実施例では、心理ストレスパラダイムとして、認知的競合課題（ストループ課題）実験を異なる課題実施文脈を設定して行った。
本実施例では、心理ストレスパラダイムとして、認知的競合課題（ストループ課題）実験を異なる課題実施文脈を設定して行った。
【0127】
（方法）
（参加者）
２０代から３０代の２６名の健常な成人被験者が、ストループ課題の実験に、異なる日に２回参加した。参加者らは、実験前にインフォームドコンセントに同意した。全参加者は、神経性および／もしくは精神性の病気、または臨床薬物治療条件下における急性および／もしくは慢性疼痛を経験していないことを自己報告した。本実施例を、大阪大学医学部付属病院の倫理委員会ならびにヘルシンキ宣言による承認下で実施した。
（参加者）
２０代から３０代の２６名の健常な成人被験者が、ストループ課題の実験に、異なる日に２回参加した。参加者らは、実験前にインフォームドコンセントに同意した。全参加者は、神経性および／もしくは精神性の病気、または臨床薬物治療条件下における急性および／もしくは慢性疼痛を経験していないことを自己報告した。本実施例を、大阪大学医学部付属病院の倫理委員会ならびにヘルシンキ宣言による承認下で実施した。
【0128】
（方法）
図１０に、実施したストループ課題の概要を示す。本実施例においては、参加者がフォントの色を回答する「フォントカラー課題」を用いた。認知的負荷が低い一致条件では、文字の意味（例．「あか」）とフォントの色が同じであった。認知的負荷が高く、回答時に競合を生じる不一致条件では、文字の意味とフォントの色が異なっていた。課題は、全部で３ブロックからなり、各ブロックは５０試行（各条件２５試行）を含んでいた。被検者は出来る限り早く回答するように求められ、間違った回答をした場合は、ブザーが鳴るエラーフィードバックを与えられた。ストレスなし条件では、課題試行中に監視者がおらず、また、実験実施者も視界に入らず、単独で実験を行った。ストレスあり条件では、未知の男性監視者が左斜め横１メートル以内に座り、終始、無言で、マスクを着用して参加者を監視した。
図１０に、実施したストループ課題の概要を示す。本実施例においては、参加者がフォントの色を回答する「フォントカラー課題」を用いた。認知的負荷が低い一致条件では、文字の意味（例．「あか」）とフォントの色が同じであった。認知的負荷が高く、回答時に競合を生じる不一致条件では、文字の意味とフォントの色が異なっていた。課題は、全部で３ブロックからなり、各ブロックは５０試行（各条件２５試行）を含んでいた。被検者は出来る限り早く回答するように求められ、間違った回答をした場合は、ブザーが鳴るエラーフィードバックを与えられた。ストレスなし条件では、課題試行中に監視者がおらず、また、実験実施者も視界に入らず、単独で実験を行った。ストレスあり条件では、未知の男性監視者が左斜め横１メートル以内に座り、終始、無言で、マスクを着用して参加者を監視した。
【0129】
（ＥＥＧデータ記録）
市販されているＢｉｏ－Ａｍｐｌｉｆｉｅｒ（ＥＥＧ１２００; Ｎｉｈｏｎ Ｋｏｄｅｎ）を使用して、７つの頭皮Ａｇ／ＡｇＣｌ電極（Ｆｐ１，Ｆｐ２，Ｆ３，Ｆ４，Ｃ３，Ｃ４，Ｐｚ）からＥＥＧを記録した。最前部の電極であるＦｐ１とＦｐ２をＥＯＧ補正のために使用した。脳波導出用の基準電極を両方の耳たぶに装着し、アース電極を額の中央部に置いた。サンプリングレートは１，０００Ｈｚで、０．３～１２０Ｈｚの範囲の帯域で増幅した。全電極のインピーダンスは１５ｋΩ未満であった。
市販されているＢｉｏ－Ａｍｐｌｉｆｉｅｒ（ＥＥＧ１２００; Ｎｉｈｏｎ Ｋｏｄｅｎ）を使用して、７つの頭皮Ａｇ／ＡｇＣｌ電極（Ｆｐ１，Ｆｐ２，Ｆ３，Ｆ４，Ｃ３，Ｃ４，Ｐｚ）からＥＥＧを記録した。最前部の電極であるＦｐ１とＦｐ２をＥＯＧ補正のために使用した。脳波導出用の基準電極を両方の耳たぶに装着し、アース電極を額の中央部に置いた。サンプリングレートは１，０００Ｈｚで、０．３～１２０Ｈｚの範囲の帯域で増幅した。全電極のインピーダンスは１５ｋΩ未満であった。
【0130】
（ＥＥＧ分析）
連続的なＥＥＧデータを、ＥＯＧノイズを減衰するため、以下の回帰フィルターを生のＥＥＧデータに適用した：
［数５］
生ＥＥＧ＝β×（Ｆｐ１＋Ｆｐ２）＋Ｃ
推定ＥＥＧ＝生ＥＥＧ－β×（Ｆｐ１＋Ｆｐ２）
β：回帰係数
Ｃ：切片
推定ＥＥＧ：推定されたＥＥＧ
Ｆｐ１とＦｐ２データは、垂直眼球運動、もしくは瞬目活動の電位を増幅するために加算して用いた。ＶＥＯＧ補正後、低周波数と高周波数成分を除くために、０．３Ｈｚ～４０Ｈｚの帯域周波数フィルターをかけた。文字提示前２００ミリ秒から提示後８００ミリ秒までの脳波データエポックを条件ごとに切り出した（５０エポック×２条件×２文脈条件）。刺激呈示前の平均振幅を用いてベースライン補正を行い、±５０μＶ以上の振幅が見られるエポックはその後の解析から除去した。生き残ったエポックを個人ごとに平均したあと、２条件（一致、不一致）と全電極（Ｆ３，Ｆ４，Ｃ３，Ｃ４，Ｐｚ）を含む全データの最大振幅絶対値を用いて標準化した。
連続的なＥＥＧデータを、ＥＯＧノイズを減衰するため、以下の回帰フィルターを生のＥＥＧデータに適用した：
［数５］
生ＥＥＧ＝β×（Ｆｐ１＋Ｆｐ２）＋Ｃ
推定ＥＥＧ＝生ＥＥＧ－β×（Ｆｐ１＋Ｆｐ２）
β：回帰係数
Ｃ：切片
推定ＥＥＧ：推定されたＥＥＧ
Ｆｐ１とＦｐ２データは、垂直眼球運動、もしくは瞬目活動の電位を増幅するために加算して用いた。ＶＥＯＧ補正後、低周波数と高周波数成分を除くために、０．３Ｈｚ～４０Ｈｚの帯域周波数フィルターをかけた。文字提示前２００ミリ秒から提示後８００ミリ秒までの脳波データエポックを条件ごとに切り出した（５０エポック×２条件×２文脈条件）。刺激呈示前の平均振幅を用いてベースライン補正を行い、±５０μＶ以上の振幅が見られるエポックはその後の解析から除去した。生き残ったエポックを個人ごとに平均したあと、２条件（一致、不一致）と全電極（Ｆ３，Ｆ４，Ｃ３，Ｃ４，Ｐｚ）を含む全データの最大振幅絶対値を用いて標準化した。
【0131】
（判別分析）
心理ストレス文脈の不快度判別器をシグモイド関数フィッティングにより作成した。特徴量として、結果で示すように、文脈の違いが、陰性電位効果の持続性に顕著に表れているので、特定の時間範囲における陰性効果の占有率（陰性電位の総時間（もしくは陰性電位データポイント数）／全体時間（もしくは、総データポイント数）×１００）を考案して使用した。フィッティング関数における変調点の数値を算出し、不快度判定値とした。判定値に基づき、異なる文脈におけるストレス度合い、すなわち不快度の判別を行った。
心理ストレス文脈の不快度判別器をシグモイド関数フィッティングにより作成した。特徴量として、結果で示すように、文脈の違いが、陰性電位効果の持続性に顕著に表れているので、特定の時間範囲における陰性効果の占有率（陰性電位の総時間（もしくは陰性電位データポイント数）／全体時間（もしくは、総データポイント数）×１００）を考案して使用した。フィッティング関数における変調点の数値を算出し、不快度判定値とした。判定値に基づき、異なる文脈におけるストレス度合い、すなわち不快度の判別を行った。
【0132】
（結果）
（行動データ）
図１１に、心理ストレスなし文脈と心理ストレスあり文脈における、一致、不一致条件の回答に要した反応時間の比較結果を示す。どちらの文脈でも、一致条件に比べ不一致条件の方が回答に要した反応時間は長くなっており、不一致条件の認知的負荷が高いことが確認できた。一方で、文脈間の反応時間には有意差が見られず、顕在的なストレスの違いは反応時間には現れなかった。これは、参加者の課題中のストレスに関する主観報告が、文脈間でほとんど違いが見られなかったこととも一致する結果であった。
（行動データ）
図１１に、心理ストレスなし文脈と心理ストレスあり文脈における、一致、不一致条件の回答に要した反応時間の比較結果を示す。どちらの文脈でも、一致条件に比べ不一致条件の方が回答に要した反応時間は長くなっており、不一致条件の認知的負荷が高いことが確認できた。一方で、文脈間の反応時間には有意差が見られず、顕在的なストレスの違いは反応時間には現れなかった。これは、参加者の課題中のストレスに関する主観報告が、文脈間でほとんど違いが見られなかったこととも一致する結果であった。
【0133】
しかしながら、図１２の前頭部の事象関連電位の一致条件と不一致条件の比較では、文脈間に顕著な違いが見られた。課題試行中の監視がないストレスなし文脈では、不一致条件の波形が、一致条件の波形に比べ下方向、すなわち陽性電位方向にシフトしているのが観察された（刺激提示後３００ミリ秒あたりから）。一方、監視があるストレスあり条件では、そのような陽性電位方向へのシフト見られず、逆に陰性方向へのシフトが見られた。この陰性方向へのシフトは刺激提示後２００ミリ秒後から持続していた。
【0134】
この陰性電位効果の持続性を占有率という特徴量として表現し、ストレスなし２６サンプル、ストレスあり２６サンプルの判別分析を行った（図１３）。判別器としてシグモイド近似関数を用い、以下のような数式で表現された。
［数６］
Ｙ＝３０．７７９＋２５．９０８／（１＋１０（２７．００３－Ｘ）×０．１２３）
変調点による心理的不快度判定値は「４３．５」であった。特に、この不快度判定値は、今回判別するストレスが、回答時間や主観評価に現れないが、前頭部優位の抑制機能関連の脳活動変動を生み出していることから、潜在的な心理的不快度を判別する閾値であることは特筆に値する。判定値４３．５以上を「ストレスあり」、４３．５未満を「ストレスなし」として、５２サンプルを判別したところ、全体の判別精度は「６８％」であり、チャンスレベル（５０％）より２０％ほど高い数値を示した。
［数６］
Ｙ＝３０．７７９＋２５．９０８／（１＋１０（２７．００３－Ｘ）×０．１２３）
変調点による心理的不快度判定値は「４３．５」であった。特に、この不快度判定値は、今回判別するストレスが、回答時間や主観評価に現れないが、前頭部優位の抑制機能関連の脳活動変動を生み出していることから、潜在的な心理的不快度を判別する閾値であることは特筆に値する。判定値４３．５以上を「ストレスあり」、４３．５未満を「ストレスなし」として、５２サンプルを判別したところ、全体の判別精度は「６８％」であり、チャンスレベル（５０％）より２０％ほど高い数値を示した。
【0135】
図１４に、本実施例の応用例を示す。職場や教育環境などにおいて、どのような環境を設定すれば、気づきにくい潜在的なストレスを軽減したり、または、適度なストレスを惹起し、作業効率をあげたり、動機付けを向上させるかは重要な問題である。そこで、図に示したように、実際の、もしくはシュミレーションした複数の環境条件を用意し、本実施例で用いた認知課題などを行い、持続性陰性電位活動の出現パターンを振幅、周波数、または本発明における占有率を用いて、条件間比較やランキングする。個人、もしくは集団のランキングを算出することで、職場や教育現場での実際の環境設定を客観的に行うことができ、産業上、環境評価や商品評価などに応用展開が期待できる。
【0136】
（実施例４：個人モデルによる快・不快痛みの判別推定：サンプル増幅法による自己複製型特徴量の利用）
【0137】
（方法）
（参加者）
２０代から７０代の１７０名の健常な成人被験者が、高温刺激パラダイムの実験に参加した。参加者らは、実験前にインフォームドコンセントに同意した。全参加者は、神経性および／もしくは精神性の病気、または臨床薬物治療条件下における急性および／もしくは慢性疼痛を経験していないことを自己報告した。本実施例を、大阪大学医学部付属病院の倫理委員会ならびにヘルシンキ宣言による承認下で実施した。
（参加者）
２０代から７０代の１７０名の健常な成人被験者が、高温刺激パラダイムの実験に参加した。参加者らは、実験前にインフォームドコンセントに同意した。全参加者は、神経性および／もしくは精神性の病気、または臨床薬物治療条件下における急性および／もしくは慢性疼痛を経験していないことを自己報告した。本実施例を、大阪大学医学部付属病院の倫理委員会ならびにヘルシンキ宣言による承認下で実施した。
【0138】
（実験刺激、ならびに手続）
温度刺激システム(Ｐａｔｈｗａｙ：Ｍｅｄｏｃ Ｃｏ．，Ｌｔｄ．，ＲａｍａｔＹｉｓｈａｉ， Ｉｓｒａｅｌ)を使用して、高温刺激を参加者の右前腕に与えた。高温刺激は、６レベルの温度強度（４０℃～５０℃の２℃ずつ上昇）を包含していた。各温度レベルは３回の刺激からなり、５秒間持続するプラトーを有し、基準温度（３５℃）からの上昇および下降待機時間は、約１０秒であった。各レベルにおける３回の刺激の後、ブロック間間隔は１００秒間持続した。参加者は、コンピュータ化された可視化アナログスケール（ＣＯＶＡＳ）上で、０から１００の範囲（０＝「無痛」；１００＝「耐えられない疼痛」）で疼痛強度を連続的に評価した。ＣＯＶＡＳデータを刺激強度の変更と同時に記録した。
温度刺激システム(Ｐａｔｈｗａｙ：Ｍｅｄｏｃ Ｃｏ．，Ｌｔｄ．，ＲａｍａｔＹｉｓｈａｉ， Ｉｓｒａｅｌ)を使用して、高温刺激を参加者の右前腕に与えた。高温刺激は、６レベルの温度強度（４０℃～５０℃の２℃ずつ上昇）を包含していた。各温度レベルは３回の刺激からなり、５秒間持続するプラトーを有し、基準温度（３５℃）からの上昇および下降待機時間は、約１０秒であった。各レベルにおける３回の刺激の後、ブロック間間隔は１００秒間持続した。参加者は、コンピュータ化された可視化アナログスケール（ＣＯＶＡＳ）上で、０から１００の範囲（０＝「無痛」；１００＝「耐えられない疼痛」）で疼痛強度を連続的に評価した。ＣＯＶＡＳデータを刺激強度の変更と同時に記録した。
【0139】
（ＥＥＧデータ記録）
前頭部4電極から脳波を実験中連続的に脳波計測装置（Ｎｅｕｒｏｆａｘ、日本光電）で記録した。導出電極は左右耳朶とし、同側の頭皮上電極に連結した。アース電極は前額中央に設置した。サンプリング周波数は１０００Ｈｚ、０．３～１２０Ｈｚで増幅した。インピーダンスは１５ｋΩ以下を保持した。
前頭部4電極から脳波を実験中連続的に脳波計測装置（Ｎｅｕｒｏｆａｘ、日本光電）で記録した。導出電極は左右耳朶とし、同側の頭皮上電極に連結した。アース電極は前額中央に設置した。サンプリング周波数は１０００Ｈｚ、０．３～１２０Ｈｚで増幅した。インピーダンスは１５ｋΩ以下を保持した。
【0140】
（ＥＥＧ分析）
（振幅の特徴量抽出のための前駆処理）
分析手順を図１９に示す。具体的には以下のとおりである。
・特徴量抽出
1. 電極抽出：Fp1,Fp2, F3, F4の4電極に、バーチャルEOG電極（主成分分析）を加える。
2. 瞬目（EOG）除去：主成分分析によりEOG成分（第一成分）抽出→元データに回帰フィルターをかける。
3. 筋電位（EMG）減衰：30Hzのハイカットフィルターをかける。
4. レベル1（40℃）とレベル6（50℃）において、刺激提示後5秒～15秒までの脳波を切り出し、「絶対値化」する（2×3刺激＝6エポック）。
5. レスト区間（刺激提示前30秒間）を用いて「ｚ値化」する。
・特徴量の増幅・複製
1. 実測サンプル１：
5秒区間を1ポイントずつずらしながら「移動平均」をかける。
各レベルにつき、「30,003サンプル(10,001サンプル×3刺激)」を作成する。
2.実測サンプル２：
2.1 1秒区間をポイントのオーバーラップなしに、各レベルにつき「30サンプル」を作成する
2.2 各電極における実測サンプルの分布特性を用いて正規乱数、もしくはピアソンシステム乱数を発生し、10,000サンプルの「自己複製型特徴量」を増幅・作成する。
＊「サンプル増幅法」は、サンプルが少ない場合でも、その分布特性を用いて、大幅にサンプル数を増加させる技法である。
・判別モデル作成＆評価
1. 機械学習により、個人の判別モデルを作成する(ＬＡＳＳＯとベイズ最適化による係数の決定)。
2. 他のサンプルを用いて、モデルの汎化能力を確認する。
・判別モデル作成の実際
参加者１名（ＩＤ１８５）のＥＥＧデータを判別モデル作成用に用いた。なお、個人判別モデルは１７０名全員で作成可能であるが、本実施例ではたった1名で作成した判別モデルがどの程度、有効に機能し、他の被検者にも汎化できるかを示すために、個別モデルを一名だけで作成した。眼球運動ノイズ（ＥＯＧ）を最初に減衰した。４電極のデータを用いて主成分分析を行い、ＥＯＧの第一成分を抽出した後、ＥＯＧ成分だけを強調して残すために。０．５から３０Ｈｚの帯域通過フィルターをかけた。その後、回帰フィルターを用いてＥＯＧ成分を除去した。ＥＯＧ除去後、３０Ｈｚの高域カットフィルターをかけ、筋電ノイズを軽減した。最後に「ノイズ逆利用」するために、主成分分析で得たＥＯＧ成分のバーチャルチャンネルを追加し、５電極に増幅した。刺激提示後５秒から１５秒までのエポックを「痛みなし」のレベル１（４０℃）と「耐えられない、不快な痛み」のレベル６（５０℃）に対して、3刺激ずつ切り出し、振幅を絶対値化した。
（振幅の特徴量抽出のための前駆処理）
分析手順を図１９に示す。具体的には以下のとおりである。
・特徴量抽出
1. 電極抽出：Fp1,Fp2, F3, F4の4電極に、バーチャルEOG電極（主成分分析）を加える。
2. 瞬目（EOG）除去：主成分分析によりEOG成分（第一成分）抽出→元データに回帰フィルターをかける。
3. 筋電位（EMG）減衰：30Hzのハイカットフィルターをかける。
4. レベル1（40℃）とレベル6（50℃）において、刺激提示後5秒～15秒までの脳波を切り出し、「絶対値化」する（2×3刺激＝6エポック）。
5. レスト区間（刺激提示前30秒間）を用いて「ｚ値化」する。
・特徴量の増幅・複製
1. 実測サンプル１：
5秒区間を1ポイントずつずらしながら「移動平均」をかける。
各レベルにつき、「30,003サンプル(10,001サンプル×3刺激)」を作成する。
2.実測サンプル２：
2.1 1秒区間をポイントのオーバーラップなしに、各レベルにつき「30サンプル」を作成する
2.2 各電極における実測サンプルの分布特性を用いて正規乱数、もしくはピアソンシステム乱数を発生し、10,000サンプルの「自己複製型特徴量」を増幅・作成する。
＊「サンプル増幅法」は、サンプルが少ない場合でも、その分布特性を用いて、大幅にサンプル数を増加させる技法である。
・判別モデル作成＆評価
1. 機械学習により、個人の判別モデルを作成する(ＬＡＳＳＯとベイズ最適化による係数の決定)。
2. 他のサンプルを用いて、モデルの汎化能力を確認する。
・判別モデル作成の実際
参加者１名（ＩＤ１８５）のＥＥＧデータを判別モデル作成用に用いた。なお、個人判別モデルは１７０名全員で作成可能であるが、本実施例ではたった1名で作成した判別モデルがどの程度、有効に機能し、他の被検者にも汎化できるかを示すために、個別モデルを一名だけで作成した。眼球運動ノイズ（ＥＯＧ）を最初に減衰した。４電極のデータを用いて主成分分析を行い、ＥＯＧの第一成分を抽出した後、ＥＯＧ成分だけを強調して残すために。０．５から３０Ｈｚの帯域通過フィルターをかけた。その後、回帰フィルターを用いてＥＯＧ成分を除去した。ＥＯＧ除去後、３０Ｈｚの高域カットフィルターをかけ、筋電ノイズを軽減した。最後に「ノイズ逆利用」するために、主成分分析で得たＥＯＧ成分のバーチャルチャンネルを追加し、５電極に増幅した。刺激提示後５秒から１５秒までのエポックを「痛みなし」のレベル１（４０℃）と「耐えられない、不快な痛み」のレベル６（５０℃）に対して、3刺激ずつ切り出し、振幅を絶対値化した。
【0141】
（自己複製型特徴量のサンプル増幅）
上記のノイズ処理をされたエポックデータから、図１９の手順に示すように、２通りの方法でサンプル増幅用の実測値サンプルを得た。第一は、各エポックに対し、５秒の区間単位を用いて、１ポイントずつずらす移動平均をかけ、各刺激に対し１０００１サンプルを得た（レベル１と６に対し、各３０００３サンプル）。第二は、本当の意味での、小さい実測サンプルから大きな複製型サンプルを得る方法であり、同じ１０秒の刺激区間を1秒ごとに重複しないように切り出し、１０サンプルの平均値を得た。これらのサンプルに対し、２種類のサンプル増幅法を行った。1つ目が、レベル１と６の平均値と標準偏差を用いて「正規乱数」を発生し（本実施例では１００００回）、２００００サンプルの「自己複製型特徴量」を得る方法である。2つ目が、平均値、標準偏差、尖度、歪度の分布特性値を用いて「ピアソンシステム乱数」を発生させて増幅特徴量を得る方法である。サンプル増幅法の背景にある理由は、１）観測データは、ランダムに発生するノイズや、粗雑な分布を示すことがあるので、モデル作成には理想的な場合でない場合がある、２）実際の痛みモニタリング時には、十分な数（例えば、１０００や１００００）のサンプルを得る時間が持てない場合が多い、ことがあり、少ないサンプルから大量サンプルを自己複製する技術は、個人に合わせた判別モデルを作成するためにも重要なモデル作成技術の一モジュールである。
上記のノイズ処理をされたエポックデータから、図１９の手順に示すように、２通りの方法でサンプル増幅用の実測値サンプルを得た。第一は、各エポックに対し、５秒の区間単位を用いて、１ポイントずつずらす移動平均をかけ、各刺激に対し１０００１サンプルを得た（レベル１と６に対し、各３０００３サンプル）。第二は、本当の意味での、小さい実測サンプルから大きな複製型サンプルを得る方法であり、同じ１０秒の刺激区間を1秒ごとに重複しないように切り出し、１０サンプルの平均値を得た。これらのサンプルに対し、２種類のサンプル増幅法を行った。1つ目が、レベル１と６の平均値と標準偏差を用いて「正規乱数」を発生し（本実施例では１００００回）、２００００サンプルの「自己複製型特徴量」を得る方法である。2つ目が、平均値、標準偏差、尖度、歪度の分布特性値を用いて「ピアソンシステム乱数」を発生させて増幅特徴量を得る方法である。サンプル増幅法の背景にある理由は、１）観測データは、ランダムに発生するノイズや、粗雑な分布を示すことがあるので、モデル作成には理想的な場合でない場合がある、２）実際の痛みモニタリング時には、十分な数（例えば、１０００や１００００）のサンプルを得る時間が持てない場合が多い、ことがあり、少ないサンプルから大量サンプルを自己複製する技術は、個人に合わせた判別モデルを作成するためにも重要なモデル作成技術の一モジュールである。
【0142】
（個人判別モデル作成、ならびに汎化テスト）
図１９に記したように、上記の自己複製型特徴量を用いて、一人のサンプルから個人判別モデルを機械学習により作成し、他の１６９名の痛みなしと不快な痛みの判別を行った。具体的な方法は、図２０に示した。
・モデル作成
１．個人（ID185）の複製型20,000サンプルを学習データとして用い、Lassoとベイズ最適化により、ロジスティック回帰モデルのハイパーパラメタ（λ）を決定し、特徴量（振幅5個）の重みづけ係数とモデル切片を決定する。
＊20,000サンプルは、正規乱数とピアソンシステム乱数を用いて増幅したものを用いる。
・モデル評価
２．決定した個人モデルを用いて、他169名のテストデータの判別推定を行う。
*テストデータは、同じ実験デザインの脳波データを刺激提示後5～15秒までの絶対平均振幅（レストでｚ値化）を用いた（3刺激×2レベル×169名）。
3. 169名の平均正答率をモデルの判別精度とする。
図１９に記したように、上記の自己複製型特徴量を用いて、一人のサンプルから個人判別モデルを機械学習により作成し、他の１６９名の痛みなしと不快な痛みの判別を行った。具体的な方法は、図２０に示した。
・モデル作成
１．個人（ID185）の複製型20,000サンプルを学習データとして用い、Lassoとベイズ最適化により、ロジスティック回帰モデルのハイパーパラメタ（λ）を決定し、特徴量（振幅5個）の重みづけ係数とモデル切片を決定する。
＊20,000サンプルは、正規乱数とピアソンシステム乱数を用いて増幅したものを用いる。
・モデル評価
２．決定した個人モデルを用いて、他169名のテストデータの判別推定を行う。
*テストデータは、同じ実験デザインの脳波データを刺激提示後5～15秒までの絶対平均振幅（レストでｚ値化）を用いた（3刺激×2レベル×169名）。
3. 169名の平均正答率をモデルの判別精度とする。
【0143】
・モデル作成・評価の実際
すなわち、被検者１名（ＩＤ１８５）の増幅された２００００サンプルを学習データとして、Ｌ１正則化とベイズ最適化を用いてハイパーパラメタ（λ）を決定し、ロジスティック回帰モデルの５特徴量係数を決定した。なお、ハイパーパラメタ決定には、他の方法、例えば、グリッドサーチを用いることもできる。両者の違いは、判別モデルである目的関数を決定する時に行うハイパーパラメタ決定過程において、ベイズ最適化が局所解（例.特定領域のパラメター値周辺、もしくはその組み合わせ）への偏向を避けつつ大局解（全域でゲインが最も大きいパラメター設定）を効率よく見つける方法であるのに対し、グリッドサーチは、全組み合わせ（例.サポートベクターマシンにおけるＣｏｓｔとγの組み合わせ）、もしくは設定範囲（本実施例におけるλの範囲）を網羅的に調べ、交差検定で検証する点に方法論上、大きな違いがある。得られた個人判別モデルを用いて、残り１６９名の痛みなしと不快な痛みのサンプルを一人ずつ判別推定した。図２０に記したように、ノイズレベルを下げるために、各レベルにおける各刺激の提示後５秒から１５秒全体の平均振幅を算出し、一人につき６サンプル（３刺激×２レベル）をテストデータとして用いた。最初に、２００００の実測サンプルの分布特性に基づき、正規乱数（平均、ＳＤの利用）とピアソンシステム乱数（平均、ＳＤ、尖度、歪度の利用）の発生法により増幅された複製サンプルを用いた結果を示す。次に、判別精度が高い乱数発生法を用いて、３０サンプルから２００００サンプルへの、文字通り「少数から多数へのサンプル増幅法」を試し、個人判別モデルを用いた痛み判別推定の可能性を検討する。
すなわち、被検者１名（ＩＤ１８５）の増幅された２００００サンプルを学習データとして、Ｌ１正則化とベイズ最適化を用いてハイパーパラメタ（λ）を決定し、ロジスティック回帰モデルの５特徴量係数を決定した。なお、ハイパーパラメタ決定には、他の方法、例えば、グリッドサーチを用いることもできる。両者の違いは、判別モデルである目的関数を決定する時に行うハイパーパラメタ決定過程において、ベイズ最適化が局所解（例.特定領域のパラメター値周辺、もしくはその組み合わせ）への偏向を避けつつ大局解（全域でゲインが最も大きいパラメター設定）を効率よく見つける方法であるのに対し、グリッドサーチは、全組み合わせ（例.サポートベクターマシンにおけるＣｏｓｔとγの組み合わせ）、もしくは設定範囲（本実施例におけるλの範囲）を網羅的に調べ、交差検定で検証する点に方法論上、大きな違いがある。得られた個人判別モデルを用いて、残り１６９名の痛みなしと不快な痛みのサンプルを一人ずつ判別推定した。図２０に記したように、ノイズレベルを下げるために、各レベルにおける各刺激の提示後５秒から１５秒全体の平均振幅を算出し、一人につき６サンプル（３刺激×２レベル）をテストデータとして用いた。最初に、２００００の実測サンプルの分布特性に基づき、正規乱数（平均、ＳＤの利用）とピアソンシステム乱数（平均、ＳＤ、尖度、歪度の利用）の発生法により増幅された複製サンプルを用いた結果を示す。次に、判別精度が高い乱数発生法を用いて、３０サンプルから２００００サンプルへの、文字通り「少数から多数へのサンプル増幅法」を試し、個人判別モデルを用いた痛み判別推定の可能性を検討する。
【0144】
（結果と考察）
個人判別モデルを作成する被検者（ＩＤ１８５）が、実験で使われた熱刺激に対し不快を感じていたかどうかを調べるために、ＣＯＶＡＳのスコアを算出した。刺激提示時間中の最大値を各レベル各刺激で算出しプロットしたものが図２１である。４０℃から４４℃までは、スコアは０に近く、ほとんど痛みを感じていないが、４８℃になると急激に痛みの不快度が増し、レベル６の５０℃では、ついに不快度は天井に到達した。そこで、個人判別モデル作成においては、レベル１と６の実測サンプルを用いてサンプル増幅を行った。各レベル３０００３サンプル（３刺激×１０００１サンプル）の分布は図２２に示した通りであり、レベル６の方が、平均値もＳＤも大きく、広がりも大きいのが視察から理解できる。
個人判別モデルを作成する被検者（ＩＤ１８５）が、実験で使われた熱刺激に対し不快を感じていたかどうかを調べるために、ＣＯＶＡＳのスコアを算出した。刺激提示時間中の最大値を各レベル各刺激で算出しプロットしたものが図２１である。４０℃から４４℃までは、スコアは０に近く、ほとんど痛みを感じていないが、４８℃になると急激に痛みの不快度が増し、レベル６の５０℃では、ついに不快度は天井に到達した。そこで、個人判別モデル作成においては、レベル１と６の実測サンプルを用いてサンプル増幅を行った。各レベル３０００３サンプル（３刺激×１０００１サンプル）の分布は図２２に示した通りであり、レベル６の方が、平均値もＳＤも大きく、広がりも大きいのが視察から理解できる。
【0145】
そこで、最初に正規乱数を用いて、両レベルの分布特性を自己複製し、各電極においてサンプルを１００００まで増幅した結果が図２３である。この増幅法により、正規分布の特性を持つ複製サンプルが得られ、個人レベルで判別モデル作成が可能となった。ロジスティック回帰モデルを使いベイズ最適化によりハイパーパラメタを決定して得られた個人判別モデルを図２４に示す。5電極特徴量の重みづけ係数は、「Ｆｐ１＝０．７９」、「Ｆｐ２＝１．２３」、「Ｆ３＝０．４７」、「Ｆ４＝２．８９」、「ＥＯＧ＝３．０１」、モデル切片は「―１．３６」であった。このロジスティック回帰モデルを使い、残り１６９名の痛み不快度をレベル１とレベル６のサンプル（2レベル×３サンプル）を使って判別推定したところ、平均判別精度は「７０％」に到達した。正答率が７０％以上の被検者は約半数、チャンスレベル以上は約7割に到達した。
【0146】
次に、同じ各レベルの３０００３サンプルを用い、ピアソンシステム乱数発生法を用いてサンプル増幅を行った。これは、平均とＳＤに加え、尖度と歪度の情報も用いて、サンプル増幅する方法である。図２５に、各レベルにおいて増幅されたサンプル（各１００００）の分布特性を示す。両レベルの分布は、レベル６がより裾野広く、オーバーラップしないエリアが良好にみられる。この複製サンプルから作成された判別モデルを図２６に示す。５特徴量の係数は、「Ｆｐ１＝１．８４」、「Ｆｐ２＝１．６１」、「Ｆ３＝４．７１」、「Ｆ４＝２．５３」、「ＥＯＧ＝５．２２」、モデル切片は「―３．３３」であった。このモデルを用いた１６９名の平均判別精度は「７２％」で、正規乱数の場合に比べ２％ほど高かった。正答率が７０％以上の被検者は約5割、チャンスレベル以上は4名ほど増え１１４名、約７割に達した。
【0147】
最後に、判別精度が比較して高かったピアソンシステム乱数発生法を用いて、文字通り、少ないサンプルから多数サンプルを増幅する方法を試した。ここでは、刺激提示後５～１５秒間を１０分割し、1秒ごとの絶対平均振幅を算出し、各レベル、各電極のサンプル数は、たった３０個にとどまった。各３０サンプルの分布特性（平均、ＳＤ、尖度、歪度）を算出し、各レベル、各電極で１００００サンプルを複製・増幅した。その結果を図２７に示す。図２５の多サンプルの場合と比較して、むしろ、自然で滑らかな曲線を持つ分布特性を示している。このサンプルを用いて、ロジスティック回帰モデルを作成したところ図２８のようなモデルが得られた（「Ｆｐ１＝０．８５」、「Ｆｐ２＝０．６２」、「Ｆ３＝１．５４」、「Ｆ４＝１．０９」、「ＥＯＧ＝０．８１」、モデル切片は「―０．９３」）。１６９名の判別精度は、わずかに上昇し、「７２％」に到達した。正答率７０％の被検者が６名増えたのが精度向上に寄与していると言える。
【0148】
本実施例は、さまざまな制約があると想定される痛みモニタリング状況において、いかに対象者から短い時間で効率よくサンプルを収集し、個人に特化した不快痛判別モデルを作成するかという課題に取り組んだ実験である。本実験は、各レベルの痛み刺激は３回だけ、しかも、１５秒間の短い提示条件を採用したものである。臨床現場で、その場で即興的に痛み判別モデルを作成する必要性が生じ、痛みなし、痛みありのリファレンス刺激を、本実施例と同じ条件で臨床群に提示したとしても、2分程度の時間しか必要としない。その限られたサンプルでも、本実施例におけるサンプル増幅法を用い、自己複製型特徴量を生成すれば、個人に特化した判別モデルが作成できることは特筆に値し、痛み判別装置における重要な要素技術の一つになりえると言える。また、実施例２、３において抽出された比較的数が少ない事象関連電位の特徴量においても、本実施例のサンプル増幅は、有効に使える意味で、拡張性が高い技術であるといえる。
【0149】
（注釈）
以上のように、本発明の好ましい実施形態を用いて本発明を例示してきたが、本発明は、請求の範囲によってのみその範囲が解釈されるべきであることが理解される。本明細書において引用した特許、特許出願および文献は、その内容自体が具体的に本明細書に記載されているのと同様にその内容が本明細書に対する参考として援用されるべきであることが理解される。本願は日本国特許出願特願２０１７－１４６５５３（２０１７年７月２８日）に対して優先権主張を行うものであり、それらの内容全体が本明細書において参照として援用される。
以上のように、本発明の好ましい実施形態を用いて本発明を例示してきたが、本発明は、請求の範囲によってのみその範囲が解釈されるべきであることが理解される。本明細書において引用した特許、特許出願および文献は、その内容自体が具体的に本明細書に記載されているのと同様にその内容が本明細書に対する参考として援用されるべきであることが理解される。本願は日本国特許出願特願２０１７－１４６５５３（２０１７年７月２８日）に対して優先権主張を行うものであり、それらの内容全体が本明細書において参照として援用される。
【産業上の利用可能性】
【0150】
本発明は、快不快を判定することができ、快不快に関する診断、治療をよりきめ細やかに行うことができる。
【符号の説明】
【0151】
１０００：刺激呈示部
１５００：対象
２０００：脳波データ取得部
２５００：脳波計
３０００：快不快判定値生成部
４０００：快不快判定部
５０９９：対象
５１００：快不快判定システム
５２００：脳波測定部
５２２０：脳波計
５２５０：脳波記録センサー
５２７０：脳波増幅部
５３００：快不快判定装置
５４００：脳波信号処理部
５５００：脳波特徴量抽出部
５６００：快不快判定部
５７００：快不快判定値生成部
５８００：判定レベル可視化部
５９００：刺激装置部
５９２０：刺激呈示部
５９６０：刺激情報可視化部
１５００：対象
２０００：脳波データ取得部
２５００：脳波計
３０００：快不快判定値生成部
４０００：快不快判定部
５０９９：対象
５１００：快不快判定システム
５２００：脳波測定部
５２２０：脳波計
５２５０：脳波記録センサー
５２７０：脳波増幅部
５３００：快不快判定装置
５４００：脳波信号処理部
５５００：脳波特徴量抽出部
５６００：快不快判定部
５７００：快不快判定値生成部
５８００：判定レベル可視化部
５９００：刺激装置部
５９２０：刺激呈示部
５９６０：刺激情報可視化部
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】
【図20】
【図21】
【図22】
【図23】
【図24】
【図25】
【図26】
【図27】
【図28】
