JP5993968B2 - 睡眠段階推定装置、方法およびプログラム - Google Patents
睡眠段階推定装置、方法およびプログラム Download PDFInfo
- Publication number
- JP5993968B2 JP5993968B2 JP2015021631A JP2015021631A JP5993968B2 JP 5993968 B2 JP5993968 B2 JP 5993968B2 JP 2015021631 A JP2015021631 A JP 2015021631A JP 2015021631 A JP2015021631 A JP 2015021631A JP 5993968 B2 JP5993968 B2 JP 5993968B2
- Authority
- JP
- Japan
- Prior art keywords
- sleep stage
- time
- sleep
- heart rate
- rate variability
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Landscapes
- Measuring Pulse, Heart Rate, Blood Pressure Or Blood Flow (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
- Measurement And Recording Of Electrical Phenomena And Electrical Characteristics Of The Living Body (AREA)
Description
その第1の態様は、前記被験者の睡眠心電図集合として与えられる、第1の単位時間ごとの睡眠段階の正解ラベルが付与された睡眠時の時間情報を含む睡眠心電図の集合を取得し、この取得された睡眠心電図の集合から、前記第1の単位時間ごとに心拍変動特徴量を算出し、当該算出された心拍変動特徴量を睡眠段階別に集計して心拍変動記憶部に格納する手段と、前記睡眠心電図集合の時間軸を前記第2の単位時間で分割した各時間帯ごとに、前記第1の単位時間ごとの睡眠段階の出現頻度を集計して出現頻度記憶部に格納する手段と、前記第2の単位時間で分割した各時間帯ごとに、前記第1の時間ごとの睡眠段階の遷移の頻度を集計して遷移頻度記憶部に格納する手段を、さらに備えるようにしたものである。
(原理)
この発明は、入眠後の経過時間に従い睡眠段階の出現頻度と遷移確率が変化することに着目し、心拍が脳活動以外の影響を強く受けた場合でも、より尤もらしい睡眠段階を推定できるようにしたものである。例えば、入眠後の深い眠りから、突然覚醒して再び深い眠りに落ちる可能性は低いが、心拍のみで推定すると0.52の確率でその状態が推定される場合がある。この発明では、このようなめったに起こらない遷移や推定の自信度が低い場合に対する推定確率を下げることで、推定精度を向上させる。
(構成)
図3は、この発明の一実施形態に係る睡眠段階推定装置の機能構成を示すブロック図である。
本実施形態における睡眠段階推定装置1は、例えばサーバコンピュータ又はパーソナルコンピュータからなり、制御ユニット10と、記憶ユニット20と、入出力インタフェースユニット30を備えている。このうち入出力インタフェースユニット30は、心電計2から出力された心電データを取り込むと共に、キーボードまたはタッチパネル等の入力部3および液晶デバイス等の表示部4との間でそれぞれ操作データおよび表示データの入出力を行う。
次に、以上のように構成された装置による睡眠段階推定動作を説明する。
(1)学習フェーズ
学習フェーズでは、被験者の心電波形を計測しながら、医療従事者等の測定者がn秒ごとに被験者の睡眠段階を目視観察してその観察結果を「睡眠段階の正解ラベル」として入力部3により入力する。そして、睡眠段階推定装置1が、上記心電波形から心拍変動の特徴量を算出して心拍変動テーブル21に格納すると共に、上記入力された睡眠段階の正解ラベルをもとに時間別の睡眠段階出現頻度と睡眠段階遷移頻度を算出してその算出値をそれぞれ出現頻度テーブル212および遷移頻度テーブル213に格納することにより、行われる。図4は、この学習フェーズにおける睡眠段階推定装置1の処理手順と処理内容を示すフローチャートである。
推定フェーズでは、推定対象者の心電図データのみを計測し、この計測された心電図データから算出される心拍変動特徴量と、心拍変動テーブル211、出現頻度テーブル212および遷移頻度テーブル213に格納されたモデルデータに基づいて、上記推定対象者の睡眠段階St がn秒ごとに推定される。図5はその処理手順と処理内容を示すフローチャートである。
(B) 時間tが、睡眠時間Tよりも小さい間では以下の処理を実行する。
(1) 睡眠段階S_{t,g} を、時間t−1の睡眠段階St-1、時間帯Ct 、時間tにおける心拍変動特徴量ht を用いてP(St =i|ht ,St-1=j,Ct =m)の値に従ってサンプリングする。
(2) 時間t=t+1とする。
(C) ギブスサンプリングカウンタgをインクリメント(g=g+1)する。
以上詳述したように一実施形態では、学習フェーズにおいて、被験者の心電図データから睡眠段階ごとの心拍変動特徴量を算出して心拍変動テーブル211に格納すると共に、時間経過に応じた睡眠段階の出現頻度と遷移頻度を求めてこれらをそれぞれ出現頻度テーブル212および遷移頻度テーブル213に格納する。そして推定フェーズにおいて、推定対象者の心電図データから算出した心拍変動特徴量と、上記各テーブル211,212,213に格納された睡眠段階のモデルデータとに基づいて、上記推定対象者の睡眠段階を推定するようにしている。
なお、この発明は上記実施形態に限定されるものではない。例えば、前記実施形態では睡眠段階推定装置1内に記憶ユニット20を設けた場合を例にとって説明したが、記憶ユニット20をクラウドに設けられたデータベースサーバなどに備え、睡眠段階推定装置1とこのデータベースサーバとの間で通信を行うことによりデータの書き込みおよび読み出しを行うようにしてもよい。
Claims (5)
- 睡眠心電図の集合から第1の単位時間ごとに算出された心拍変動特徴量を睡眠段階別に集計した心拍変動情報を記憶した心拍変動記憶部と、
前記睡眠心電図の時間軸を前記第1の単位時間より長い第2の単位時間で分割した各時間帯ごとに、前記第1の単位時間ごとの睡眠段階の出現頻度を集計した出現頻度情報を記憶した出現頻度記憶部と、
前記分割した各時間帯ごとに、前記第1の時間ごとの睡眠段階の遷移の頻度を集計した遷移頻度情報を記憶した遷移頻度記憶部と、
推定対象の被験者から取得した睡眠心電図から、前記第2の単位時間で分割した各時間帯ごとに、前記第1の時間ごとの心拍変動の特徴量を算出する心拍変動特徴量算出手段と、
前記算出された心拍変動の特徴量の時間的な変化と、前記心拍変動記憶部、出現頻度記憶部および遷移頻度記憶部に記憶された情報とに基づいて、前記被験者の前記第1の単位時間ごとの睡眠段階の確率分布を推定する推定手段と
を具備することを特徴とする睡眠段階推定装置。 - 前記被験者の睡眠心電図集合として与えられる、前記第1の単位時間ごとの睡眠段階の正解ラベルが付与された睡眠時の時間情報を含む睡眠心電図の集合を取得する手段と、
前記取得された睡眠心電図の集合から、前記第1の単位時間ごとに心拍変動特徴量を算出し、当該算出された心拍変動特徴量を睡眠段階別に集計して前記心拍変動記憶部に格納する手段と、
前記睡眠心電図集合の時間軸を前記第2の単位時間で分割した各時間帯ごとに、前記第1の単位時間ごとの睡眠段階の出現頻度を集計して、集計された出現頻度情報を前記出現頻度記憶部に格納する手段と、
前記分割した各時間帯ごとに、前記第1の時間ごとの睡眠段階の遷移の頻度を集計して、集計された遷移頻度情報を前記遷移頻度記憶部に格納する手段と
を、さらに具備することを特徴とする請求項1記載の睡眠段階推定装置。 - 前記心拍変動特徴量算出手段は、心拍変動特徴量として、心拍変動RRI(R-R Interval)の平均値、RRIの標準偏差、前後するRRI間の差の二乗平均値の平方根(RMSSD)、前後するRRI間の時間的差分が予め設定した時間値以上となった回数の発生割合(pNN50)、前後するRRIのローレンツプロットが描く形状の長辺の長さL、RRIのローレンツプロットが描く形状の短辺の長さT、LおよびTから求められる交感神経活動指標(CSI: Cardiac Sympathetic Index)、LおよびTから求められる副交感神経の活動指標(CVI: Cardiac Vagal Index)、RRIのパワースペクトル波形の高周波成分(HF)、または低周波成分(LF)のうちの少なくとも1つ、或いはこれらの特徴量から算出されるその他の特徴量を算出することを特徴とする請求項1または2記載の睡眠段階推定装置。
- 睡眠心電図の集合から第1の単位時間ごとに算出された心拍変動特徴量を睡眠段階別に集計した心拍変動情報を記憶した心拍変動記憶部と、前記睡眠心電図の時間軸を前記第1の単位時間より長い第2の単位時間で分割した各時間帯ごとに、前記第1の単位時間ごとの睡眠段階の出現頻度を集計した出現頻度情報を記憶した出現頻度記憶部と、前記分割した各時間帯ごとに、前記第1の時間ごとの睡眠段階の遷移の頻度を集計した遷移頻度情報を記憶した遷移頻度記憶部とを有するモデルデータベースを備えた睡眠段階推定装置が実行する睡眠方法推定動作方法であって、
前記睡眠段階推定装置が、推定対象の被験者から取得した睡眠心電図から、前記第2の単位時間で分割した各時間帯ごとに、前記第1の時間ごとの心拍変動の特徴量を算出する過程と、
前記睡眠段階推定装置が、前記算出された心拍変動の特徴量の時間的な変化と、前記心拍変動記憶部、出現頻度記憶部および遷移頻度記憶部に記憶された情報とに基づいて、前記被験者の前記第1の単位時間ごとの睡眠段階の確率分布を推定する過程と
を具備することを特徴とする睡眠段階推定動作方法。 - 請求項1乃至3のいずれかに記載の睡眠段階推定装置が具備する各手段が行う処理を、前記睡眠段階推定装置が備えるコンピュータに実行させるプログラム。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2015021631A JP5993968B2 (ja) | 2015-02-05 | 2015-02-05 | 睡眠段階推定装置、方法およびプログラム |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2015021631A JP5993968B2 (ja) | 2015-02-05 | 2015-02-05 | 睡眠段階推定装置、方法およびプログラム |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2016140729A JP2016140729A (ja) | 2016-08-08 |
JP5993968B2 true JP5993968B2 (ja) | 2016-09-21 |
Family
ID=56568030
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2015021631A Active JP5993968B2 (ja) | 2015-02-05 | 2015-02-05 | 睡眠段階推定装置、方法およびプログラム |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5993968B2 (ja) |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7024187B2 (ja) * | 2017-02-13 | 2022-02-24 | 東京電力ホールディングス株式会社 | 学習装置及び属性推定装置 |
BR112020013317A2 (pt) * | 2018-03-22 | 2020-12-01 | Universidad De La Sabana | processo para classificar a profundidade anestésica em intervenções com anestesia total intravenosa |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2006008743A2 (en) * | 2004-07-21 | 2006-01-26 | Widemed Ltd. | Sleep quality indicators |
JP6040874B2 (ja) * | 2013-06-25 | 2016-12-07 | トヨタ自動車株式会社 | 睡眠段階推定装置 |
JP6102581B2 (ja) * | 2013-07-05 | 2017-03-29 | 富士通株式会社 | 睡眠支援装置、睡眠支援方法及びプログラム |
-
2015
- 2015-02-05 JP JP2015021631A patent/JP5993968B2/ja active Active
Also Published As
Publication number | Publication date |
---|---|
JP2016140729A (ja) | 2016-08-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Pham et al. | Heart rate variability in psychology: A review of HRV indices and an analysis tutorial | |
US11877861B2 (en) | Methods and systems for labeling sleep states | |
EP3427653B1 (en) | Biological information analysis device, system, and program | |
JP6122884B2 (ja) | 作業覚醒度推定装置、方法およびプログラム | |
Ho et al. | Predicting survival in heart failure case and control subjects by use of fully automated methods for deriving nonlinear and conventional indices of heart rate dynamics | |
JP5365199B2 (ja) | ライフスタイルを評価するための方法、情報処理装置及びプログラム | |
CN106256316B (zh) | 用于评估生理老化水平的方法和设备 | |
Keissar et al. | Coherence analysis between respiration and heart rate variability using continuous wavelet transform | |
CN107427267B (zh) | 用于导出主体的精神状态的方法和装置 | |
Cicone et al. | How nonlinear-type time-frequency analysis can help in sensing instantaneous heart rate and instantaneous respiratory rate from photoplethysmography in a reliable way | |
JP6513005B2 (ja) | 疲労度計 | |
Valenza et al. | A nonlinear heartbeat dynamics model approach for personalized emotion recognition | |
JP6579890B2 (ja) | 疲労度計 | |
JP6600732B2 (ja) | 睡眠段階推定装置、方法およびプログラム | |
Valenza et al. | Instantaneous transfer entropy for the study of cardiovascular and cardiorespiratory nonstationary dynamics | |
JP6608592B2 (ja) | 概日時間の差の生成 | |
CN103892812A (zh) | 匹配度分析技术 | |
WO2017027232A1 (en) | System and method for sympathetic and parasympathetic activity monitoring by heartbeat | |
JP6559095B2 (ja) | 睡眠モデル作成装置、睡眠段階推定装置、方法およびプログラム | |
JP6518570B2 (ja) | 睡眠段階推定装置、方法およびプログラム | |
Alves et al. | Linear and complex measures of heart rate variability during exposure to traffic noise in healthy women | |
JP5993968B2 (ja) | 睡眠段階推定装置、方法およびプログラム | |
US20210308415A1 (en) | Massage apparatus, system and method capable of deriving a parameter of an individual | |
JP6778657B2 (ja) | モデル選択方法、モデル選択装置およびプログラム | |
Aboamer et al. | Nonlinear features of heart rate variability in paranoid schizophrenic |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20160531 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20160722 |
|
TRDD | Decision of grant or rejection written | ||
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20160816 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20160822 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5993968 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |