(57)【特許請求の範囲】
【請求項１】
睡眠心電図の集合から第１の単位時間ごとに算出された心拍変動特徴量を睡眠段階別に集計した心拍変動情報を記憶した心拍変動記憶部と、
前記睡眠心電図の時間軸を前記第１の単位時間より長い第２の単位時間で分割した各時間帯ごとに、前記第１の単位時間ごとの睡眠段階の出現頻度を集計した出現頻度情報を記憶した出現頻度記憶部と、
前記分割した各時間帯ごとに、前記第１の時間ごとの睡眠段階の遷移の頻度を集計した遷移頻度情報を記憶した遷移頻度記憶部と、
推定対象の被験者から取得した睡眠心電図から、前記第２の単位時間で分割した各時間帯ごとに、前記第１の時間ごとの心拍変動の特徴量を算出する心拍変動特徴量算出手段と、
前記算出された心拍変動の特徴量の時間的な変化と、前記心拍変動記憶部、出現頻度記憶部および遷移頻度記憶部に記憶された情報とに基づいて、前記被験者の前記第１の単位時間ごとの睡眠段階の確率分布を推定する推定手段と
を具備することを特徴とする睡眠段階推定装置。
【請求項２】
前記被験者の睡眠心電図集合として与えられる、前記第１の単位時間ごとの睡眠段階の正解ラベルが付与された睡眠時の時間情報を含む睡眠心電図の集合を取得する手段と、
前記取得された睡眠心電図の集合から、前記第１の単位時間ごとに心拍変動特徴量を算出し、当該算出された心拍変動特徴量を睡眠段階別に集計して前記心拍変動記憶部に格納する手段と、
前記睡眠心電図集合の時間軸を前記第２の単位時間で分割した各時間帯ごとに、前記第１の単位時間ごとの睡眠段階の出現頻度を集計して、集計された出現頻度情報を前記出現頻度記憶部に格納する手段と、
前記分割した各時間帯ごとに、前記第１の時間ごとの睡眠段階の遷移の頻度を集計して、集計された遷移頻度情報を前記遷移頻度記憶部に格納する手段と
を、さらに具備することを特徴とする請求項１記載の睡眠段階推定装置。
【請求項３】
前記心拍変動特徴量算出手段は、心拍変動特徴量として、心拍変動RRI（R-R Interval）の平均値、RRIの標準偏差、前後するRRI間の差の二乗平均値の平方根（RMSSD）、前後するRRI間の時間的差分が予め設定した時間値以上となった回数の発生割合（pNN50）、前後するRRIのローレンツプロットが描く形状の長辺の長さＬ、RRIのローレンツプロットが描く形状の短辺の長さＴ、ＬおよびＴから求められる交感神経活動指標（CSI: Cardiac Sympathetic Index）、ＬおよびＴから求められる副交感神経の活動指標（CVI: Cardiac Vagal Index）、RRIのパワースペクトル波形の高周波成分（HF）、または低周波成分（LF）のうちの少なくとも１つ、或いはこれらの特徴量から算出されるその他の特徴量を算出することを特徴とする請求項１または２記載の睡眠段階推定装置。
【請求項４】
睡眠心電図の集合から第１の単位時間ごとに算出された心拍変動特徴量を睡眠段階別に集計した心拍変動情報を記憶した心拍変動記憶部と、前記睡眠心電図の時間軸を前記第１の単位時間より長い第２の単位時間で分割した各時間帯ごとに、前記第１の単位時間ごとの睡眠段階の出現頻度を集計した出現頻度情報を記憶した出現頻度記憶部と、前記分割した各時間帯ごとに、前記第１の時間ごとの睡眠段階の遷移の頻度を集計した遷移頻度情報を記憶した遷移頻度記憶部とを有するモデルデータベースを備えた睡眠段階推定装置が実行する睡眠方法推定動作方法であって、
前記睡眠段階推定装置が、推定対象の被験者から取得した睡眠心電図から、前記第２の単位時間で分割した各時間帯ごとに、前記第１の時間ごとの心拍変動の特徴量を算出する過程と、
前記睡眠段階推定装置が、前記算出された心拍変動の特徴量の時間的な変化と、前記心拍変動記憶部、出現頻度記憶部および遷移頻度記憶部に記憶された情報とに基づいて、前記被験者の前記第１の単位時間ごとの睡眠段階の確率分布を推定する過程と
を具備することを特徴とする睡眠段階推定動作方法。
【請求項５】
請求項１乃至３のいずれかに記載の睡眠段階推定装置が具備する各手段が行う処理を、前記睡眠段階推定装置が備えるコンピュータに実行させるプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
この発明は、人の睡眠段階を推定するための睡眠段階推定装置、方法およびプログラムに関する。
【背景技術】
【０００２】
睡眠段階とは、眠りの深さを表す指標であり、国際的基準では２〜６段階で定義されている。睡眠段階を時系列的に表した睡眠経過図は、医療やヘルスケアの分野において睡眠の状態を把握する上で、重要な情報である。睡眠段階は、ポリグラフ検査にて得られる脳波データまたは筋電データを利用して、専門の技師が視察判定ルールであるＲ＆Ｋ法に基づいて時間区間（約３０秒）ごとに判定される。しかし、このポリグラフ検査は、ユーザにとって身体的負荷が大きいため、近年では身体的負荷の少ない計測装置で取得可能な心拍データを用いた睡眠段階推定技術に関する研究が盛んに行われている。
【０００３】
心拍データを用いた従来の睡眠段階推定手法は、例えば心拍データから睡眠段階に関連する特徴量を複数生成し、機械学習（サポートベクターマシーンやランダムフォレスト等）を利用して推定するものとなっている（例えば非特許文献１を参照）。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】Xiao,M et al., “Sleep stages classification based on heart rate variability and random forest” Biomedical Signal Processing and Control, p624-633, (2013).
【発明の概要】
【発明が解決しようとする課題】
【０００５】
ところが、心拍に反映される自律神経活動は、睡眠時の脳活動以外の影響、例えば血圧調整や体温調整、身体の動きや呼吸等の影響も受けて変化する。このため、非特許文献１に記載された技術では、脳活動から得られる睡眠段階に対して心拍から完全に一致して推定することは困難である。
【０００６】
この発明は上記事情に着目してなされたもので、その目的とするところは、心拍が脳活動以外の影響を強く受けた場合でも、睡眠段階をより正確に推定できるようにした睡眠段階推定装置、方法およびプログラムを提供することにある。
【課題を解決するための手段】
【０００７】
上記目的を達成するためにこの発明の第１の観点は、睡眠心電図の集合から第１の単位時間ごとに算出された心拍変動特徴量を睡眠段階別に集計した心拍変動情報を記憶した心拍変動記憶部と、前記睡眠心電図の時間軸を前記第１の単位時間より長い第２の単位時間で分割した各時間帯ごとに、前記第１の単位時間ごとの睡眠段階の出現頻度を集計した出現頻度情報を記憶した出現頻度記憶部と、前記分割した各時間帯ごとに、前記第１の時間ごとの睡眠段階の遷移の頻度を集計した遷移頻度情報を記憶した遷移頻度記憶部とを、モデルデータベースとして備える。そして、推定対象の被験者から取得した睡眠心電図から、前記第２の単位時間で分割した各時間帯ごとに、前記第１の時間ごとの心拍変動の特徴量を算出し、この算出された心拍変動の特徴量の時間的な変化と、前記心拍変動記憶部、出現頻度記憶部および遷移頻度記憶部に記憶された情報とに基づいて、前記被験者の前記第１の単位時間ごとの睡眠段階の確率分布を推定するようにしたものである。
【０００８】
この発明の第１の観点には以下のような態様が考えられる。  
その第１の態様は、前記被験者の睡眠心電図集合として与えられる、第１の単位時間ごとの睡眠段階の正解ラベルが付与された睡眠時の時間情報を含む睡眠心電図の集合を取得し、この取得された睡眠心電図の集合から、前記第１の単位時間ごとに心拍変動特徴量を算出し、当該算出された心拍変動特徴量を睡眠段階別に集計して心拍変動記憶部に格納する手段と、前記睡眠心電図集合の時間軸を前記第２の単位時間で分割した各時間帯ごとに、前記第１の単位時間ごとの睡眠段階の出現頻度を集計して出現頻度記憶部に格納する手段と、前記第２の単位時間で分割した各時間帯ごとに、前記第１の時間ごとの睡眠段階の遷移の頻度を集計して遷移頻度記憶部に格納する手段を、さらに備えるようにしたものである。
【０００９】
第２の態様は、心拍変動の特徴量を算出する際に、心拍変動特徴量として、心拍変動RRI（R-R Interval）の平均値、RRIの標準偏差、前後するRRI間の差の二乗平均値の平方根（RMSSD）、前後するRRI間の時間的差分が予め設定した時間値以上となった回数の発生割合（pNN50）、前後するRRIのローレンツプロットが描く形状の長辺の長さＬ、RRIのローレンツプロットが描く形状の短辺の長さＴ、ＬおよびＴから求められる交感神経活動指標（CSI: Cardiac Sympathetic Index）、ＬおよびＴから求められる副交感神経の活動指標（CVI: Cardiac Vagal Index）、RRIのパワースペクトル波形の高周波成分（HF）、または低周波成分（LF）のうちの少なくとも１つ、或いはこれらの特徴量から算出されるその他の特徴量を算出するようにしたものである。
【発明の効果】
【００１１】
この発明の第１の観点によれば、推定対象者の心電図データから算出した心拍変動特徴量と、モデルデータベースに予め記憶された心拍変動情報、出現頻度情報および遷移頻度情報とに基づいて、上記推定対象者の睡眠段階が推定される。このため、心拍が睡眠時の脳活動以外の影響、例えば血圧調整や体温調整、身体の動きや呼吸等の影響も受けて変化したとしても、睡眠段階を正確に推定することが可能となる。
【００１２】
この発明の第１の観点に係る第１の態様によれば、学習フェーズにおいて心拍変動情報、出現頻度情報および遷移頻度情報を生成してモデルデータベースに格納する手段を備えたことにより、睡眠段階推定装置において自立的に最適な心拍変動情報、出現頻度情報および遷移頻度情報を生成し格納することが可能となる。
【００１３】
この発明の第１の観点に係る第２の態様によれば、心拍変動の特徴量として多種類の特徴量の中から適当なものを選択することができる。
【００１５】
すなわちこの発明によれば、心拍が脳活動以外の影響を強く受けた場合でも、睡眠段階をより正確に推定することが可能な睡眠段階推定装置、方法およびプログラムを提供することができる。
【図面の簡単な説明】
【００１６】
【図１】睡眠段階の出現頻度の一例を時間帯別に示した図。
【図２】睡眠段階遷移パターンの出現確率の一例を示す図。
【図３】この発明の一実施形態に係る睡眠段階推定装置の機能構成を示すブロック図。
【図４】図３に示した睡眠段階推定装置の学習フェーズにおける処理手順と処理内容を示すフローチャート。
【図５】図３に示した睡眠段階推定装置の推定フェーズにおける処理手順と処理内容を示すフローチャート。
【図６】図３に示した睡眠段階推定装置が備える睡眠段階別心拍変動テーブルの一例を示す図。
【図７】図３に示した睡眠段階推定装置が備える時間別睡眠段階出現頻度テーブルの一例を示す図。
【図８】図３に示した睡眠段階推定装置が備える時間別睡眠段階遷移頻度テーブルの一例を示す図。
【図９】時系列心拍データの観測手法と睡眠段階推定値の出力手法の一例を示す図。
【図１０】ローレンツプロットの一例を示す図。
【図１１】心拍変動特徴量の種類を一覧表示した図。
【発明を実施するための形態】
【００１７】
以下、図面を参照してこの発明に係わる実施形態を説明する。  
（原理）
この発明は、入眠後の経過時間に従い睡眠段階の出現頻度と遷移確率が変化することに着目し、心拍が脳活動以外の影響を強く受けた場合でも、より尤もらしい睡眠段階を推定できるようにしたものである。例えば、入眠後の深い眠りから、突然覚醒して再び深い眠りに落ちる可能性は低いが、心拍のみで推定すると０．５２の確率でその状態が推定される場合がある。この発明では、このようなめったに起こらない遷移や推定の自信度が低い場合に対する推定確率を下げることで、推定精度を向上させる。
【００１８】
図１は、４５名分の睡眠段階データを利用して、３０分ごとの各睡眠段階の出現回数を算出した結果の一例を示している。この例から明らかなように、入眠開始から３０分間におけるREM（レム睡眠）の出現回数（図１（ａ））は、図１（ｂ），（ｃ）にそれぞれ示す入眠から１８０分および３６０分経過後の３０分間におけるREMの出現回数より少なくなる。
【００１９】
また、図２は同じく４５名分の睡眠段階データを利用して、３０分ごとの各睡眠段階の遷移パターンの出現確率を算出した結果の一例を示している。同図に示すように、入眠開始から３０分間におけるREMからNREM（ノンレム睡眠）への遷移確率（図２（ａ））は、図２（ｂ），（ｃ）にそれぞれ示す入眠から１５０分経過後および３００分経過後の３０分間におけるREMからNREMへの遷移確率よりも高くなることが分かる。
【００２０】
そこでこの発明では、経過時間に応じた睡眠段階の出現頻度と遷移確率を学習フェーズにより睡眠段階の生成モデルに組み込み、推定フェーズにおいて、ユーザの心拍データの特徴量を求め、この心拍データの特徴量と上記睡眠段階の生成モデルに基づいて睡眠段階を推定するようにしている。このようにすることで、心拍が脳活動以外の影響を強く受けた場合でも、睡眠段階をより正確に推定することが可能となる。
【００２１】
［第１の実施形態］
（構成）
図３は、この発明の一実施形態に係る睡眠段階推定装置の機能構成を示すブロック図である。  
本実施形態における睡眠段階推定装置１は、例えばサーバコンピュータ又はパーソナルコンピュータからなり、制御ユニット１０と、記憶ユニット２０と、入出力インタフェースユニット３０を備えている。このうち入出力インタフェースユニット３０は、心電計２から出力された心電データを取り込むと共に、キーボードまたはタッチパネル等の入力部３および液晶デバイス等の表示部４との間でそれぞれ操作データおよび表示データの入出力を行う。
【００２２】
記憶ユニット２０は、記憶媒体としてＨＤＤ（Hard Disk Drive）やＳＳＤ（Solid State Drive）等の書き込みおよび読み出しが可能な不揮発性メモリを用いたもので、本実施形態を実施する上で必要な記憶部として、モデルデータベース２１と、推定データ記憶部２２を備えている。
【００２３】
モデルデータベース２０は、心拍変動テーブル２１１、出現頻度テーブル２１２および遷移頻度テーブル２１３を有する。心拍変動テーブル２１１は、睡眠段階別に、心拍特徴ベクトルの識別情報（心拍特徴ベクトルＩＤ）に関連付けて心拍特徴量レベルと出現頻度を格納するために用いられる。図６にその一例を示す。出現頻度テーブル２１２は、時間帯別に当該時間帯に出現する睡眠段階と出現頻度を格納するために用いられる。図７にその一例を示す。遷移頻度テーブル２１３は、時間帯別に当該時間帯に発生する睡眠段階遷移パターンとその出現頻度を格納するために用いられる。図８はその一例を示すものである。推定データ記憶部２２は、睡眠段階の推定結果を表すデータを格納するために用いられる。
【００２４】
制御ユニット１０は、中央処理ユニット（ＣＰＵ；Central Processing Unit）とメモリを備えたもので、本実施形態を実施する上で必要な制御処理機能として、心拍変動特徴量算出部１１と、睡眠段階別心拍変動解析部１２と、時間別睡眠段階出現頻度解析部１３と、時間別睡眠段階遷移頻度解析部１４と、睡眠段階確率推定部１５と、推定データ出力制御部１６を有している。なお、これらの制御処理機能はいずれも図示しないプログラムメモリに格納されたプログラムを上記ＣＰＵに実行させることにより実現される。
【００２５】
心拍変動特徴量算出部１１は、心電計２から入出力インタフェースユニット３０を介して取り込んだ心電図データをもとに、心拍変動の特徴量を算出する。具体的には、心電波形からピーク間隔を表すRRI（R-RInterval）を検出し、この検出されたRRIをｎ秒ごとに分割して、フレーム長Ｎ[sec]、スライド幅ｎ[sec]で、心拍変動の特徴ベクトルｈt ＝｛ｈt1，ｈt2，…，ｈtk｝を算出する。ここで、ｔは経過時間を示し、ｔ＝｛ｎ，２*ｎ、３*ｎ，…｝のようにｎ秒ずつカウントされる。
【００２６】
睡眠段階別心拍変動解析部１２は、学習フェーズにおいて、上記心拍変動特徴量算出部１１により算出されたｈt を、時刻情報ｔをもとに、入力部３において入力されたｎ秒ごとの睡眠段階の正解ラベルと対応付け、各睡眠段階の各心拍特徴量の分布を算出する。そして、この睡眠段階ごとに算出された各心拍特徴量の分布を、モデルデータベース２１内の心拍変動テーブル２１１に格納する処理を行う。
【００２７】
時間別睡眠段階出現頻度解析部１３は、学習フェーズにおいて、入力部３において入力されたｎ秒ごとの睡眠段階正解ラベルの時系列データをｑ分（例えば３０分）ごとに分割し、この分割された各時間帯ｍにおける各睡眠段階ｉの出現回数を計算する。そして、この算出された各睡眠段階ｉの出現回数を、モデルデータベース２１内の出現頻度テーブル２１２に格納する処理を行う。
【００２８】
時間別睡眠段階遷移頻度解析部１４は、学習フェーズにおいて、上記時間別睡眠段階出現頻度解析部１３と同様に、ｑ分ごとに分割された各時間帯ｍにおいて、ｎ秒ごとの各睡眠段階遷移の出現回数を算出する。そして、この算出された各睡眠段階遷移の出現回数を、モデルデータベース２１内の上記遷移頻度テーブル２１３に格納する処理を行う。
【００２９】
睡眠段階確率推定部１５は、推定フェーズにおいて、心電計２から入力された心電図データをもとに上記心拍変動特徴量算出部１１により算出された心拍変動特徴ベクトルと、上記モデルデータベース２１の各テーブル２１１，２１２，２１３に記憶されている情報とを用いて、ｎ秒ごとの睡眠段階Ｓt を推定する。そして、その推定結果を表すデータを推定データ記憶部２２に格納する処理を行う。
【００３０】
推定データ出力制御部１６は、入力部３において入力された出力指示に応じ、上記推定データ記憶部２２から該当するユーザの該当する時間帯における睡眠段階推定結果を表すデータを読み出す。そして、この読み出された推定データを表示させるための表示データを生成し、当該生成された表示データを入出力インタフェースユニット３０を介して表示部４へ出力し表示させる処理を行う。
【００３１】
（動作）
次に、以上のように構成された装置による睡眠段階推定動作を説明する。  
（１）学習フェーズ
学習フェーズでは、被験者の心電波形を計測しながら、医療従事者等の測定者がｎ秒ごとに被験者の睡眠段階を目視観察してその観察結果を「睡眠段階の正解ラベル」として入力部３により入力する。そして、睡眠段階推定装置１が、上記心電波形から心拍変動の特徴量を算出して心拍変動テーブル２１に格納すると共に、上記入力された睡眠段階の正解ラベルをもとに時間別の睡眠段階出現頻度と睡眠段階遷移頻度を算出してその算出値をそれぞれ出現頻度テーブル２１２および遷移頻度テーブル２１３に格納することにより、行われる。図４は、この学習フェーズにおける睡眠段階推定装置１の処理手順と処理内容を示すフローチャートである。
【００３２】
学習フェーズによる計測が開始されると、心拍変動特徴量算出部１１の制御の下、先ずステップＳ１１により心電計２から出力された被験者の心電図データが入出力インタフェースユニット３０を介して取り込まれる。そして、ステップＳ１２により上記取り込まれた心電図データをもとに心拍変動の特徴量が以下のように算出される。
【００３３】
すなわち、先ず心電図データから心電波形のピーク間隔を表すRRI（R-RInterval）が検出され、この検出されたRRIがｎ秒ごとに分割される。そして、フレーム長Ｎ[sec]、スライド幅ｎ[sec]として、心拍変動の特徴ベクトルｈt ＝｛ｈt1，ｈt2，…，ｈtk｝が算出される。ここで、ｔは経過時間を示し、ｔ＝｛ｎ，２*ｎ、３*ｎ，…｝のようにｎ秒ずつカウントされる。具体例を図９に示す。
【００３４】
心拍変動特徴量としては、例えば心拍変動RRI（RRInterval）の平均値、RRIの標準偏差、前後するRRI間の差の二乗平均値の平方根（RMSSD）、前後するRRI間の時間的差分が５０ミリ秒以上となった回数の発生割合（pNN50）、前後するRRIのローレンツプロットが描く形状の長辺の長さＬ、RRIのローレンツプロットが描く形状の短辺の長さＴ、ＬおよびＴから求められる交感神経活動指標（CSI: Cardiac Sympathetic Index）、ＬおよびＴから求められる副交感神経の活動指標（CVI: Cardiac Vagal Index）、RRIのパワースペクトル波形の高周波成分（HF）、または低周波成分（LF）のうちの少なくとも１つ、或いはこれらの特徴量から算出されるその他の特徴量が用いられる。なお、RRIのローレンツブロットの一例を図１０に、また上記した各心拍変動特徴量の定義を図１１にそれぞれ示す。なお、心拍変動特徴量の詳細な算出方法については、例えば非特許文献１に記載されている。
【００３５】
上記心拍変動特徴ベクトルｈt の各特徴量の値は、例えば０〜１の範囲に正規化した上でＬレベルで離散値化しておくとよい。例えば、ｔ＝ｎ（最初のｎ秒区間）のときの特徴ベクトルの１要素目はｈn1＝｛０，１，…，Ｌ｝のように表される。
【００３６】
続いて睡眠段階別心拍変動解析部１２において、上記心拍変動特徴量算出部１１により算出された心拍変動特徴ベクトルｈt をもとに、睡眠段階別の心拍変動が以下のように解析される。すなわち、上記算出された心拍変動特徴ベクトルｈt は、時刻情報tをもとに、入力部３で入力されるｎ秒ごとの睡眠段階の正解ラベルと対応付けられる。そして、各睡眠段階の各心拍特徴量の分布が算出される。
【００３７】
次にステップＳ１３において、時間別睡眠段階出現頻度解析部１３により時間帯（ｍ分）ごとの各睡眠段階の出現頻度が以下のように集計される。すなわち、先ず入力部３において入力されたｎ秒ごとの睡眠段階正解ラベルの時系列データがｑ分（例えば３０分）ごとに分割され、この分割された各時間帯ｍにおける各睡眠段階ｉの出現回数が計算される。このときｍ＝１は０分〜（ｑ×１）分の区間、ｍ＝２は（ｑ×１）分〜（ｑ×２）分の区間と表される。そして、この算出された各睡眠段階ｉの出現回数は、例えば図７に示すように、時間帯に対応付けられて睡眠段階を表す情報と共に出現頻度テーブル２１２に格納される。
【００３８】
続いてステップＳ１４において、時間別睡眠段階遷移頻度解析部１４により、時間帯（ｍ分）ごとの各睡眠段階遷移パターンの出現頻度が以下のように集計される。すなわち、先ずｑ分ごとに分割された各時間帯ｍにおいて、ｎ秒ごとの各睡眠段階遷移の出現回数が算出される。
【００３９】
例えば、ある時間帯ｍにおいて、最初の０〜ｎ秒区間の睡眠段階が「覚醒（WAKE）」、ｎ〜２ｎ秒区間の睡眠段階が「レム睡眠（REM）」の場合、ｎ〜２ｎ秒区間の睡眠段階の遷移は「覚醒（WAKE）→レム睡眠（REM）」の遷移パターンとしてカウントされる。なお、最初の０〜ｎ秒区間の睡眠段階の遷移パターンはカウントしなくてよいものとする。そして、以上のように算出された各睡眠段階遷移の出現回数が、例えば図８に示すように、時間帯に対応付けられて睡眠段階遷移パターンを表す情報と共に遷移頻度テーブル２１３に格納される。
【００４０】
最後にステップＳ１５において、睡眠段階別心拍変動解析部１２により、睡眠段階ごとの心拍変動特徴量が集計され、その算出結果である出現頻度が、例えば図６に示すように心拍特徴ベクトルＩＤおよび心拍特徴量レベルと共に、各睡眠段階に対応付けられて睡眠段階別心拍変動テーブル２１１に格納される。例えば、睡眠段階の正解ラベルが「覚醒（WAKE）」のときの全時刻をt_wake ＝｛ｎ，５*ｎ，８*ｎ｝としたとき、全時刻に対してｈt_wake,1＝ｌとなるカウントが記憶される。このとき、心拍変動特徴量の離散値は０≦ｌ≦Ｌなので、各値についてカウントされる。
【００４１】
（２）推定フェーズ
推定フェーズでは、推定対象者の心電図データのみを計測し、この計測された心電図データから算出される心拍変動特徴量と、心拍変動テーブル２１１、出現頻度テーブル２１２および遷移頻度テーブル２１３に格納されたモデルデータに基づいて、上記推定対象者の睡眠段階Ｓt がｎ秒ごとに推定される。図５はその処理手順と処理内容を示すフローチャートである。
【００４２】
先ずステップＳ２１において、心拍変動特徴量算出部１１により、心電計２から推定対象者の１回の睡眠における心電図データが取り込まれ、また睡眠段階別心拍変動解析部１２により、入力部３により測定者が入力したｎ秒ごとの睡眠段階の正解ラベルが取り込まれる。そしてステップＳ２２において、心拍変動特徴量算出部１１により、上記取り込まれた心電図データからｎ秒ごとの心拍変動特徴量が算出される。その算出手法は、前記学習フェーズの場合と同じである。
【００４３】
次にステップＳ２３において、睡眠段階確率推定部１５の制御の下、推定対象者の睡眠段階の推定が以下のように行われる。すなわち、先ずステップＳ２３によりｎ秒ごとの睡眠段階変数がランダムに初期化される。例えば、ギブスサンプリングカウンタｇが１で初期化される。
【００４４】
次にステップＳ２４において、ｎ秒ごとの睡眠段階変数がギブスサンプリングにより更新され、ステップＳ２５によりギブスサンプリングカウンタｇがギブスサンプリング繰り返し数になったと判定されるまで、以下の処理が繰り返し実施される。
【００４５】
(A) 初期値として、時間ｔ＝０、睡眠段階Ｓ0 ＝”Wake”に設定する。  
(B) 時間ｔが、睡眠時間Ｔよりも小さい間では以下の処理を実行する。  
(1) 睡眠段階Ｓ_{t,g} を、時間ｔ−１の睡眠段階Ｓt-1、時間帯Ｃt 、時間ｔにおける心拍変動特徴量ｈt を用いてＰ（Ｓt ＝ｉ｜ｈt ，Ｓt-1＝ｊ，Ｃt ＝ｍ）の値に従ってサンプリングする。  
(2) 時間ｔ＝ｔ＋１とする。  
(C) ギブスサンプリングカウンタｇをインクリメント（ｇ＝ｇ＋１）する。
【００４６】
続いてステップＳ２６において、各時間ｔについて、ｇ個の睡眠段階サンプリング履歴集合｛Ｓ_{t,g} ｝から、各睡眠段階の出現分布を集計し、時間ｔにおける睡眠段階確率分布として出力する。睡眠段階を出力する場合は、サンプリング履歴集合｛Ｓ_{t,g} ｝における最頻値となる睡眠段階を、推定結果として推定データ記憶部２２に格納する。
【００４７】
Ｐ（Ｓt ＝ｉ｜ｈt ，Ｓt-1＝ｊ，Ｃt ＝ｍ）は以下の式によって得られる。  
【数１】
0005993968.tif 000002
【００４８】
ここで、θmiは時間帯mにおいて睡眠段階がｉになる確率を表し、以下の式にて算出される。なお、Ｍmiは時間帯ｍにて睡眠段階ｉが出現する回数を表す。
【００４９】
【数２】
0005993968.tif 000003
【００５０】
また、λmjiは時間帯ｍにおいて、睡眠段階がｊからｉへ遷移する確率を表し、以下の式により算出される。なお、Ｍmjiは時間帯ｍにて睡眠段階がｊからｉへ遷移する状態が出現する回数を表す。  
【数３】
0005993968.tif 000004
【００５１】
さらにφikhtkは、睡眠段階ｉにおける各心拍特徴量の出現確率を表すもので、以下のように算出される。  
【数４】
0005993968.tif 000005
【００５２】
睡眠段階ｉであるときの時間ｔにて新たに得られる心拍変動特徴ベクトルｋ番目の要素ｈtkがｌであったとき、Ｍiklはモデルデータベースを参照して、ｈtkが値ｌを満たす時の出現頻度の回数を表している。
【００５３】
なお、上記推定データ記憶部２２に記憶された推定対象者の睡眠段階の推定結果を表すデータは、推定データ出力制御部１６の制御の下、推定データ記憶部２２から読み出される。そして、この読み出された推定データをもとに表示データが生成され、この表示データは入出力インタフェースユニット３０を介して表示部４に供給され表示される。図９の下段にその表示結果の一例を示す。
【００５４】
（実施形態の効果）
以上詳述したように一実施形態では、学習フェーズにおいて、被験者の心電図データから睡眠段階ごとの心拍変動特徴量を算出して心拍変動テーブル２１１に格納すると共に、時間経過に応じた睡眠段階の出現頻度と遷移頻度を求めてこれらをそれぞれ出現頻度テーブル２１２および遷移頻度テーブル２１３に格納する。そして推定フェーズにおいて、推定対象者の心電図データから算出した心拍変動特徴量と、上記各テーブル２１１，２１２，２１３に格納された睡眠段階のモデルデータとに基づいて、上記推定対象者の睡眠段階を推定するようにしている。
【００５５】
したがって、心拍が睡眠時の脳活動以外の影響、例えば血圧調整や体温調整、身体の動きや呼吸等の影響も受けて変化したとしても、睡眠段階を正確に推定することが可能となる。また、学習フェーズを備えることで、最適な心拍変動テーブル２１１、出現頻度テーブル２１２および遷移頻度テーブル２１３を構築することができる。
【００５６】
［他の実施形態］
なお、この発明は上記実施形態に限定されるものではない。例えば、前記実施形態では睡眠段階推定装置１内に記憶ユニット２０を設けた場合を例にとって説明したが、記憶ユニット２０をクラウドに設けられたデータベースサーバなどに備え、睡眠段階推定装置１とこのデータベースサーバとの間で通信を行うことによりデータの書き込みおよび読み出しを行うようにしてもよい。
【００５７】
また、前記実施形態では心電計により得られる心電図データをもとに心拍変動特徴量を算出するようにしたが、脈波を計測してその計測データをもとに心拍変動特徴量を算出するようにしてもよい。
【００５８】
さらに、推定結果を表すデータは、通信回線を介して他の端末またはサーバに送信するようにしてもよく、その他制御ユニットが備える制御機能、当該制御機能による処理手順と処理内容、推定結果の表示フォーマットなどについても、この発明の要旨を逸脱しない範囲で種々変形して実施可能である。
【００５９】
要するにこの発明は、上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合せにより種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態に亘る構成要素を適宜組み合せてもよい。
【符号の説明】
【００６０】
１…睡眠段階推定装置、２…心電計、３…入力部、４…表示部、１０…制御ユニット、１１…心拍変動特徴量算出部、１２…睡眠段階別心拍変動解析部、１３…時間別睡眠段階出現頻度解析部、１４…時間別睡眠段階遷移頻度解析部、１５…睡眠段階確率推定部、１６…推定データ出力制御部、２０…記憶ユニット、２１…モデルデータベース、２２…推定データ記憶部、２１１…心拍変動テーブル、２１２…出現頻度テーブル、２１３…遷移頻度テーブル、３０…入出力インタフェースユニット。
【図１】
0005993968.tif 000006
【図２】
0005993968.tif 000007
【図３】
0005993968.tif 000008
【図４】
0005993968.tif 000009
【図５】
0005993968.tif 000010
【図６】
0005993968.tif 000011
【図７】
0005993968.tif 000012
【図８】
0005993968.tif 000013
【図９】
0005993968.tif 000014
【図１０】
0005993968.tif 000015
【図１１】
0005993968.tif 000016