JP6932351B2 - Sleep stage determination device, sleep stage determination method and program - Google Patents
Sleep stage determination device, sleep stage determination method and program Download PDFInfo
- Publication number
- JP6932351B2 JP6932351B2 JP2017173563A JP2017173563A JP6932351B2 JP 6932351 B2 JP6932351 B2 JP 6932351B2 JP 2017173563 A JP2017173563 A JP 2017173563A JP 2017173563 A JP2017173563 A JP 2017173563A JP 6932351 B2 JP6932351 B2 JP 6932351B2
- Authority
- JP
- Japan
- Prior art keywords
- period
- rem sleep
- determined
- sleep
- rem
- 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
- 238000000034 method Methods 0.000 title claims description 130
- 230000008667 sleep stage Effects 0.000 title claims description 86
- 230000036385 rapid eye movement (rem) sleep Effects 0.000 claims description 187
- 230000008569 process Effects 0.000 claims description 39
- 238000003066 decision tree Methods 0.000 claims description 16
- 230000007423 decrease Effects 0.000 claims description 14
- 238000010801 machine learning Methods 0.000 claims description 13
- 230000007958 sleep Effects 0.000 description 26
- 230000008452 non REM sleep Effects 0.000 description 10
- 238000010586 diagram Methods 0.000 description 6
- 238000007637 random forest analysis Methods 0.000 description 5
- 238000004422 calculation algorithm Methods 0.000 description 3
- 238000001514 detection method Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 230000029058 respiratory gaseous exchange Effects 0.000 description 3
- 210000004556 brain Anatomy 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000012790 confirmation Methods 0.000 description 2
- 230000002068 genetic effect Effects 0.000 description 2
- 238000010606 normalization Methods 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- 206010062519 Poor quality sleep Diseases 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000037007 arousal Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 238000003745 diagnosis Methods 0.000 description 1
- 230000004424 eye movement Effects 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000003340 mental effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 208000019116 sleep disease Diseases 0.000 description 1
- 230000003860 sleep quality Effects 0.000 description 1
- 239000007787 solid Substances 0.000 description 1
Images
Landscapes
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Description
本発明は、レム睡眠状態を判定する睡眠段階判定装置及び睡眠段階判定方法、並びに睡眠段階の判定を実行するプログラムに関する。 The present invention relates to a sleep stage determination device for determining a REM sleep state, a sleep stage determination method, and a program for executing sleep stage determination.
医療現場では、睡眠障害などによる睡眠の質の低下を診断するために、被測定者の睡眠段階を測定することが行われている。人間の睡眠段階は、睡眠の深さの観点で6段階に分類したものが知られており、その6つの睡眠段階は、眠りが浅い段階から順に、覚醒、レム睡眠、ノンレム睡眠(ステージ1〜4)と呼ばれている。これらの6段階の睡眠段階の判定は、従来、例えば被計測者の顔や頭部に多数の電極を装着して、その多数の電極から脳波、眼球運動、及び顎筋電を測定し、測定結果の解析により行われていた。 In the medical field, the sleep stage of the subject is measured in order to diagnose the deterioration of sleep quality due to sleep disorders and the like. Human sleep stages are known to be classified into 6 stages from the viewpoint of sleep depth, and the 6 sleep stages are awakening, REM sleep, and non-REM sleep (stages 1 to 1) in order from the lightest sleep stage. It is called 4). Conventionally, these six stages of sleep are determined by attaching a large number of electrodes to the face or head of the person to be measured and measuring brain waves, eye movements, and jaw myoelectricity from the large number of electrodes. It was done by analyzing the results.
このような顔や頭部に多数の電極を装着した状態での睡眠段階の検査は、通常、医療機関に宿泊して、長時間連続して電極を身体に装着して行う検査であり、日常と異なる状況での睡眠であるため、被測定者(患者)に精神的な負担と肉体的な負担を強いることになる。また、電極を使って取得したデータは、専門知識と経験を持つ医師が解析して判定する必要があり、睡眠状態を簡単に判定できるものではないという問題があった。 Such an examination at the sleep stage with a large number of electrodes attached to the face or head is usually an examination performed by staying at a medical institution and continuously attaching the electrodes to the body for a long period of time. Because it sleeps in a different situation from the above, it imposes a mental burden and a physical burden on the subject (patient). Further, the data acquired by using the electrodes needs to be analyzed and judged by a doctor having specialized knowledge and experience, and there is a problem that the sleep state cannot be easily judged.
これらの問題を解決するために、専門医師による診断を不要とする睡眠段階推定手法が従来から数多く提案されている。
例えば、特許文献1には、遺伝的アルゴリズムによる学習手法を改良したDatabase-based Compact Genetic Algorithmと称される手法で、マットレス型圧力センサの検出データから睡眠段階を推定する技術が記載されている。この特許文献1に記載された技術は、マットレス型圧力センサが検出した被測定者の体動と心拍に基づいて、睡眠段階を推定するものである。このようなマットレス型圧力センサを使って睡眠段階を推定することで、被測定者に負担を強いることなく、被測定者の睡眠状態を推定することができる。
In order to solve these problems, many sleep stage estimation methods that do not require diagnosis by a specialist doctor have been proposed.
For example, Patent Document 1 describes a technique for estimating a sleep stage from detection data of a mattress type pressure sensor by a method called Database-based Compact Genetic Algorithm, which is an improvement of a learning method using a genetic algorithm. The technique described in Patent Document 1 estimates the sleep stage based on the body movement and heartbeat of the subject detected by the mattress type pressure sensor. By estimating the sleep stage using such a mattress type pressure sensor, it is possible to estimate the sleep state of the person to be measured without imposing a burden on the person to be measured.
ところで、上述した6段階の睡眠段階を考えた場合に、眠りが浅い段階から順に、覚醒、レム睡眠、4つのノンレム睡眠(ステージ1〜4)が存在するが、これら6つの段階は、検出データを6つの範囲に区切って得るような単純な処理では判定できないという問題がある。特に、レム睡眠は、ノンレム睡眠や覚醒とは全く異なる睡眠状態であり、従来、レム睡眠を精度良く判定することは非常に困難であった。
特許文献1に記載されるように、マットレス型圧力センサのような圧力センサを使って睡眠段階を判定することができれば、被測定者に負担を強いることなく睡眠状態の診断が可能になる。しかし、従来から提案されている手法では、レム睡眠の判定精度が高くなく、結果的に睡眠段階の推定の信頼性がそれほど良くないという問題があった。
By the way, when considering the above-mentioned 6 stages of sleep, there are awakening, REM sleep, and 4 non-REM sleeps (stages 1 to 4) in order from the light sleep stage, and these 6 stages are detection data. There is a problem that it cannot be determined by a simple process such as dividing the above into six ranges. In particular, REM sleep is a sleep state completely different from non-REM sleep and awakening, and it has been very difficult to accurately determine REM sleep in the past.
As described in Patent Document 1, if the sleep stage can be determined by using a pressure sensor such as a mattress type pressure sensor, it is possible to diagnose the sleep state without imposing a burden on the person to be measured. However, the conventionally proposed method has a problem that the determination accuracy of REM sleep is not high, and as a result, the reliability of estimation of the sleep stage is not so good.
本発明の目的は、被測定者に負担をかけずに、被測定者のレム睡眠を精度良く判定できるようにすることにある。 An object of the present invention is to make it possible to accurately determine the REM sleep of a person to be measured without imposing a burden on the person to be measured.
上記課題を解決するため、本発明の睡眠段階判定装置は、データ処理部と睡眠段階判定部とを備える。
睡眠段階判定部は、被測定者の第1の期間毎の心拍の中央値の増加率が、所定の閾値以上であるとき、第1のレム睡眠期間の開始と判定し、第1の期間毎の心拍の中央値の減少率が減少したとき、第1のレム睡眠期間の終了を判定すると共に、被測定者の体動を第2の期間毎に周波数解析したデータの分析で、レム睡眠の判定を得る決定木を機械学習により複数生成し、生成した複数の決定木で得られたレム睡眠か否かの判定の多数決から、第2の期間が第2のレム睡眠期間であると判定する。
睡眠段階判定部は、第1のレム睡眠期間と判定した期間と第2のレム睡眠期間と判定した期間とが重なる期間を含み、かつ一定間隔以内に第2のレム睡眠期間と判定した期間が存在する場合に、その一定間隔以内に第2のレム睡眠期間と判定した期間を、レム睡眠期間と確定する。
In order to solve the above problems, the sleep stage determination device of the present invention includes a data processing unit and a sleep stage determination unit.
The sleep stage determination unit determines that the start of the first REM sleep period is started when the rate of increase in the median heart rate of the subject for each first period is equal to or higher than a predetermined threshold, and every first period. When the rate of decrease in the median heart rate of REM sleep decreases, the end of the first REM sleep period is determined, and the body movements of the subject are frequency-analyzed for each second period. A plurality of decision trees for which a judgment is obtained are generated by machine learning, and it is determined that the second period is the second REM sleep period from the majority decision of whether or not the determination tree is REM sleep obtained from the generated plurality of decision trees. ..
The sleep stage determination unit includes a period in which the period determined to be the first REM sleep period and the period determined to be the second REM sleep period overlap, and the period determined to be the second REM sleep period within a certain interval is If present, the period determined to be the second REM sleep period within the fixed interval is determined as the REM sleep period.
また、本発明の睡眠段階判定方法は、第1のレム睡眠期間判定処理と、第2のレム睡眠期間判定処理と、睡眠段階判定処理と、を含む。
第1のレム睡眠期間判定処理は、被測定者の第1の期間毎の心拍の中央値の増加率が、所定の閾値以上であるとき、第1のレム睡眠期間の開始と判定し、第1の期間毎の心拍の中央値の減少率が減少したとき、第1のレム睡眠期間の終了を判定する。
第2のレム睡眠期間判定処理は、被測定者の体動を第2の期間毎に周波数解析したデータの分析で、レム睡眠の確率値を得る決定木を機械学習により複数生成し、生成した複数の決定木で得られたレム睡眠か否かの判定の多数決から、第2の期間が第2のレム睡眠期間であると判定する。
睡眠段階判定処理は、第1のレム睡眠期間と判定した期間と第2のレム睡眠期間と判定した期間とが重なる期間を含み、かつ一定間隔以内に第2のレム睡眠期間と判定した期間が存在する場合に、その一定間隔以内に第2のレム睡眠期間と判定した期間を、レム睡眠期間と確定する。
In addition, the sleep stage determination method of the present invention includes a first REM sleep period determination process, a second REM sleep period determination process, and a sleep stage determination process.
The first REM sleep period determination process determines that the start of the first REM sleep period is started when the rate of increase in the median heart rate of the subject for each first period is equal to or greater than a predetermined threshold, and the first REM sleep period is determined. When the rate of decrease in the median heart rate for each period decreases, the end of the first REM sleep period is determined.
The second REM sleep period determination process is an analysis of data obtained by frequency-analyzing the body movements of the subject for each second period, and generates a plurality of determination trees for obtaining the probability value of REM sleep by machine learning. From the majority decision of whether or not it is REM sleep obtained from a plurality of determination trees, it is determined that the second period is the second REM sleep period.
The sleep stage determination process includes a period in which the period determined to be the first REM sleep period and the period determined to be the second REM sleep period overlap, and the period determined to be the second REM sleep period within a certain interval is If present, the period determined to be the second REM sleep period within the fixed interval is determined as the REM sleep period.
また、本発明のプログラムは、上記睡眠段階判定方法の第1のレム睡眠期間判定処理と第2のレム睡眠期間判定処理と睡眠段階判定処理とを実行する手順を、コンピュータに実行させるものである。 Further, the program of the present invention causes a computer to execute a procedure for executing the first REM sleep period determination process, the second REM sleep period determination process, and the sleep stage determination process of the sleep stage determination method. ..
本発明によれば、圧力センサから得た被測定者の睡眠中の圧力データに基づいて、レム睡眠の期間を正確に判定できるようになる。 According to the present invention, the duration of REM sleep can be accurately determined based on the pressure data during sleep of the subject obtained from the pressure sensor.
以下、本発明の一実施の形態例(以下、「本例」と称する)について、添付図面を参照して説明する。
[1.睡眠段階判定装置の構成]
図1は、本例の睡眠段階判定装置10の構成を示すブロック図である。図2は、本例の睡眠段階判定装置10を使って睡眠段階の判定を行う状態の例を示す図である。
本例の睡眠段階判定装置10は、被測定者の体動をマットレスセンサ2で圧力データとして取得する。マットレスセンサ2は、例えば図2に示すように、被測定者Aが睡眠を行うベッド1のマットレスの上に敷いて、被測定者Aの睡眠中の上半身の体動を圧力の変化として検出する。被測定者Aの下側になるマットレスの上にマットレスセンサ2を配置するのは一例であり、例えばマットレスの中にマットレスセンサ2を内蔵させてもよい。
図2では、ベッド1の脇に睡眠段階判定装置10を設置し、マットレスセンサ2と睡眠段階判定装置10をケーブルで接続した例を示すが、例えばマットレスセンサ2が取得した圧力データを、無線伝送で別の部屋の睡眠段階判定装置10に伝送するようにしてもよい。
Hereinafter, an example of an embodiment of the present invention (hereinafter, referred to as “this example”) will be described with reference to the accompanying drawings.
[1. Configuration of sleep stage determination device]
FIG. 1 is a block diagram showing a configuration of the sleep
The sleep
FIG. 2 shows an example in which the sleep
睡眠段階判定装置10は、図1に示すように、生体データ取得部11と、生体データ処理部12と、睡眠段階判定部13と、出力部14とを備える。
生体データ取得部11は、マットレスセンサ2が出力する圧力データを取得する。生体データ取得部11が取得した圧力データは、データ処理部12に供給される。生体データ処理部12は、供給される圧力データの解析処理を行う。具体的には、圧力データをサンプリングしてデジタルデータ化し、そのデジタルデータ化された圧力データの高速フーリエ変換処理(以下、「FFT処理」と称する)により、各周波数の成分を取得する。そして、取得したそれぞれの周波数成分の状態から、心拍と呼吸の状態を判別する。すなわち、圧力データをFFT処理で得た周波数成分の内で、心拍及び呼吸に相当する周波数成分から、心拍数と呼吸数を判別する。
さらに、生体データ処理部12は、FFT処理で得た周波数成分を使って、睡眠段階の解析に必要なデータ処理を行う。
As shown in FIG. 1, the sleep
The biological
Further, the biological
そして、生体データ処理部12で解析した結果が睡眠段階判定部13に供給される。睡眠段階判定部13は、生体データ処理部12で解析した結果に基づいて、被測定者Aの睡眠状態を判定する。睡眠段階判定部13は、レム睡眠の睡眠状態と、レム睡眠以外の睡眠状態の判定を行う。さらに、睡眠段階判定部13は、レム睡眠以外の睡眠状態として、覚醒とノンレム睡眠(ステージ1〜4)の5段階の判定を行う。なお、以下の説明では、レム睡眠の睡眠状態と、レム睡眠以外の睡眠状態との判定を行う処理を中心に説明し、レム睡眠以外の睡眠状態を5段階(覚醒及びノンレム睡眠のステージ1〜4)に判定する処理の説明は省略する。なお、ノンレム睡眠の睡眠段階は、睡眠の深さと対応付けて考えられており、ステージ1から4の順に睡眠が深くなることを示す。例えば、ステージ2はステージ1よりも睡眠が深く、「ステージ2以下」と記載した場合には、ステージ2〜4を示す。
Then, the result analyzed by the biological
出力部14は、睡眠段階判定部13が判定した睡眠段階を出力する。出力部14は、例えば表示装置により構成され、睡眠段階を表示する。あるいは、出力部14を記録装置として構成して、一晩の睡眠状態を記録するようにしてもよい。
The
[2.睡眠段階判定装置のハードウェア構成例]
図3は、睡眠段階判定装置10をコンピュータ装置で構成した場合のハードウェア構成例を示す。
コンピュータ装置Cは、バスC8に接続されたCPU(Central Processing Unit:中央処理装置)C1、ROM(Read Only Memory)C2、及びRAM(Random Access Memory)C3を備える。さらに、コンピュータ装置Cは、不揮発性ストレージC4、ネットワークインターフェイスC5、入力装置C6、及び表示装置C7を備える。
[2. Hardware configuration example of sleep stage determination device]
FIG. 3 shows an example of hardware configuration when the sleep
The computer device C includes a CPU (Central Processing Unit) C1, a ROM (Read Only Memory) C2, and a RAM (Random Access Memory) C3 connected to the bus C8. Further, the computer device C includes a non-volatile storage C4, a network interface C5, an input device C6, and a display device C7.
CPU C1は、睡眠段階判定装置10の生体データ処理部12や睡眠段階判定部13が備える各機能を実現するソフトウェアのプログラムコードをROM C2から読み出して実行する。圧力データを周波数解析するFFT処理についても、該当する処理を実行するプログラムをROM C2から読み出して、CPU C1が実行する。RAM C3には、演算処理の途中に発生した変数やパラメータ等が一時的に書き込まれる。
The CPU C1 reads out from the ROM C2 the program code of the software that realizes each function included in the biological
不揮発性ストレージC4としては、例えば、HDD(Hard disk drive)、SSD(Solid State Drive)、フレキシブルディスク、光ディスク、光磁気ディスク、CD−ROM、CD−R、磁気テープ、不揮発性のメモリ等が用いられる。この不揮発性ストレージC4には、OS(Operating System)、各種のパラメータの他に、コンピュータ装置Cを睡眠段階判定装置10として機能させるためのプログラムが記録されている。また、睡眠段階判定部13が判定した睡眠段階についてのデータが、不揮発性ストレージC4に記録される。
As the non-volatile storage C4, for example, HDD (Hard disk drive), SSD (Solid State Drive), flexible disk, optical disk, optical magnetic disk, CD-ROM, CD-R, magnetic tape, non-volatile memory and the like are used. Be done. In this non-volatile storage C4, in addition to the OS (Operating System) and various parameters, a program for making the computer device C function as the sleep
ネットワークインターフェイスC5には、例えば、NIC(Network Interface Card)等が用いられ、端子が接続されたLAN(Local Area Network)、専用線等を介して各種のデータを送受信することが可能である。例えば、コンピュータ装置Cは、マットレスセンサ2が出力する圧力データを、ネットワークインターフェイスC5を介して取得する。入力装置C6は、例えばキーボード等の機器で構成され、この入力装置C6により、睡眠段階判定装置10で睡眠段階を判定する期間の設定や、睡眠段階の表示形態の指示等が行われる。表示装置C7は、睡眠段階判定装置10で判定された睡眠段階を表示する。
なお、睡眠段階判定装置10をコンピュータ装置から構成するのは一例であり、例えばFFT処理などのデータ処理を、専用のハードウェアを用意して行うようにしてもよい。
For the network interface C5, for example, a NIC (Network Interface Card) or the like is used, and various data can be transmitted / received via a LAN (Local Area Network) to which terminals are connected, a dedicated line, or the like. For example, the computer device C acquires the pressure data output by the
It should be noted that the sleep
[3.レム睡眠の判定処理]
次に、本例の睡眠段階判定装置10が行うレム睡眠の判定処理について説明する。
本例の睡眠段階判定装置10は、2種類の手法によりレム睡眠の判定処理(第1のレム睡眠期間判定処理及び第2のレム睡眠期間判定処理)を行い、その2種類の手法によるレム睡眠の判定結果を総合的に判断して、最終的にレム睡眠の期間の判定結果を得る。ここでは、2種類の手法を手法1及び手法2と称する。
2種類の手法の詳細は後述するが、手法1は、心拍数変動に着目したレム睡眠判定手法であり、手法2は、全周波数成分を用いた機械学習によるレム睡眠判定手法である。睡眠段階判定装置10は、これらの2つの手法を使ったハイブリッド構造でレム睡眠判定を行うものである。
[3. REM sleep judgment processing]
Next, the REM sleep determination process performed by the sleep
The sleep
The details of the two types of methods will be described later, but the method 1 is a REM sleep determination method focusing on heart rate fluctuation, and the
まず、手法1によるレム睡眠の判定処理について説明する。
この手法1では、一定時間の心拍情報から心拍数変動の中周波成分を三角関数回帰モデル(Trigonometric Function Model:TFM)を用いてリアルタイムに推定し、算出した心拍変動から睡眠深度を推定する。ここでは、レム睡眠時に心拍数が揺らぐという生理学的な観点を利用して、レム睡眠を推定する。
First, the REM sleep determination process according to Method 1 will be described.
In this method 1, the medium frequency component of the heart rate variability is estimated in real time from the heart rate information for a certain period of time using a trigonometric function model (TFM), and the sleep depth is estimated from the calculated heart rate variability. Here, REM sleep is estimated using the physiological viewpoint that the heart rate fluctuates during REM sleep.
図4は、手法1によるレム睡眠の判定処理を示すフローチャートである。
まず、生体データ処理部12は、FFT処理で周波数解析したデータから、心拍の値を取得する。生体データ処理部12は、検出された心拍について、一定時間(ここでは5分間)毎の中央値HRを算出する(ステップS111)。
そして、生体データ処理部12は、検出した心拍の中央値HRの5分毎の変動量ΔHRを算出し、その5分毎の変動量ΔHRを閾値TH1と比較する(ステップS112)。
ここでは、推定する時刻から直前5分間の心拍の中央値をΔHRcurr、5分前から10分前までの心拍の中央値をΔHRprevとしたとき、次の[数1]式により、変動量ΔHRを得る。
FIG. 4 is a flowchart showing the determination process of REM sleep by the method 1.
First, the biological
Then, the biological
Here, when the median heart rate for 5 minutes immediately before the estimated time is ΔHR curr and the median heart rate from 5 minutes to 10 minutes before is ΔHR prev , the amount of fluctuation is determined by the following equation [Equation 1]. Obtain ΔHR.
ステップS112では、このようにして得た5分毎の変動量ΔHRと閾値TH1とを比較する。このステップS112の比較では、閾値TH1を0.04とし、ΔHR>0.04の条件を満たすか否かを判断する。
ステップS112で、ΔHR>0.04の条件を満たさない場合には(ステップS112のNO)、生体データ処理部12は、該当する条件を満たすまで待機する。
そして、ΔHR>0.04の条件を満たすと判断したとき(ステップS112のYES)、生体データ処理部12は、レム睡眠期間(手法1によるレム睡眠期間)の開始と判断する(ステップS113)。その後、生体データ処理部12は、変動量ΔHRの減少率が減少したか否かを判断する(ステップS114)。なお、「減少率が減少」するとは、例えば、ΔHRの値が0.07から0.06に変わることをいう。ここで、変動量ΔHRの減少率が減少しないとき(ステップS114のNO)、生体データ処理部12は、変動量ΔHRの減少率が減少するまで待機する。また、変動量ΔHRの減少率が減少したとき(ステップS114のYES)、生体データ処理部12は、レム睡眠期間(手法1によるレム睡眠期間)の終了を判断し(ステップS115)、ステップS112の変動量ΔHRと閾値TH1との比較に戻る。
このように、手法1では、心拍数の変動量に基づいてレム睡眠の期間を判断することができる。
In step S112, the variation amount ΔHR obtained in this way every 5 minutes is compared with the threshold value TH1. In the comparison in step S112, the threshold value TH1 is set to 0.04, and it is determined whether or not the condition of ΔHR> 0.04 is satisfied.
If the condition of ΔHR> 0.04 is not satisfied in step S112 (NO in step S112), the biometric
Then, when it is determined that the condition of ΔHR> 0.04 is satisfied (YES in step S112), the biological
As described above, in the method 1, the period of REM sleep can be determined based on the fluctuation amount of the heart rate.
次に、手法2によるレム睡眠の判定処理について説明する。
この手法2は、アンサンブル学習による機械学習手法であるランダムフォレスト(Random Forests)を用いて、レム睡眠とレム睡眠以外とを判定する手法である。この手法2では、一定期間(ここでは30秒間)のマットレスセンサ2が取得した圧力データに対してFFT処理を施すことで周波数解析を行って、レム睡眠を判定するものである。
Next, the REM sleep determination process according to the
This
図5は、手法2によるレム睡眠の判定処理を示すフローチャートである。
まず、生体データ処理部12は、生体データ取得部11が取得した生体データを一定期間(1エポック区間:30秒間)ごとに分割し、一定期間ごとの生体データに対してFFT処理で周波数解析を行うことで周波数表現された生体データ生成する(ステップS121)。そして、FFT処理された生体データの全周波数成分について正規化を行い、全周波数成分の総和が1になるようにする。なお、ここで正規化を行うのは一例であり、正規化を行わないデータを扱うようにしてもよい。
そして、生体データ処理部12は、その一定期間毎の圧力データを使って、ランダムフォレストの機械学習手法による2値判定を行う(ステップS122)。ランダムフォレストの機械学習による2値判定の詳細については後述する。
FIG. 5 is a flowchart showing the determination process of REM sleep by the
First, the biometric
Then, the biological
その後、生体データ処理部12は、ステップS122による2値判定で、レム睡眠の条件を満たすか否かを判断する(ステップS123)。この判断でレム睡眠の条件を満たさない場合には(ステップS123のNO)、ステップS121の処理に戻り、生体データ処理部12は、次の一定期間の処理に移る。
そして、この判断でレム睡眠の条件を満たした場合には(ステップS123のYES)、生体データ処理部12は、該当する一定期間(30秒間)をレム睡眠期間(手法2によるレム睡眠期間)に設定し(ステップS124)、その後ステップS121の処理に戻り、次の一定期間の処理に移る。
After that, the biological
Then, when the condition of REM sleep is satisfied by this determination (YES in step S123), the biological
図6は、ランダムフォレストの機械学習による2値判定状態を説明する図である。
まず、図6の上側に示すように、生体データ処理部12は、30秒間の圧力データD1を取得する。圧力データD1のグラフの縦軸はセンサ値であり、横軸は時間(30秒間)を示す。
この圧力データD1について、生体データ処理部12がFFT処理を施すことで、周波数成分毎の強度値を示したFFTデータD2を得る。FFTデータD2のグラフの縦軸は強度値であり、横軸は周波数(0.016Hz〜5.0Hz)を示す。ここでは、FFTデータD2として、5Hz以下の周波数成分のみの300次元ベクトルのデータとする。
FIG. 6 is a diagram for explaining a binary determination state by machine learning of a random forest.
First, as shown on the upper side of FIG. 6, the biological
The biometric
生体データ処理部12は、この300次元ベクトルのデータから、ランダムフォレストと称される機械学習のアルゴリズムを使って、予測モデルである決定木を複数生成する。すなわち、図6に示すように、複数の決定木T1,T2,・・・,TN(Nは任意の整数)を生成する。それぞれの決定木は、各周波数成分の状態がどのようであるかに基づいて分岐して、各睡眠段階(ここではレム睡眠段階)の判定を得るものである。例えば、図6の左下に示すように、決定木T1のある段階Saでは、周波数成分0.3Hzの正規化値が0.25より小さいか否かを判断し、その判断で、正規化値が0.25より小さいとき未確定と判定し、0.25以上のとき、レム睡眠と判定する。そして、その未確定の場合(0.3Hzの成分の正規化値が0.25より小さいとき)に、次の段階Sbで、周波数成分1.5Hzの正規化値が0.05より大きいか否かを判断し、その判断で、正規化値が0.05より大きいとき、ノンレム睡眠と判定し、正規化値が0.05以下のとき、レム睡眠と判定する。
The biological
このようにして、決定木T1に基づいてFFTデータD2を解析することで、レム睡眠か否かの判定q1を得る。また、決定木T1とは別の条件の決定木T2に基づいてFFTデータD2を解析することで、レム睡眠か否かの判定q2を得る。このようにして、それぞれ異なる条件の決定木T1,T2,・・・,TN毎に、レム睡眠か否かの判定q1,q2,・・・,qNを得る。
そして、生体データ処理部12は、複数の決定木T1,T2,・・・,TN毎に得たレム睡眠か否かの判定q1,q2,・・・,qNの多数決で、現在の30秒間の期間がレム睡眠か否かを判定する。
In this way, the decision tree T by analyzing the FFT data D2 based on 1, to obtain a judgment q 1 whether REM sleep. Further, by analyzing the FFT data D2 based on the decision tree T 2 under a condition different from that of the decision tree T 1 , the determination q 2 of whether or not it is REM sleep is obtained. In this way, for each decision tree T 1 , T 2 , ..., TN under different conditions, the determination q 1 , q 2 , ..., Q N of whether or not it is REM sleep is obtained.
Then, biometric
図7は、手法1でレム睡眠と判定した区間と、手法2でレム睡眠と判定した区間とを使って、レム睡眠期間を確定する処理例を示すフローチャートである。このレム睡眠期間を確定する処理は、睡眠段階判定部13が判定する。
まず、睡眠段階判定部13は、手法1でレム睡眠と判定した区間があるか否かを判断する(ステップS131)。この判断で、手法1によりレム睡眠と判定した区間がない場合には(ステップS131のNO)、睡眠段階判定部13は、探索した期間内にレム睡眠の期間がないと判断して、レム睡眠期間を探索する処理を終了する。
そして、手法1によりレム睡眠と判定した区間がある場合には(ステップS131のYES)、睡眠段階判定部13は、手法1によりレム睡眠と判定した区間と、手法2によりレム睡眠と判定した区間とが重なる区間があるか否かを判断する(ステップS132)。
FIG. 7 is a flowchart showing a processing example of determining the REM sleep period by using the section determined to be REM sleep by the method 1 and the section determined to be REM sleep by the
First, the sleep
Then, when there is a section determined to be REM sleep by the method 1 (YES in step S131), the sleep
さらに、ステップS132の判断で、手法1によりレム睡眠と判定した区間と、手法2によりレム睡眠と判定した区間とが重なる区間がある場合には(ステップS132のYES)、手法1によりレム睡眠と判定した区間を含み、一定期間以内(ここでは2分以内)に手法2でレム睡眠と判定した区間が存在する全ての期間を、レム睡眠期間と確定する(ステップS133)。また、ステップS132の判断で、手法1によりレム睡眠と判定した区間と、手法2によりレム睡眠と判定した区間とが重なっていない場合には(ステップS132のNO)、手法1によりレム睡眠と判定した区間を、レム睡眠期間と確定する(ステップS134)。
Further, if there is a section in which the section determined to be REM sleep by the method 1 and the section determined to be REM sleep by the
図8は、図7のフローチャートに示した手法1での判定結果と手法2での判定結果から、レム睡眠の期間を確定する処理状態の例を示す説明図である。
図8Aの例は、ある期間に、手法1による判定でレム睡眠と判定した区間があり、その区間と重なる期間に、手法2による判定でレム睡眠と判定した区間がない場合を示す。この図8Aに示す状態の場合には、睡眠段階判定部13は、手法1でレム睡眠と判定した区間を、レム睡眠区間と確定する。この図8Aの確定処理は、図7のフローチャートのステップS134でレム睡眠と確定する処理に相当する。
FIG. 8 is an explanatory diagram showing an example of a processing state for determining the period of REM sleep from the determination result in the method 1 and the determination result in the
The example of FIG. 8A shows a case where there is a section determined to be REM sleep by the determination by the method 1 in a certain period, and there is no section determined to be REM sleep by the determination by the
図8Bの例は、ある期間に、手法1による判定でレム睡眠と判定した区間があり、その区間と重なる期間に、手法2による判定でレム睡眠と判定した区間がある場合を示す。この場合には、睡眠段階判定部13は、手法1による判定でレム睡眠と判定した区間を全てレム睡眠期間と確定すると共に、手法2による判定でレム睡眠と判定した区間についても、レム睡眠期間と確定する。なお、ここでの手法2による判定でレム睡眠と判定した区間は、手法2によって、2分以内にレム睡眠と判定した区間がある期間を、連続したレム睡眠期間とする。この図8Bの確定処理は、図7のフローチャートのステップS133でレム睡眠と確定する処理に相当する。
The example of FIG. 8B shows a case where there is a section determined to be REM sleep by the determination by the method 1 in a certain period, and there is a section determined to be REM sleep by the determination by the
図8Cの例は、ある期間に、手法2による判定でレム睡眠と判定した区間があり、手法1では該当する区間をレム睡眠と判定しない場合を示す。この場合には、睡眠段階判定部13は、該当する区間をレム睡眠とは確定しない。
The example of FIG. 8C shows a case where there is a section determined to be REM sleep by the determination by the
図9は、被測定者の睡眠状態を判定した例を示す。図9の特性BW及び特性HRの横軸は時間を示し、この図9では数時間の睡眠中の特性BW及び特性HRを示す。
図9の上側に示す特性BWは、被測定者の睡眠状態(睡眠段階)を、脳波計により測定した例(つまり本例の睡眠段階判定装置10とは別の装置により測定した例)を示す。この特性BWのグラフの縦軸において、Aと示す位置は覚醒を示し、1,2,3,4の数字で示す位置はノンレム睡眠のステージ1〜4を示す。図9に示す特性BWでは、レム睡眠の期間の判定を行っていないが、ノンレム睡眠のステージ1と覚醒との間に、レム睡眠が存在することになる。
FIG. 9 shows an example of determining the sleep state of the person to be measured. The horizontal axis of the characteristic BW and the characteristic HR in FIG. 9 indicates the time, and in this FIG. 9, the characteristic BW and the characteristic HR during sleep for several hours are shown.
The characteristic BW shown on the upper side of FIG. 9 shows an example in which the sleep state (sleep stage) of the subject is measured by an electroencephalograph (that is, an example measured by a device different from the sleep
ここで、図9の特性BWに重ねて示す区間R1,R2,R3,R4,R5は、本例の睡眠段階判定装置10が手法1と手法2の双方の判定結果から確定したレム睡眠期間である。
Here, the sections R1, R2, R3, R4, and R5 shown overlaid on the characteristic BW of FIG. 9 are REM sleep periods determined by the sleep
図9の下側に示す特性HRは、心拍数の変化を示す。この心拍数を示す特性HRは、図2に示すマットレスセンサ2が検出した圧力データに基づいて、睡眠段階判定装置10が取得したものである。この心拍数の変動を示す特性HRの内で、前後の心拍数よりも一時的に高くなるピークの位置ra1,ra2,ra3,ra4,ra5(図9で破線の円で示す部分のほぼ中央部)は、手法1でレム睡眠と判定した区間に相当する。
手法1でレム睡眠と判定した区間ra1,ra2,ra3,ra4,ra5は、特性HRの心拍数がピークとなる位置であり、比較的短時間の区間である。
The characteristic HR shown at the bottom of FIG. 9 shows the change in heart rate. The characteristic HR indicating the heart rate is acquired by the sleep
The sections ra1, ra2, ra3, ra4, and ra5 determined to be REM sleep in the method 1 are positions where the heart rate of the characteristic HR peaks, and are relatively short sections.
また、図9に示す複数の区間rbは、それぞれが手法2でレム睡眠と判定した区間である。図9では、レム睡眠期間R4及びR5の付近のみ、手法2でレム睡眠と判定した区間rbを示し、その他の区間は図示を省略する。
本例の場合、手法2でレム睡眠と判定する1つの区間rbは30秒間であり、その30秒間の区間rbが2分以内に存在したとき、手法2で判定したレム睡眠が連続していると見なして、1つのレム睡眠期間(例えば期間R4やR5)が確定する。但し、本例の場合には、その30秒間の区間rbが2分以内に存在した期間が、手法1でレム睡眠と判定した区間ra1〜ra5と重なっていない箇所については、レム睡眠でないと判定する。
Further, each of the plurality of sections rb shown in FIG. 9 is a section determined to be REM sleep by the
In the case of this example, one section rb determined to be REM sleep by the
このようにして、心拍に基づいてレム睡眠を判定する手法1と、マットレスセンサ2が検出した圧力データの周波数分布の機械学習に基づいてレム睡眠を判定する手法2とを組み合わせて、レム睡眠の期間を判定することで、非常に精度の高いレム睡眠の判定ができるようになる。
具体的には、心拍に基づいてレム睡眠を判定する手法1では、図9の区間ra1〜ra5で分かるように、比較的短時間の区間をレム睡眠と判定しており、実際のレム睡眠の長さとは一致していない可能性が高い。
一方、検出データの周波数分布の機械学習からレム睡眠を判定する手法2では、例えば図9の区間R4と区間R5の間の、最終的にレム睡眠と確定しない期間であっても、レム睡眠と判定される区間rbが存在し、レム睡眠と判定する精度がそれ程高くない可能性があった。
In this way, the method 1 for determining REM sleep based on the heartbeat and the
Specifically, in the method 1 for determining REM sleep based on the heartbeat, as can be seen from the sections ra1 to ra5 in FIG. 9, a relatively short period is determined as REM sleep, and the actual REM sleep It is likely that it does not match the length.
On the other hand, in the
ここで、本例の場合には、この手法1と手法2を組み合わせて、少なくとも手法1でレム睡眠と判定された区間があり、かつ手法2でレム睡眠と判定される区間rbが一定期間(2分間)以内に存在した場合に、レム期間R1〜R5と確定するようにしたことで、非常に精度の高いレム睡眠期間の判定が行えるようになる。
しかも本例の場合には、手法1の判定を行うための心拍データと、手法2の判定を行うための体動の周波数分布のデータのいずれについても、マットレスセンサ2が取得した被測定者の体動に基づいた圧力データであり、脳波を測定する場合と異なり、被測定者に負担をかけずに、的確なレム睡眠の判定ができるようになる。
Here, in the case of this example, by combining the method 1 and the
Moreover, in the case of this example, both the heartbeat data for determining the method 1 and the frequency distribution data of the body movement for determining the
なお、上述した実施の形態例では、手法1の判定を行うための心拍データと、手法2の判定を行うための体動の周波数分布のデータの双方を、同じセンサ(マットレスセンサ2)の出力データから得るようにしたが、例えば心拍データについては、別のセンサから取得するようにして、手法2の判定を行うための体動の周波数分布のデータのみを、マットレスセンサ2から得るようにしてもよい。
In the above-described embodiment, both the heartbeat data for determining the method 1 and the frequency distribution data of the body movement for determining the
また、上述した実施の形態例で説明した、手法1の判定を行うための心拍の判断期間(ステップS111での5分間)や、手法2の判定を行うための周波数分布の集計期間(ステップS121での30秒間)の長さは、好適な一例を示したものであり、これらの値に限定されるものではない。さらに、上述した実施の形態例で説明した閾値の値についても、上述した値は一例であり、その他の値を設定してもよい。 Further, the heartbeat determination period (5 minutes in step S111) for determining the method 1 and the frequency distribution aggregation period for determining the method 2 (step S121) described in the above-described embodiment example. The length of (30 seconds) in) shows a suitable example and is not limited to these values. Further, as for the threshold value described in the above-described embodiment, the above-mentioned value is an example, and other values may be set.
1…ベッド、2…マットレスセンサ、10…睡眠段階判定装置、11…生体データ取得部、12…生体データ処理部、13…睡眠段階判定部、14…出力部、A…被測定者、C…コンピュータ装置、C1…CPU、C2…ROM、C3…RAM、C4…不揮発性ストレージ、C5…ネットワークインターフェイス表示部、C6…入力装置、C7…表示装置、C8…バス
1 ... bed, 2 ... mattress sensor, 10 ... sleep stage determination device, 11 ... biometric data acquisition unit, 12 ... biometric data processing unit, 13 ... sleep stage determination unit, 14 ... output unit, A ... subject, C ... Computer device, C1 ... CPU, C2 ... ROM, C3 ... RAM, C4 ... Non-volatile storage, C5 ... Network interface display, C6 ... Input device, C7 ... Display device, C8 ... Bus
Claims (5)
前記被測定者の体動を第2の期間毎に周波数解析し、その周波数解析したデータの分析で、レム睡眠の判定を得る決定木を機械学習により複数生成し、生成した複数の決定木で得られたレム睡眠か否かの判定の多数決から、前記第2の期間が第2のレム睡眠期間であると判定するデータ処理部と、
前記第1のレム睡眠期間と判定した期間と前記第2のレム睡眠期間と判定した期間とが重なる期間を含み、かつ一定間隔以内に第2のレム睡眠期間と判定した期間が存在する場合に、その一定間隔以内に第2のレム睡眠期間と判定した期間を、レム睡眠期間と確定する睡眠段階判定部と、を備える
睡眠段階判定装置。 When the rate of increase in the median heart rate for each first period of the subject is equal to or greater than a predetermined threshold, it is determined that the first REM sleep period has started, and the median heart rate for each first period is determined. When the rate of decrease is reduced, the end of the first REM sleep period is determined and
The body movement of the person to be measured is frequency-analyzed every second period, and a plurality of decision trees for obtaining a determination of REM sleep are generated by machine learning by analyzing the frequency-analyzed data, and the generated multiple decision trees are used. From the obtained majority decision on whether or not it is REM sleep, the data processing unit that determines that the second period is the second REM sleep period, and
When the period determined to be the first REM sleep period and the period determined to be the second REM sleep period overlap, and there is a period determined to be the second REM sleep period within a certain interval. , A sleep stage determination device including a sleep stage determination unit that determines a period determined as a second REM sleep period within a certain interval as a REM sleep period.
請求項1に記載の睡眠段階判定装置。 The sleep stage determination device according to claim 1, wherein the sleep stage determination unit determines that the first REM sleep period that does not overlap with the second REM sleep period is also a true REM sleep period.
前記心拍は、前記生体データ取得部が取得した前記圧力センサの出力データから得ると共に、前記第2のレム睡眠期間を判定する体動についても、前記生体データ取得部が取得した前記圧力センサの出力データから得るようにした
請求項1又は2に記載の睡眠段階判定装置。 It is provided with a biological data acquisition unit that acquires output data of a pressure sensor that detects the body movement of the person to be measured.
The heartbeat is obtained from the output data of the pressure sensor acquired by the biometric data acquisition unit, and the output of the pressure sensor acquired by the biometric data acquisition unit is also obtained for the body movement for determining the second REM sleep period. The sleep stage determination device according to claim 1 or 2, which is obtained from the data.
前記被測定者の体動を第2の期間毎に周波数解析したデータの分析で、レム睡眠の判定を得る決定木を機械学習により複数生成し、生成した複数の決定木で得られたレム睡眠か否かの判定の多数決から、前記第2の期間が第2のレム睡眠期間であると判定する第2のレム睡眠期間判定処理と、
前記第1のレム睡眠期間と判定した期間と前記第2のレム睡眠期間と判定した期間とが重なる期間を含み、かつ一定間隔以内に第2のレム睡眠期間と判定した期間が存在する場合に、その一定間隔以内に第2のレム睡眠期間と判定した期間を、レム睡眠期間と確定する睡眠段階判定処理と、を含む
睡眠段階判定方法。 When the rate of increase in the median heart rate for each first period of the subject is equal to or greater than a predetermined threshold value, it is determined that the first REM sleep period has started, and the median heart rate for each first period is determined. When the decrease rate of the first REM sleep period is determined, the first REM sleep period determination process for determining the end of the first REM sleep period and
In the analysis of the data obtained by frequency-analyzing the body movement of the person to be measured for each second period, a plurality of decision trees for obtaining the determination of REM sleep were generated by machine learning, and the REM sleep obtained from the generated plurality of decision trees was generated. From the majority decision of whether or not, the second REM sleep period determination process for determining that the second period is the second REM sleep period, and
When the period determined to be the first REM sleep period and the period determined to be the second REM sleep period overlap, and there is a period determined to be the second REM sleep period within a certain interval. , A sleep stage determination method including a sleep stage determination process for determining a period determined to be a second REM sleep period within a certain interval as a REM sleep period.
前記被測定者の体動を第2の期間毎に周波数解析し、その周波数解析したデータの分析で、レム睡眠の判定を得る決定木を機械学習により複数生成し、生成した複数の決定木で得られたレム睡眠か否かの判定の多数決から、前記第2の期間が第2のレム睡眠期間であると判定する第2のレム睡眠期間判定手順と、
前記第1のレム睡眠期間と判定した期間と前記第2のレム睡眠期間と判定した期間とが重なる期間を含み、かつ一定間隔以内に第2のレム睡眠期間と判定した期間が存在する場合に、その一定間隔以内に第2のレム睡眠期間と判定した期間を、真のレム睡眠期間と判定する睡眠段階判定手順と、
をコンピュータに実行させるプログラム。
When the rate of increase in the median heart rate for each first period of the subject is equal to or greater than a predetermined threshold value, it is determined that the first REM sleep period has started, and the median heart rate for each first period is determined. The first REM sleep period determination procedure for determining the end of the first REM sleep period and the procedure for determining the end of the first REM sleep period when the decrease rate of
The body movement of the person to be measured is frequency-analyzed every second period, and a plurality of decision trees for obtaining a judgment of REM sleep are generated by machine learning by analyzing the frequency-analyzed data, and the generated multiple decision trees are used. From the obtained majority decision on whether or not it is REM sleep, the second REM sleep period determination procedure for determining that the second period is the second REM sleep period, and the second REM sleep period determination procedure.
When the period determined to be the first REM sleep period and the period determined to be the second REM sleep period overlap, and there is a period determined to be the second REM sleep period within a certain interval. , The sleep stage determination procedure for determining the period determined to be the second REM sleep period within the fixed interval as the true REM sleep period,
A program that causes a computer to run.
Applications Claiming Priority (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2017066199 | 2017-03-10 | ||
JP2017066199 | 2017-03-10 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2018149260A JP2018149260A (en) | 2018-09-27 |
JP6932351B2 true JP6932351B2 (en) | 2021-09-08 |
Family
ID=63681154
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2017173563A Active JP6932351B2 (en) | 2017-03-10 | 2017-09-08 | Sleep stage determination device, sleep stage determination method and program |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP6932351B2 (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7216408B2 (en) * | 2019-03-22 | 2023-02-01 | 国立大学法人電気通信大学 | Device for determining sleep apnea syndrome, method for determining sleep apnea syndrome, and program for determining sleep apnea syndrome |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP4551148B2 (en) * | 2004-07-21 | 2010-09-22 | 株式会社ジェピコ | Sleep analyzer |
WO2006121455A1 (en) * | 2005-05-10 | 2006-11-16 | The Salk Institute For Biological Studies | Dynamic signal processing |
JP6213983B2 (en) * | 2013-06-11 | 2017-10-18 | 国立大学法人電気通信大学 | Sleep stage estimation apparatus and method, and program |
JP5583257B2 (en) * | 2013-09-20 | 2014-09-03 | 三菱電機株式会社 | Biological state acquisition device, biological state acquisition program, device provided with biological state acquisition device, and air conditioner |
WO2015119726A2 (en) * | 2014-01-02 | 2015-08-13 | Intel Corporation (A Corporation Of Delaware) | Identifying and characterizing nocturnal motion and stages of sleep |
-
2017
- 2017-09-08 JP JP2017173563A patent/JP6932351B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
JP2018149260A (en) | 2018-09-27 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP6124140B2 (en) | Assessment of patient cognitive function | |
JP5005539B2 (en) | Assessment of sleep quality and sleep-disordered breathing based on cardiopulmonary coupling | |
EP1105812B1 (en) | System and method for facilitating clinical decision making | |
US20140206946A1 (en) | Apparatus and method for measuring stress based on behavior of a user | |
CN106413541B (en) | System and method for diagnosing sleep | |
US20200178875A1 (en) | Predictive neuromarkers of alzheimer's disease | |
JP2016508382A (en) | Detection of sleep apnea using respiratory signals | |
JP6515670B2 (en) | Sleep depth estimation device, sleep depth estimation method, and program | |
US9591972B1 (en) | System and method for evaluating a state of health of a coronary artery | |
KR102236791B1 (en) | System and method for supporting diagnostic for patient based on eeg analysis | |
JP6957011B2 (en) | Sleep stage determination device, sleep stage determination method and program | |
JP6845404B2 (en) | Sleep stage determination method, sleep stage determination device, and sleep stage determination program | |
González et al. | Poincaré plot analysis of cerebral blood flow signals: Feature extraction and classification methods for apnea detection | |
CA2991001A1 (en) | Method and system for predicting of acute hypotensive episodes | |
JP7216408B2 (en) | Device for determining sleep apnea syndrome, method for determining sleep apnea syndrome, and program for determining sleep apnea syndrome | |
US20210315507A1 (en) | System and methods for consciousness evaluation in non-communicating subjects | |
EP3203900A1 (en) | Weaning readiness indicator, sleeping status recording device, and air providing system applying nonlinear time-frequency analysis | |
JP6932351B2 (en) | Sleep stage determination device, sleep stage determination method and program | |
US20230107435A1 (en) | Brain Wave Analysis of Human Brain Cortical Function | |
US10070812B2 (en) | Method for improved seizure detection | |
EP4391902A1 (en) | Fetal heart rate signal processing | |
JP6925056B2 (en) | Sleep stage determination method, sleep stage determination device, and sleep stage determination program | |
KR20200078431A (en) | Method and apparatus for obtaining heart and lung sounds | |
JP7125087B2 (en) | Sleep stage estimation device, sleep stage estimation method and sleep stage estimation program | |
Melcer et al. | Wavelet representation of the corneal pulse for detecting ocular dicrotism |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
AA64 | Notification of invalidation of claim of internal priority (with term) |
Free format text: JAPANESE INTERMEDIATE CODE: A241764 Effective date: 20171017 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20171012 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20200903 |
|
A711 | Notification of change in applicant |
Free format text: JAPANESE INTERMEDIATE CODE: A711 Effective date: 20200909 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A821 Effective date: 20200909 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20210625 |
|
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: 20210720 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20210805 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 6932351 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S111 | Request for change of ownership or part of ownership |
Free format text: JAPANESE INTERMEDIATE CODE: R313117 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |