JP6925057B2 - 睡眠段階判定方法、睡眠段階判定装置、及び睡眠段階判定プログラム - Google Patents

Description

本発明は、睡眠段階の判定技術に関する。

睡眠段階判定の国際基準であるRechtschaffen & Kales（以下、「Ｒ＆Ｋ法」と略称する）法では、図１に示すように睡眠状態が６段階に分類されている。睡眠の浅い方から順に覚醒段階（「Ｗ」と表記）、レム（ＲＥＭ：Rapid Eye Movement）睡眠段階（「Ｒ」と表記）、ノンレム睡眠の段階１〜４（「１」、「２」、「３」、「４」で標記）に分類し、ノンレム段階４を最も深い睡眠としている。睡眠段階の推定は、認知症等の医療診断のみならず、介護の現場でも有用である。たとえば、レム睡眠の周期を図ることで、高齢者の行動管理や、睡眠障害等の発現を予測することができる。

図１の睡眠段階データは、診断対象者に特殊な器具を装着させ、脳波（ＥＥＧ：Electroencephalogram）や眼球運動（ＥＯＧ：Electro-oculogram）等を測定して専門医の経験に基づいて作成されたものである。この睡眠段階データは、主として精神疾患等の医療判定に用いられる精密なデータである。睡眠段階の判定を介護の現場に適用する場合、高齢者や身障者に電極等の器具を一晩中装着させてデータを取ることは非現実的である。また、被介護者の睡眠障害の発現を予測するわけではないので、精密な睡眠分析までは要しないが、睡眠段階はなるべく正確に推定されるのが望ましい。

一方、無拘束で身体データを取得して睡眠段階を判定する手法が提案されている（たとえば、非特許文献１参照）。この方法は、圧電フィルム等を適用した無拘束のエアマットレス型センサを用いて、心拍、体動等の身体データを計測し、得られたデータから睡眠段階を推定する。被験者に電極等の器具を装着させないので、被験者にストレスを与えずに身体データを取得することができる。

図２は、従来の無拘束型の睡眠段階判定方法を示す図である。無拘束型センサで得られた身体情報から心拍数データを取り出して睡眠段階を判定する。心拍数が低いほど睡眠が深く、心拍数が高いほど覚醒に近づくという一般的見地に基づくものである。この従来方法は、心拍数データをＦＦＴ（高速フーリエ変換）で周波数領域のデータに変換し、バンドパスフィルタで、パワースペクトルのピーク値の周期の半分から１３５分までの周期の中周波成分を取り出す。取り出した中周波成分にＩＦＦＴ（逆フーリエ変換）を施して時間領域の心拍数データに戻した後に、６段階に離散化する。図２に示す方法を「従来手法」と称する。

また、心拍数から睡眠段階を推定するために、生体データから適切な周波数成分のデータを取り出すためのフィルタが提案されている（たとえば、特許文献１参照）。

特開２０１４−２３９７８９号公報

T. Watanabe and K. Watanabe, "Non-contact Method for Sleep Stage Estimation" IEEE Transaction on Biomedical Engineering, No. 10, Vol. 51, pp. 1735-1748, 2004

図２の従来方法では睡眠の大きなリズム（周期）を把握することができるが、突然の覚醒状態への遷移や、レム睡眠の特徴である心拍の細かい振動を判断することができない。

また、睡眠の周期（心拍の周期）自体が変動するにもかかわらず、平均的な周期に固定して睡眠段階を推定しているので、実際の睡眠段階からのズレが生じる可能性がある。

睡眠段階判定を高齢者等の介護の現場に適用する場合には、無拘束で取得したデータからレム睡眠や覚醒状態への遷移をできるだけ正確に判定できることが望ましい。

そこで、本発明は取得した身体データから、簡便かつより高い精度で睡眠段階を推定することのできる手法と構成を提供することを課題とする。

上記課題を解決するために、本発明では、睡眠段階に応じてその睡眠段階の推定に最適な身体データを用いて睡眠段階を推定することで、推定精度を向上する。

本発明の第１の側面では、睡眠段階判定装置は、
被験者の心拍データを取得するデータ取得部と、
前記心拍データからレム睡眠段階を推定する判定部であって、前記心拍データを周波数領域の成分に変換し、前記周波数領域のデータから、レム睡眠に関係のある２．５秒〜１５０分の周期にわたる範囲の周期成分を抽出して時間領域のデータに変換し、前記時間領域のデータのうち、所定の閾値を超える領域をレム睡眠段階と判定する判定部、とを有する。

本発明の第２の側面では、睡眠段階判定装置は、
被験者の体動データを取得するデータ取得部と、
前記体動データの全体平均と標準偏差を算出し、前記標準偏差に所定の係数を乗算した値を前記全体平均に足し合わせた値を閾値として設定し、前記閾値を超える領域を覚醒段階と判定する判定部と、
を有する。

本発明の第３の側面では、睡眠段階判定装置は、
被験者の呼吸データを取得するデータ取得部と、
前記呼吸データに基づいて、単位時間当たりの呼吸数の平均値と標準偏差を算出し、前記標準偏差を前記平均値で除算した変動係数に所定の係数を乗算した値を閾値として設定し、前記閾値以下の領域をノンレム睡眠段階と判定する判定部と、
を有する。

本発明の第４の側面では、睡眠段階判定装置は、
被験者の身体データを取得するデータ取得部と、
前記身体データから心拍データと、呼吸データと、体動データをそれぞれ分離し、前記呼吸データを用いたノンレム睡眠段階の判定と、前記心拍データを用いたレム睡眠段階の判定と、前記体動データを用いた覚醒段階の判定を個別または順次に行って睡眠段階を判定する判定部と、
を有する。

測定した身体データから、簡便かつより高い精度で睡眠段階を推定することができる。また、睡眠の段階ごとに個別の推定が可能になる。

Ｒ＆Ｋ国際分類による睡眠段階を示す図である。 心拍数データを用いた従来の睡眠段階判定手法を示す図である。 本発明の睡眠段階判定法の基本概念を説明する図である。 実施形態の睡眠段階判定方法のフローチャートである。 ノンレム睡眠段階の判定を説明する図である。 レム睡眠段階の判定を説明する図である。 覚醒段階の判定を説明する図である。 各段階の判定をトータルした判定結果を説明する図である。 実施形態のレム判定結果を従来手法と比較して示す図である。 実施形態の覚醒判定結果を従来手法と比較して示す図である。 実施形態のトータルの判定結果とＲ＆Ｋ法による判定との一致度を示す図である。 実施形態の睡眠段階判定結果とＲ＆Ｋ法との一致率を、従来手法とＲ＆Ｋ法との一致率と比較して示す図である。 実施形態の睡眠段階判定結果とＲ＆Ｋ法との一致率を、従来手法とＲ＆Ｋ法との一致率と比較して示す図である。 実施形態の睡眠段階判定装置の概略構成図である。

図３は、本発明の睡眠段階判定の基本概念を説明する図である。実施形態では、睡眠の段階に応じて、その睡眠段階の推定に最も適した身体データを用いることで推定精度を向上する。具体的には、覚醒段階の推定には体動データを、レム睡眠段階の推定には心拍データを、ノンレム睡眠段階の推定には呼吸データを用いて、各睡眠段階を判定する。これは以下の知見に基づく。

心拍数は、レム睡眠段階ではノンレム睡眠段階と比べて同等かやや高く、覚醒段階よりも低い。心拍の揺らぎは、レム睡眠段階ではノンレム睡眠段階と比較して極めて大きくなり、覚醒段階と同等かやや大きい。すなわち、レム睡眠では心拍数は覚醒時と同じ程度にまで大きくなり、かつ心拍の変動が最も顕著になる。そこで、レム睡眠段階の判定に、心拍数と心拍の揺らぎを用いる。

体動は、睡眠が深いほど小さく、睡眠が浅いほど大きくかつ頻繁になる。そこで、覚醒段階の判定に体動データを用いる。

呼吸は、急速な眼球運動を伴わない深い睡眠（ノンレム睡眠）では安定するため，呼吸の標準偏差（分散）が小さくなる。そこで、呼吸データを用いてノンレム睡眠段階を推定する。

このようなアプローチにより、個々の睡眠段階を個別に推定することもできるし、任意の睡眠段階の推定を組み合わせた総合的な推定も可能になる。また、本発明の測定は拘束型、無拘束型の測定を問わず、脳波情報等を用いるＲ＆Ｋ法の推定結果に近づけることができる。

図４は、実施形態の睡眠段階判定方法のフローチャートである。まず、生体センサで、被験者の心拍データ、体動データ、及び呼吸データを含む身体データを取得する（Ｓ１０１）。取得した身体データから呼吸データ、心拍データ、及び体動データを分離する（Ｓ１０２）。生体センサの出力には多数の種類の波形データが含まれ、たとえば周波数分離等の手法により各データに分離する。

次に、各睡眠段階を判定する。たとえば、呼吸データに１つ以上の第１閾値を設定し、閾値比較によりノンレム睡眠段階を判定する（Ｓ１０３）。呼吸データに基づくノンレム睡眠の推定の詳細は後述する。また、心拍データに第２閾値を設定し、閾値比較によりレム睡眠段階を判定する（Ｓ１０４）。さらに、体動データに第３閾値を設定し、閾値判定により覚醒段階を判定する（Ｓ１０５）。ステップＳ１０３〜Ｓ１０５はこの順で順次実行してもよいし、また、覚醒段階の判定、レム睡眠段階の判定、ノンレム睡眠段階の判定の順に優先して判定結果に反映するのであれば、判定自体は個別に並列または順不同に実施してもよい。例えば、覚醒段階の判定とレム睡眠段階の判定とが競合する場合に、覚醒段階の判定をレム睡眠段階の判定に優先して判定結果に反映するのであれば、判定順を問わず同一の結果を得ることができる。同様に、覚醒段階の判定とノンレム睡眠段階の判定の場合は、覚醒段階の判定をノンレム睡眠段階の判定に優先して反映する。レム睡眠段階の判定とノンレム睡眠段階の判定が競合する場合は、レム睡眠段階の判定をノンレム睡眠段階の判定に優先して反映する。これにより、判定の順序を問わず同一の結果を得ることができる。

判定結果をデータ出力する（Ｓ１０６）。データの出力形式としては、図１のような時系列に記述された全睡眠段階のデータを出力してもよいし、覚醒判定がなされたときに優先的にアラーム等の通知を出力してもよい。また、レム判定がなされたときに、覚醒判定と異なる種類のアラームを発生させてもよい。

図５は、ノンレム睡眠段階の判定を説明する図である。実施形態では、呼吸の分散（または標準偏差）を平均値で除算した値を変動係数として用いてノンレム睡眠段階を推定する（変動係数＝標準偏差／平均）。まず、生体センサで取得した身体データから、呼吸データを分離する（手順１）。分離した呼吸データから、たとえば５分ごとに、単位時間あたりの呼吸数の標準偏差と平均値を計算して、５分間変動係数を求める。変動係数は平均呼吸数からのズレを示す。標準偏差を平均値で除算することで、被験者間のばらつきを低減する。なお、変動係数を求めるための時間は５分に限定されない。変動係数の算出のために１００回程度の呼吸数があればよいので、被験者に応じて４分ごとに変動係数を算出してもよいし、８分ごとに算出してもよい。

次に、変動係数を用いて、ノンレム判定用の第１閾値を少なくともひとつ設定する。図２の例では、ノンレム判定用に、ＴＨ１−ａとＴＨ１−ｂの２つの閾値を設定する。ＴＨ１−ａは、変動係数の全体平均に０．６を乗算した値であり、ＴＨ１−ｂは変動係数の全体平均に０．２を乗算した値である。閾値の設定のために変動係数の全体平均に乗算される係数０．６と０．２は実験的に得られた値であり、正答率を最も高くする値である。

５分間変動係数がＴＨ１−ａを超える領域をノンレム睡眠段階１（浅い睡眠）と判定し、５分間変動係数がＴＨ１−ａ以下の領域をノンレム睡眠段階２〜４（深い睡眠）と判定する。ＴＨ１−ａを超える領域には、ノンレム睡眠段階１の他に、レム睡眠段階と覚醒段階も含まれているが、図５の処理はノンレム睡眠段階の推定を適正に行う処理なので、ノンレム睡眠段階１以上をノンレム１として扱う。

２つ目の閾値ＴＨ１−ｂを用いて、ノンレム睡眠段階２とノンレム睡眠段階３以下を区別してもよい。ＴＨ１−ｂを超え、かつＴＨ１−ａ以下の領域をノンレム睡眠段階２と推定する。ＴＨ１−ｂ以下の領域を、ノンレム睡眠段階３または４と推定する。

この方法により、少なくとも被験者がノンレム睡眠段階１以上の浅い睡眠状態にあるのか、ノンレム睡眠段階２以下の深い睡眠状態にあるのかを、精度よく推定することができる。また、ノンレム判定用の閾値を２以上設定することで、被験者がノンレム睡眠段階２にあるのか、３以下にあるのかを区別することができる。ノンレム睡眠段階３と４はともに深い眠りであり、これら２つを同じ段階として取り扱う分類法もある。したがって、「ノンレム睡眠段階３以下」という推定でも十分であるが、必要に応じて、３つめの閾値、たとえば全体の変動係数に０．１を乗算したＴＨ１−ｃを設定してノンレム睡眠段階３と４を区別してもよい。

呼吸データを用いた上述の判別法では、ノンレム睡眠段階を精度良く推定することはできるが、ノンレム睡眠段階とレム睡眠段階の区別はされていない。そこで別途心拍データを用いて、レム睡眠段階とノンレム睡眠段階を区別する。レム睡眠段階は、心拍データから個別に判定することもできるが、上述したノンレム睡眠の推定結果を利用して推定する場合は、「ノンレム睡眠段階１以上」という広い枠の中からレム睡眠段階を特定するという意味で「ＲＥＭ補正」と称してもよい。

図６は、レム睡眠段階の判定を説明する図である。実施形態では、心拍の揺らぎを利用してレム睡眠段階を推定する。まず、生体センサで取得した身体データから、心拍データを分離する（手順１）。心拍データは９０分前後の大きな周期で緩やかに変化するが、大きな周期の中に細かな振動成分が多数含まれている。次に、分離した心拍データにＦＦＴを適用して周波数領域の成分に変換し、バンドバスフィルタで所定の周期成分を抽出した後に、ＩＦＦＴにより時間領域のデータに戻す（手順２）。

実施形態の特徴として、ＦＦＴで得られた周波数成分のうち、２．５秒〜１３５分までの周期成分を抽出する。この抽出法は、特許文献１の従来法や他の一般的な睡眠段階推定で２６分から１３５分の周期成分を抽出する点と異なる。従来法は、睡眠の周期自体に含まれる細かな変動を考慮していない。これに対し、実施形態では、心拍に生じる２．５秒から２５分の細かい周期の心拍の変動（揺らぎ）を考慮するため、ＦＦＴ成分の中から、２．５秒から１３５分までの周期成分を抽出してＩＦＦＴを適用する。これにより、時間領域で、心拍の揺らぎを表わす波を復元する（手順３）。周期成分の抽出は、２．５秒から１３５分までの周期に限定されない。心拍揺らぎには個人差があるので、被験者に応じて２秒から１５０分の周期の中で抽出範囲を調整してもよい。２５分以下の細かい振動周期の成分から１２０分以上の長い変化周期の成分までが含まれればよいので、２秒から１２０分までの周期成分を抽出してもよいし、５秒〜１５０分までの周期成分を抽出してもよい。

次に、ＩＦＦＴ後の波形データに対して、５分間の平均を計算し、レム睡眠段階判定用の閾値ＴＨ２を設定する（手順４）。ここで設定される閾値ＴＨ２は、ＩＦＦＴ後のデータ全体の平均に、標準偏差に０．２を乗じた値を足し算した値である（全体平均＋標準偏差×０．２）。係数０．２は実験的に得られた値であり、最も正答率を高くする値である。

閾値ＴＨ２を超える区間をレム睡眠段階と判定する（手順５）。太線で示される区間がレム睡眠段階として判定された区間である。図５で取得されたノンレム判定データを下地データとして用いる場合は、ノンレム睡眠段階１と推定された区間の中からレム睡眠段階の区間が特定される。この意味で、図６の手順５で得られる判定結果を「ＲＥＭ補正」後の睡眠段階判定結果と称してもよい。

図７は、覚醒段階の判定を説明する図である。覚醒段階は体動データに基づいて推定される。まず、生体センサで取得した身体データから、体動データを分離する（手順１）。分離した体動データから、１分ごとの平均値を算出し、閾値ＴＨ３を設定する（手順２）。閾値ＴＨ３は、一晩の睡眠の全体動データの平均に、標準偏差に０．２を乗算した値を足し算した値である（全体平均＋標準偏差×０．２）。平均値の算出区間は１分に限定されず、２分ごと、３分ごと等であってもよい。ただし、覚醒段階では体動の大きさ、頻度ともに増大するので、心拍データや呼吸データの平均化時間よりも短い平均化時間に設定するのが望ましい。

閾値ＴＨ３を超える区間を覚醒段階と判定する（手順３）。黒丸で示される区間が覚醒段階として判定された区間である。図６のレム判定されたデータを下地データとして用いる場合は、ノンレム睡眠段階１と推定された区間の中でレム睡眠段階と推定された区間も含め覚醒段階の区間が特定される。この意味で、図７の手順３で得られる判定結果を「ＷＡＫＥ補正」と称してもよい。図５のノンレム睡眠段階の判定データと、図６のレム判定されたデータの両方を下地データとして用いる場合は、図７の最下段の図のように、睡眠段階を総合的に推定したデータを生成することができる。

図８は、実施形態のノンレム判定、レム判定、及び覚醒判定の利用態様を説明する図である。図４を参照して説明したように、睡眠段階に応じて最適な身体データを用いた推定を行うので、覚醒判定、レム判定、ノンレム判定の順に判定を優先して判定結果に反映するのであれば、ノンレム睡眠段階の判定、レム睡眠段階の判定、及び覚醒段階の判定は、個別に行ってもよいし、順次行ってもよい。判定順を問わず同一の結果を得ることができる。

図８の例では、まず閾値ＴＨ１−ａとＴＨ１−ｂを用いてノンレム睡眠段階を推定する（処理（ａ））。この推定により、ノンレム睡眠段階１、レム睡眠段階、及び覚醒段階を含む浅い睡眠区間と、ノンレム睡眠段階２〜４の深い睡眠区間とが区別される。例えば、被介護者が夜中に起きて徘徊する可能性がある場合に、ＴＨ１−ａを超える区間を重点的にチェックすることで、徘徊を未然に防止することができる。また、介護施設等で排泄ケアのタイミングを、より正確に特定することもできる。

処理（ａ）の結果に基づいて、レム睡眠段階をより正確に絞り込んでもよい（処理（ｂ））。実施形態では心拍の揺らぎに基づいてレム睡眠段階を推定しているので、従来例と比較してより正確にレム睡眠段階を推定することができる。

処理（ａ）の結果に基づいて、覚醒段階を推定してもよい（処理（ｃ））。さらに、処理（ｂ）の後に処理（ｃ）を実施して、処理（ｂ）と処理（ｃ）で得られたデータを合算して、最終的に処理（ｄ）のデータを生成してもよい。

図９は、実施形態の方法によるレム判定の効果を示す図である。グレーのラインはＲ＆Ｋ法による精密な睡眠判定結果である。実施形態の狙いは、センサで取得した身体データからの睡眠段階の推定結果をできるだけＲ＆Ｋ法による推定結果に近づけることである。比較例として、特許文献１の従来法（図２参照）によるレム判定結果を点線で示す。図中、太線で示す部分が実施例のレム判定による推定結果である。

従来法では、睡眠周期に含まれる細かい変動成分が考慮されていないので、レム睡眠段階が正確に推定されていない。これに対し、実施形態の方法は心拍変化の大きな周期とともに、心拍変化の周期に含まれる細かい周波数変動成分を取り込むので、Ｒ＆Ｋ法によるレム睡眠段階の推定結果に近づけることができる。

図１０は、実施形態の方法による覚醒判定の効果を示す図である。グレーのラインはＲ＆Ｋ法による睡眠判定結果である。比較例として、特許文献１の従来法（図２参照）による覚醒判定結果を点線で示す。図中、太線で示す部分が実施例の覚醒判定（ＷＡＫＥ補正とも称する）による推定結果である。従来法では、主として心拍データを用いて睡眠段階を判定しているので、レム睡眠段階と覚醒段階が明確に区別されていない。これに対し、覚醒判定に体動データを用いる実施形態の方法では、Ｒ＆Ｋ法による判定結果との近似性が高い。

図１１は、全睡眠段階における実施形態の手法のＲ＆Ｋ法に対する一致率を示す図である。このデータは、一人の被験者に対して、Ｒ＆Ｋ法による睡眠段階の推定と、実施形態の方法による睡眠段階推定を行って、その比較結果から得られたものである。図中、点線がＲ＆Ｋ法による推定結果、実線が実施形態の手法による推定結果である。太い実線で示す部分は、Ｒ＆Ｋ法との完全一致が得られた区間である。斜め線は、実施形態の推定結果がＲ＆Ｋ法による睡眠段階の推定から１段階だけズレた区間である。

同一の被験者での比較によると、完全一致は５７．８％、レム睡眠段階の推定に着目した一致率は６０．９％である。１段階のズレを許容すると、一致率は８６．４％になる。

図１２は、２０代から７０代までの男女２６人の平均データにおけるＲ＆Ｋ法との一致率を示す図である。図９〜１１が一人の被験者についての推定データを用いて比較しているのに対し、図１２では年齢、性別が広い範囲にわたって異なる複数の被験者の平均推定データを用いている。比較例として、特許文献１の従来法とＲ＆Ｋ法との一致率を示している。睡眠段階の推定結果の完全一致率は、従来法で３２％であるのに対し、実施形態の方法では４２％に向上している。

図１３は、図１２と同じ２６人の平均データにおける睡眠段階ごとの一致率を示す。実施形態の特徴として、覚醒段階の一致率が２７％も向上している。これは覚醒段階の判定に体動データを用いているためである。また、レム睡眠段階の一致率が６％も向上している。従来からレム睡眠段階の推定は困難とされており、６％の上昇効果は極めて大きい。さらに、ノンレム段階２の一致率も７％向上している。ノンレム段階２の推定精度が向上することによって、たとえば、高齢者の深い睡眠中に排尿ケアを行うことで、睡眠を優先させた適切なケアを行うことができる。

図１４は、睡眠段階推定装置１０の概略構成図である。睡眠段階推定装置１０は、プロセッサ１１、ユーザインタフェース１６、メモリ１７、通信インタフェース１８、及びドライブ１９を有する。ユーザインタフェース１６は、キーボード、マウス、タッチパネル等の入力操作部と、モニタディスプレイ、スピーカ等の表示・出力部を含む。メモリ１７は、ＲＯＭ（Read Only Memory）、ＲＡＭ（Random Access Memory）、ハードディスクドライブ等を含む。通信インタフェース１８は、ネットワークを介して他の装置やセンサと通信を行う。ドライブ１９は、図示しないリムーバブル記憶媒体が挿入されたときにそれらを駆動する。

プロセッサ１１は、データ取得部１２と、判定部２０と、合成部２１を含み、判定部２０は、ノンレム判定部１３、レム判定部１４、及び覚醒判定部１５を含む。データ取得部１２は、通信インタフェース１８を介して外部の生体センサから入力される身体データを取得する。ノンレム判定部１３は、データ取得部１２で取得された身体データから呼吸データを分離し、所定時間ごとに呼吸数の変動係数を計算する。変動係数は、単位時間当たりの呼吸数の標準偏差を平均値で除算して算出される（変動係数＝標準偏差／平均）。さらに、変動係数の全体の平均値に所定の係数（たとえば０．６）を乗算して、ノンレム判定用の第１閾値を設定し、第１閾値以下の領域をノンレム睡眠段階２〜４と推定する。

レム判定部１４は、データ取得部１２で取得された身体データから、心拍データを分離して周波数領域のデータに変換し、２５分以下の振動周期から１２０分以上の周期にわたる周期成分（たとえば２．５秒から１３５分までの周期成分）を抽出する。抽出した周期成分を時間領域の心拍データに戻して、所定の時間ごとに平均化する。次に、全体の平均値と標準偏差を求め、標準偏差に所定の係数（たとえば０．２）を乗算した値を全体の平均値に足し算して、レム判定用の第２閾値を設定する。第２閾値を超える領域をレム睡眠段階と推定する。

覚醒判定部１５は、データ取得部１２で取得された身体データから、体動データを分離し、所定時間ごと（たとえば１分ごと）に体動を平均化する。次に体動の全体の平均値と標準偏差を求め、標準偏差に所定の係数（たとえば０．２）を乗算した値を全体の平均値に足し算して、覚醒判定用の第３閾値を設定する。第３閾値を超える領域を覚醒段階と推定する。

合成部２１は、ノンレム判定部１３と、レム判定部１４と、覚醒判定部１５の判定結果を合成して、ひとつの睡眠段階判定データとして出力してもよい。レム判定部１４は、ノンレム判定部１３の判定結果を下地データとして使用して、レム睡眠段階の判定結果を付加してもよい（ＲＥＭ補正）。覚醒判定部１５は、ノンレム判定部１３の判定結果、またはノンレム判定部１３とレム判定部１４の判定結果に、覚醒段階の判定結果を付加してもよい（ＷＡＫＥ補正）。

ノンレム判定部１３と、レム判定部１４と、覚醒判定部１５の判定結果は個別にユーザインタフェース１６から出力されてもよいし、合成部２１で合成された後に出力されてもよい。

ノンレム判定部１３、レム判定部１４、覚醒判定部１５、及び合成部２１の動作をコンピュータプログラムで実現する場合は、メモリ１７に格納された睡眠段階判定プログラムをプロセッサ１１で実行する。睡眠段階判定プログラムは通信インタフェース１８を介してインストールされてもよいし、ドライブ１９を介してリムーバブル記憶媒体からインストールされてもよい。この場合、プログラムは、プロセッサに、
被験者の心拍データを取得する手順と、
前記心拍データを周波数領域の成分に変換する手順と、
前記周波数領域のデータから、２．５秒〜１５０分の周期にわたる範囲の周期成分を抽出して時間領域のデータに変換する手順と、
前記時間領域のデータのうち、所定の閾値を超える領域をレム睡眠段階と判定する手順とを実行させる。

あるいは、コンピュータに、
被験者の体動データを取得する手順と、
前記体動データの全体平均と標準偏差を算出する手順と、
前記標準偏差に所定の係数を乗算した値を前記全体平均に足し合わせた値を閾値として設定する手順と、
前記閾値を超える領域を覚醒段階と判定する手順と、
を実行させてもよい。

また、コンピュータに、
被験者の呼吸データを取得する手順と、
前記呼吸データに基づいて単位時間当たりの呼吸数の全体平均と標準偏差を算出する手順と、
前記標準偏差に所定の係数を乗算した値を前記全体平均に足し合わせた値を閾値として設定する手順と、
前記閾値以下の領域をノンレム睡眠段階と判定する手順と、
を実行させてもよい。

上述した構成と手法により、測定された身体データから、簡単な手法で、より正確に睡眠段階を推定することができる。

１０ 睡眠段階判定装置
１１ プロセッサ
１２ データ取得部
１３ ノンレム判定部
１４ レム判定部
１５ 覚醒判定部
１６ ユーザインタフェース
１７ メモリ
１８ 通信インタフェース
１９ ドライブ
２０ 判定部
２１ 合成部

Claims (6)

1. 被験者の身体データから呼吸データと、心拍データと、体動データを取得するデータ取得部と、
   前記呼吸データに基づいて、単位時間当たりの呼吸数の平均値と標準偏差を算出し、前記標準偏差を前記平均値で除算して単位時間ごとに変動係数を求め、前記変動係数の全体平均に所定の係数を乗算した値を閾値として設定し、前記単位時間ごとの変動係数を前記閾値と比較してノンレム睡眠段階１〜４を推定する判定部と、
   を有し、
   前記判定部は、前記心拍データを用いてレム睡眠段階を判定するときに、前記ノンレム睡眠段階と推定された区間を含めて前記レム睡眠段階を判定することを特徴とする睡眠段階判定装置。
2. 前記データ取得部は、無拘束で測定された前記被験者の前記身体データを取得し、
   前記判定部は、前記呼吸データを用いた前記ノンレム睡眠段階の判定と、前記心拍データを用いた前記レム睡眠段階の判定と、前記体動データを用いた覚醒段階の判定を個別または順次に行って睡眠段階を判定することを特徴とする請求項１に記載の睡眠段階判定装置。
3. 前記判定部は、前記ノンレム睡眠段階と判定されたデータを下地データに用い、前記ノンレム睡眠段階と推定された区間を含めて、前記心拍データを用いて前記レム睡眠段階を判定し、前記レム睡眠段階と推定された区間を含めて、前記体動データを用いて前記覚醒段階を判定する、ことを特徴とする請求項２に記載の睡眠段階判定装置。
4. 前記判定部は、前記覚醒段階の判定結果、前記レム睡眠段階の判定結果、前記ノンレム睡眠段階の判定結果の順に判定を優先することを特徴とする請求項３に記載の睡眠段階判定装置。
5. 被験者の身体データから呼吸データと、心拍データと、体動データを取得し、
   前記呼吸データに基づいて単位時間当たりの呼吸数の平均値と標準偏差を算出し、
   前記標準偏差を前記平均値で除算して単位時間ごとに変動係数を求め、前記変動係数の全体平均に所定の係数を乗算した値を閾値として設定し、
   前記単位時間ごとの変動係数を前記閾値と比較してノンレム睡眠段階１〜４を推定し、
   前記心拍データを用いてレム睡眠段階を判定するときに、前記ノンレム睡眠段階と推定された区間を含めて前記レム睡眠段階を判定する、
   ことを特徴とする睡眠段階判定方法。
6. コンピュータに、
   被験者の身体データから呼吸データ、心拍データと、体動データを取得する手順と、
   前記呼吸データに基づいて単位時間当たりの呼吸数の平均値と標準偏差を算出する手順と、
   前記標準偏差を前記平均値で除算して単位時間ごとに変動係数を求め、前記変動係数の全体平均に所定の係数を乗算した値を閾値として設定する手順と、
   前記単位時間ごとの変動係数を前記閾値と比較してノンレム睡眠段階１〜４を推定する手順と、
   前記心拍データを用いてレム睡眠段階を判定するときに、前記ノンレム睡眠段階と推定された区間を含めて前記レム睡眠段階を判定する手順と、
   を実行させる睡眠段階判定プログラム。
   

Images