JP6932351B2 - Sleep stage determination device, sleep stage determination method and program - Google Patents

Sleep stage determination device, sleep stage determination method and program Download PDF

Info

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
Application number
JP2017173563A
Other languages
Japanese (ja)
Other versions
JP2018149260A (en
Inventor
圭樹 ▲高▼玉
圭樹 ▲高▼玉
原田 智広
智広 原田
知里 上原
知里 上原
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
THE UNIVERSITY OF ELECTRO-COMUNICATINS
Tokyo Metropolitan Public University Corp
Original Assignee
THE UNIVERSITY OF ELECTRO-COMUNICATINS
Tokyo Metropolitan Public University Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by THE UNIVERSITY OF ELECTRO-COMUNICATINS, Tokyo Metropolitan Public University Corp filed Critical THE UNIVERSITY OF ELECTRO-COMUNICATINS
Publication of JP2018149260A publication Critical patent/JP2018149260A/en
Application granted granted Critical
Publication of JP6932351B2 publication Critical patent/JP6932351B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

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.

特開2014−239789号公報Japanese Unexamined Patent Publication No. 2014-239789

ところで、上述した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.

本発明の一実施の形態例による睡眠段階判定装置の構成例を示すブロック図である。It is a block diagram which shows the structural example of the sleep stage determination apparatus by one Embodiment of this invention. 本発明の一実施の形態例による睡眠段階の判定状態の例を示す図である。It is a figure which shows the example of the determination state of the sleep stage by the example of one Embodiment of this invention. 本発明の一実施の形態例の睡眠段階判定装置のハードウェア構成例を示すブロック図である。It is a block diagram which shows the hardware configuration example of the sleep stage determination apparatus of one Embodiment of this invention. 本発明の一実施の形態例による心拍成分の判定処理例を示すフローチャートである。It is a flowchart which shows the determination processing example of the heartbeat component by one Embodiment of this invention. 本発明の一実施の形態例による全周波数成分の判定処理例を示すフローチャートである。It is a flowchart which shows the determination processing example of all frequency components by one Embodiment of this invention. 本発明の一実施の形態例によるレム睡眠の判定処理例を示す図である。It is a figure which shows the determination processing example of REM sleep by one Embodiment of this invention. 本発明の一実施の形態例によるレム睡眠の判定処理例を示すフローチャートである。It is a flowchart which shows the determination processing example of REM sleep by one Embodiment of this invention. 本発明の一実施の形態例によるレム睡眠の判定例を示す図である。It is a figure which shows the determination example of REM sleep by one Embodiment of this invention. 本発明の一実施の形態例による実際の検出データからレム睡眠を判定した例を示す図である。It is a figure which shows the example which determined the REM sleep from the actual detection data by one Embodiment of this invention.

以下、本発明の一実施の形態例(以下、「本例」と称する)について、添付図面を参照して説明する。
[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 stage determination device 10 of this example. FIG. 2 is a diagram showing an example of a state in which a sleep stage is determined using the sleep stage determination device 10 of this example.
The sleep stage determination device 10 of this example acquires the body movement of the person to be measured as pressure data by the mattress sensor 2. As shown in FIG. 2, for example, the mattress sensor 2 is laid on the mattress of the bed 1 in which the subject A sleeps, and detects the body movement of the upper body of the subject A during sleep as a change in pressure. .. The mattress sensor 2 is placed on the mattress below the person A to be measured as an example. For example, the mattress sensor 2 may be built in the mattress.
FIG. 2 shows an example in which the sleep stage determination device 10 is installed beside the bed 1 and the mattress sensor 2 and the sleep stage determination device 10 are connected by a cable. For example, the pressure data acquired by the mattress sensor 2 is wirelessly transmitted. It may be transmitted to the sleep stage determination device 10 in another room.

睡眠段階判定装置10は、図1に示すように、生体データ取得部11と、生体データ処理部12と、睡眠段階判定部13と、出力部14とを備える。
生体データ取得部11は、マットレスセンサ2が出力する圧力データを取得する。生体データ取得部11が取得した圧力データは、データ処理部12に供給される。生体データ処理部12は、供給される圧力データの解析処理を行う。具体的には、圧力データをサンプリングしてデジタルデータ化し、そのデジタルデータ化された圧力データの高速フーリエ変換処理(以下、「FFT処理」と称する)により、各周波数の成分を取得する。そして、取得したそれぞれの周波数成分の状態から、心拍と呼吸の状態を判別する。すなわち、圧力データをFFT処理で得た周波数成分の内で、心拍及び呼吸に相当する周波数成分から、心拍数と呼吸数を判別する。
さらに、生体データ処理部12は、FFT処理で得た周波数成分を使って、睡眠段階の解析に必要なデータ処理を行う。
As shown in FIG. 1, the sleep stage determination device 10 includes a biological data acquisition unit 11, a biological data processing unit 12, a sleep stage determination unit 13, and an output unit 14.
The biological data acquisition unit 11 acquires the pressure data output by the mattress sensor 2. The pressure data acquired by the biological data acquisition unit 11 is supplied to the data processing unit 12. The biological data processing unit 12 analyzes the supplied pressure data. Specifically, the pressure data is sampled and converted into digital data, and the components of each frequency are acquired by a high-speed Fourier transform process (hereinafter, referred to as “FFT process”) of the digitalized pressure data. Then, the state of heartbeat and the state of respiration are discriminated from the state of each acquired frequency component. That is, the heart rate and the respiration rate are discriminated from the frequency components corresponding to the heartbeat and the respiration among the frequency components obtained by FFT processing the pressure data.
Further, the biological data processing unit 12 uses the frequency component obtained by the FFT processing to perform data processing necessary for the analysis of the sleep stage.

そして、生体データ処理部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 data processing unit 12 is supplied to the sleep stage determination unit 13. The sleep stage determination unit 13 determines the sleep state of the person to be measured A based on the result analyzed by the biological data processing unit 12. The sleep stage determination unit 13 determines the sleep state of REM sleep and the sleep state other than REM sleep. Further, the sleep stage determination unit 13 determines five stages of awakening and non-REM sleep (stages 1 to 4) as sleep states other than REM sleep. In the following description, the process of determining the sleep state of REM sleep and the sleep state other than REM sleep will be mainly described, and the sleep states other than REM sleep will be divided into 5 stages (awakening and non-REM sleep stages 1 to 1). The description of the process for determining in 4) will be omitted. The sleep stage of non-rem sleep is considered in association with the depth of sleep, and indicates that sleep becomes deeper in the order of stages 1 to 4. For example, stage 2 sleeps deeper than stage 1, and when it is described as "stage 2 or less", it indicates stages 2 to 4.

出力部14は、睡眠段階判定部13が判定した睡眠段階を出力する。出力部14は、例えば表示装置により構成され、睡眠段階を表示する。あるいは、出力部14を記録装置として構成して、一晩の睡眠状態を記録するようにしてもよい。 The output unit 14 outputs the sleep stage determined by the sleep stage determination unit 13. The output unit 14 is composed of, for example, a display device, and displays the sleep stage. Alternatively, the output unit 14 may be configured as a recording device to record the sleep state overnight.

[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 stage determination device 10 is configured by a computer device.
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 data processing unit 12 and the sleep stage determination unit 13 of the sleep stage determination device 10 and executes the program code. Regarding the FFT process for frequency analysis of pressure data, the CPU C1 executes the FFT process by reading the program for executing the corresponding process from the ROM C2. Variables, parameters, etc. generated during the arithmetic processing are temporarily written in the RAM C3.

不揮発性ストレージ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 stage determination device 10 is recorded. Further, the data about the sleep stage determined by the sleep stage determination unit 13 is recorded in the non-volatile storage C4.

ネットワークインターフェイス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 mattress sensor 2 via the network interface C5. The input device C6 is composed of a device such as a keyboard, for example, and the input device C6 sets a period for determining the sleep stage by the sleep stage determination device 10, and instructs the display form of the sleep stage. The display device C7 displays the sleep stage determined by the sleep stage determination device 10.
It should be noted that the sleep stage determination device 10 is configured from a computer device as an example, and data processing such as FFT processing may be performed by preparing dedicated hardware.

[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 stage determination device 10 of this example will be described.
The sleep stage determination device 10 of this example performs REM sleep determination processing (first REM sleep period determination processing and second REM sleep period determination processing) by two types of methods, and REM sleep by the two types of methods. The judgment result of the above is comprehensively judged, and finally the judgment result of the period of REM sleep is obtained. Here, the two types of methods are referred to as method 1 and method 2.
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 method 2 is a REM sleep determination method by machine learning using all frequency components. The sleep stage determination device 10 determines REM sleep with a hybrid structure using these two methods.

まず、手法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 data processing unit 12 acquires a heartbeat value from the data frequency-analyzed by the FFT process. The biological data processing unit 12 calculates the median HR at regular time intervals (here, 5 minutes) for the detected heartbeat (step S111).
Then, the biological data processing unit 12 calculates the fluctuation amount ΔHR every 5 minutes of the detected median heart rate HR, and compares the fluctuation amount ΔHR every 5 minutes with the threshold value TH1 (step S112).
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.

Figure 0006932351
Figure 0006932351

ステップ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 data processing unit 12 waits until the corresponding condition is satisfied.
Then, when it is determined that the condition of ΔHR> 0.04 is satisfied (YES in step S112), the biological data processing unit 12 determines that the REM sleep period (REM sleep period according to the method 1) has started (step S113). After that, the biological data processing unit 12 determines whether or not the reduction rate of the fluctuation amount ΔHR has decreased (step S114). In addition, "decrease rate" means, for example, that the value of ΔHR changes from 0.07 to 0.06. Here, when the decrease rate of the fluctuation amount ΔHR does not decrease (NO in step S114), the biological data processing unit 12 waits until the decrease rate of the fluctuation amount ΔHR decreases. Further, when the decrease rate of the fluctuation amount ΔHR decreases (YES in step S114), the biological data processing unit 12 determines the end of the REM sleep period (REM sleep period according to the method 1) (step S115), and in step S112. Returning to the comparison between the fluctuation amount ΔHR and the threshold value TH1.
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 method 2 will be described.
This method 2 is a method of determining REM sleep and non-REM sleep by using Random Forests, which is a machine learning method based on ensemble learning. In this method 2, the pressure data acquired by the mattress sensor 2 for a certain period (here, 30 seconds) is subjected to FFT processing to perform frequency analysis to determine REM sleep.

図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 method 2.
First, the biometric data processing unit 12 divides the biometric data acquired by the biometric data acquisition unit 11 at regular intervals (1 epoch section: 30 seconds), and performs frequency analysis on the biometric data at regular intervals by FFT processing. By doing so, biometric data expressed in frequency is generated (step S121). Then, all frequency components of the FFT-processed biological data are normalized so that the sum of all frequency components becomes 1. Note that normalization is performed here as an example, and data that is not normalized may be handled.
Then, the biological data processing unit 12 uses the pressure data at regular intervals to perform binary determination by a machine learning method of random forest (step S122). The details of binary judgment by machine learning of random forest will be described later.

その後、生体データ処理部12は、ステップS122による2値判定で、レム睡眠の条件を満たすか否かを判断する(ステップS123)。この判断でレム睡眠の条件を満たさない場合には(ステップS123のNO)、ステップS121の処理に戻り、生体データ処理部12は、次の一定期間の処理に移る。
そして、この判断でレム睡眠の条件を満たした場合には(ステップS123のYES)、生体データ処理部12は、該当する一定期間(30秒間)をレム睡眠期間(手法2によるレム睡眠期間)に設定し(ステップS124)、その後ステップS121の処理に戻り、次の一定期間の処理に移る。
After that, the biological data processing unit 12 determines whether or not the condition of REM sleep is satisfied by the binary determination in step S122 (step S123). If the condition of REM sleep is not satisfied by this determination (NO in step S123), the process returns to the process of step S121, and the biometric data processing unit 12 moves to the process for the next fixed period.
Then, when the condition of REM sleep is satisfied by this determination (YES in step S123), the biological data processing unit 12 sets the corresponding fixed period (30 seconds) into the REM sleep period (REM sleep period according to the method 2). The setting (step S124) is performed, and then the process returns to the process of step S121, and the process proceeds to the next fixed period process.

図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 data processing unit 12 acquires the pressure data D1 for 30 seconds. The vertical axis of the graph of the pressure data D1 is the sensor value, and the horizontal axis is the time (30 seconds).
The biometric data processing unit 12 performs FFT processing on the pressure data D1 to obtain FFT data D2 showing the intensity value for each frequency component. The vertical axis of the graph of FFT data D2 is the intensity value, and the horizontal axis is the frequency (0.016 Hz to 5.0 Hz). Here, the FFT data D2 is 300-dimensional vector data containing only frequency components of 5 Hz or less.

生体データ処理部12は、この300次元ベクトルのデータから、ランダムフォレストと称される機械学習のアルゴリズムを使って、予測モデルである決定木を複数生成する。すなわち、図6に示すように、複数の決定木T,T,・・・,T(Nは任意の整数)を生成する。それぞれの決定木は、各周波数成分の状態がどのようであるかに基づいて分岐して、各睡眠段階(ここではレム睡眠段階)の判定を得るものである。例えば、図6の左下に示すように、決定木Tのある段階Saでは、周波数成分0.3Hzの正規化値が0.25より小さいか否かを判断し、その判断で、正規化値が0.25より小さいとき未確定と判定し、0.25以上のとき、レム睡眠と判定する。そして、その未確定の場合(0.3Hzの成分の正規化値が0.25より小さいとき)に、次の段階Sbで、周波数成分1.5Hzの正規化値が0.05より大きいか否かを判断し、その判断で、正規化値が0.05より大きいとき、ノンレム睡眠と判定し、正規化値が0.05以下のとき、レム睡眠と判定する。 The biological data processing unit 12 generates a plurality of decision trees, which are prediction models, from the data of the 300-dimensional vector by using a machine learning algorithm called a random forest. That is, as shown in FIG. 6, a plurality of decision trees T 1 , T 2 , ..., TN (N is an arbitrary integer) are generated. Each decision tree branches based on what the state of each frequency component is, and obtains the judgment of each sleep stage (here, the REM sleep stage). For example, as shown in the lower left of FIG. 6, determining the tree T 1 of a certain stage Sa, normalized values of the frequency components 0.3Hz will determine 0.25 or less, in the determination, normalization value When is less than 0.25, it is determined to be unconfirmed, and when it is 0.25 or more, it is determined to be REM sleep. Then, in the undetermined case (when the normalized value of the 0.3 Hz component is smaller than 0.25), whether or not the normalized value of the frequency component 1.5 Hz is greater than 0.05 in the next step Sb. When the normalized value is greater than 0.05, it is determined to be non-REM sleep, and when the normalized value is 0.05 or less, it is determined to be REM sleep.

このようにして、決定木Tに基づいてFFTデータD2を解析することで、レム睡眠か否かの判定qを得る。また、決定木Tとは別の条件の決定木Tに基づいてFFTデータD2を解析することで、レム睡眠か否かの判定qを得る。このようにして、それぞれ異なる条件の決定木T,T,・・・,T毎に、レム睡眠か否かの判定q,q,・・・,qを得る。
そして、生体データ処理部12は、複数の決定木T,T,・・・,T毎に得たレム睡眠か否かの判定q,q,・・・,qの多数決で、現在の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 data processing unit 12, a plurality of decision trees T 1, T 2, · · ·, determines whether REM sleep obtained for each T N q 1, q 2, ···, the q N majority Then, it is determined whether or not the current period of 30 seconds is REM sleep.

図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 method 2. The process of determining the REM sleep period is determined by the sleep stage determination unit 13.
First, the sleep stage determination unit 13 determines whether or not there is a section determined to be REM sleep by the method 1 (step S131). In this determination, if there is no section determined to be REM sleep by the method 1 (NO in step S131), the sleep stage determination unit 13 determines that there is no REM sleep period within the searched period, and determines that there is no REM sleep period. The process of searching for a period ends.
Then, when there is a section determined to be REM sleep by the method 1 (YES in step S131), the sleep stage determination unit 13 determines the section determined to be REM sleep by the method 1 and the section determined to be REM sleep by the method 2. It is determined whether or not there is a section that overlaps with (step S132).

さらに、ステップ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 method 2 overlap in the judgment of step S132 (YES in step S132), the REM sleep is determined by the method 1. All the periods including the determined section and the section determined to be REM sleep by the method 2 exist within a certain period (here, within 2 minutes) are determined as the REM sleep period (step S133). If the section determined to be REM sleep by method 1 and the section determined to be REM sleep by method 2 do not overlap in the determination of step S132 (NO in step S132), it is determined to be REM sleep by method 1. The section is determined as the REM sleep period (step S134).

図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 method 2 shown in the flowchart of FIG. 7.
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 method 2 in the period overlapping the section. In the case of the state shown in FIG. 8A, the sleep stage determination unit 13 determines the section determined to be REM sleep by the method 1 as the REM sleep section. The confirmation process of FIG. 8A corresponds to the process of determining REM sleep in step S134 of the flowchart of FIG. 7.

図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 method 2 in the period overlapping the section. In this case, the sleep stage determination unit 13 determines all the sections determined to be REM sleep by the determination by the method 1 as the REM sleep period, and also determines the REM sleep period for the sections determined to be REM sleep by the determination by the method 2. Is confirmed. The section determined to be REM sleep by the determination by the method 2 here is defined as a continuous REM sleep period during which there is a section determined to be REM sleep within 2 minutes by the method 2. The confirmation process of FIG. 8B corresponds to the process of determining REM sleep in step S133 of the flowchart of FIG. 7.

図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 method 2 in a certain period, and the corresponding section is not determined to be REM sleep in the method 1. In this case, the sleep stage determination unit 13 does not determine the corresponding section as REM sleep.

図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 stage determination device 10 of this example). .. On the vertical axis of the graph of this characteristic BW, the positions indicated by A indicate arousal, and the positions indicated by the numbers 1, 2, 3 and 4 indicate stages 1 to 4 of non-REM sleep. In the characteristic BW shown in FIG. 9, the period of REM sleep is not determined, but REM sleep exists between stage 1 of non-REM sleep and awakening.

ここで、図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 stage determination device 10 of this example from the determination results of both method 1 and method 2. be.

図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 stage determination device 10 based on the pressure data detected by the mattress sensor 2 shown in FIG. Within the characteristic HR showing the fluctuation of the heart rate, the positions of the peaks that are temporarily higher than the heart rate before and after the heart rate ra1, ra2, ra3, ra4, ra5 (almost the central part of the part indicated by the broken line circle in FIG. 9). ) Corresponds to the section determined to be REM sleep by the method 1.
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 method 2. In FIG. 9, the section rb determined to be REM sleep by the method 2 is shown only in the vicinity of the REM sleep periods R4 and R5, and the other sections are not shown.
In the case of this example, one section rb determined to be REM sleep by the method 2 is 30 seconds, and when the 30-second section rb exists within 2 minutes, the REM sleep determined by the method 2 is continuous. Assuming that, one REM sleep period (for example, period R4 or R5) is determined. However, in the case of this example, a portion where the period in which the 30-second section rb exists within 2 minutes does not overlap with the sections ra1 to ra5 determined to be REM sleep by the method 1 is determined not to be REM sleep. do.

このようにして、心拍に基づいてレム睡眠を判定する手法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 method 2 for determining REM sleep based on the machine learning of the frequency distribution of the pressure data detected by the mattress sensor 2 are combined to obtain REM sleep. By determining the period, it becomes possible to determine REM sleep with extremely high accuracy.
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 method 2 for determining REM sleep from the machine learning of the frequency distribution of the detected data, for example, even in the period between the section R4 and the section R5 in FIG. There was a section rb to be determined, and there was a possibility that the accuracy of determining REM sleep was not so high.

ここで、本例の場合には、この手法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 method 2, there is at least a section determined to be REM sleep by the method 1, and the section rb determined to be REM sleep by the method 2 is a certain period ( By determining the REM period R1 to R5 when it exists within (2 minutes), it becomes possible to determine the REM sleep period with extremely high accuracy.
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 method 2 are obtained by the mattress sensor 2 of the person to be measured. It is pressure data based on body movement, and unlike the case of measuring brain waves, it becomes possible to accurately determine REM sleep without imposing a burden on the person to be measured.

なお、上述した実施の形態例では、手法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 method 2 are output from the same sensor (mattress sensor 2). Although it was obtained from the data, for example, the heartbeat data was acquired from another sensor so that only the data of the frequency distribution of the body movement for determining the method 2 was obtained from the mattress sensor 2. May be good.

また、上述した実施の形態例で説明した、手法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)

被測定者の第1の期間毎の心拍の中央値の増加率が、所定の閾値以上であるとき、第1のレム睡眠期間の開始と判定し、前記第1の期間毎の心拍の中央値の減少率が減少したとき、前記第1のレム睡眠期間の終了を判定すると共に、
前記被測定者の体動を第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.
さらに、前記睡眠段階判定部は、前記第2のレム睡眠期間と重なっていない前記第1のレム睡眠期間についても、真のレム睡眠期間と判定する
請求項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.
被測定者の第1の期間毎の心拍の中央値の増加率が、所定の閾値以上であるとき、第1のレム睡眠期間の開始と判定し、前記第1の期間毎の心拍の中央値の減少率が減少したとき、前記第1のレム睡眠期間の終了を判定する第1のレム睡眠期間判定処理と、
前記被測定者の体動を第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.
被測定者の第1の期間毎の心拍の中央値の増加率が、所定の閾値以上であるとき、第1のレム睡眠期間の開始と判定し、前記第1の期間毎の心拍の中央値の減少率が減少したとき、前記第1のレム睡眠期間の終了を判定する第1のレム睡眠期間判定手順と、
前記被測定者の体動を第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.
JP2017173563A 2017-03-10 2017-09-08 Sleep stage determination device, sleep stage determination method and program Active JP6932351B2 (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

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