JP7419719B2 - 睡眠ステージ推定装置、睡眠ステージ推定方法及びプログラム - Google Patents

Description

本開示は、睡眠状態モニタリング技術に関する。

睡眠状態を判定する睡眠状態モニタリングシステムとして、脳波による睡眠ステージ（浅い睡眠、深い睡眠、レム睡眠、ノンレム睡眠等）の判定が広範に利用されている（例えば、特許文献１）。また、脳波に比較して取得が容易な体動や脈波を利用した睡眠モニタリングシステムもまた提案されている（例えば、特許文献２）。

特開２０１１－８３３９３号公報 特開２０１８－１６１４３２号公報

上記特許文献１及び特許文献２においては、ある時点の睡眠ステージを判定する際に、時系列的に遷移する睡眠ステージ遷移の発生確率が考慮されていない、即ち、睡眠ステージの判定に際し、時系列が考慮されていないという課題があった。

本発明は上記課題に鑑みてなされたもので、その目的とするところは、時系列を考慮した睡眠ステージ推定技術を提供することである。

上記課題を解決するため、本開示の一態様は、被検者の拍動データ及び体動データを取得する被検者データ取得部と、前記拍動データから特徴量系列を算出し、学習済み睡眠ステージ確率推定モデルを利用して、前記算出した特徴量系列から前記被検者の睡眠ステージ確率系列を推定する睡眠ステージ確率推定部と、前記体動データから体動量系列を算出し、学習済み睡眠ステージ遷移確率推定モデルを利用して、前記算出した体動量系列から前記被検者の睡眠ステージ遷移確率系列を推定する睡眠ステージ遷移確率推定部と、前記睡眠ステージ確率系列を観測素性とし、前記睡眠ステージ遷移確率系列を遷移素性とする素性関数から構成される学習済み条件付き確率場モデルを利用して、前記睡眠ステージ確率系列及び前記睡眠ステージ遷移確率系列から前記被検者の睡眠ステージ系列を推定する睡眠ステージ推定部と、を有し、前記学習済み条件付き確率場モデルは、前記睡眠ステージ確率系列または前記睡眠ステージ遷移確率系列の学習系列データと正解ラベル列の確率の交差エントロピーを最小化するパラメータを求めることにより学習される睡眠ステージ推定装置に関する。

本開示によると、時系列を考慮した睡眠ステージ推定技術を提供することができる。

本開示の一実施例による睡眠ステージ推定装置を示す概略図である。 本開示の一実施例による睡眠ステージ推定装置のハードウェア構成を示すブロック図である。 本開示の一実施例による睡眠ステージ推定装置の機能構成を示すブロック図である。 本開示の一実施例による睡眠ステージ推定処理を示すフローチャートである。

以下の実施例では、被検者の睡眠中の時系列データから睡眠ステージ系列を推定する睡眠ステージ推定装置が開示される。

［本開示の概要］
後述される実施例を概略すると、図１に示されるように、睡眠ステージ推定装置１００は、被検者の睡眠中の脈波や心拍等の拍動データ及び加速度や角速度等の体動データを取得すると、拍動データから一定時間毎の特徴量（例えば、心拍間隔の平均、標準偏差、ＬＦ／ＨＦ比など）を算出し、時系列に得られた特徴量から特徴量系列を生成する。そして、睡眠ステージ推定装置１００は、学習済み睡眠ステージ確率推定モデルを利用して、各時刻における特徴量から、被検者が各睡眠ステージにいると推定される確率を示す睡眠ステージ確率を推定する。例えば、睡眠ステージ確率推定モデルは、ロジスティック回帰モデルによって実現されてもよく、３つ以上の睡眠ステージ（例えば、浅い睡眠、深い睡眠、レム睡眠、ノンレム睡眠等）に分類する場合、多クラスロジスティック回帰モデルによって実現されてもよい。

他方、睡眠ステージ推定装置１００は、各時刻における体動データから各時刻の回転の角速度を算出して体動の情報を算出する。そして、睡眠ステージ確率と体動の情報から、学習済み睡眠ステージ遷移確率推定モデルを利用して、被検者がある睡眠ステージから他の睡眠ステージに遷移したと推定される確率を示す時系列を考慮した睡眠ステージ遷移確率を推定する。例えば、睡眠ステージ遷移確率推定モデルは、シグモイド関数によって実現されてもよい。

そして、睡眠ステージ推定装置１００は、学習済み条件付き確率場（Ｃｏｎｄｉｔｉｏｎａｌ Ｒａｎｄｏｍ Ｆｉｅｌｄ：ＣＲＦ）モデルを利用して、推定された各時刻の睡眠ステージ確率と睡眠ステージ遷移確率とから睡眠ステージ系列を推定する。例えば、条件付き確率場モデルは、睡眠ステージ確率を観測素性とし、睡眠ステージ遷移確率を遷移素性とする素性関数によって実現される。

このようにして、睡眠ステージ推定装置１００は、被検者に装着されたセンサによって睡眠中の被検者から時系列に収集された拍動データ及び体動データから、被検者の睡眠ステージ系列を推定することができる。

［ハードウェア構成］
ここで、睡眠ステージ推定装置１００は、例えば、サーバなどの計算装置であってもよく、図２に示されるようなハードウェア構成を有してもよい。すなわち、睡眠ステージ推定装置１００は、バスＢを介し相互接続されるドライブ装置１０１、補助記憶装置１０２、メモリ装置１０３、ＣＰＵ（Ｃｅｎｔｒａｌ Ｐｒｏｃｅｓｓｉｎｇ Ｕｎｉｔ）１０４、インタフェース装置１０５及び通信装置１０６を有する。

睡眠ステージ推定装置１００における後述される各種機能及び処理を実現するプログラム又は命令を含む各種コンピュータプログラムは、ＣＤ－ＲＯＭ（Ｃｏｍｐａｃｔ Ｄｉｓｋ－Ｒｅａｄ Ｏｎｌｙ Ｍｅｍｏｒｙ）などの記録媒体１０７によって提供されてもよい。プログラムを記憶した記録媒体１０７がドライブ装置１０１にセットされると、プログラムが記録媒体１０７からドライブ装置１０１を介して補助記憶装置１０２にインストールされる。但し、プログラムのインストールは必ずしも記録媒体１０７により行う必要はなく、ネットワークなどを介し何れかの外部装置からダウンロードするようにしてもよい。補助記憶装置１０２は、インストールされたプログラムを格納すると共に、必要なファイルやデータなどを格納する。メモリ装置１０３は、プログラムの起動指示があった場合に、補助記憶装置１０２からプログラムやデータを読み出して格納する。補助記憶装置１０２及びメモリ装置１０３は、プログラム又は命令を格納する非一時的なコンピュータ可読記憶媒体として実現される。プロセッサとして機能するＣＰＵ１０４は、メモリ装置１０３に格納されたプログラムやプログラムを実行するのに必要なパラメータなどの各種データに従って、睡眠ステージ推定装置１００の各種機能及び処理を実行する。インタフェース装置１０５は、ネットワーク又は外部装置に接続するための通信インタフェースとして用いられる。通信装置１０６は、外部装置と通信するための各種通信処理を実行する。

しかしながら、睡眠ステージ推定装置１００は、上述したハードウェア構成に限定されるものでなく、例えば、睡眠ステージ推定装置１００による機能及び処理の１つ以上を実現する１つ以上の回路などの他の何れか適切なハードウェア構成により実現されてもよく、例えば、被験者の腕に装着される腕時計型のウェアラブルタイプや、被験者の耳に挿入するイヤホン型のヒアラブルタイプの装置として実現されてもよい。

［睡眠ステージ推定装置］
次に、図３を参照して、本開示の一実施例による睡眠ステージ推定装置１００を説明する。図３は、本開示の一実施例による睡眠ステージ推定装置１００の機能構成を示すブロック図である。

図３に示されるように、睡眠ステージ推定装置１００は、被検者データ取得部１１０、睡眠ステージ確率推定部１２０、睡眠ステージ遷移確率推定部１３０及び睡眠ステージ推定部１４０を有する。

被検者データ取得部１１０は、被検者の拍動データ及び体動データを取得する。例えば、拍動データは、睡眠中の被検者の腕や耳などに装着された脈波センサ又は心拍センサによって収集されてもよく、体動データは、拍動データと同じく、睡眠中の被検者の腕や耳などに装着された体動を検知する加速度センサ又は体動センサなどによって収集されてもよい。

睡眠ステージ推定装置１００は、センサに有線又は無線接続され、被検者データ取得部１１０は、収集された拍動データ及び体動データを通信回線を介し取得してもよい。あるいは、収集された拍動データ及び体動データは、センサに有線又は無線接続された計算装置（図示せず）に送信及び格納され、被検者データ取得部１１０は、当該計算装置から拍動データ及び体動データを取得してもよい。被検者データ取得部１１０は、取得した拍動データ及び／又は体動データを睡眠ステージ確率推定部１２０及び睡眠ステージ遷移確率推定部１３０に提供する。

睡眠ステージ確率推定部１２０は、拍動データから特徴量系列を算出し、学習済み睡眠ステージ確率推定モデルを利用して、算出した特徴量系列から被検者の睡眠ステージ確率系列を推定する。具体的には、被検者データ取得部１１０から拍動データを取得すると、睡眠ステージ確率推定部１２０は、拍動データの心拍間隔系列から所定時間毎（例えば、３０秒毎、１分毎、５分毎など）に心拍間隔の平均、標準偏差、ＬＦ／ＨＦ比などの何れかを示す特徴量ｆ_ｔｉ（ｉ＝０，１，・・・，ｎ）を算出し、算出した各特徴量ｆ_ｔｉから特徴量系列ｆ_ｔ０，ｆ_ｔ１，・・・，ｆ_ｔｎを生成する。そして、睡眠ステージ確率推定部１２０は、各特徴量ｆ_ｔｉを学習済み睡眠ステージ確率推定モデルに入力し、時刻ｔｉの特徴量ｆ_ｔｉに対する各睡眠ステージの睡眠ステージ確率Ｐ_{ｓｔａｇｅ＿ｔｉ}を推定し、時刻ｔｉの特徴量ｆ_ｔｉに対する確率Ｐ_{ｓｔａｇｅ＿ｔｉ}が最大となる睡眠ステージ確率Ｐ_{ｍａｘ＿ｓｔａｇｅ＿ｔｉ}の睡眠ステージ確率系列Ｐ_{ｍａｘ＿ｓｔａｇｅ＿ｔ０}，Ｐ_{ｍａｘ＿ｓｔａｇｅ＿ｔ１}，・・・，Ｐ_{ｍａｘ＿ｓｔａｇｅ＿ｔｎ}を睡眠ステージ推定部１４０に提供する。

一実施例では、学習済み睡眠ステージ確率推定モデルは、特徴量系列の特徴量から各睡眠ステージの確率を推定するロジスティック回帰モデルから構成されてもよい。ロジスティック回帰モデルは、例えば、睡眠ステージ集合Ｓが｛ｓ_０（＝浅い睡眠），ｓ_１（＝深い睡眠），ｓ_２（＝レム睡眠），ｓ_３（＝ノンレム睡眠）｝である場合、多クラスロジスティック回帰モデルが利用されてもよい。すなわち、各時刻の心拍間隔の平均、標準偏差、ＬＦ／ＨＦ比などの特徴量と、当該時刻の脳波データから決定された正解の睡眠ステージとのペアから構成される訓練データが予め用意される。何れか適切なロジスティック曲線などの多クラスロジスティック回帰モデルに対して訓練データを利用した教師有り学習が実行される。多クラスロジスティック回帰モデルは、多クラスロジスティック回帰モデルの学習に関する何れか周知の学習手法に従って学習可能であり、訓練データの特徴量に対して正解の睡眠ステージが正しく予測できるように、多クラスロジスティック回帰モデルの重みベクトルなどのパラメータが更新される。例えば、訓練データと正解データの確率の交差エントロピーを求め、その値が最小となるパラメータを求めてもよい。具体的には、確率的勾配降下法などの手法を用いて、交差エントロピーの各パラメータの偏微分値を使ってパラメータを繰り返し更新し、交差エントロピーが最小となるパラメータを求めてもよい。なお、実際上は、何れか適切なライブラリの関数を利用して学習可能である。

睡眠ステージ遷移確率推定部１３０は、体動データから体動量系列を算出し、学習済み睡眠ステージ遷移確率推定モデルを利用して、算出した体動量系列から被検者の睡眠ステージ遷移確率系列を推定する。具体的には、被検者データ取得部１１０から体動データを取得すると、睡眠ステージ遷移確率推定部１３０は、体動データから所定時間毎（例えば、３０秒毎、１分毎、５分毎など）の角速度などの体動量ｍ_ｔｉ（ｉ＝０，１，・・・，ｎ）を算出し、算出した各体動量ｍ_ｔｉから体動量系列ｍ_ｔ０，ｍ_ｔ１，・・・，ｍ_ｔｎを生成する。そして、睡眠ステージ遷移確率推定部１３０は、各体動量ｍ_ｔｉを学習済み睡眠ステージ遷移確率推定モデルに入力し、時刻ｔｉの体動量ｍ_ｔｉに対する睡眠ステージ遷移確率Ｐ_{ｔｒａｎｓｉｔｉｏｎ＿ｔｉ}を推定し、推定した睡眠ステージ遷移確率系列Ｐ_{ｔｒａｎｓｉｔｉｏｎ＿ｔ０}，Ｐ_{ｔｒａｎｓｉｔｉｏｎ＿ｔ１}，・・・，Ｐ_{ｔｒａｎｓｉｔｉｏｎ＿ｔｎ}を睡眠ステージ推定部１４０に提供する。

一実施例では、学習済み睡眠ステージ遷移確率推定モデルは、体動量系列の体動量から睡眠ステージ間の遷移確率を推定するシグモイド関数から構成されてもよい。例えば、シグモイド関数は、
Ｐｒ（ｍ_ｔｉ）＝１／（１＋ｅｘｐ（－ａｘ＋ｂ））
として定式化されてもよい。なお、パラメータａ，ｂは、後述される条件付き確率場（ＣＲＦ）モデルの学習時に決定されてもよい。なお、睡眠ステージ遷移確率推定モデルのパラメータは、以下の条件付き確率場（ＣＲＦ）の学習と同時に行ってもよく、条件付き確率場（ＣＲＦ）のパラメータの一部として学習してもよい。

睡眠ステージ推定部１４０は、学習済み条件付き確率場モデルを利用して、睡眠ステージ確率系列及び睡眠ステージ遷移確率系列から被検者の睡眠ステージ系列を推定する。具体的には、睡眠ステージ確率推定部１２０及び睡眠ステージ遷移確率推定部１３０から睡眠ステージ確率系列Ｐ_{ｍａｘ＿ｓｔａｇｅ＿ｔ０}，Ｐ_{ｍａｘ＿ｓｔａｇｅ＿ｔ１}，・・・，Ｐ_{ｍａｘ＿ｓｔａｇｅ＿ｔｎ}及び睡眠ステージ遷移確率系列Ｐ_{ｔｒａｎｓｉｔｉｏｎ＿ｔ０}，Ｐ_{ｔｒａｎｓｉｔｉｏｎ＿ｔ１}，・・・，Ｐ_{ｔｒａｎｓｉｔｉｏｎ＿ｔｎ}をそれぞれ取得すると、睡眠ステージ推定部１４０は、各時点の睡眠ステージ確率Ｐ_{ｍａｘ＿ｓｔａｇｅ＿ｔｉ}及び睡眠ステージ遷移確率Ｐ_{ｔｒａｎｓｉｔｉｏｎ＿ｔｉ}をそれぞれ観測素性及び遷移素性として学習済み条件付き確率場（ＣＲＦ）モデルに入力し、各時点の睡眠ステージＳ_ｔｉを推定し、睡眠ステージ系列Ｓ_ｔ０，Ｓ_ｔ１，・・・，Ｓ_ｔｎを生成する。

一実施例では、学習済み条件付き確率場モデルは、睡眠ステージ確率を観測素性とし、睡眠ステージ遷移確率を遷移素性とする素性関数から構成されてもよい。睡眠ステージ推定部１４０は、当該素性関数を利用して、各時点の睡眠ステージ確率Ｐ_{ｍａｘ＿ｓｔａｇｅ＿ｔｉ}及び睡眠ステージ遷移確率Ｐ_{ｔｒａｎｓｉｔｉｏｎ＿ｔｉ}から睡眠ステージＳ_ｔｉを推定する。例えば、学習系列データと正解ラベル列の確率の交差エントロピーを求め、その値が最小となるパラメータを求めてもよい。具体的には、確率的勾配降下法や動的計画法などの手法を用いて、交差エントロピーの各パラメータの偏微分値を使ってパラメータを繰り返し更新し、交差エントロピーが最小となるパラメータを求めてもよい。なお、実際上は、何れか適切なライブラリの関数を利用して学習可能である。

具体的には、λを条件付き確率場のパラメータとし、ａ，ｂを睡眠ステージ遷移確率のパラメータとした場合、交差エントロピーＥ（λ，ａ，ｂ）を最小化するパラメータλ，ａ，ｂが同時に学習される。すなわち、交差エントロピーＥを各パラメータλ，ａ，ｂで偏微分した値を用いて、確率的勾配降下法などの手法によってパラメータを繰り返し更新し、交差エントロピーＥが最小となるパラメータλ，ａ，ｂを求める。

例えば、Ｎ個の学習データ

に対して、交差エントロピーＥは、

である。ただし、対数尤度ｌ（ｘ，ｙ）は、時系列の長さをＭとして、

であり、

である。

ここで、素性関数

は、例えば、睡眠ステージ確率（Ｋ個）、定数１、睡眠ステージ遷移確率ｐ_ｒで構成され、

である。ただし、体動量ｑの睡眠ステージ遷移確率ｐ_ｒは、

である。
素性関数の例として、ｘ_ｍ＝［ｐ_１，ｐ_２，ｐ_３，ｐ_４，ｑ］，ｙ_ｍ－１＝３，ｙ_ｍ＝２である場合、

である。

［睡眠ステージ推定処理］
次に、図４を参照して、本開示の一実施例による睡眠ステージ推定処理を説明する。当該睡眠ステージ推定処理は、上述した睡眠ステージ推定装置１００によって実現され、例えば、睡眠ステージ推定装置１００のプロセッサがプログラム又は命令を実行することによって実現されてもよい。図４は、本開示の一実施例による睡眠ステージ推定処理を示すフローチャートである。

図４に示されるように、ステップＳ１０１において、睡眠ステージ推定装置１００は、拍動データ及び体動データを取得する。拍動データ及び体動データは上述した脈波センサ又は心拍センサ、及び加速度センサ又は体動センサにより収集される。

ステップＳ１０２において、睡眠ステージ推定装置１００は、特徴量系列から睡眠ステージ確率系列を推定する。具体的には、睡眠ステージ推定装置１００は、所定時間毎の心拍間隔の平均、標準偏差、ＬＦ／ＨＦ比などの特徴量を算出し、算出した特徴量を学習済み睡眠ステージ確率推定モデルに入力する。睡眠ステージ推定装置１００は、各時刻について学習済み睡眠ステージ確率推定モデルから睡眠ステージ確率を取得し、睡眠ステージ確率系列を生成する。

ステップＳ１０３において、睡眠ステージ推定装置１００は、体動量系列から睡眠ステージ遷移確率系列を推定する。具体的には、睡眠ステージ推定装置１００は、所定時間毎の角速度などの体動量を算出し、算出した体動量を学習済み睡眠ステージ遷移確率推定モデルに入力する。睡眠ステージ推定装置１００は、各時刻について学習済み睡眠ステージ遷移確率推定モデルから睡眠ステージ遷移確率を取得し、睡眠ステージ遷移確率系列を生成する。

ステップＳ１０４において、睡眠ステージ推定装置１００は、睡眠ステージ系列を推定する。具体的には、睡眠ステージ推定装置１００は、ステップＳ１０２において生成された睡眠ステージ確率系列と、ステップＳ１０３において生成された睡眠ステージ遷移確率系列とから、各時刻における睡眠ステージ確率及び睡眠ステージ遷移確率をそれぞれ観測素性及び遷移素性として学習済み睡眠ステージ推定モデルに入力し、各時刻における睡眠ステージを推定する。全ての時刻について睡眠ステージを推定すると、睡眠ステージ推定装置１００は、睡眠ステージ系列を生成する。

本願発明は、上記実施形態に限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で種々に変形することが可能である。また、各実施形態は可能な限り適宜組み合わせて実施してもよく、その場合組み合わせた効果が得られる。さらに、上記実施形態には種々の段階の発明が含まれており、開示される複数の構成要件における適当な組み合わせにより種々の発明が抽出され得る。例えば、実施形態に示される全構成要件からいくつかの構成要件が削除されても、発明が解決しようとする課題の欄で述べた課題が解決でき、発明の効果の欄で述べられている効果が得られる場合には、この構成要件が削除された構成が発明として抽出され得る。

以下に、本願出願の当初の特許請求の範囲に記載された発明を付記する。
［付記］
本開示の一態様では、被検者の拍動データ及び体動データを取得する被検者データ取得部と、
前記拍動データから特徴量系列を算出し、学習済み睡眠ステージ確率推定モデルを利用して、前記算出した特徴量系列から前記被検者の睡眠ステージ確率系列を推定する睡眠ステージ確率推定部と、
前記体動データから体動量系列を算出し、学習済み睡眠ステージ遷移確率推定モデルを利用して、前記算出した体動量系列から前記被検者の睡眠ステージ遷移確率系列を推定する睡眠ステージ遷移確率推定部と、
学習済み条件付き確率場モデルを利用して、前記睡眠ステージ確率系列及び前記睡眠ステージ遷移確率系列から前記被検者の睡眠ステージ系列を推定する睡眠ステージ推定部と、
を有する睡眠ステージ推定装置が提供される。

一実施例では、前記特徴量系列は、所定時間毎の心拍間隔の平均、標準偏差及びＬＦ／ＨＦ比の１つ以上の時系列データを含んでもよい。

一実施例では、前記学習済み睡眠ステージ確率推定モデルは、前記特徴量系列の特徴量から各睡眠ステージの確率を推定するロジスティック回帰モデルから構成されてもよい。

一実施例では、前記学習済み睡眠ステージ確率推定モデルは、前記拍動データの特徴量と、対応する正解の睡眠ステージとのペアから構成される訓練データを利用した教師有り学習によって取得されてもよい。

一実施例では、前記体動量系列は、角速度の時系列データを含んでもよい。

一実施例では、前記学習済み睡眠ステージ遷移確率推定モデルは、前記体動量系列の体動量から睡眠ステージ間の遷移確率を推定するシグモイド関数から構成されてもよい。

一実施例では、前記学習済み条件付き確率場モデルは、睡眠ステージ確率を観測素性とし、睡眠ステージ遷移確率を遷移素性とする素性関数から構成されてもよい。

本開示の他の態様では、
被検者の拍動データ及び体動データを取得するステップと、
前記拍動データから特徴量系列を算出し、学習済み睡眠ステージ確率推定モデルを利用して、前記算出した特徴量系列から前記被検者の睡眠ステージ確率系列を推定するステップと、
前記体動データから体動量系列を算出し、学習済み睡眠ステージ遷移確率推定モデルを利用して、前記算出した体動量系列から前記被検者の睡眠ステージ遷移確率系列を推定するステップと、
学習済み条件付き確率場モデルを利用して、前記睡眠ステージ確率系列及び前記睡眠ステージ遷移確率系列から前記被検者の睡眠ステージ系列を推定するステップと、
を有する睡眠ステージ推定方法が提供される。

本開示の他の態様では、
睡眠ステージ推定装置のコンピュータが実行するプログラムであって、前記コンピュータに、
被検者の拍動データ及び体動データを取得する処理と、
前記拍動データから特徴量系列を算出し、学習済み睡眠ステージ確率推定モデルを利用して、前記算出した特徴量系列から前記被検者の睡眠ステージ確率系列を推定する処理と、
前記体動データから体動量系列を算出し、学習済み睡眠ステージ遷移確率推定モデルを利用して、前記算出した体動量系列から前記被検者の睡眠ステージ遷移確率系列を推定する処理と、
学習済み条件付き確率場モデルを利用して、前記睡眠ステージ確率系列及び前記睡眠ステージ遷移確率系列から前記被検者の睡眠ステージ系列を推定する処理と、
を実行させるプログラムが提供される。

１００ 睡眠ステージ推定装置
１１０ 被検者データ取得部
１２０ 睡眠ステージ確率推定部
１３０ 睡眠ステージ遷移確率推定部
１４０ 睡眠ステージ推定部

Claims (9)

1. 被検者の拍動データ及び体動データを取得する被検者データ取得部と、
   前記拍動データから特徴量系列を算出し、学習済み睡眠ステージ確率推定モデルを利用して、前記算出した特徴量系列から前記被検者の睡眠ステージ確率系列を推定する睡眠ステージ確率推定部と、
   前記体動データから体動量系列を算出し、学習済み睡眠ステージ遷移確率推定モデルを利用して、前記算出した体動量系列から前記被検者の睡眠ステージ遷移確率系列を推定する睡眠ステージ遷移確率推定部と、
   前記睡眠ステージ確率系列を観測素性とし、前記睡眠ステージ遷移確率系列を遷移素性とする素性関数から構成される学習済み条件付き確率場モデルを利用して、前記睡眠ステージ確率系列及び前記睡眠ステージ遷移確率系列から前記被検者の睡眠ステージ系列を推定する睡眠ステージ推定部と、を有し、
   前記学習済み条件付き確率場モデルは、前記睡眠ステージ確率系列または前記睡眠ステージ遷移確率系列の学習系列データと正解ラベル列の確率の交差エントロピーを最小化するパラメータを求めることにより学習される、
   睡眠ステージ推定装置。
2. 前記特徴量系列は、所定時間毎の心拍間隔の平均、標準偏差及びＬＦ／ＨＦ比の１つ以上の時系列データを含む、請求項１記載の睡眠ステージ推定装置。
3. 前記学習済み睡眠ステージ確率推定モデルは、前記特徴量系列の特徴量から各睡眠ステージの確率を推定するロジスティック回帰モデルから構成される、請求項１又は２記載の睡眠ステージ推定装置。
4. 前記学習済み睡眠ステージ確率推定モデルは、前記拍動データの特徴量と、対応する正解の睡眠ステージとのペアから構成される訓練データを利用した教師有り学習によって取得される、請求項１乃至３何れか一項記載の睡眠ステージ推定装置。
5. 前記体動量系列は、角速度の時系列データを含む、請求項１乃至４何れか一項記載の睡眠ステージ推定装置。
6. 前記学習済み睡眠ステージ遷移確率推定モデルは、前記体動量系列の体動量から睡眠ステージ間の遷移確率を推定するシグモイド関数から構成される、請求項１乃至５何れか一項記載の睡眠ステージ推定装置。
7. 前記学習済み条件付き確率場モデルは、前記交差エントロピーを前記パラメータで偏微分した値を用いて、確率的勾配降下法または動的計画法によって前記パラメータを繰り返し更新することにより、前記最小化するパラメータを求める、請求項１乃至６何れか一項記載の睡眠ステージ推定装置。
8. 睡眠ステージ推定装置における睡眠ステージ推定方法であって、
   被検者の拍動データ及び体動データを取得するステップと、
   プロセッサが、前記拍動データから特徴量系列を算出し、学習済み睡眠ステージ確率推定モデルを利用して、前記算出した特徴量系列から前記被検者の睡眠ステージ確率系列を推定するステップと、
   前記プロセッサが、前記体動データから体動量系列を算出し、学習済み睡眠ステージ遷移確率推定モデルを利用して、前記算出した体動量系列から前記被検者の睡眠ステージ遷移確率系列を推定するステップと、
   前記プロセッサが、前記睡眠ステージ確率系列を観測素性とし、前記睡眠ステージ遷移確率系列を遷移素性とする素性関数から構成される学習済み条件付き確率場モデルを利用して、前記睡眠ステージ確率系列及び前記睡眠ステージ遷移確率系列から前記被検者の睡眠ステージ系列を推定するステップと、を有し、
   前記学習済み条件付き確率場モデルは、前記睡眠ステージ確率系列または前記睡眠ステージ遷移確率系列の学習系列データと正解ラベル列の確率の交差エントロピーを最小化するパラメータを求めることにより学習される、
   睡眠ステージ推定方法。
9. 睡眠ステージ推定装置のコンピュータが実行するプログラムであって、前記コンピュータに、
   被検者の拍動データ及び体動データを取得する処理と、
   前記拍動データから特徴量系列を算出し、学習済み睡眠ステージ確率推定モデルを利用して、前記算出した特徴量系列から前記被検者の睡眠ステージ確率系列を推定する処理と、
   前記体動データから体動量系列を算出し、学習済み睡眠ステージ遷移確率推定モデルを利用して、前記算出した体動量系列から前記被検者の睡眠ステージ遷移確率系列を推定する処理と、
   前記睡眠ステージ確率系列を観測素性とし、前記睡眠ステージ遷移確率系列を遷移素性とする素性関数から構成される学習済み条件付き確率場モデルを利用して、前記睡眠ステージ確率系列及び前記睡眠ステージ遷移確率系列から前記被検者の睡眠ステージ系列を推定する処理として実行させ、
   前記学習済み条件付き確率場モデルは、前記睡眠ステージ確率系列または前記睡眠ステージ遷移確率系列の学習系列データと正解ラベル列の確率の交差エントロピーを最小化するパラメータを求めることにより学習される、
   プログラム。

Images