JP6318047B2

JP6318047B2 - 疾病予測装置、及びプログラム

Info

Publication number: JP6318047B2
Application number: JP2014164703A
Authority: JP
Inventors: 徳人今野
Original assignee: Nihon Kohden Corp
Current assignee: Nihon Kohden Corp
Priority date: 2014-08-13
Filing date: 2014-08-13
Publication date: 2018-04-25
Anticipated expiration: 2034-08-13
Also published as: JP2016040693A; US20160048649A1

Description

本発明は疾病予測装置、疾病予測方法、及びプログラムに関する。

近年、世界的に高齢化社会の問題が深刻になってきている。特に日本は、世界の中でも高齢化の問題が深刻である。日本の社会保障については、生産年齢人口３〜４人で高齢者１人を支える騎馬戦型から、生産年齢人口１人で高齢者１人を支えるおんぶ型にシフトするといわれている。このような社会構造の中で以下の事項を特に考慮する必要がある。

第１に、高齢者割合の増加に伴い、医療費の増大が顕著に進む可能性がある。そのため、疾病患者をすぐに治療して病院から早く出す（退院させる）ことが重要となる。第２に、高齢者の活用の観点から、病院に入れない（病気にならない）ということも重要となる。この２点への対応として、疾病が悪化する前に疾病の危険性を即座に把握することが重要となる。

ところで最近では、コンピュータの処理能力の向上に伴い、大量（Ｖｏｌｕｍｅ）のデータを高速（Ｖｅｌｏｃｉｔｙ）に、様々な形式（Ｖａｒｉｅｔｙ）で扱うことが可能となってきている。このような背景から様々な解析手法や機械学習、データマイニングといった技術が様々な分野で活用されてきており、医療の分野においても疾病の予測等にこれらの技術を用いることが検討されてきている。

そこで、統計解析等を用いた疾病の危険性を予測する先行技術について説明する。特許文献１には、被験者から取得した唾液データと予め記憶していた相関データを比較することにより生活習慣病の判定を行う装置が開示されている。特許文献２には、被験者の体重と健診データ（総コレステロール等）との相関を解析し、当該解析結果から健康状態を推定する手法が開示されている。

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

上述の特許文献１に記載の装置は、唾液データと相関データの比較を行うことにより、唾液の取得時における生活習慣病の危険度を判定している。また特許文献２に記載の手法は、体重を基に検診時の健康状態の判定を行うものである。すなわち特許文献１及び２の技術は、複数の生体パラメータの関係性（相関）を基に、ある時点での疾病や健康状態のリスクを分析するものである。換言すると特許文献１及び２の技術は、将来的な疾病の危険性等を事前に予測することができない。なお一般的な生体情報モニタは、各生体パラメータの変動状況から将来の予測を行う（例えば心拍数が急激に減少している場合に心臓異常を素早く察知する）ことはできるが、複数のパラメータが関連する疾病に関する予測を行うことはできない。

そのため、複数のパラメータが関連する疾病の発生予測を行うことが出来る技術の確立が望まれている。ここでパラメータとは生体パラメータ（血圧、呼吸数、体温、脈波、心拍数等）のみならず、環境因子（温度、湿度、照度、騒音等）、被験者の属性（性別、年齢、居住地等）等を含む概念である。

本発明は上記した課題を鑑みてなされたものであり、複数のパラメータが関連するような疾病を事前に予測することが出来る疾病予測装置、疾病予測方法、及びプログラムを提供することを主たる目的とする。

本発明にかかる疾病予測装置の一態様は、
被験者の生体状態を示す第１パラメータを取得する第１パラメータ取得部と、
前記被験者の生体状態を示す第２パラメータを取得する第２パラメータ取得部と、
前記第１パラメータと前記第２パラメータとの関係を示す統計値を算出する統計値算出部と、
前記第１パラメータと前記第２パラメータのから疾病の予兆を定義した定義情報を記憶する記憶部と、
前記統計値の経時的変化と、前記定義情報と、を基に前記被験者の疾病の予兆を解析する解析部と、を備える、ものである。

また本発明にかかる疾病予測装置の別の一態様は、
被験者の生体状態を示す第１パラメータを取得する第１パラメータ取得部と、
前記被験者を取り巻く環境因子または属性に関する第２パラメータを取得する第２パラメータ取得部と、
前記第１パラメータと前記第２パラメータとの関係を示す統計値を算出する統計算出部と、
前記第１パラメータと前記第２パラメータの医学的関係から導き出される疾病の予兆を定義した定義情報を記憶する記憶部と、
前記統計値の経時的変化と、前記疾病情報と、を基に前記被験者の疾病の予兆を解析する解析部と、を備える、ものである。

定義情報は、各生体パラメータ間の関係から疾病の予兆を定義している。解析部は、第１パラメータと第２パラメータから算出した統計値の変化を用いて解析を行う。当該統計値の変化は、被験者の身体の変化を示す有効な指標である。解析部は、この統計値の変化を定義情報と比較することにより、被験者の将来的な疾病発生の可能性を予測することが出来る。

本発明は、複数のパラメータが関連するような疾病を事前に予測することが出来る疾病予測装置、疾病予測方法、及びプログラムを提供することができる。

実施の形態１にかかる疾病予測装置１の構成を示すブロック図である。実施の形態１にかかる心拍数（ＨＲ）と脈拍数（ＰＲ）の関係を示す散布図である。実施の形態１にかかる定義情報の例を示す図である。実施の形態１にかかる心拍数（ＨＲ）と脈拍数（ＰＲ）から求めた傾きの変化と疾病予兆の関係を示した概念図である実施の形態１にかかる心拍数（ＨＲ）と呼吸数（ＲＲ）の関係を示す散布図である。

＜実施の形態１＞
以下、図面を参照して本発明の実施の形態について説明する。図１は、本実施の形態にかかる疾病予測装置の構成を示すブロック図である。疾病予測装置１は、好適には複数の生体パラメータを取得する医療機器であり、例えば生体情報モニタである。

疾病予測装置１は、第１パラメータ取得部１１、第２パラメータ取得部１２、統計値算出部１３、記憶部１４、解析部１５、及び出力部１６を備える。

第１パラメータ取得部１１は、被験者の生体から各種の生体パラメータを取得する。第１パラメータ取得部１１は、例えば被験者の生体に取り付けられる電極（図示せず）等と接続し、当該電極から生体信号を検出することにより生体パラメータを取得する。第１パラメータ取得部１１は、取得値を統計値算出部１３に供給する。ここで第１パラメータ取得部１１が取得する生体パラメータは、例えば呼吸数（ＲＲ：Respiration Rate）、心電図（ＥＣＧ）、脈拍数（ＰＲ：Pulse Rate）、心拍数（ＨＲ：Heart Rate）などである。

第２パラメータ取得部１２は、第１パラメータ取得部１１と同様に被験者の生体から各種の生体パラメータを取得する。ここで第２パラメータ取得部１２は、第１パラメータ取得部１１とは異なる生体パラメータを取得して、統計値算出部１３に取得値を供給する。

なお以下の説明では、第１パラメータ取得部１１が取得した生体パラメータを第１パラメータと呼称し、第２パラメータ取得部１２が取得した生体パラメータを第２パラメータと呼称する。

統計値算出部１３は、第１パラメータ取得部１１の取得値と、第２パラメータ取得部１２の取得値と、を基に統計値を算出する。ここで統計値とは、例えば相関係数、標準偏差（ＳＤ：Standard Deviation）、両パラメータの取得値を座標配置した際の近似式の傾き等である。図２を参照して、統計値の意味について詳細に説明する。図２において第１パラメータが心拍数（ＨＲ）であり、第２パラメータが脈拍数（ＰＲ）であるものとする。

図２に示すように各時点における心拍数（ＨＲ）と脈拍数（ＰＲ）の値を組とした座標系を定義する。例えば図２の（ア）は、ある時点での心拍数（ＨＲ）が６０であり、脈拍数（ＰＲ）が７６であることを意味している。そして、統計値算出部１３は、定期的に相関係数、標準偏差、近似式（ＨＲ／ＰＲ）の傾きを算出する。ここで統計値算出部１３は、解析部１５の解析内容に応じて算出する統計値の種別を切り替える。例えば統計値算出部１３は、標準偏差が基準を越えればノイズと判定して処理を終了し（図３（Ａ）（ウ）を用いて後述する。）、そうでなければ近似式(ＨＲ／ＰＲ)の傾きも算出する。統計値算出部１３は、算出した統計値を解析部１５に供給する。

記憶部１４は、第１パラメータと第２パラメータの関係から疾病の予兆を定義した定義情報を記憶する。記憶部１４は、好適には疾病予測装置１内の２次記憶装置（例えばハードディスク等）である。なお記憶部１４は、疾病予測装置１に着脱可能な単体装置（例えばＵＳＢ（Universal Serial Bus）メモリ等）であってもよい。

記憶部１４が記憶する定義情報は、例えば第１パラメータと第２パラメータから算出された統計値の経時的変化の傾向と疾病発生のリスクを定義したものである。具体例は解析部１５の処理と共に後述する。

解析部１５は、記憶部１４に記憶された定義情報と、統計値算出部１３が算出した統計値（相関係数、標準偏差、近似式の傾き）と、に基づいて被験者の疾病の予兆を解析する。換言すると解析部１５は、表面的には異常値ではないものの、将来的に疾病（例えば心疾患）の危険性があることを予測するものである。当該解析処理の詳細は、具体的なデータ例を用いて後述する。

出力部１６は、音または表示により使用者（医師、看護師、被験者等）に疾病の予兆があることを通知する。出力部１６は、一般的な生体情報モニタの液晶モニタやスピーカ及びその周辺回路等である。なお出力部１６は、疾病予測装置１内の通信機能と連携して疾病の予兆を検出したことを他の端末（例えば担当看護師の携帯端末）に送信するような出力を行ってもよい。

続いて解析部１５による解析処理の具体例を説明する。

（解析部１５の解析例（心疾患の予測））
はじめに記憶部１４に記憶された定義情報の例を図３（Ａ）及び図３（Ｂ）を参照して説明する。図３（Ａ）は、第１パラメータが心拍数（ＨＲ）であり、第２パラメータが脈拍数（ＰＲ）である場合の両パラメータの関係性を定める定義情報である。心拍数（ＨＲ）は、心電図（ＥＣＧ）のＱＲＳ波を基に算出した心臓の拍動数（１分間の拍動回数）である。また脈拍数（ＰＲ）は、動脈の拍動の変化を末梢部位（例えば指先）に取り付けたプローブでカウントした１分間の数値である。両パラメータは心臓の拍動に起因する数値であるため、同一被験者から算出した心拍数（ＨＲ）と脈拍数（ＰＲ）は原理的には同数となる。しかし、様々な理由から両者が一致しない場合がある。

換言すると、理想的には心拍数（ＨＲ）と脈拍数（ＰＲ）から算出した相関係数は１に限りなく近い。また理想的には、心拍数（ＨＲ）と脈拍数（ＰＲ）から算出した標準偏差は、ばらつきが小さいために限りなく小さい値となる。ここで相関係数も減少し、脈拍数（ＰＲ）の平均値（以下、ＰＲａｖｅとも記載する。）に０．１（１０％）を乗じた値よりも脈拍数（ＰＲ）の標準偏差（ＳＤ）が大きくなっている場合（｜ＳＤ｜＞＝０．１＊ＰＲａｖｅ）、指先から測定することが多い脈拍数（ＰＲ）の値が体動ノイズによって変動していることが想定される（図３（Ａ）（ウ））。そのため、この場合には疾病発生の予兆とは判断しないことが定義されている。なお図３（Ａ）のケースでは体動ノイズに関する定義のみが記載されているが、必ずしもこれに限られず、ＳｐＯ２プローブの接触不良や外来光等に起因するノイズを考慮した定義が行われても勿論良い。以下、ノイズ以外のケース（図３（Ａ）（ア）（イ））について検討する。

相関係数の値が減少している場合、何らかの要因により両パラメータの値の差が大きくなってきたことが想定される。この場合において、両パラメータの最新データ及び最新データから遡った所定数のデータを用いて求めた近似式の傾き（例えば最小二乗法で算出した１次式の傾きである。以下の例では、「近似１次式（ＨＲ／ＰＲ）の傾き」と表現する。）が大きくなってきた場合、１つの脈拍がカウントされる間に、心臓が異常に速く駆動して複数回の心拍がカウントされたことが想定される。すなわち、頻拍等の心臓異常の危険が想定される（図３（Ａ）（ア））。また、一つの心拍に対して脈拍が非常に遅いとも捉えることができ、このような観点からも心肺機能の循環不全等の心臓異常が疑われる（図３（Ａ）（ア））。

一方、心拍数（ＨＲ）と脈拍数（ＰＲ）の相関係数の値が減少している場合において、近似１次式（ＨＲ／ＰＲ）の傾きが大きくなってきた場合を検討する。この場合、１つの脈拍に対して、心拍数としてカウントされない拍動があることが想定される。そのため、不整脈の疑いがあることを意味している（図３（Ａ）（イ））。

定義情報（図３（Ａ））は、このように第１パラメータと第２パラメータの関係から、統計値がどのように変化した場合に、どのような疾病の危険性があるかを定義したものである。換言すると第１パラメータと第２パラメータから算出された統計値（相関係数、標準偏差、近似１次式（ＨＲ／ＰＲ）の傾き）の経時的変化と疾病の予兆の関係を定めるものである。

図４は、心拍数（ＨＲ）と脈拍数（ＰＲ）から求めた近似１次式（ＨＲ／ＰＲ）の傾きの変化と疾病予兆の関係を示した概念図である。ノイズ等の影響を受けておらず、かつ身体の状態に異常がない場合には、心拍数（ＨＲ）と脈拍数（ＰＲ）はほぼ同じ値となる。これに対して近似１次式（ＨＲ／ＰＲ）の傾きが増加した場合、上述のように心臓異常の予兆と捉えることが出来る。一方、近似１次式（ＨＲ／ＰＲ）の傾きが減少した場合、上述のように不整脈等の予兆と捉えることが出来る。

解析部１５は、この定義情報（図３（Ａ））を基に疾病の予兆の予測を行う。前提として統計値算出部１３は、心拍数（ＨＲ）と脈拍数（ＰＲ）の相関係数、標準偏差、及び近似１次式（ＨＲ／ＰＲ）の傾きを逐次算出する。解析部１５は、心拍数（ＨＲ）と脈拍数（ＰＲ）のデータを蓄積していき、蓄積データを用いて解析を行う。

統計値算出部１３による統計値の算出処理の詳細について説明する。統計値算出部１３は、一般的な相関係数算出式を用いて心拍数（ＨＲ）と脈拍数（ＰＲ）の相関係数を算出する。

また統計値算出部１３は、一般的な標準偏差算出式に従って脈拍数（ＰＲ）の標準偏差を算出する。

更に統計値算出部１３は、最小二乗法等を用いて心拍数（ＨＲ）と脈拍数（ＰＲ）を１次関数で表した場合の傾きを算出する。例えば統計値算出部１３は、図２において最新のデータから１０個のデータを取り出し、当該１０個のデータを用いて最小二乗法を行って近似１次式（ＨＲ／ＰＲ）の傾きを算出する。

統計値算出部１３は、このように最新のデータから遡って所定数以上のデータを用いて統計値（相関係数、標準偏差、近似１次式（ＨＲ／ＰＲ）の傾き）を算出することが好ましい。例えば最新のデータと１つ前のデータのみを用いて近似１次式（ＨＲ／ＰＲ）の傾き等を求めた場合、体動の影響等を強く受けてしまい、誤った値を算出してしまう恐れがある。しかし所定数以上のデータを用いることにより、統計値算出部１３はこのような外れ値の影響を小さくすることができ、精度の良い統計値を算出することが出来る。

解析部１５は、定義情報（図３（Ａ））を参照し、統計値算出部１３が算出した統計値の経時的変化と定義情報（図３（Ａ））を比較して疾病予兆の解析を行う。例えば相関係数が減少し、かつ脈拍数（ＰＲ）の標準偏差が平均脈拍数（ＰＲ）の１０％未満であり、近似１次式（ＨＲ／ＰＲ）の傾きが増加している場合、解析部１５は心臓異常の恐れがあると解析する。

なお解析部１５は、定義情報に合致した状態となった場合に即座に疾病の予兆があると判定するのではなく、一定時間を超えて定義情報に合致した状態となった場合に被験者に疾病の予兆があると判定してもよい。これにより解析部１５は、ノイズ等の影響によって一時的に測定値が変化したような場合を適切にキャンセルできるになり、より精度の高い疾病予兆の解析を実現できる。

上述の処理では、統計値として（相関係数、標準偏差、近似１次式（ＨＲ／ＰＲ）の傾き）を全て用いたが、必ずしもこれに限らない。例えば解析部１５は、相関係数と標準偏差のみを求めて近似式を使用せずに解析処理を行ってもよい。

心拍数（ＨＲ）と脈拍数（ＰＲ）は同じ値になることが一般的であるため、近似１次式（ＨＲ／ＰＲ）の傾きが変化することにより、何らかの異常があることがわかる。また解析部１５は、近似１次式（ＨＲ／ＰＲ）の傾きの変化の方向により、どのような異常が生じているかを判断することが出来る。何らかの身体的な異常が生じてきた場合には心拍数（ＨＲ）と脈拍数（ＰＲ）の値にバラつきが生じ始めるため、相関係数の値が下がることが想定される。このため解析部１５は、標準偏差や相関係数の値も考慮したうえで判断することにより更に精度良く疾病（主に心疾患）の予測を行うことが出来る。なお、近似１次式（ＨＲ／ＰＲ）の傾きを用いることなく、相関係数や標準偏差のみを用いて解析を行うことも可能である。上述のように心拍数（ＨＲ）と脈拍数（ＰＲ）は同じ値になることが理想であるため、両パラメータのバラツキが大きくなってきた場合、原因は明らかではないものの、何らかの異常が生じたと解析することが出来るためである。

（解析部１５の解析例（心拍数と呼吸を用いた予測））
図３（Ａ）では心拍数（ＨＲ）と脈拍数（ＰＲ）の関係から疾病予測を行う定義情報を用いたが、その他の生体パラメータを用いて同様の解析を行うこともできる。図３（Ｂ）は、心拍数（ＨＲ）と呼吸数（ＲＲ）と疾病の発生の関係を定義した定義情報である。図５は、ある被験者の心拍数（ＨＲ）と呼吸数（ＲＲ）の関係をプロットした図である。図５に示すように一般的に心拍数（ＨＲ）と呼吸数（ＲＲ）の相関は小さい。換言すると両パラメータの関連は小さい。

しかしながら、相関係数が増加し（相関のある方向に値が変化し）、近似１次式（ＲＲ／ＨＲ）の傾きが増加している場合、一つの心拍数が測定される間に呼吸数が異常に多く検出されていることを意味している。すなわちこのケースでは過呼吸等の恐れがある（図３（Ｂ）（ア））。そのため、リハビリテーションの実行時にこの様な状況を検出した場合、運動過多の恐れがあり、運動停止命令を音出力することやリハビリプログラムを見直す等の対応を取るべきである。

一方、相関係数が増加し（相関のある方向に値が変化し）、近似１次式（ＲＲ／ＨＲ）の傾きが減少している場合、一つの呼吸が測定される間に心拍が異常に多く検出されていることを意味している。この場合には、心拍数のダブルカウント等の恐れ、すなわち運動性の脚ブロック等の不整脈の恐れがある。そのため、例えばリハビリテーションの実行中等にこのような状況を検出した場合、不整脈の予兆があるため、運動停止命令を音出力することやリハビリプログラムを見直す等の対応を取るべきである。

解析部１５は、この定義情報（図３（Ｂ））と統計値算出部１３が算出した統計値を比較して疾病予兆の解析を行う。例えば相関係数が増加し、近似１次式（ＲＲ／ＨＲ）の傾きが増加している場合、解析部１５は運動過多の恐れがあると解析する。なお本例においても解析部１５は、近似１次式（ＲＲ／ＨＲ）の傾きのみを用いて簡易的に疾病の予兆の解析を行うことも可能である。

以上が解析部１５による解析処理の具体例である。なお定義情報は図３に示すものに限られず、他の生体パラメータを用いて定義することも勿論可能である。また疾病予測装置１の使用開始後に、ユーザが図示しない入力装置（マウス、キーボード等）を用いて新たな定義情報を定義することも可能である。

続いて本実施の形態にかかる疾病予測装置１の効果について説明する。定義情報は、各生体パラメータ間の関係から疾病の予兆を定義している。解析部１５は、第１パラメータと第２パラメータから算出した統計値の変化を用いて解析を行う。当該統計値の変化は、被験者の身体の変化を示す有効な指標である。解析部１５は、この統計値の変化を定義情報と比較することにより、被験者の将来的な疾病発生の可能性を予測することが出来る。疾病発生の可能性を事前に予測できることにより、各生体パラメータが異常値（例えば心電図がＶＦ状態）となる前であってもアラームによる通知等をすることが出来る。

より詳細には、解析部１５は第１パラメータと第２パラメータの近似式の傾きの変化に応じて疾病の予兆を解析している（図３（Ａ）等）。当該近似式の傾きの算出は計算量の少ない処理であり、疾病予測装置１が生体情報モニタ等の装置である場合にも容易に組み込むことが出来る処理である。

更に解析部１５は、第１パラメータと第２パラメータの相関係数や標準偏差の変化も考慮して疾病の予兆を解析している。相関係数は両パラメータの相関関係を示し、標準偏差はデータのバラツキを定義する。解析部１５は、相関係数の変化や標準偏差の変化を扱うことにより、身体状態が変化してきたことを客観的に検出できるようになる。

以上、本発明者によってなされた発明を実施の形態に基づき具体的に説明したが、本発明は既に述べた実施の形態に限定されるものではなく、その要旨を逸脱しない範囲において種々の変更が可能であることはいうまでもない。

例えば上述の例では第１パラメータ及び第２パラメータが被験者の生体から取得する生体パラメータであるものとして説明したが、必ずしもこれに限られない。すなわち第１パラメータ及び第２パラメータの一方が環境因子（気温、湿度、照度、騒音）や属性情報（性別、年齢、居住地）等であってもよい。この場合であっても、第１パラメータと第２パラメータの関係を医学的な見地から定義情報として定義することにより、疾病予兆の解析を行うことが出来る。

上記の統計値算出部１３及び解析部１５の各処理は、疾病予測装置１内で動作するコンピュータプログラムとして実現することができる。すなわち疾病予測装置１は、一般的なコンピュータが備えるＣＰＵ（Central Processing Unit）、ハードディスクドライブ、キャッシュメモリ等も備えているものとする。

プログラムは、様々なタイプの非一時的なコンピュータ可読媒体（non-transitory computer readable medium）を用いて格納され、コンピュータに供給することができる。非一時的なコンピュータ可読媒体は、様々なタイプの実体のある記録媒体（tangible storage medium）を含む。非一時的なコンピュータ可読媒体の例は、磁気記録媒体（例えばフレキシブルディスク、磁気テープ、ハードディスクドライブ）、光磁気記録媒体（例えば光磁気ディスク）、ＣＤ−ＲＯＭ（Read Only Memory）、ＣＤ−Ｒ、ＣＤ−Ｒ／Ｗ、半導体メモリ（例えば、マスクＲＯＭ、ＰＲＯＭ（Programmable ROM）、ＥＰＲＯＭ（Erasable PROM）、フラッシュＲＯＭ、ＲＡＭ（random access memory））を含む。また、プログラムは、様々なタイプの一時的なコンピュータ可読媒体（transitory computer readable medium）によってコンピュータに供給されてもよい。一時的なコンピュータ可読媒体の例は、電気信号、光信号、及び電磁波を含む。一時的なコンピュータ可読媒体は、電線及び光ファイバ等の有線通信路、又は無線通信路を介して、プログラムをコンピュータに供給できる。

１疾病予測装置
１１第１パラメータ取得部
１２第２パラメータ取得部
１３統計値算出部
１４記憶部
１５解析部
１６出力部

Claims

被験者の生体状態を示す第１パラメータとして心拍数（ＨＲ）を取得する第１パラメータ取得部と、
前記被験者の生体状態を示す第２パラメータとして脈拍数（ＰＲ）を取得する第２パラメータ取得部と、
前記第１パラメータである心拍数（ＨＲ）と前記第２パラメータである脈拍数（ＰＲ）との関係を示す近似１次式（ＨＲ／ＰＲ）を算出する統計値算出部と、
前記近似１次式（ＨＲ／ＰＲ）の傾きから疾病の予兆を定義した定義情報を記憶する記憶部と、
前記近似１次式（ＨＲ／ＰＲ）の傾きが増加した場合に前記被験者に心臓異常の恐れがあると解析し、前記近似１次式（ＨＲ／ＰＲ）の傾きが減少した場合に前記被験者に不整脈の恐れがあると解析する解析部と、を備える疾病予測装置。
前記統計値算出部は、前記第１パラメータ及び前記第２パラメータの最新データから順に所定数以上のデータを抽出し、抽出したデータを用いて前記近似１次式（ＨＲ／ＰＲ）を算出する、ことを特徴とする請求項１に記載の疾病予測装置。
前記解析部は、前記近似１次式（ＨＲ／ＰＲ）の変化が前記定義情報に定義する状態となってから一定時間経過した場合に、前記被験者に疾病の予兆があると判断する、ことを特徴とする請求項１または請求項２に記載の疾病予測装置。
コンピュータに、
被験者の生体状態を示す第１パラメータとして心拍数（ＨＲ）の値と、前記被験者の生体状態を示す第２パラメータとして脈拍数（ＰＲ）の値と、の関係を示す近似１次式（ＨＲ／ＰＲ）を算出する統計算出ステップと、
前記近似１次式（ＨＲ／ＰＲ）の傾きが増加した場合に前記被験者に心臓異常の恐れがあると解析し、前記近似１次式（ＨＲ／ＰＲ）の傾きが減少した場合に前記被験者に不整脈の恐れがあると解析する解析ステップと、
を実行させる、プログラム。