ホーム > 特許ランキング > 国立大学法人電気通信大学
特許7216408睡眠時無呼吸症候群判定装置、睡眠時無呼吸症候群判定方法および睡眠時無呼吸症候群判定プログラム
(19)【発行国】日本国特許庁(JP)
(12)【公報種別】特許公報(B2)
(11)【特許番号】
(24)【登録日】2023-01-24
(45)【発行日】2023-02-01
(54)【発明の名称】睡眠時無呼吸症候群判定装置、睡眠時無呼吸症候群判定方法および睡眠時無呼吸症候群判定プログラム
(51)【国際特許分類】
A61B 5/16 20060101AFI20230125BHJP
A61B 5/11 20060101ALI20230125BHJP
A61B 5/113 20060101ALI20230125BHJP
A61B 5/08 20060101ALI20230125BHJP
【ＦＩ】
A61B5/16 130
A61B5/11 100
A61B5/113
A61B5/08
【請求項の数】
8
(21)【出願番号】P 2019055132
(22)【出願日】2019-03-22
(65)【公開番号】P2020151358
(43)【公開日】2020-09-24
【審査請求日】2022-03-22
(73)【特許権者】
【識別番号】504133110
【氏名又は名称】国立大学法人電気通信大学
(74)【代理人】
【識別番号】110000925
【氏名又は名称】弁理士法人信友国際特許事務所
(72)【発明者】
【氏名】▲高▼玉 圭樹
(72)【発明者】
【氏名】中理 怡恒
【審査官】湊 和也
(56)【参考文献】
【文献】国際公開第２０１８／２２１７５０（ＷＯ，Ａ１）
【文献】特開２０１８－１４９２６０（ＪＰ，Ａ）
【文献】特開２０１８－１４９２５９（ＪＰ，Ａ）
(58)【調査した分野】(Int.Cl.，ＤＢ名)
Ａ６１Ｂ ５／１６－５／１８
Ａ６１Ｂ ５／１１－５／１１３
Ａ６１Ｂ ５／０８－５／０９７
(57)【特許請求の範囲】
【請求項1】
被測定者の睡眠中の体動または圧力変化に基づいた生体振動データを得る生体データ取得部と、前記生体データ取得部が取得した前記生体振動データから、一定時間ごとに複数の特徴量を算出する生体データ処理部と、
前記生体データ処理部が算出した複数の特徴量から、睡眠段階を判別するための分岐条件を持った複数の決定木を用いて睡眠段階を判別した後、その判別に使用された決定木の分類基準の組み合わせで得られる判別結果に基づいて、睡眠時無呼吸症候群の判定を行う判定部と、を備える
睡眠時無呼吸症候群判定装置。
【請求項2】
前記複数の決定木を学習する学習部を備え、前記学習部による複数の決定木の組み合わせの学習は、正解となる睡眠段階の情報を取得し、正解と一致度が高い睡眠段階の判別ができる決定木を選ぶようにした
請求項１に記載の睡眠時無呼吸症候群判定装置。
【請求項3】
前記判定部が睡眠時無呼吸症候群の判定を行う際には、前記複数の決定木で分岐する際の、特徴量の組み合わせの種類数と、各特徴量の出現回数を使用する請求項２に記載の睡眠時無呼吸症候群判定装置。
【請求項4】
被測定者の睡眠中の体動または圧力変化に基づいた生体振動データを得る生体データ取得部と、前記生体データ取得部が取得した前記生体振動データから、一定時間ごとに前記生体振動データに含まれる周波数成分のパワースペクトルを算出する生体データ処理部と、
前記生体データ処理部が算出した周波数成分のパワースペクトルから、睡眠段階を判定するための、それぞれ異なる分岐条件を持った複数の決定木を用いて睡眠段階を判別した後、その判別に使用された周波数成分のパワースペクトルの寄与度を算出し、算出した寄与度に基づいて睡眠時無呼吸症候群の判定を行う判定部と、を備える
睡眠時無呼吸症候群判定装置。
【請求項5】
被測定者の睡眠中の体動または圧力変化に基づいた生体振動データを得る生体データ取得処理と、前記生体データ取得処理により取得した前記生体振動データから、一定時間ごとに複数の特徴量を算出する生体データ処理と、
前記生体データ処理で算出した複数の特徴量から、睡眠段階を判別するための分岐条件を持った複数の決定木を用いて睡眠段階を判別した後、その判断に使用された決定木の分類基準の組み合わせで得られる判別結果に基づいて、睡眠時無呼吸症候群の判定を行う判定処理と、を含む
睡眠時無呼吸症候群判定方法。
【請求項6】
被測定者の睡眠中の体動または圧力変化に基づいた生体振動データを得る生体データ取得処理と、前記生体データ取得処理により取得した前記生体振動データから、一定時間ごとに前記生体振動データに含まれる周波数成分のパワースペクトルを算出する生体データ処理と、
前記生体データ処理により算出した周波数成分のパワースペクトルから、睡眠段階を判定するための、それぞれ異なる分岐条件を持った複数の決定木を用いて睡眠段階を判別した後、その判別に使用された周波数成分のパワースペクトルの寄与度を算出し、算出した寄与度に基づいて睡眠時無呼吸症候群の判定を行う判定処理と、を含む
睡眠時無呼吸症候群判定方法。
【請求項7】
被測定者の睡眠中の体動または圧力変化に基づいた生体振動データを得る生体データ取得手順と、前記生体データ取得手順により取得した前記生体振動データから、一定時間ごとに複数の特徴量を算出する生体データ処理手順と、
前記生体データ処理手順で算出した複数の特徴量から、睡眠段階を判別するための分岐条件を持った複数の決定木を用いて睡眠段階を判別した後、その判断に使用された決定木の分類基準の組み合わせで得られる判別結果に基づいて、睡眠時無呼吸症候群の判定を行う判定手順と、
をコンピュータに実行させる睡眠時無呼吸症候群判定プログラム。
【請求項8】
被測定者の睡眠中の体動または圧力変化に基づいた生体振動データを得る生体データ取得手順と、前記生体データ取得手順により取得した前記生体振動データから、一定時間ごとに前記生体振動データに含まれる周波数成分のパワースペクトルを算出する生体データ処理手順と、
前記生体データ処理手順により算出した周波数成分のパワースペクトルから、睡眠段階を判定するための、それぞれ異なる分岐条件を持った複数の決定木を用いて睡眠段階を判別した後、その判別に使用された周波数成分のパワースペクトルの寄与度を算出し、算出した寄与度に基づいて睡眠時無呼吸症候群の判定を行う判定手順と、
をコンピュータに実行させる睡眠時無呼吸症候群判定プログラム。
【発明の詳細な説明】
【技術分野】
【0001】
本発明は、被験者の睡眠時無呼吸症候群を判定する睡眠時無呼吸症候群判定装置、睡眠時無呼吸症候群判定方法および睡眠時無呼吸症候群判定プログラムに関する。
【背景技術】
【0002】
医療現場では、睡眠障害や睡眠時無呼吸症候群を診断するために、被測定者の睡眠状態を測定することが行われている。人間の睡眠段階は、睡眠の深さの観点で６段階に分類したものが知られており、その６つの睡眠段階は、眠りが浅い段階から順に、覚醒、レム睡眠、ノンレム睡眠（ステージ１～４）と呼ばれている。これらの６段階の睡眠段階の判定は、従来、例えば被計測者の顔や頭部に多数の電極を装着して、その多数の電極から脳波、眼球運動、および顎筋電を測定し、測定結果の解析により行われていた。
【0003】
また、睡眠時無呼吸症候群の患者が睡眠時に無呼吸になるのは、睡眠段階が覚醒の状態であることが多く、無呼吸症候群の診断を行う上でも、睡眠段階を測定する必要がある。
但し、睡眠時無呼吸症候群の判定を行うためには、睡眠段階の測定の他に、口と鼻の気流などの呼吸に伴った空気の流れの測定や、胸部および腹部の換気運動などの様々な測定を同時に行う必要がある。そして、睡眠段階の解析結果と呼吸状態の測定結果などに基づいて、医師が無呼吸症候群であるか否かを診断している。
但し、睡眠時無呼吸症候群の判定を行うためには、睡眠段階の測定の他に、口と鼻の気流などの呼吸に伴った空気の流れの測定や、胸部および腹部の換気運動などの様々な測定を同時に行う必要がある。そして、睡眠段階の解析結果と呼吸状態の測定結果などに基づいて、医師が無呼吸症候群であるか否かを診断している。
【0004】
このような診断を行うために必要な、顔や頭部に多数の電極を装着した状態での睡眠の検査は、通常、医療機関に宿泊して、長時間連続して電極を身体に装着して行う検査であり、被測定者（患者）に精神的な負担と肉体的な負担を強いることになる。また、取得したデータは、専門知識と経験を持つ医師が解析して判定する必要があり、睡眠時無呼吸症候群を簡単に判定できるものではないという問題があった。
【0005】
睡眠段階の測定に関する問題を解決するために、専門医師による診断を不要とする睡眠段階推定手法は、従来から数多く提案されている。
例えば、特許文献１には、遺伝的アルゴリズムによる学習手法を改良したDatabase-based Compact Genetic Algorithmと称される手法で、マットレス型圧力センサの検出データから睡眠段階を推定する技術が記載されている。この特許文献１に記載された技術は、マットレス型圧力センサが検出した被測定者の体動と心拍に基づいて、睡眠段階を推定するものである。このようなマットレス型圧力センサを使って睡眠段階を推定することで、被測定者に負担を強いることなく、被測定者の睡眠状態を推定することができる。
例えば、特許文献１には、遺伝的アルゴリズムによる学習手法を改良したDatabase-based Compact Genetic Algorithmと称される手法で、マットレス型圧力センサの検出データから睡眠段階を推定する技術が記載されている。この特許文献１に記載された技術は、マットレス型圧力センサが検出した被測定者の体動と心拍に基づいて、睡眠段階を推定するものである。このようなマットレス型圧力センサを使って睡眠段階を推定することで、被測定者に負担を強いることなく、被測定者の睡眠状態を推定することができる。
【先行技術文献】
【特許文献】
【0006】
【文献】特開２０１４－２３９７８９号公報
【発明の概要】
【発明が解決しようとする課題】
【0007】
上述したように、無呼吸症候群の患者が睡眠時に無呼吸になる状態は、睡眠段階が覚醒の状態であることが多いため、マットレス型圧力センサの検出データ、つまり体動の検出データを使って睡眠段階を検出することは、無呼吸症候群の判定の一つの指標になる。
【0008】
マットレス型圧力センサの検出データから睡眠段階を検出する際には、例えばマットレス型圧力センサの検出データから、特定の閾値を超える体動を検出したとき、睡眠段階が覚醒状態であると判別する処理を行うことが考えられる。正常者の睡眠の場合、このような体動の閾値超えの判別を行うことで、睡眠段階が覚醒か否かを比較的高い精度で判別することができる。
【0009】
しかしながら、睡眠時無呼吸症候群の患者の場合、体動の変化パターンが正常者の睡眠時とは全く異なる場合があり、正常者の睡眠段階と同じ判別処理を行っても、覚醒と判別できないケースが多々発生する。具体的には、睡眠時無呼吸症候群の患者が無呼吸状態の場合には、睡眠段階が覚醒状態であっても、体幹が動くような大きな体動が殆どない状態で無呼吸となることがあり、閾値を用いた体動の判定では覚醒と検出できないケースが発生してしまう。
したがって、従来提案されているマットレス型圧力センサの検出データから睡眠段階を検出する処理を、睡眠時無呼吸症候群の判定用にそのまま適用することは困難であった。
したがって、従来提案されているマットレス型圧力センサの検出データから睡眠段階を検出する処理を、睡眠時無呼吸症候群の判定用にそのまま適用することは困難であった。
【0010】
本発明は、被測定者への負担をかけずに睡眠時無呼吸症候群の判定を行うことが可能な無呼吸症候群判定装置、無呼吸症候群判定方法および無呼吸症候群判定プログラムを提供することを目的とする。
【課題を解決するための手段】
【0011】
第１の発明の睡眠時無呼吸症候群判定装置は、被測定者の睡眠中の体動または圧力変化に基づいた生体振動データを得る生体データ取得部と、生体データ取得部が取得した生体振動データから、一定時間ごとに複数の特徴量を算出する生体データ処理部と、生体データ処理部が算出した複数の特徴量から、睡眠段階を判別するための分岐条件を持った複数の決定木を用いて睡眠段階を判別した後、その判断に使用された決定木の分類基準の組み合わせで得られる判別結果に基づいて、睡眠時無呼吸症候群の判定を行う判定部と、を備える。
【0012】
また、第１の発明の睡眠時無呼吸症候群判定方法は、被測定者の睡眠中の体動または圧力変化に基づいた生体振動データを得る生体データ取得処理と、生体データ取得処理により取得した生体振動データから、一定時間ごとに複数の特徴量を算出する生体データ処理と、生体データ処理で算出した複数の特徴量から、睡眠段階を判別するための分岐条件を持った複数の決定木を用いて睡眠段階を判別した後、その判断に使用された決定木の分類基準の組み合わせで得られる判別結果に基づいて、睡眠時無呼吸症候群の判定を行う判定処理と、を含む。
【0013】
また、第１の発明の睡眠時無呼吸症候群判定プログラムは、上記の第１の睡眠時無呼吸症候群判定方法が行う各処理を手順としてコンピュータに実行させるものである。
【0014】
第２の発明の睡眠時無呼吸症候群判定装置は、被測定者の睡眠中の体動または圧力変化に基づいた生体振動データを得る生体データ取得部と、生体データ取得部が取得した生体振動データから、一定時間ごとに生体振動データに含まれる周波数成分のパワースペクトルを算出する生体データ処理部と、生体データ処理部が算出した周波数成分のパワースペクトルから、睡眠段階を判定するための、それぞれ異なる分岐条件を持った複数の決定木を用いて睡眠段階を判別した後、その判別に使用された周波数成分のパワースペクトルの寄与度を算出し、算出した寄与度に基づいて睡眠時無呼吸症候群の判定を行う判定部と、を備える。
【0015】
また、第２の発明の睡眠時無呼吸症候群判定方法は、被測定者の睡眠中の体動または圧力変化に基づいた生体振動データを得る生体データ取得処理と、生体データ取得処理により取得した生体振動データから、一定時間ごとに生体振動データに含まれる周波数成分のパワースペクトルを算出する生体データ処理と、生体データ処理により算出した周波数成分のパワースペクトルから、睡眠段階を判定するための、それぞれ異なる分岐条件を持った複数の決定木を用いて睡眠段階を判別した後、その判別に使用された周波数成分のパワースペクトルの寄与度を算出し、算出した寄与度に基づいて睡眠時無呼吸症候群の判定を行う判定処理と、を含む。
【0016】
また、第２の発明の睡眠時無呼吸症候群判定プログラムは、上記の第２の発明の睡眠時無呼吸症候群判定方法が行う各処理を手順としてコンピュータに実行させるものである。
【発明の効果】
【0017】
本発明によると、被測定者の睡眠中の体動または圧力変化に基づいた生体振動データから、精度の高い睡眠時無呼吸症候群の判定ができるようになる。したがって、マットレス型などの圧力センサを敷いた状態で睡眠するだけで、簡単に精度の高い睡眠時無呼吸症候群の判定ができるようになる。
【図面の簡単な説明】
【0018】
【図1】本発明の第１の実施の形態例による睡眠時無呼吸症候群判定装置の構成例を示すブロック図である。
【図2】本発明の第１の実施の形態例による睡眠時無呼吸症候群の判定状態の例を示す図である。
【図3】本発明の第１の実施の形態例の睡眠時無呼吸症候群判定装置のハードウェア構成例を示すブロック図である。
【図4】睡眠段階の変化例を、健常者（ａ）と睡眠時無呼吸症候群の患者（ｂ）とで比較した図である。
【図5】の覚醒（ＷＡＫＥ）時の特性例を、健常者（ａ）と睡眠時無呼吸症候群の患者（ｂ）とで比較した図である。
【図6】本発明の第１の実施の形態例による睡眠時無呼吸症候群を判定するための学習フェーズ（ａ）および判定フェーズ（ｂ）の流れを示すフローチャートである。
【図7】睡眠時無呼吸症候群の患者の生体振動データの例（ａ）と、その一部を拡大したもの（ｂ）を示す図である。
【図8】本発明の第１の実施の形態例による決定木の例を示す図である。
【図9】本発明の第１の実施の形態例による生体振動データから特徴量を示すデータを得るまでの処理状態の例を示す図である。
【図10】本発明の第１の実施の形態例による特徴量を示すデータからランダムフォレストモデルを使った処理までの例を示す図である。
【図11】本発明の第１の実施の形態例によるランダムフォレストモデルの判別特徴量のカウント処理の例を示す図である。
【図12】本発明の第１の実施の形態例による判定結果の例を示す図である。
【図13】本発明の第２の実施の形態例による睡眠時無呼吸症候群判定処理を判定するための学習フェーズ（ａ）および判定フェーズ（ｂ）の流れを示すフローチャートである。
【図14】本発明の第２の実施の形態例による判定処理の原理を説明するための、睡眠時無呼吸症候群の患者（ａ）と健常者（ｂ）とで周波数成分の寄与度を比較した図である。
【図15】本発明の第２の実施の形態例による生体振動データの取得からランダムフォレストモデルの構築までの例を示す図である。
【図16】本発明の第２の実施の形態例によるランダムフォレストモデルから各周波数の寄与度を算出する処理の例を示す図である。
【図17】本発明の第２の実施の形態例による周波数分布を分割する例を示す図である。
【図18】本発明の第２の実施の形態例による指標の定義を示す図である。
【図19】本発明の第２の実施の形態例による判定結果の例を示す図である。
【発明を実施するための形態】
【0019】
＜１．第１の実施の形態例＞
以下、本発明の第１の実施の形態例について、図１～図１２を参照して説明する。
［１－１．睡眠時無呼吸症候群判定装置の構成］
図１は、本実施の形態例の睡眠時無呼吸症候群判定装置１０の構成を示すブロック図である。
図２は、本実施の形態例の睡眠時無呼吸症候群判定装置１０を使って睡眠時無呼吸症候群の判定を行う状態の例を示す図である。
以下、本発明の第１の実施の形態例について、図１～図１２を参照して説明する。
［１－１．睡眠時無呼吸症候群判定装置の構成］
図１は、本実施の形態例の睡眠時無呼吸症候群判定装置１０の構成を示すブロック図である。
図２は、本実施の形態例の睡眠時無呼吸症候群判定装置１０を使って睡眠時無呼吸症候群の判定を行う状態の例を示す図である。
【0020】
本実施の形態例の睡眠時無呼吸症候群判定装置１０は、被測定者の生体振動（体動）をマットレスセンサ２で圧力データとして取得する。マットレスセンサ２は、例えば図２に示すように、被測定者Ａが睡眠を行うベッド１のマットレスの上あるいは下に敷いて、被測定者Ａの睡眠中の上半身の生体振動を圧力の変化として検出する。被測定者Ａの下側になるマットレスの上にマットレスセンサ２を配置するのは一例であり、例えばマットレスの中にマットレスセンサ２を内蔵させてもよい。
図２では、ベッド１の脇に睡眠時無呼吸症候群判定装置１０を設置し、マットレスセンサ２と睡眠時無呼吸症候群判定装置１０をケーブルで接続した例を示すが、例えばマットレスセンサ２が取得した圧力データ（生体振動データ）を、無線伝送で別の部屋の睡眠時無呼吸症候群判定装置１０に伝送するようにしてもよい。
なお、以下の説明では、マットレスセンサ２が出力する圧力データを、生体振動データと称する。
図２では、ベッド１の脇に睡眠時無呼吸症候群判定装置１０を設置し、マットレスセンサ２と睡眠時無呼吸症候群判定装置１０をケーブルで接続した例を示すが、例えばマットレスセンサ２が取得した圧力データ（生体振動データ）を、無線伝送で別の部屋の睡眠時無呼吸症候群判定装置１０に伝送するようにしてもよい。
なお、以下の説明では、マットレスセンサ２が出力する圧力データを、生体振動データと称する。
【0021】
睡眠時無呼吸症候群判定装置１０は、図１に示すように、生体データ取得部１１と、生体データ処理部１２と、判定部１３と、出力部１４と、学習部１５とを備える。
生体データ取得部１１は、マットレスセンサ２が出力する生体振動データを取得する生体データ取得処理を行う。生体データ取得部１１が取得した生体振動データは、生体データ処理部１２に供給される。
生体データ取得部１１は、マットレスセンサ２が出力する生体振動データを取得する生体データ取得処理を行う。生体データ取得部１１が取得した生体振動データは、生体データ処理部１２に供給される。
【0022】
生体データ処理部１２は、供給される生体振動データをサンプリングしてデジタルデータ化し、そのデジタルデータ化された生体振動データの所定時間ごとの平均値を算出する。例えば、生体データ処理部１２は、生体振動データの平均値を１秒ごとに取得する。それぞれの平均値は、同じ期間のデータに対して平均化したものである。
【0023】
そして、生体データ処理部１２は、所定時間ごとに得られた平均値を演算処理して、複数種類の特徴量を算出する。例えば、１秒ごとの生体振動データの平均値ｘから、以下に示す複数の特徴量ＳＤ、Range、ＳＵＭ、Square、ＬＣ、ＲＭＳを、３０秒ごとに算出する。以下の説明では、生体データ処理部１２は、１秒ごとの生体振動データの平均値ｘから、３０秒ごとに特徴量を求めるものとする。
【0024】
各特徴量について順に説明する。
特徴量ＳＤは、１秒ごとの生体振動データの平均値ｘの３０秒間の分散を示す。
特徴量Rangeは、３０秒間の平均値ｘの最大値と最小値の差を示す。
特徴量ＳＵＭは、３０秒間の平均値ｘの積算値である。
特徴量Squareは、３０秒間の平均値ｘの二乗値の積算値である。
特徴量ＬＣは、１秒ごとの生体振動データの平均値ｘｎと、その平均値ｘｎの１つ前（１秒前）の平均値ｘｎ－１との差分の、３０秒間の積算値である。
特徴量ＲＭＳは、３０秒間の平均値ｘの二乗値の積算値を３０秒間で平均してルートを求めたものである。
特徴量ＳＤは、１秒ごとの生体振動データの平均値ｘの３０秒間の分散を示す。
特徴量Rangeは、３０秒間の平均値ｘの最大値と最小値の差を示す。
特徴量ＳＵＭは、３０秒間の平均値ｘの積算値である。
特徴量Squareは、３０秒間の平均値ｘの二乗値の積算値である。
特徴量ＬＣは、１秒ごとの生体振動データの平均値ｘｎと、その平均値ｘｎの１つ前（１秒前）の平均値ｘｎ－１との差分の、３０秒間の積算値である。
特徴量ＲＭＳは、３０秒間の平均値ｘの二乗値の積算値を３０秒間で平均してルートを求めたものである。
【0025】
生体データ処理部１２で３０秒ごとに得られた、これらの特徴量ＳＤ、Range、ＳＵＭ、Square、ＬＣ、ＲＭＳは、判定部１３および学習部１５に供給される。
ここで、生体データ処理部１２では、３０秒ごとの複数の特徴量ＳＤ、Range、ＳＵＭ、Square、ＬＣ、ＲＭＳを示すデータに対して、睡眠段階を付与するラベリング処理が行われ、そのラベリング処理が行われた特徴量のデータが、判定部１３および学習部１５に供給される。
ここで付与される睡眠段階は、既に知られた睡眠段階推定手法を適用して、生体振動データから取得したものである。例えば、生体振動データの特定の周波数成分のレベルと閾値とを比較して、判断した睡眠段階を、その期間の特徴量のデータに対して付与する。なお、ここで生体振動データから睡眠段階を判断する処理は、睡眠時無呼吸症候群の患者であることを考慮しない、いわゆる健常者用の睡眠段階の判断処理でよい。
ここで、生体データ処理部１２では、３０秒ごとの複数の特徴量ＳＤ、Range、ＳＵＭ、Square、ＬＣ、ＲＭＳを示すデータに対して、睡眠段階を付与するラベリング処理が行われ、そのラベリング処理が行われた特徴量のデータが、判定部１３および学習部１５に供給される。
ここで付与される睡眠段階は、既に知られた睡眠段階推定手法を適用して、生体振動データから取得したものである。例えば、生体振動データの特定の周波数成分のレベルと閾値とを比較して、判断した睡眠段階を、その期間の特徴量のデータに対して付与する。なお、ここで生体振動データから睡眠段階を判断する処理は、睡眠時無呼吸症候群の患者であることを考慮しない、いわゆる健常者用の睡眠段階の判断処理でよい。
【0026】
そして、学習部１５は、生体データ処理部１２で得た生体振動データの特徴量を使って、学習部１５が、機械学習の１つであるランダムフォレストを用いて、特徴量のデータから覚醒（ＷＡＫＥ）か否かを分類する複数の決定木を生成する。
学習部１５での複数の決定木の生成は、学習フェーズで行われる。学習部１５で生成された複数の決定木は、判定フェーズで利用され、判定部１３が、複数の決定木を使った判定結果から、睡眠段階を判定して、睡眠時無呼吸症候群か否かを判定する。
学習部１５での複数の決定木の生成は、学習フェーズで行われる。学習部１５で生成された複数の決定木は、判定フェーズで利用され、判定部１３が、複数の決定木を使った判定結果から、睡眠段階を判定して、睡眠時無呼吸症候群か否かを判定する。
【0027】
また、上述したラベリング処理で、生体データ処理部１２が３０秒ごとの複数の特徴量に睡眠段階を付与する際には、生体振動データとは全く別のデータの解析で、被測定者の睡眠段階を判断して、対応した期間の特徴量のデータに対して付与するようにしてもよい。例えば、被測定者に取り付けた電極から得たデータや呼吸などを検出するセンサのデータに基づいて、被測定者の睡眠段階を判断してもよい。
【0028】
そして、判定フェーズでは、学習部１５が学習フェーズで得た複数の決定木を使って、判定部１３が、生体データ処理部１２で得た生体振動データから、睡眠時無呼吸症候群か否かを判別する。
このときには、生成された複数の決定木で、覚醒（ＷＡＫＥ）と判定されたときのルール、つまり決定木の分類経路を学習部１５が抽出し、抽出したルールの種類数によって、判定部１３が睡眠時無呼吸症候群か否かを判定する。
このときには、生成された複数の決定木で、覚醒（ＷＡＫＥ）と判定されたときのルール、つまり決定木の分類経路を学習部１５が抽出し、抽出したルールの種類数によって、判定部１３が睡眠時無呼吸症候群か否かを判定する。
【0029】
出力部１４は、判定部１３が判定した睡眠時無呼吸症候群か否かの結果を出力する。出力部１４は、例えば表示装置により構成され、睡眠時無呼吸症候群の判定結果を表示する。あるいは、出力部１４を記録装置として構成して、一晩の睡眠状態などと共に睡眠時無呼吸症候群の判定結果を記録するようにしてもよい。また、出力部１４が表示または記録を行う際には、睡眠時無呼吸症候群の判定結果だけでなく、睡眠段階の判定結果の表示または記録を同時に行うようにしてもよい。
【0030】
［１－２．睡眠時無呼吸症候群判定装置のハードウェア構成例］
図３は、睡眠時無呼吸症候群判定装置１０をコンピュータ装置で構成した場合のハードウェア構成例を示す。
コンピュータ装置Ｃは、バスＣ８に接続されたＣＰＵ（Central Processing Unit：中央処理装置）Ｃ１、ＲＯＭ（Read Only Memory）Ｃ２、およびＲＡＭ（Random Access Memory）Ｃ３を備える。さらに、コンピュータ装置Ｃは、不揮発性ストレージＣ４、ネットワークインターフェイスＣ５、入力装置Ｃ６、および表示装置Ｃ７を備える。
図３は、睡眠時無呼吸症候群判定装置１０をコンピュータ装置で構成した場合のハードウェア構成例を示す。
コンピュータ装置Ｃは、バスＣ８に接続されたＣＰＵ（Central Processing Unit：中央処理装置）Ｃ１、ＲＯＭ（Read Only Memory）Ｃ２、およびＲＡＭ（Random Access Memory）Ｃ３を備える。さらに、コンピュータ装置Ｃは、不揮発性ストレージＣ４、ネットワークインターフェイスＣ５、入力装置Ｃ６、および表示装置Ｃ７を備える。
【0031】
ＣＰＵ Ｃ１は、睡眠時無呼吸症候群判定装置１０の生体データ処理部１２や判定部１３が備える各機能を実現するソフトウェアのプログラムコードをＲＯＭ Ｃ２から読み出して実行する。圧力データを周波数解析するＦＦＴ処理についても、該当する処理を実行するプログラムをＲＯＭ Ｃ２から読み出して、ＣＰＵ Ｃ１が実行する。ＲＡＭ Ｃ３には、演算処理の途中に発生した変数やパラメータ等が一時的に書き込まれる。
【0032】
不揮発性ストレージＣ４としては、例えば、ＨＤＤ（Hard disk drive）、ＳＳＤ（Solid State Drive）、フレキシブルディスク、光ディスク、光磁気ディスク、ＣＤ－ＲＯＭ、ＣＤ－Ｒ、磁気テープ、不揮発性のメモリ等が用いられる。この不揮発性ストレージＣ４には、ＯＳ（Operating System）、各種のパラメータの他に、コンピュータ装置Ｃを睡眠時無呼吸症候群判定装置１０として機能させるためのプログラムが記録されている。また、判定部１３が判定した睡眠段階についてのデータが、不揮発性ストレージＣ４に記録される。
【0033】
ネットワークインターフェイスＣ５には、例えば、ＮＩＣ（Network Interface Card）等が用いられ、端子が接続されたＬＡＮ（Local Area Network）、専用線等を介して各種のデータを送受信することが可能である。例えば、コンピュータ装置Ｃは、マットレスセンサ２が出力する圧力データを、ネットワークインターフェイスＣ５を介して取得する。入力装置Ｃ６は、例えばキーボード等の機器で構成され、この入力装置Ｃ６により、睡眠時無呼吸症候群判定装置１０で睡眠時無呼吸症候群を判定する期間の設定や、判定結果の表示形態の指示等が行われる。表示装置Ｃ７には、睡眠時無呼吸症候群判定装置１０で睡眠時無呼吸症候群の判定結果が表示される。
なお、睡眠時無呼吸症候群判定装置１０を、プログラム（ソフトウェア）の実行で判定装置として機能するコンピュータ装置から構成するのは一例であり、睡眠時無呼吸症候群判定装置１０の一部または全ての処理を実行する専用のハードウェアを用意してもよい。
なお、睡眠時無呼吸症候群判定装置１０を、プログラム（ソフトウェア）の実行で判定装置として機能するコンピュータ装置から構成するのは一例であり、睡眠時無呼吸症候群判定装置１０の一部または全ての処理を実行する専用のハードウェアを用意してもよい。
【0034】
［１－３．睡眠時無呼吸症候群を判定する原理］
次に、睡眠時無呼吸症候群判定装置１０が睡眠時無呼吸症候群を判定する原理について、図４および図５を参照して説明する。
図４は、ある時間の間、健常者が睡眠段階を測定した例（ａ）と、睡眠時無呼吸症候群の患者が睡眠段階を測定した例（ｂ）を比較した図である。なお、以下の説明では、睡眠時無呼吸症候群の患者は、ＳＡＳ患者とも称する。
次に、睡眠時無呼吸症候群判定装置１０が睡眠時無呼吸症候群を判定する原理について、図４および図５を参照して説明する。
図４は、ある時間の間、健常者が睡眠段階を測定した例（ａ）と、睡眠時無呼吸症候群の患者が睡眠段階を測定した例（ｂ）を比較した図である。なお、以下の説明では、睡眠時無呼吸症候群の患者は、ＳＡＳ患者とも称する。
【0035】
図４のそれぞれの例は、横軸が睡眠時間を示し、縦軸が睡眠段階を示す。縦軸の睡眠段階は、最も上側が、最も睡眠段階が浅い覚醒（ＷＡＫＥ）を示し、下側に行くほど、レム睡眠（ＲＥＭ）、ステージ１のノンレム睡眠（ＮＲＥＭ１）、ステージ２のノンレム睡眠（ＮＲＥＭ２）、ステージ３のノンレム睡眠（ＮＲＥＭ３）、ステージ４のノンレム睡眠（ＮＲＥＭ４）と順に変化して、深い睡眠段階になることを示す。
但し、図４の例では、（ａ）に示す健常者と、（ｂ）に示すＳＡＳ患者のいずれの場合も、最も深い睡眠段階は、ステージ３のノンレム睡眠（ＮＲＥＭ３）であり、ステージ４のノンレム睡眠（ＮＲＥＭ４）には到達していない。
図４の健常者（ａ）とＳＡＳ患者（ｂ）とを比較すると判るように、ＳＡＳ患者は無呼吸が原因での覚醒（ＷＡＫＥ）が多発している。
但し、図４の例では、（ａ）に示す健常者と、（ｂ）に示すＳＡＳ患者のいずれの場合も、最も深い睡眠段階は、ステージ３のノンレム睡眠（ＮＲＥＭ３）であり、ステージ４のノンレム睡眠（ＮＲＥＭ４）には到達していない。
図４の健常者（ａ）とＳＡＳ患者（ｂ）とを比較すると判るように、ＳＡＳ患者は無呼吸が原因での覚醒（ＷＡＫＥ）が多発している。
【0036】
図５は、健常者とＳＡＳ患者とで、マットレスセンサ２が出力する生体振動データの相違を示すものである。図５の（ａ）は、睡眠段階が覚醒（ＷＡＫＥ）時に、健常者に多く現れる生体振動データの例である。図５の（ｂ）は、睡眠段階が覚醒（ＷＡＫＥ）時に、ＳＡＳ患者に多く現れる生体振動データの例を示す。図５に示す各特性のグラフの横軸は時間、縦軸は生体振動のレベルを示す。
【0037】
健常者の場合、図５の（ａ）に示すように、睡眠段階が覚醒（ＷＡＫＥ）時には、寝返り等の体動に連動して、一時的に比較的高いレベルの生体振動が検出されることが多い。例えば、図５の（ａ）の左側に示す特性は、大きな寝返りがあった場合であり、図５の（ａ）の右側に示す特性は、比較的小さな体の動きがあった場合である。
一方、ＳＡＳ患者の場合、睡眠段階が覚醒（ＷＡＫＥ）時には、様々な変化状態の生体振動が検出される。すなわち、図５の（ｂ）の左端の特性で示すように、健常者と同様な一時的に高いレベルの生体振動となることもあるが、図５の（ｂ）の右側の３つの特性で示すように、無呼吸による小さな体動が検出される。
ＳＡＳ患者の場合、覚醒（ＷＡＫＥ）時には、舌が気道を塞ぐことにより無呼吸が起こり、その影響が拍動（心拍）の変化に現れる。そして、健常者では検出されない無呼吸による小さな体動が周期的に検出される。但し、図５の（ｂ）の右端に示すように、健常者と同等の動きが覚醒（ＷＡＫＥ）時に検出される場合もある。
一方、ＳＡＳ患者の場合、睡眠段階が覚醒（ＷＡＫＥ）時には、様々な変化状態の生体振動が検出される。すなわち、図５の（ｂ）の左端の特性で示すように、健常者と同様な一時的に高いレベルの生体振動となることもあるが、図５の（ｂ）の右側の３つの特性で示すように、無呼吸による小さな体動が検出される。
ＳＡＳ患者の場合、覚醒（ＷＡＫＥ）時には、舌が気道を塞ぐことにより無呼吸が起こり、その影響が拍動（心拍）の変化に現れる。そして、健常者では検出されない無呼吸による小さな体動が周期的に検出される。但し、図５の（ｂ）の右端に示すように、健常者と同等の動きが覚醒（ＷＡＫＥ）時に検出される場合もある。
【0038】
本実施の形態例では、このような図５の（ａ）に示す健常者の覚醒（ＷＡＫＥ）時における生体振動と、図５の（ｂ）に示すようなＳＡＳ患者の覚醒（ＷＡＫＥ）時における様々な生体振動とを、ランダムフォレストを利用した学習処理で区別することで、睡眠時無呼吸症候群の判定を行うものである。
【0039】
［１－４．睡眠時無呼吸症候群の判定処理の流れ］
図６は、本実施の形態例の睡眠時無呼吸症候群判定装置１０が睡眠時無呼吸症候群を判定する処理の流れを示すフローチャートである。
睡眠時無呼吸症候群を判定する際には、最初に、図６の（ａ）に示す学習フェーズで学習を行って、ランダムフォレストモデルを構築する処理が行われる。その後、図６の（ｂ）に示す判定フェーズで、被測定者の生体振動を取得して、ランダムフォレストモデルを使って、睡眠時無呼吸症候群を判定する処理が行われる。
まず、図６の（ａ）に示す学習フェーズについて説明すると、生体データ取得部１１が被測定者の睡眠中の生体振動データの取得を開始し、生体データ処理部１２は、一定時間ｔａ秒ごとに、マットレスセンサ２が出力する生体振動データ（圧力データ）の平均値を算出し、算出した平均値を記録する（ステップＳ１１）。ここでのｔａ秒は、例えば１秒とする。
図６は、本実施の形態例の睡眠時無呼吸症候群判定装置１０が睡眠時無呼吸症候群を判定する処理の流れを示すフローチャートである。
睡眠時無呼吸症候群を判定する際には、最初に、図６の（ａ）に示す学習フェーズで学習を行って、ランダムフォレストモデルを構築する処理が行われる。その後、図６の（ｂ）に示す判定フェーズで、被測定者の生体振動を取得して、ランダムフォレストモデルを使って、睡眠時無呼吸症候群を判定する処理が行われる。
まず、図６の（ａ）に示す学習フェーズについて説明すると、生体データ取得部１１が被測定者の睡眠中の生体振動データの取得を開始し、生体データ処理部１２は、一定時間ｔａ秒ごとに、マットレスセンサ２が出力する生体振動データ（圧力データ）の平均値を算出し、算出した平均値を記録する（ステップＳ１１）。ここでのｔａ秒は、例えば１秒とする。
【0040】
そして、生体データ処理部１２は、ｔａ秒ごとの平均値から、ｔｂ秒ごとに複数の特徴量を算出する（ステップＳ１２）。ｔｂ秒は、ｔａ秒よりも十分に長い時間であり、例えば３０秒とする。ここでの複数の特徴量の一例としては、先に説明した分散量を示す特徴量ＳＤ、最大値と最小値の差の特徴量Range、積算値の特徴量ＳＵＭ、積算値の二乗の特徴量Square、ｔａ秒ごとの生体振動データの平均値ｘｎとその平均値ｘの１つ前の平均値ｘｎ－１との差分の積算値ＬＣ、ｔａ秒間ごとの平均値の二乗値の積算値をｔｂ秒間で平均してルートを求めた値であるＲＭＳがある。
【0041】
次に、生体データ処理部１２は、ｔｂ秒ごとの特徴量のデータに対して、該当する特徴量の生体振動データを取得した期間での被測定者の睡眠段階のデータを取る。ここで、被測定者の睡眠段階のデータは、生体データ処理部１２による生体振動データの解析から取得される。あるいは、睡眠時無呼吸症候群判定装置１０とは別の測定装置を使って、被測定者の睡眠段階のデータを取得してもよい。
このようにして被測定者の睡眠段階のデータを取得すると、生体データ処理部１２は、ｔｂ秒ごとの特徴量のデータに、睡眠段階を付与するラベリング処理を行う（ステップＳ１３）。ここでのラベリング処理で付与する睡眠段階は、少なくとも睡眠段階が、覚醒（ＷＡＫＥ）か、または覚醒以外（Ｎｏｎ－ＷＡＫＥ）かの２段階である。
生体データ処理部１２で得られた睡眠段階がラベリングされた特徴量のデータは、学習部１５に送られる。
このようにして被測定者の睡眠段階のデータを取得すると、生体データ処理部１２は、ｔｂ秒ごとの特徴量のデータに、睡眠段階を付与するラベリング処理を行う（ステップＳ１３）。ここでのラベリング処理で付与する睡眠段階は、少なくとも睡眠段階が、覚醒（ＷＡＫＥ）か、または覚醒以外（Ｎｏｎ－ＷＡＫＥ）かの２段階である。
生体データ処理部１２で得られた睡眠段階がラベリングされた特徴量のデータは、学習部１５に送られる。
【0042】
そして、学習部１５は、機械学習の１つであるランダムフォレストモデル（ＲＦモデル）を用いた学習で、複数の特徴量から睡眠段階を判定するための、それぞれ異なる分岐条件を持った複数の決定木を構築する（ステップＳ１４）。ここで、複数の決定木を構築する際には、それぞれの決定木の深さ（分岐条件の段階数）と、構築する決定木の数を予め決めておく。例えば、決定木の深さ（分岐数）を１０段階、構築する決定木の数を３００本とする。
また、分岐条件を決める際には、２つに分岐する処理を、例えばジニ係数と称される指標を使って、覚醒（ＷＡＫＥ）と覚醒以外（Ｎｏｎ－ＷＡＫＥ）の２つに大きく分かれるような条件とする。
また、学習部１５が、複数の決定木の組み合わせの学習を行う際には、複数の特徴量の組み合わせの種類数と、特徴量の重要度と、各特徴量の出現回数とを判断して行う。例えば、睡眠時無呼吸症候群を判定する上で、どの特徴量を使うのが重要かを学習する。
ここまでが学習フェーズで行われ、学習部１５は、構築したランダムフォレストモデル（ＲＦモデル）を記憶する。
また、分岐条件を決める際には、２つに分岐する処理を、例えばジニ係数と称される指標を使って、覚醒（ＷＡＫＥ）と覚醒以外（Ｎｏｎ－ＷＡＫＥ）の２つに大きく分かれるような条件とする。
また、学習部１５が、複数の決定木の組み合わせの学習を行う際には、複数の特徴量の組み合わせの種類数と、特徴量の重要度と、各特徴量の出現回数とを判断して行う。例えば、睡眠時無呼吸症候群を判定する上で、どの特徴量を使うのが重要かを学習する。
ここまでが学習フェーズで行われ、学習部１５は、構築したランダムフォレストモデル（ＲＦモデル）を記憶する。
【0043】
次に、図６の（ｂ）に示す判定フェーズでの、被測定者の生体振動から睡眠時無呼吸症候群を判定する処理の流れを説明する。
判定フェーズでは、学習フェーズと同様に、ステップＳ１１でのｔａ秒ごとの被測定者の生体振動の取得処理と、ステップＳ１２でのｔｂ秒間ごとの特徴量の算出処理が行われる。
その後、判定部１３は、学習フェーズでのステップＳ１４で学習部１５が生成したランダムフォレストモデルの複数の決定木に、覚醒（ＷＡＫＥ）とラベリング処理が行われた特徴量のデータを入力する（ステップＳ１５）。
この複数の決定木に入力されたデータに基づいて、判定部１３は、各特徴量のデータの分類処理が行われ、分岐結果が覚醒（ＷＡＫＥ）となる際に使用された経路(根ノードから葉ノードまでの一連の特徴量)を抽出する（ステップＳ１６）。このとき、判定部１３は、決定木の各分岐時においてどの種類の特徴量のデータを用いていずれに分岐するかを判断したのかに基づいて、特徴量の種類ごとに判断回数をカウントする。
判定フェーズでは、学習フェーズと同様に、ステップＳ１１でのｔａ秒ごとの被測定者の生体振動の取得処理と、ステップＳ１２でのｔｂ秒間ごとの特徴量の算出処理が行われる。
その後、判定部１３は、学習フェーズでのステップＳ１４で学習部１５が生成したランダムフォレストモデルの複数の決定木に、覚醒（ＷＡＫＥ）とラベリング処理が行われた特徴量のデータを入力する（ステップＳ１５）。
この複数の決定木に入力されたデータに基づいて、判定部１３は、各特徴量のデータの分類処理が行われ、分岐結果が覚醒（ＷＡＫＥ）となる際に使用された経路(根ノードから葉ノードまでの一連の特徴量)を抽出する（ステップＳ１６）。このとき、判定部１３は、決定木の各分岐時においてどの種類の特徴量のデータを用いていずれに分岐するかを判断したのかに基づいて、特徴量の種類ごとに判断回数をカウントする。
【0044】
さらに、判定部１３は、ステップＳ１６で抽出された経路の各分岐で使用した特徴量を比較して、特徴量の種類と数で、被測定者がＳＡＳ患者か否かを判定する（ステップＳ１７）。ここでの特徴の比較とは、例えば上述した複数の特徴量ＳＤ、Range、ＳＵＭ、Square、ＬＣ、ＲＭＳの各々を分岐時に判断した回数（出現回数）、分岐時に判断した特徴量の種類の合計の数、などがある。出現回数については、全ての特徴量の出現回数を判断してもよいが、複数の種類の特徴量の内で重要と思われる特徴量を判断した回数を、ＳＡＳ患者か否かを判定する上で優先して扱う、などの様々な条件を付けて、判別に使用された決定木を使って分類した際の分類基準の組み合わせを設定することが考えられる。
特徴の数の多数決で判定する具体例については後述する。そして、出力部１４が、ステップＳ１７での判定結果を出力する。
特徴の数の多数決で判定する具体例については後述する。そして、出力部１４が、ステップＳ１７での判定結果を出力する。
【0045】
［１－５．判別処理の詳細］
次に、図６のフローチャートに示す流れの各処理の詳細について、図７～図１２を参照して説明する。
図７は、睡眠時無呼吸症候群の患者（ＳＡＳ患者）の生体振動データの例を示す。図７の横軸は時間であり、図７の（ａ）は長時間の測定データを示し、図７の（ｂ）は、その長時間の測定データの一部の区間を拡大して示す。
次に、図６のフローチャートに示す流れの各処理の詳細について、図７～図１２を参照して説明する。
図７は、睡眠時無呼吸症候群の患者（ＳＡＳ患者）の生体振動データの例を示す。図７の横軸は時間であり、図７の（ａ）は長時間の測定データを示し、図７の（ｂ）は、その長時間の測定データの一部の区間を拡大して示す。
【0046】
図７の（ａ），（ｂ）に示すグラフの横軸は時間であり、縦軸は、６段階の睡眠段階の変化ＰＳＧと、生体振動データのレベルＭを示す。６段階の睡眠段階の変化ＰＳＧは、本実施の形態例の睡眠時無呼吸症候群判定装置１０とは別の測定装置で、睡眠ポリグラフ検査法（polysomnography）によって、被測定者の睡眠段階を計測したものである。ここでは、この６段階の睡眠段階の変化ＰＳＧを、正解の睡眠段階とも称する。
図７に示す６段階の睡眠段階は、上側から順に、覚醒（Ｗ）、レム睡眠（Ｒ）、ノンレム睡眠のステージ１（ＮＲ１）、ノンレム睡眠のステージ２（ＮＲ２）、ノンレム睡眠のステージ３（ＮＲ３）、ノンレム睡眠のステージ４（ＮＲ４）である。
図７に示す６段階の睡眠段階は、上側から順に、覚醒（Ｗ）、レム睡眠（Ｒ）、ノンレム睡眠のステージ１（ＮＲ１）、ノンレム睡眠のステージ２（ＮＲ２）、ノンレム睡眠のステージ３（ＮＲ３）、ノンレム睡眠のステージ４（ＮＲ４）である。
【0047】
図７の（ａ），（ｂ）に示す生体振動データのレベルＭは、体動などで高くなるが、例えば健常者の場合には、閾値ＴＨ１をレベルＭが超えたとき、睡眠段階が覚醒であると判別することができる。
ところが、図７の例は、ＳＡＳ患者の生体振動データであり、図５で説明したように、生体振動データが低い状態であっても、睡眠段階が覚醒となることがある。
ところが、図７の例は、ＳＡＳ患者の生体振動データであり、図５で説明したように、生体振動データが低い状態であっても、睡眠段階が覚醒となることがある。
【0048】
図７の（ａ）の覚醒（Ｗ）の箇所に、丸の印（○）を付けた箇所は、正解の睡眠段階ＰＳＧで覚醒（Ｗ）が検出された箇所が、生体振動データのレベルＭと閾値ＴＨ１との比較でも、覚醒（Ｗ）と判別された箇所である。
一方、図７の（ａ）の覚醒（Ｗ）の箇所に、バツの印（×）を付けた箇所は、正解の睡眠段階ＰＳＧで覚醒（Ｗ）が検出された箇所が、生体振動データのレベルＭと閾値ＴＨ１との比較で、覚醒（Ｗ）と判別されなかった箇所である。
この図７の（ａ）に示すように、ＳＡＳ患者の場合、睡眠段階が覚醒（Ｗ）であっても、多様な生体振動状態が存在する。
本実施の形態では、ランダムフォレストモデルでの学習で、この多様な生体振動状態を捉えるようにして、ＳＡＳ患者か否かを判別するようにしたものである。
一方、図７の（ａ）の覚醒（Ｗ）の箇所に、バツの印（×）を付けた箇所は、正解の睡眠段階ＰＳＧで覚醒（Ｗ）が検出された箇所が、生体振動データのレベルＭと閾値ＴＨ１との比較で、覚醒（Ｗ）と判別されなかった箇所である。
この図７の（ａ）に示すように、ＳＡＳ患者の場合、睡眠段階が覚醒（Ｗ）であっても、多様な生体振動状態が存在する。
本実施の形態では、ランダムフォレストモデルでの学習で、この多様な生体振動状態を捉えるようにして、ＳＡＳ患者か否かを判別するようにしたものである。
【0049】
図８は、ランダムフォレストモデルの決定木の概要を示す。
ランダムフォレストモデルの決定木は、学習データ集合Ｓから、部分集合Ｓ１，Ｓ２，・・・，ＳＮtreeをランダムに抽出することで構築される。ＳＮtreeは、構築される決定木の数である。決定木の分岐箇所に示すｙｉは、ｉ番目の特徴量を示す。すなわち、分岐箇所では、複数個の特徴量の内のどの特徴量を判断するのかを示す。
複数の決定木は、データを部分的に学習することで、決定木ごとに異なる学習を行うデータとなっており、複数のルールを生成したことになる。つまり、ＳＡＳ患者の場合における多様な覚醒時の生体振動データの波形を捉えることになる。
ランダムフォレストモデルの決定木は、学習データ集合Ｓから、部分集合Ｓ１，Ｓ２，・・・，ＳＮtreeをランダムに抽出することで構築される。ＳＮtreeは、構築される決定木の数である。決定木の分岐箇所に示すｙｉは、ｉ番目の特徴量を示す。すなわち、分岐箇所では、複数個の特徴量の内のどの特徴量を判断するのかを示す。
複数の決定木は、データを部分的に学習することで、決定木ごとに異なる学習を行うデータとなっており、複数のルールを生成したことになる。つまり、ＳＡＳ患者の場合における多様な覚醒時の生体振動データの波形を捉えることになる。
【0050】
図８の例について説明すると、例えば図８の左端に示す部分集合Ｓ１の決定木としては、最初の分岐で、特徴量ｙ８が１未満か否かが判断され、特徴量ｙ８が１未満である（True）とき、次の分岐で、特徴量ｙ７が６未満か否かが判断される。そして、この判断で、特徴量ｙ７が６未満のものが、睡眠段階が覚醒（Ｗ）のもの（○印）と、睡眠段階が覚醒（Ｗ）以外のもの（×印）に分離される。
また、特徴量ｙ８が１未満でない場合（False）には、次の分岐で、特徴量ｙ２が８未満か否かが判断される。そして、この判断で、特徴量ｙ2が８未満のものが、睡眠段階が覚醒（Ｗ）のもの（○印）と、睡眠段階が覚醒（Ｗ）以外のもの（×印）に分離される。
このようにして、それぞれの決定木の最終的な分岐で、睡眠段階が覚醒（Ｗ）のもの（○印）と、睡眠段階が覚醒（Ｗ）以外のもの（×印）に分離されて行く。なお、図８の例では、説明を簡単にするために分岐箇所が２つのものを示すが、これは一例であり、１つの決定木がより多くの分岐箇所を持ってもよい。
また、特徴量ｙ８が１未満でない場合（False）には、次の分岐で、特徴量ｙ２が８未満か否かが判断される。そして、この判断で、特徴量ｙ2が８未満のものが、睡眠段階が覚醒（Ｗ）のもの（○印）と、睡眠段階が覚醒（Ｗ）以外のもの（×印）に分離される。
このようにして、それぞれの決定木の最終的な分岐で、睡眠段階が覚醒（Ｗ）のもの（○印）と、睡眠段階が覚醒（Ｗ）以外のもの（×印）に分離されて行く。なお、図８の例では、説明を簡単にするために分岐箇所が２つのものを示すが、これは一例であり、１つの決定木がより多くの分岐箇所を持ってもよい。
【0051】
ここで、このように複数の決定木が構築されるとき、決定木で分岐する際の条件に使用した特徴量の寄与度を算出する。
特徴量の寄与度Imp(y)は、例えば次の式で算出される。
特徴量の寄与度Imp(y)は、例えば次の式で算出される。
【0052】
【数1】
【0053】
［数１］式において、ｙｉはｉ番目の特徴量、Ｎtreeは決定木の数、Δ（Ｔｊ（ｙｉ））は、決定木Ｔｊにおける特徴量ｙｉを用いたノードとそのノードの分岐後のジニ係数の差を示す。
ここでのジニ係数は、図８に示すように、該当する分岐点で分岐したことで、睡眠段階が覚醒のもの（○印）と、睡眠段階が覚醒（Ｗ）以外のもの（×印）に明確に分離されるとき高くし、そうでない場合に低くする係数値である。ここでの睡眠段階は、図６のステップＳ１３でラベリングした睡眠段階を利用する。
ここでのジニ係数は、図８に示すように、該当する分岐点で分岐したことで、睡眠段階が覚醒のもの（○印）と、睡眠段階が覚醒（Ｗ）以外のもの（×印）に明確に分離されるとき高くし、そうでない場合に低くする係数値である。ここでの睡眠段階は、図６のステップＳ１３でラベリングした睡眠段階を利用する。
【0054】
ジニ係数の差が大きいとき、図８の左下に示すように、該当する分岐点での分岐で、睡眠段階が覚醒のもの（○印）と、睡眠段階が覚醒（Ｗ）以外のもの（×印）に明確に分離される状態を示す。
一方、ジニ係数の差が小さいとき、図８の左下に示すように、該当する分岐点での分岐で、睡眠段階が覚醒のもの（○印）と、睡眠段階が覚醒（Ｗ）以外のもの（×印）に明確に分離される状態を示す。
一方、ジニ係数の差が小さいとき、図８の左下に示すように、該当する分岐点での分岐で、睡眠段階が覚醒のもの（○印）と、睡眠段階が覚醒（Ｗ）以外のもの（×印）に明確に分離される状態を示す。
【0055】
このようにして、ジニ係数の差に基づいて、睡眠段階が覚醒のものと覚醒以外のものに判別するために必要な特徴量の寄与度を、それぞれの特徴量ごとに算出することで、どの特徴量を使った分岐を行う決定木を利用することが、睡眠段階が覚醒か否かを判別できるのか判るようになる。
【0056】
次に、生体振動データから特徴量の寄与度の判別を使った判定の流れを、図９～図１１を参照して説明する。
図９は、１人の被測定者の睡眠開始（入眠）から起床までの数時間の生体振動データの波形の例を示す。
本実施の形態例の場合、図９に示すような数時間の生体振動データから、３０秒間ごとに特徴量を算出する処理が行われる。例えば、特徴量ＳＤ、Range、ＳＵＭ、Square、ＬＣ、ＲＭＳを、３０秒ごとに算出する。
そして、３０秒ごとの各特徴量を示すデータに、その時点での睡眠段階（ＷＡＫＥまたはＮｏｎ－ＷＡＬＫＥ）をラベリングする処理が行われる。
図９は、１人の被測定者の睡眠開始（入眠）から起床までの数時間の生体振動データの波形の例を示す。
本実施の形態例の場合、図９に示すような数時間の生体振動データから、３０秒間ごとに特徴量を算出する処理が行われる。例えば、特徴量ＳＤ、Range、ＳＵＭ、Square、ＬＣ、ＲＭＳを、３０秒ごとに算出する。
そして、３０秒ごとの各特徴量を示すデータに、その時点での睡眠段階（ＷＡＫＥまたはＮｏｎ－ＷＡＬＫＥ）をラベリングする処理が行われる。
【0057】
次に、図１０に示すように、睡眠段階がラベリングされた３０秒ごとの各特徴量を示すデータから、ランダムフォレストモデル（ＲＦモデル）の決定木を構築する処理が行われる。
このとき、各決定木の深さ（分岐数）と、構築する決定木の合計数を予め決めておく。例えば、決定木の深さを１０段階、決定木の合計数を３００とする。
このとき、各決定木の深さ（分岐数）と、構築する決定木の合計数を予め決めておく。例えば、決定木の深さを１０段階、決定木の合計数を３００とする。
【0058】
そして、構築したランダムフォレストモデル（ＲＦモデル）の複数の決定木に、睡眠段階が覚醒（ＷＡＫＥ）にラベリングされた特徴量のデータを入力する。
【0059】
次に、図１１に示すように、各決定木で、覚醒（ＷＡＫＥ）と分類する際に通過した決定木の経路を抽出し、１つの決定木の経路の途中で分岐を判断した特徴量の組み合わせを、種類ごとにカウントする。
図１１の下側に示す特徴量ごとに示す値は、そのカウント値の一例である。
例えば、図１１の例では、特徴量ＳＤ，ＳＵＭの組み合わせが８０回、特徴量ＬＣ，ＲＭＳの組み合わせが１００回、特徴量ＳＤ，Rangeの組み合わせが８回、特徴量Square，ＳＵＭの組み合わせが７９回、・・・とカウントする。なお、ここで、１０回を超えた組み合わせを１個とカウントする。そして、例えば、経路の特徴量の組み合わせの種類のカウント数が閾値を超えるとき、ＳＡＳ患者であると判別する。
図１１の下側に示す特徴量ごとに示す値は、そのカウント値の一例である。
例えば、図１１の例では、特徴量ＳＤ，ＳＵＭの組み合わせが８０回、特徴量ＬＣ，ＲＭＳの組み合わせが１００回、特徴量ＳＤ，Rangeの組み合わせが８回、特徴量Square，ＳＵＭの組み合わせが７９回、・・・とカウントする。なお、ここで、１０回を超えた組み合わせを１個とカウントする。そして、例えば、経路の特徴量の組み合わせの種類のカウント数が閾値を超えるとき、ＳＡＳ患者であると判別する。
【0060】
図１２は、本実施の形態例によるＳＡＳ患者か否かの判別結果の一例を示す。図１２の縦軸は、経路の特徴量の組み合わせの種類のカウント数を示している。
図１２の例は、９名のＳＡＳ患者Ａ，Ｂ，Ｃ，・・・，Ｉと、９名の健常者（ＳＡＳ患者でない者）ａ，ｂ，ｃ，・・・，ｉとを判定した結果である。図１２に示す１８名の内で、×印を付与した者は、正しく判定できなかった場合を示している。
図１２では、経路の特徴量の組み合わせの種類のカウント数が１００を閾値とし、カウント数が１００以上であるとき、ＳＡＳ患者であると判別し、カウント数が１００未満のとき、ＳＡＳ患者でないと判別している。
図１２の例は、９名のＳＡＳ患者Ａ，Ｂ，Ｃ，・・・，Ｉと、９名の健常者（ＳＡＳ患者でない者）ａ，ｂ，ｃ，・・・，ｉとを判定した結果である。図１２に示す１８名の内で、×印を付与した者は、正しく判定できなかった場合を示している。
図１２では、経路の特徴量の組み合わせの種類のカウント数が１００を閾値とし、カウント数が１００以上であるとき、ＳＡＳ患者であると判別し、カウント数が１００未満のとき、ＳＡＳ患者でないと判別している。
【0061】
本実施の形態例の場合、図１２から判るように、ＳＡＳ患者９名の内で、６名をＳＡＳ患者であると正しく判別することができ、３名をＳＡＳ患者でないと誤判定した。一方、健常者については、７名を健常者と正しく判定し、２名をＳＡＳ患者と誤判定した。
この図１２から判るように、高い精度でＳＡＳ患者か否かを判別することができる。なお、図１２に示す結果は判定についての１つのシミュレーション結果であり、決定木による学習を繰り返すことで、判定の精度をより高くすることができるようになる。
この図１２から判るように、高い精度でＳＡＳ患者か否かを判別することができる。なお、図１２に示す結果は判定についての１つのシミュレーション結果であり、決定木による学習を繰り返すことで、判定の精度をより高くすることができるようになる。
【0062】
＜２．第２の実施の形態例＞
以下、本発明の第２の実施の形態例について、図１３～図１９を参照して説明する。
第２の実施の形態例においても、睡眠時無呼吸症候群判定装置１０の構成については、第１の実施の形態例で図１～図３で説明した構成をそのまま適用する。
本実施の形態例では、睡眠時無呼吸症候群判定装置１０の生体データ処理部１２での生体データ処理、判定部１３での判定処理、及び学習部１５での学習処理が、第１の実施の形態例と異なる。
以下、本発明の第２の実施の形態例について、図１３～図１９を参照して説明する。
第２の実施の形態例においても、睡眠時無呼吸症候群判定装置１０の構成については、第１の実施の形態例で図１～図３で説明した構成をそのまま適用する。
本実施の形態例では、睡眠時無呼吸症候群判定装置１０の生体データ処理部１２での生体データ処理、判定部１３での判定処理、及び学習部１５での学習処理が、第１の実施の形態例と異なる。
【0063】
［２－１．睡眠段階判定処理の流れ］
図１３は、本実施の形態例の睡眠時無呼吸症候群判定装置１０が睡眠時無呼吸症候群を判定する処理の流れを示すフローチャートである。
睡眠時無呼吸症候群を判定する際には、最初に、図１３の（ａ）に示す学習フェーズで学習を行って、ランダムフォレストモデルを構築する処理が行われる。その後、図１３の（ｂ）に示す判定フェーズで、被測定者の生体振動を取得して、ランダムフォレストモデルを使って、睡眠時無呼吸症候群を判定する処理が行われる。
まず、図１３の（ａ）に示す学習フェーズについて説明すると、生体データ取得部１１が被測定者の睡眠中の生体振動データの取得を開始し、生体データ処理部１２は、一定時間ｔｃ秒ごとに、マットレスセンサ２が出力する生体振動データ（圧力データ）を取得し、取得した生体振動データを記録する。ここでのｔｃ秒は、例えば１秒とする。
ここで、生体振動データをｔｃ秒ごと（１秒ごと）に算出する際には、過去の所定時間（例えば６４秒間）の生体振動データをフーリエ変換して、所定時間での周波数ごとのパワースペクトルを算出する（ステップＳ２１）。
図１３は、本実施の形態例の睡眠時無呼吸症候群判定装置１０が睡眠時無呼吸症候群を判定する処理の流れを示すフローチャートである。
睡眠時無呼吸症候群を判定する際には、最初に、図１３の（ａ）に示す学習フェーズで学習を行って、ランダムフォレストモデルを構築する処理が行われる。その後、図１３の（ｂ）に示す判定フェーズで、被測定者の生体振動を取得して、ランダムフォレストモデルを使って、睡眠時無呼吸症候群を判定する処理が行われる。
まず、図１３の（ａ）に示す学習フェーズについて説明すると、生体データ取得部１１が被測定者の睡眠中の生体振動データの取得を開始し、生体データ処理部１２は、一定時間ｔｃ秒ごとに、マットレスセンサ２が出力する生体振動データ（圧力データ）を取得し、取得した生体振動データを記録する。ここでのｔｃ秒は、例えば１秒とする。
ここで、生体振動データをｔｃ秒ごと（１秒ごと）に算出する際には、過去の所定時間（例えば６４秒間）の生体振動データをフーリエ変換して、所定時間での周波数ごとのパワースペクトルを算出する（ステップＳ２１）。
【0064】
そして、生体データ処理部１２は、フーリエ変換した周波数ごとのパワースペクトルのデータに、睡眠段階をラベリングする（ステップＳ２２）。ここでのラベリング処理で付与する睡眠段階は、少なくとも睡眠段階が、覚醒（ＷＡＫＥ）か、または覚醒以外（Ｎｏｎ－ＷＡＫＥ）かの２段階である。このラベリングする睡眠段階は、第１の実施の形態例でラベリングした睡眠段階と同様に、生体データ処理部１２による生体振動データの解析から取得される。あるいは、睡眠時無呼吸症候群判定装置１０とは別の測定装置を使って、被測定者の睡眠段階のデータを取得してもよい。
生体データ処理部１２で得られた睡眠段階がラベリングされたパワースペクトルのデータは、学習部１５に送られる。
生体データ処理部１２で得られた睡眠段階がラベリングされたパワースペクトルのデータは、学習部１５に送られる。
【0065】
そして、学習部１５は、機械学習の１つであるランダムフォレストモデル（ＲＦモデル）を用いた学習で、複数のパワースペクトルから睡眠段階を判定するための、それぞれ異なる分岐条件を持った複数の決定木を構築する（ステップＳ２３）。ここで、複数の決定木を構築する際には、それぞれの決定木の深さ（分岐条件の段階数）と、構築する決定木の数を予め決めておく。
ここまでが学習フェーズで行われ、学習部１５は、構築したランダムフォレストモデル（ＲＦモデル）を記憶する。
ここまでが学習フェーズで行われ、学習部１５は、構築したランダムフォレストモデル（ＲＦモデル）を記憶する。
【0066】
次に、図１３の（ｂ）に示す判定フェーズでの、被測定者の生体振動から睡眠時無呼吸症候群を判定する処理の流れを説明する。
判定フェーズでは、学習フェーズと同様に、ステップＳ２１での所定時間での周波数ごとのパワースペクトルの算出処理が行われる。
その後、判定部１３は、睡眠段階が覚醒（ＷＡＫＥ）とラベリングされたデータについて、学習フェーズでのステップＳ２３で学習部１５が生成したランダムフォレストモデル（ＲＦモデル）から、各周波数成分の寄与度を算出する（ステップＳ２４）。
その後、判定部１３は、ジニ係数を利用して、各周波数成分の寄与度の分布を大小に分割する（ステップＳ２５）。
判定フェーズでは、学習フェーズと同様に、ステップＳ２１での所定時間での周波数ごとのパワースペクトルの算出処理が行われる。
その後、判定部１３は、睡眠段階が覚醒（ＷＡＫＥ）とラベリングされたデータについて、学習フェーズでのステップＳ２３で学習部１５が生成したランダムフォレストモデル（ＲＦモデル）から、各周波数成分の寄与度を算出する（ステップＳ２４）。
その後、判定部１３は、ジニ係数を利用して、各周波数成分の寄与度の分布を大小に分割する（ステップＳ２５）。
【0067】
さらに、判定部１３は、各周波数成分の寄与度の分布を大小に分割した関係(厳密にはｘ軸で大小，ｙ軸で上下の２つの関係)から睡眠時無呼吸症候群を判定して、睡眠時無呼吸症候群か否かを判別する（ステップＳ２６）。
判定部１３が判別した結果は、出力部１４が出力する。
判定部１３が判別した結果は、出力部１４が出力する。
【0068】
［２－２．周波数ごとのパワースペクトルから判別する原理］
次に、図１４を参照して、睡眠時無呼吸症候群の患者（ＳＡＳ患者）と、健常者との周波数ごとのパワースペクトルの寄与度の相違を説明する。
図１４の（ａ）は、睡眠段階が覚醒（ＷＡＫＥ）とラベリングされたデータで、ＳＡＳ患者の場合のパワースペクトルの寄与度を示し、図１４の（ｂ）は、睡眠段階が覚醒（ＷＡＫＥ）とラベリングされたデータで、ＳＡＳ患者でない健常者の場合のパワースペクトルの寄与度を示す。
次に、図１４を参照して、睡眠時無呼吸症候群の患者（ＳＡＳ患者）と、健常者との周波数ごとのパワースペクトルの寄与度の相違を説明する。
図１４の（ａ）は、睡眠段階が覚醒（ＷＡＫＥ）とラベリングされたデータで、ＳＡＳ患者の場合のパワースペクトルの寄与度を示し、図１４の（ｂ）は、睡眠段階が覚醒（ＷＡＫＥ）とラベリングされたデータで、ＳＡＳ患者でない健常者の場合のパワースペクトルの寄与度を示す。
【0069】
ここでは、生体データ処理部１２は、１６Ｈｚで検出される生体振動データの６４秒間のセンサ値をフーリエ変換して、［１６×６４ｓ＝１０２４］のパワースペクトルを算出し、その１０２４のパワースペクトルの内の０～８Ｈｚの成分を、５１２次元のベクトル成分に変換する処理を行う。ここで、０番目から１１番目の次元のベクトル成分は０とする。
【0070】
図１４の縦軸は、この５１２次元のベクトル成分を示し、０～８Ｈｚの周波数位置を示す。例えば、図１４の縦軸の１００の位置は、［１００÷６４］＝約１．５６Ｈｚの周波数位置を示す。横軸は、それぞれの周波数成分のパワースペクトルの寄与度を示す。
なお、図１４の例の場合、生体振動の内で、例えば呼吸に伴った生体振動は、主として１０～４０の位置に現れ、心拍に伴った生体振動は、主として４０～１２０の位置に現れる。
なお、図１４の例の場合、生体振動の内で、例えば呼吸に伴った生体振動は、主として１０～４０の位置に現れ、心拍に伴った生体振動は、主として４０～１２０の位置に現れる。
【0071】
ここで、図１４の（ａ）に示すＳＡＳ患者の場合には、図中の周波数位置１００（周波数約１．５６Ｈｚ）よりも低い周波数成分の寄与度が大きく、それよりも高い周波数成分の寄与度が小さいことが判る。
一方、図１４の（ｂ）の健常者の場合には、図中の周波数位置１００（周波数約１．５６Ｈｚ）よりも高い周波数成分の寄与度が大きく、それよりも低い周波数成分の寄与度が小さいことが判る。
本実施の形態例の場合には、この図１４に示すＳＡＳ患者と健常者との周波数成分のパワースペクトルの寄与度の相違に着目して、ランダムフォレストモデルを使った学習で、ＳＡＳ患者と健常者を判別するようにしたものである。
一方、図１４の（ｂ）の健常者の場合には、図中の周波数位置１００（周波数約１．５６Ｈｚ）よりも高い周波数成分の寄与度が大きく、それよりも低い周波数成分の寄与度が小さいことが判る。
本実施の形態例の場合には、この図１４に示すＳＡＳ患者と健常者との周波数成分のパワースペクトルの寄与度の相違に着目して、ランダムフォレストモデルを使った学習で、ＳＡＳ患者と健常者を判別するようにしたものである。
【0072】
［２－３．判別処理の詳細］
次に、図１５～図１８を参照して、ＳＡＳ患者を判別する処理の詳細について説明する。
図１５は、１人の被測定者の睡眠開始（入眠）から起床までの数時間の生体振動データの波形の例を示す。
本実施の形態例の場合、図１５に示すような数時間の生体振動データから、６４秒ごとに生体振動データ（圧力データ）の周波数を算出し、パワースペクトルを得る。但し、６４秒単位でパワースペクトルを得る際に、そのパワースペクトルを算出する幅は、１秒ごとに動かすようにしている。周波数成分のパワースペクトルを得る処理の詳細は、図１４の説明で述べた通りである。
次に、図１５～図１８を参照して、ＳＡＳ患者を判別する処理の詳細について説明する。
図１５は、１人の被測定者の睡眠開始（入眠）から起床までの数時間の生体振動データの波形の例を示す。
本実施の形態例の場合、図１５に示すような数時間の生体振動データから、６４秒ごとに生体振動データ（圧力データ）の周波数を算出し、パワースペクトルを得る。但し、６４秒単位でパワースペクトルを得る際に、そのパワースペクトルを算出する幅は、１秒ごとに動かすようにしている。周波数成分のパワースペクトルを得る処理の詳細は、図１４の説明で述べた通りである。
【0073】
生体データ処理部１２は、周波数成分のパワースペクトルのそれぞれについて、覚醒（ＷＡＫＥ）か覚醒以外（Ｎｏｎ－ＷＡＫＥ）かにラベリングする。
そして、ラベリングされた周波数成分のパワースペクトルから、ランダムフォレストモデル（ＲＦモデル）を構築する。
そして、ラベリングされた周波数成分のパワースペクトルから、ランダムフォレストモデル（ＲＦモデル）を構築する。
【0074】
図１６に示すように、生体データ処理部１２は、ランダムフォレストモデルから、各周波数の寄与度を算出する。ランダムフォレストモデルは、第１の実施の形態例で説明したように、複数の決定木を持ち、それぞれの決定木が１つ又は複数の分岐を持つ。但し、第２の実施の形態例の場合には、分岐条件が周波数になる。
ランダムフォレストモデルから、各周波数の寄与度を算出する際には、例えば次の［数２］式から、寄与度Ｉｍｐ（ｙｉ）を求める。
ランダムフォレストモデルから、各周波数の寄与度を算出する際には、例えば次の［数２］式から、寄与度Ｉｍｐ（ｙｉ）を求める。
【0075】
【数2】
【0076】
［数２］式において、ｙｉはｉ番目の周波数（図１４の縦軸の数値）、Ｎｔｒｅｅは決定木の数、Δ（Ｔｊ（ｙｉ））は、決定木Ｔｊにおける周波数ｙｉを用いたノードとそのノードの分岐後のジニ係数の差である。
【0077】
例えば、図１６にジニ係数の差が大の場合を示すように、ｙ１０の周波数位置では、決定木での分岐で、ｙ１０＜６の判別を行うことで、○印で示す覚醒（ＷＡＫＥ）か、×印データ示す覚醒以外（Ｎｏｎ－ＷＡＫＥ）かを明確に区別できている。
【0078】
一方、図１６にジニ係数の差が小の場合を示すように、ｙ５１２の周波数位置では、決定木での分岐で、ｙ５１２＜３の判別を行うことで、○印で示す覚醒（ＷＡＫＥ）と、×印データ示す覚醒以外（Ｎｏｎ－ＷＡＫＥ）とが、決定木では区別できない状態である。
【0079】
なお、それぞれのノードｒでのジニ係数ＩＧ（ｒ）は、例えば次の［数３］式から求める。
【0080】
【数3】
【0081】
［数３］式において、ｃは、クラスの数（分類数）であり、ｐ（ｉ|ｒ）は、ノードｒ内のクラスｉに属する要素の割合を示す。
また、ジニ係数の差は、例えば次の［数４］式から求める。
また、ジニ係数の差は、例えば次の［数４］式から求める。
【0082】
【数4】
【0083】
［数４］式において、ｒＬはノードｒを分割したときの左側であり、ｒＲはノードｒを分割したときの右側であり、ｐ（ｒｘ）は、分割後のノードｒｘに属する要素数の割合である。
【0084】
図１７は、周波数分布を、ジニ係数を利用して大小（上下）に分割する例を示す。
図１７（ａ）は、ＳＡＳ患者の場合であり、図１７（ｂ）は健常者の場合である。
図１７の例の場合には、横軸の寄与度ｘが、特定位置ａよりも少ないか多いかで（つまりｘ＝ａの寄与度を境界として）、大小の要素で分割する。さらに、縦軸の周波数ｙが、特定位置ｂよりも大きい要素か小さい要素かで、上下（高い周波数か低い周波数か）に分割する。
図１７（ａ）は、ＳＡＳ患者の場合であり、図１７（ｂ）は健常者の場合である。
図１７の例の場合には、横軸の寄与度ｘが、特定位置ａよりも少ないか多いかで（つまりｘ＝ａの寄与度を境界として）、大小の要素で分割する。さらに、縦軸の周波数ｙが、特定位置ｂよりも大きい要素か小さい要素かで、上下（高い周波数か低い周波数か）に分割する。
【0085】
ここで、本実施の形態例の場合には、ランダムフォレストモデルの複数の決定木による学習で、寄与度ｘと周波数ｙが、最もきれいに分割される特定位置ａ，ｂの値を探索する処理を行う。
このようにして特定位置ａ，ｂの値が定まることで、図１４で説明したＳＡＳ患者の周波数特性と、健常者の周波数特性とを区別できるようになり、被測定者の生体振動データがいずれに該当するか判別できるようになる。
このようにして特定位置ａ，ｂの値が定まることで、図１４で説明したＳＡＳ患者の周波数特性と、健常者の周波数特性とを区別できるようになり、被測定者の生体振動データがいずれに該当するか判別できるようになる。
【0086】
図１８は、ＳＡＳ患者か否かを判定する指標ＳＩＦの例を示す。
ＳＡＳ患者判定指標ＳＩＦ（Spectrum Importance Feature）は、次の［数５］式で定義される。
ＳＡＳ患者判定指標ＳＩＦ（Spectrum Importance Feature）は、次の［数５］式で定義される。
【0087】
【数5】
【0088】
［数５］式において、Overaveは上側の平均であり、Overwidhtは上側の幅であり、Undermaxは、下側の最大値である。
【0089】
図１９は、本実施の形態例によるＳＡＳ患者か否かの判別結果の一例を示す。図１９の縦軸は指標ＳＩＦの値を示す。
ここでは、指標ＳＩＦの閾値を６０００とし、閾値６０００以上であるとき、ＳＡＳ患者と判別し、閾値６０００未満としたとき、ＳＡＳ患者でないと判別したものである。
図１９の例では、９名のＳＡＳ患者Ａ，Ｂ，Ｃ，・・・，Ｉと、９名の健常者（ＳＡＳ患者でない者）ａ，ｂ，ｃ，・・・，ｉとを判定した結果を示す。この図１９の場合には、ＳＡＳ患者Ａ～Ｉと健常者ａ～ｉを、指標ＳＩＦの数値順に示している。
ここでは、指標ＳＩＦの閾値を６０００とし、閾値６０００以上であるとき、ＳＡＳ患者と判別し、閾値６０００未満としたとき、ＳＡＳ患者でないと判別したものである。
図１９の例では、９名のＳＡＳ患者Ａ，Ｂ，Ｃ，・・・，Ｉと、９名の健常者（ＳＡＳ患者でない者）ａ，ｂ，ｃ，・・・，ｉとを判定した結果を示す。この図１９の場合には、ＳＡＳ患者Ａ～Ｉと健常者ａ～ｉを、指標ＳＩＦの数値順に示している。
【0090】
図１９に示す例では、９名のＳＡＳ患者Ａ～Ｉを全てＳＡＳ患者と判別し、９名の健常者ａ～ｉを非ＳＡＳ患者と判別している。
この図１９から判るように、本実施の形態例によると、非常に精度の高い睡眠時無呼吸症候群の患者の判別でできるようになる。
この図１９から判るように、本実施の形態例によると、非常に精度の高い睡眠時無呼吸症候群の患者の判別でできるようになる。
【0091】
＜３．変形例＞
なお、上述した各実施の形態例で説明した特徴量や周波数の数値などは、好適な一例を示したものであり、これらの特徴量や周波数の数値に限定されるものではない。
なお、上述した各実施の形態例で説明した特徴量や周波数の数値などは、好適な一例を示したものであり、これらの特徴量や周波数の数値に限定されるものではない。
【符号の説明】
【0092】
１…ベッド、２…マットレスセンサ、１０…睡眠時無呼吸症候群判定装置、１１…生体データ取得部、１２…生体データ処理部、１３…判定部、１４…出力部、１５…学習部、Ａ…被測定者、Ｃ…コンピュータ装置、Ｃ１…ＣＰＵ、Ｃ２…ＲＯＭ、Ｃ３…ＲＡＭ、Ｃ４…不揮発性ストレージ、Ｃ５…ネットワークインターフェイス表示部、Ｃ６…入力装置、Ｃ７…表示装置、Ｃ８…バス
【図1】
【図2】
【図3】
【図4】
【図5】
【図6】
【図7】
【図8】
【図9】
【図10】
【図11】
【図12】
【図13】
【図14】
【図15】
【図16】
【図17】
【図18】
【図19】