(57)【要約】  本発明は、ＥＣＧ信号データのみから無呼吸または低呼吸の発生を確定する方法および装置を提供する。本方法は、睡眠中の被検者からＥＣＧ信号を取得し、それら信号をデータに変換し、呼吸努力を推定するのに適切なＥＣＧ特徴を抽出して、無呼吸および低呼吸の呼吸イベント特徴を確定するように構成された装置によって実行される。抽出されたＥＣＧ特徴は、呼吸努力の相関であり、呼吸イベントの代用測度として用いられる。本方法は、ＡＨＩすなわち無呼吸／低呼吸指数を計算することを含むことができる。本方法は、無呼吸を閉塞型無呼吸または中枢型無呼吸に分類することができる。
2010540124.tif 000002
【特許請求の範囲】
【請求項１】
ＥＣＧ信号から呼吸イベントを確定する方法であって、
ＥＣＧ信号データを取得するステップと、
前記ＥＣＧ信号データから波形形態データを抽出するステップと、
前記波形形態データから呼吸努力を推定するステップと、
を含む方法。
【請求項２】
前記呼吸努力のピークを交差検定するステップをさらに含む、請求項１に記載の方法。
【請求項３】
前記呼吸努力から呼吸パターンを特徴付けるステップと、
観測するために前記呼吸イベントを検出するステップと、
をさらに含む、請求項１または２に記載の方法。
【請求項４】
呼吸イベントを、無呼吸または低呼吸のうちの一方として分類するステップをさらに含む、請求項１〜３のいずれか一項に記載の方法。
【請求項５】
呼吸イベントを、閉塞型無呼吸または中枢型無呼吸のうちの一方として分類するステップをさらに含む、請求項１〜４のいずれか一項に記載の方法。
【請求項６】
無呼吸／低呼吸指数を計算するステップをさらに含む、請求項４または５に記載の方法。
【請求項７】
前記波形データが、ＥＣＧから導出された呼吸信号、心拍数、あるいはＲ信号の面積または振幅または周期性を含む、ＥＣＧから導出された生理学的予測因子のうちの任意の１つまたは組合せに変換される、請求項１に記載の方法。
【請求項８】
睡眠状態にある被検者に対して実施される、請求項１に記載の方法。
【請求項９】
睡眠呼吸障害を解析する方法であって、
被検者に隣接するＥＣＧセンサから生体信号を取得するステップと、
前記生体信号をデータとしてコンピュータファイルに格納するステップと、
前記生体信号データから心拍数および波形形態データを抽出するステップと、
抽出された心拍数およびＥＣＧ波形形態データから生理学的予測因子を導出するステップと、
前記導出された生理学的予測因子から呼吸イベントの開始および終了を確定するステップと、
を含む方法。
【請求項１０】
呼吸イベントを無呼吸または低呼吸として分類するステップをさらに含む、請求項９に記載の方法。
【請求項１１】
無呼吸／低呼吸指数を計算するステップと、
前記無呼吸／低呼吸指数を表示するステップと、
をさらに含む、請求項９に記載の方法。
【請求項１２】
ＥＣＧ信号から呼吸努力を確定する装置であって、
被検者から少なくとも２つの直交する信号を取得する少なくとも２つのＥＣＧ誘導と、
前記信号をデジタルデータに変換する手段と、
電子データ記憶手段と、
デジタルＥＣＧデータから波形を抽出し、呼吸イベントおよび呼吸努力を確定するようにプログラムされたマイクロプロセッサと、
を具備する装置。
【請求項１３】
無呼吸／低呼吸指数を計算するようにプログラムされたマイクロプロセッサをさらに具備する、請求項１２に記載の呼吸努力を確定する装置。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、患者から、詳細には睡眠呼吸障害、より詳細には睡眠時無呼吸、さらに詳細には閉塞型睡眠時無呼吸または中枢型睡眠時無呼吸を有する患者から、心電図データを取得しそれを分析する方法および装置に関する。
【背景技術】
【０００２】
呼吸が停止する、閉塞型睡眠時無呼吸（ＯＳＡ）および／または中枢型睡眠時無呼吸（ＣＳＡ）を含む睡眠時無呼吸（ＳＡ）は、人口のかなりの割合に影響を及ぼしている。ＯＳＡの影響は、気道が虚脱し空気が通ることができなくなる場合に出る。ＣＳＡの影響は、患者が、中枢神経系の機能に関連する原因により空気を吸い込まずかつ吐き出さない場合に出る。ＯＳＡおよびＣＳＡはともに、非効率な睡眠、さらには医学的な問題をもたらす可能性がある。
【０００３】
ＯＳＡおよびＣＳＡはともに、目下、睡眠検査室において終夜ポリグラフ（ＰＳＧ）検査によって診断される。そこでは、被検者は、さまざまな生理学的パラメータを測定するために身体に多数の電極が取り付けられた状態で睡眠をとる。こうしたＰＳＧ検査は、被検者が診療所で終夜睡眠をとる必要があるため、受けるのに費用がかかる。
【０００４】
ＳＡの発生およびタイプを測定する装置および方法を開発するために、いくつかの取組みがなされてきた。ＳＡを診断する標準ＰＳＧ基準では、個々の呼吸イベントの実際の開始時点、持続時間および部類を検出することが必要である。そして、この情報を、無呼吸低呼吸指数（ＡＨＩ）を生成することによって睡眠時無呼吸の存在および深刻さを評価するために用いる［１］。目下、最もうまくいっている、報告されている方法は、被検者が睡眠時無呼吸を有しているかを（ＡＨＩ値の確実な推定なしに）検出するか、または任意の瞬間の間に睡眠時無呼吸イベントがあるかを検出することができる［２〜５］。たとえば、米国特許第５，７６９，０８４号明細書は、ＳＡの存在を確定するために、心電図（ＥＣＧ）を測定することを含む、さまざまな睡眠パラメータのカオス処理を用いる方法を教示している。しかしながら、その方法は、ＣＳＡからＯＳＡを識別することまでは及んでいない。被検者の種々の面に配置されたセンサからＳＡを測定する方法が、米国特許第６，４１５，１７４号明細書に教示されている。この文書は、ＣＳＡとＯＳＡとを識別することができる方法は教示していない。参照により本明細書に援用される国際公開第２００５／０６７７９０号パンフレットにおいて、Ｂｕｒｔｏｎは、ＥＣＧトレースを用いて、ＳＡの症例を特定しＳＡをＯＳＡかまたはＣＳＡのいずれかとして分類することができることを指摘している。しかしながら、そのＣＳＡ信号は、背景雑音に比較して小さく、すなわち信号対雑音比が小さい。
【発明の概要】
【発明が解決しようとする課題】
【０００５】
心臓病患者の間ではＯＳＡおよびＣＳＡを有する人が多く、ＡＨＩを確定するためのＥＣＧデータの使用は、ＰＳＧデータの使用よりさらに魅力のある選択肢になっている。必要なのは、ＰＳＧ検査より簡単な装置を用いて、睡眠時無呼吸およびその種類を測定するより便利かつ確実な方法である。その方法は、個々の呼吸イベントを検出し分類することができなければならない。こうした方法が、家庭で行われることが好ましい［６〜８］。
【課題を解決するための手段】
【０００６】
本発明は、ＥＣＧから導出されたパラメータが呼吸イベントと相関され、ＥＣＧから導出された相関を用いて、特に睡眠中の被検者において呼吸イベントを特定することができる、という意外な観測結果を活用する。特定された呼吸イベントを用いて、呼吸努力を確定することができる。呼吸努力の確定を、無呼吸／低呼吸指数の計算に用いることができる。ＥＣＧから導出されたデータを用いて、ＯＳＡとＣＳＡとを識別することも可能である。本発明の目的は、ＥＣＧ信号データを用いて呼吸イベントを特定する方法を提供することである。本発明のさらなる目的は、ＥＣＧ信号データを用いて呼吸努力を確定する方法を提供することである。本発明のさらなる目的は、ＥＣＧ信号データから無呼吸／低呼吸指数を計算する方法を提供することである。
【０００７】
一態様では、本発明は、ＥＣＧ信号から呼吸イベントを確定する方法であって、ＥＣＧ信号データを取得するステップと、前記ＥＣＧ信号データから波形形態データを抽出するステップと、前記波形形態データから呼吸努力を推定するステップとを含む方法を提供する。本方法は、呼吸努力のピークを交差検定するステップを含むことができる。本方法は、前記呼吸努力から呼吸パターンを特徴付けるステップと、観測するために呼吸イベントを検出するステップと、をさらに含むことができる。本方法は、呼吸イベントを、無呼吸または低呼吸のうちの一方として分類するステップを含むことができる。本方法は、呼吸イベントを、閉塞型無呼吸または中枢型無呼吸のうちの一方として分類するステップを含むことができる。本方法は、無呼吸／低呼吸指数を計算するステップを含むことができる。波形データを、ＥＣＧから導出された呼吸信号、心拍数、あるいはＲ信号の面積または振幅または周期性を含む、ＥＣＧから導出された生理学的予測因子のうちの任意の１つまたは組合せに変換することができる。本方法は、睡眠状態にある被検者に対して実施されることが好ましい。
【０００８】
さらなる態様では、本発明は、睡眠呼吸障害を解析する方法であって、ＥＣＧセンサから生体信号を取得するステップと、生体信号をデータとしてコンピュータファイルに格納するステップと、前記生体信号データから心拍数および波形形態データを抽出するステップと、抽出された心拍数およびＥＣＧ波形形態データから生理学的予測因子を導出するステップと、前記導出された生理学的予測因子から呼吸イベントの開始および終了を確定するステップと、を含む方法を提供する。本方法は、１時間当りの１０秒間より長い呼吸イベントの平均解数から無呼吸／低呼吸指数を計算するステップを含むことができる。本方法は、無呼吸／低呼吸指数を表示するステップを含むことができる。
【０００９】
さらなる態様では、本発明は、ＥＣＧ信号から呼吸努力を確定する装置であって、被検者から少なくとも２つの直交する信号を取得する少なくとも２つのＥＣＧ誘導と、前記信号をデジタルデータに変換する手段と、電子データ記憶手段と、デジタルＥＣＧデータから波形を抽出し、呼吸イベントおよび呼吸努力を確定するようにプログラムされたマイクロプロセッサと、を備える装置を提供する。本装置は、無呼吸／低呼吸指数を計算するようにプログラムされたマイクロプロセッサをさらに備えることができる。
【図面の簡単な説明】
【００１０】
【図１】ＥＣＧデータ群から無呼吸／低呼吸指数を計算するステップの概略フローチャートを示す。
【図２】ＥＣＧデータ群から無呼吸／低呼吸指数を計算する詳細なステップのフローチャートを示す。
【図３】グラフの各部分に名称が付されている典型的なＥＣＧ誘導Ｉ群の出力のグラフを示す。
【図４】睡眠中の規則的な安定した呼吸を示す、被検者のＱＲＳ面積、Ｒ振幅および呼吸流量の時間に対するＥＣＧ信号のグラフを示す。
【図５】呼吸イベントを示す、被検者のＱＲＳ面積、Ｒ振幅および呼吸流量の時間に対するＥＣＧ信号のグラフを示す。
【図６】睡眠中の呼吸イベントにおけるＥＣＧ信号のＱＲＳ群数に対するＲ波のｍＶのグラフを示す。
【図７】睡眠中の一続きの呼吸イベントにわたる、ＥＣＧ信号におけるＱＲＳ群数に対する、時間に対するＲ−Ｒ間隔および呼吸の流量のｍＶのグラフを示す。
【図８】睡眠中の無呼吸の期間および正常な呼吸の期間中の、ＥＣＧ信号におけるＱＲＳ面積、Ｒ振幅および流量の時間に対する信号のグラフを示す。
【図９】呼吸イベント中のＥＣＧ信号におけるＱＲＳ群数に対するＲ波の振幅のグラフを示す。
【図１０】呼吸イベント中のＥＣＧ信号におけるＱＲＳ群数に対するＱＲＳ領域の面積のグラフを示す。
【図１１】呼吸イベントにわたるＥＣＧ信号におけるＱＲＳ群数に対するＲ−Ｒ間隔のグラフを示す。
【図１２】ＥＣＧ信号における呼吸イベントを通してのＱＲＳ群数に対するＱＲＳ面積のグラフを示す。
【図１３】図１３ａは第１ＥＣＧ誘導における呼吸イベントを通してのＲ波を示す出力信号を示し、図１３ｂは図１２ａに示す呼吸イベントを通しての第２ＥＣＧ誘導からのＲ波を示す出力信号を示す。
【図１４】呼吸イベントを通してのＥＣＧ信号のＱＲＳ数に対するＲ−Ｒ間隔のグラフを示す。
【図１５】呼吸イベントを通してのＥＣＧ信号におけるＱＲＳ群数に対するＥＤＲ信号のグラフを示す。
【図１６】呼吸イベントを通してのＥＣＧ信号におけるＱＲＳ群数に対するＲ波の振幅のグラフを示す。
【図１７ａ】ＯＳＡイベントの場合の時間の経過によるＥＣＧ信号のＱＲＳ面積およびＲ波振幅のグラフを示す。
【図１７ｂ】ＣＳＡイベントの場合の時間の経過によるＥＣＧ信号のＱＲＳ面積およびＲ波振幅のグラフを示す。
【発明を実施するための形態】
【００１１】
本発明は、被検者から取得されたＥＣＧ信号データから、無呼吸／低呼吸指数を確定する以前は未知であった方法を提供する。本発明の方法は、睡眠期間の無呼吸／低呼吸指数を計算するために、睡眠中に取得されたＥＣＧ信号の変動を有利に活用する。
【００１２】
種々のＥＣＧ誘導が、身体の種々の部位における電位の差を測定する。呼吸は、電極間のインピーダンス、したがって電位およびそれらの間の差に影響を与える。胸部および腹部の動きは、呼吸イベント中に変化し、したがって、電極間のインピーダンスの呼吸努力による変化の仕方を変化させる。これにより、ＥＣＧから導出された信号の間の位相変化を観察することができるようになる。
【００１３】
本方法は、呼吸イベントに対するマーカとして位相変化を用いる。本方法は、交感神経活動の亢進を伴う心拍数、収縮期血圧および拡張期血圧が上昇する短い期間がある、呼吸イベント後かつ覚醒中の短い期間を有利に活用する。無呼吸イベントは、通常、覚醒で終了し、覚醒のマーカを無呼吸終了のマーカとして用いることが可能となる。急激な心拍数の上昇は、ＥＣＧから容易に検出可能である。後述するようにＲ−Ｒ間隔の急な低下を、呼吸イベントの終了のマーカとして用いる。
【００１４】
睡眠時無呼吸は、本質的に周期的であることが多く、この結果、抽出されたＥＣＧ信号の多くもまた同様に周期的パターンを示す可能性があり、これを特に、図８において呼気の流量２の最下トレースに示す。ＱＲＳ面積およびＲ波振幅もまた、同じ周波数の周期的パターンを有することを観測することができる。（ＥＣＧ信号のパラメータを、図３においてトレース１として示す。）呼吸信号がより長期間の周期的パターンを展開するに従い、それに伴う、抽出されたＥＣＧ信号もまた、同じ周波数で周期的パターンを展開する。これを、本発明の方法において呼吸信号の推定される代用物として用いることができる。本方法では、これらパターンを認識するために、フーリエ解析技法およびデジタルフィルタリングを用いる平滑化技法を用いることができる。これら周期的パターンを含むことが分かった包絡線を用いて、本明細書で説明するようにより詳細な検査のためにデータの領域を強調する。
【００１５】
本方法は、ＥＣＧ信号データにおいて測定可能な特徴を有する呼吸イベントの開始点および終了点を特定する。最も好都合には、ＥＣＧ信号データは、被検者の睡眠中に、被検者に隣接するセンサによって取得される。データは、収集されコンピュータファイルに格納され、マイクロプロセッサおよびコンピュータプログラムによって、処理され、呼吸イベントが特定され、呼吸努力が確定され、または無呼吸／低呼吸指数が計算される。
【００１６】
表１は、ＥＣＧ信号データにおいて測定可能なパラメータのいくつかと本明細書で用いるそれらの頭字語との定義を提供する。パラメータは、ＥＣＧ信号を取得し解析する技術分野における当業者には周知であろう。本方法は、表１のパラメータに限定されず、呼吸努力との相関関係がある、ＥＣＧから導出されるいかなる信号を含んでもよい。
【００１７】
本発明の方法のステップの概略を図１のフローチャートに示す。本方法をさらに詳述するより詳細なフローチャートを図２に示す。典型的なＥＣＧトレース１を、図３において、時間に対する取得された信号の大きさのグラフとして示す。本方法の計算はすべて、マイクロプロセッサを用いて行われる。図１を参照すると、第１ステップにおいて、ＥＣＧホルター等の標準的なＥＣＧ監視機器を用いて、被検者からＥＣＧ信号を記録する。デジタルＥＣＧ信号を記録することができるいかなる好適な機器を使用してもよい。取得されたデータは、解析のためにコンピュータ媒体上のデジタルファイルに格納されることが好都合である。データを、既知の技法を用いてＥＣＧデータに対する適切な帯域幅に従ってフィルタリングしてもよい。
【００１８】
ＥＣＧ信号１は、時間の経過により、図３に示す特徴的な波形形態を組み込んだセグメントを有し、それらは、本方法の第２ステップで特定される。特に、点Ｑで開始し、Ｒで最大に到達し、Ｓで最小に到達するピークは、「ＱＲＳ」群として周知である。本方法で用いかつ図に示すさらなるパラメータは、各ＥＣＧ群のＴ波の振幅およびＴ波のグラフの下の領域と、Ｒにおける信号の振幅と、連続するＲ波間の期間と、ＱＲＳセグメントのグラフの下の面積とである。心拍数は、Ｒ−Ｒ間隔（ＲＲＩ）として知られる、図３に示すＲ波ピーク間の時間として測定される。
【００１９】
第３ステップにおいて、ＥＣＧから導出された信号を用いて、呼吸イベントの複数のＥＣＧから導出された相関を生成する。これら相関を用いて、無呼吸イベントまたは低呼吸イベントを含む、呼吸イベントの期間中にＥＣＧにおいて発生する変化を確認する。本方法で用いることができる相関の各々を本明細書で説明する。本発明の範囲は、呼吸イベントと相関するＥＣＧから導出されたいかなるパラメータの使用をも含む。そして、特定されたＥＣＧパラメータを、正常なまたは異常な心機能を表すものとして分類する。この解析では、線形近似を用いて、ＥＣＧ群のさまざまな領域を特定する。しかしながら、他の適切な解析を用いてもよい。分類は、第２ステップで特定されたＥＣＧパラメータに対する変化を検査することによって行う。
【００２０】
正常であるかまたは異常であるというパラメータの分類は、既知の相関関係を用いる。心拍数の分類に対し、本方法は、以下の相関関係を組み込む。呼吸性洞性不整脈（ＲＳＡ）は、哺乳動物において、ＥＣＧの心拍数（ＨＲ）またはＲ−Ｒ間隔（ＲＲＩ）において観察される自然に発生するリズムである。このリズムは、呼吸と心血管系との相互作用の直接の結果である。ＲＳＡは、呼吸数と同様のＲＲＩの最大および最小を示す周期的信号であることを特徴とする。通常、ＨＲは、吸気によって増大し、呼気によって低減する。本方法では、ＲＳＡ信号の変化を、生理学的予測因子として用いることができる。
【００２１】
呼吸イベント後に、ＲＲＩが減少することは周知である。これは、覚醒中の交感神経活動の亢進に関連している。覚醒中、交感神経活動の亢進が伴う、心拍数、収縮期血圧および拡張期血圧が上昇する短い期間がある。心拍数の変動もまた、交感神経活動の亢進によって低減する。本方法では、このＲＲＩの低下を生理学的予測因子として用いることができる。
【００２２】
ＥＣＧ波形の形態を用いて、変化する胸郭内インピーダンスと相関する電極に対する心臓の動きによる呼吸による変化に従って、呼吸が正常であるか否かを分類する。呼吸によってもたらされる、電極に対する心臓の回転は、Ｒ波の振幅とＲ波の面積（これは、図３に示すＱＲＳセグメント全体の面積に含まれる）との両方の振動において明らかである。イベント中に呼吸が低減するか、または完全に停止すると、信号のピーク間振動の低減が予測される。本方法では、振幅および面積の変化を生理学的予測因子として用いることができる。
【００２３】
心臓の向きが急にずれることにより、Ｒ波の振幅が急に変化する。ＥＣＧ誘導が直交している場合、ずれは反対方向であるはずである。本方法では、ＥＣＧ誘導のＲ波の振幅の反対の変化を物理的予測因子として用いることができる。
【００２４】
Ｒ波またはＴ波の振幅の漸進的な増大を、物理的予測因子として用いることができる。ＥＣＧ信号における増大するＲ波の一例を図１６に示す。
【００２５】
第４ステップにおいて、各ＥＣＧ群に対し生理学的予測因子の変化を検査することにより、本方法は、呼吸イベントの開始時にＥＣＧ群の各々が発生するか否かを判断する。そして、本方法の第４ステップでは、各呼吸イベントの開始および終了を確定する。図２は、データファイルに格納された連続するＥＣＧ群から無呼吸／低呼吸指数を計算する方法のステップを示している。
【００２６】
第５ステップにおいて、検出された呼吸イベントの数およびそれらの持続時間を用いて、無呼吸／低呼吸指数を計算する。無呼吸／低呼吸指数を、１時間当りの１０秒間より長い呼吸イベントの平均回数として計算する。
【００２７】
本方法は、有利にＥＣＧ振動データのみを用いて呼吸パターンを特徴付ける。ＥＣＧから導出される呼吸努力および抽出されたＥＣＧ信号を検査することにより、呼吸イベントに関連する呼吸パターンの変化を含むデータの領域を特定する。本明細書では、あり得る呼吸イベントの特定に役立つように調査することができる多数の変化について説明する。
【００２８】
図５は、ＥＣＧ信号パラメータ、すなわちＱＲＳ面積およびＲ振幅とともに、閉塞型睡眠時無呼吸イベントにわたって測定された気流２のグラフを示す。イベントの前後に比較して、イベント中、呼吸の大きさが低減することが予測され、それは図５において明白に示されている。あり得る呼吸イベントは、呼吸努力を検査し、呼吸努力が低減する領域の位置を特定することによって特定される。本方法は、あり得る呼吸イベントをマークするのに役立つ呼吸数の変化、すなわちイベント開始時、イベント中および覚醒時に発生する呼吸数の変化の特定を用いる。本方法は、呼吸イベントの開始および終了に発生する、ＱＲＳ面積およびＲ波振幅の急激な変化を用いることができる。心臓の向きが急にずれることにより、Ｒ波の振幅が急に変化する。ＥＣＧから導出される信号の急なずれを用いて、呼吸イベントのあり得る開始および終了をマークする。これを図６に示し、そこでは、呼吸イベント中に発生する急なずれが、Ｒ波振幅信号のこのセグメントにおいて明白である。ＡおよびＢの線は、それぞれ呼吸イベントの開始および終了をマークしている。図７は、流れにおいて明白な各閉塞型無呼吸イベント後に、ＲＲＩ信号の一時的減少を観察することができることを示す。
【００２９】
本方法は、ＥＣＧから導出された呼吸努力の測度から呼吸パターンを確定することを含むことができる。これは、図５に示すように呼吸努力の大きさが低減するパターンか、または図７に示すように呼吸数が急峻に変動するパターンか、または図６に示す呼吸努力マーカの基本となるレベルが急峻に変動するパターンを含んでもよい。本方法は、呼吸努力の複数の測度の間の位相関係が変化するパターン、一続きの呼吸の過程にわたり呼吸努力推定値の大きさが一貫して変動するパターン、または図８に示すように呼吸努力推定値の包絡線の周期的なパターンを使用してもよい。
【００３０】
より詳細には、図２を参照すると、２誘導（またはそれより多い）ＥＣＧ記録を介して、被検者から生ＥＣＧデータを取得する。これは、直交する少なくとも２つのＥＣＧ誘導を使用しなければならない。このデータを、デジタル記憶装置に保存する。ＡＨＩの計算は、格納されているＥＣＧ信号データにアクセスしそれを解析し（図２、ボックス１）、必要な場合は、適切な帯域幅に対してフィルタリングする（図２、ボックス２）ようにプログラムされたマイクロプロセッサによって開始する。取得される大量の信号データと必要な変換は、マイクロプロセッサによって行われ、マイクロプロセッサは、本技術分野において必要であることが既知である処理速度および出力を有するが、また、本明細書に記載されている方法に従って信号およびデータを取得し、格納し、処理するようにプログラムもされている。
【００３１】
格納されたデータの解析は、従来のホルター解析プログラムを使用する検査で開始し、ＥＣＧのＱＲＳ波セグメント、Ｐ波セグメントおよびＴ波セグメントを特定する（図２、ボックス３）。そして、心拍の各セグメントを、正常、心室性、動脈性またはアーチファクトとして分類する。収集されたデータの一例を図４にグラフ化している。図４は、正常に呼吸している睡眠中の被検者から抽出されたＥＣＧ特徴の表現を示す。図４において時間に対して表示されている信号は、本技術分野において既知であるＥＣＧ誘導構成の誘導Ｉから取得された信号である。特定された波の特性を検査することにより、ソフトウェアによってＥＣＧから導出された信号を生成する。生成されたＥＣＧから導出された信号は、Ｒ波振幅、ＱＲＳ面積、Ｒ−Ｒ間隔（ＲＲＩ）、Ｔ波振幅およびＴ波面積である。図４には、最初の２つの信号とともに測定された呼吸気流を示す。ソフトウェアプログラムは、ＥＣＧ信号の形態解析を用いて、Ｐ波、ＱＲＳ群およびＴ波の開始時点、終了時点およびピーク時点を含む、ＥＣＧ信号の形態的特徴を抽出する（図２、ボックス４）。これには、不整脈分類基準を適用することにより、ＱＲＳ群を正常、心室性、上室性またはアーチファクトとして分類することが含まれる。Ｒ−Ｒ間隔、Ｒ波振幅、ＱＲＳ群面積、Ｔ波振幅およびＴ波面積を含む、ＥＣＧ信号の定量的パラメータを生成する（図２、ボックス５〜ボックス９）。
【００３２】
ＥＣＧから導出された信号の多くは、呼吸によって発生する振動パターンに類似する振動パターンを示す（図４）。呼吸努力は、呼吸とともに振動するＱＲＳ面積およびＲ波振幅から推定される（図２、ボックス１０〜ボックス１１）。呼吸信号、すなわちＥＣＧから導出された呼吸信号（ＥＣＧ−Ｄｅｒｉｖｅｄ  Ｒｅｓｐｉｒａｔｉｏｎ）（ＥＤＲ）を、直交するＥＣＧ誘導からのＱＲＳセグメントの面積の比から計算することができる（図２、ボックス２６）。
【００３３】
ＥＣＧから導出された信号の各々において、呼吸の方向の変化（呼気から吸気への変化またはその逆）に関連する局所ピークを特定する（図２、ボックス１０）。これらピークを、ＥＣＧから導出された信号内の呼吸ピーク候補に対し振動パターン、期間の類似性および信号の大きさを検査することにより特定する。
【００３４】
ＥＣＧから導出された信号の各々においてこれらのピークが特定されると、それらを信号間で比較する。これにより、各信号における呼吸ピークの交差検定（クロスバリデーション）が可能になる（図２、ボックス１０）。各信号のピークの各々を検査し比較することにより、呼吸ピークの単一セットを特定する。これを、ピークの各々の間の期間と信号間のピークの振幅の変化とを比較することによって達成する（図２、ボックス１０）。上述した方法を用いて、ＥＣＧから導出された信号の変化と推定された呼吸努力とを検査することにより、あり得るＥＣＧイベントの開始および終了のマーカを特定する。イベントのあり得る開始および終了のこれらマーカを用いて、呼吸イベント候補のセットを生成する。呼吸イベントを、ルックアップテーブルおよび人工ニューラルネットワークの適用を含む、推定された呼吸およびＲ−Ｒ間隔の検出されたパターンに対する予め定義された基準または調整可能な基準と比較する。
【００３５】
本方法は、以下のパラメータを計算する。図９に、Ｒ波振幅の振動の平均のピーク間の大きさが、あり得る呼吸イベントの過程にわたり、イベントの前後のそれぞれの値に対して低減することを示す。図１０に、ＱＲＳ群の下の面積の振動の平均のピーク間の大きさが、あり得る呼吸イベントの過程にわたり、イベントの前後のそれぞれの値に対して低減することを示す。図１１に、あり得る呼吸イベントの後の平均Ｒ−Ｒ間隔が、イベントの過程にわたる平均Ｒ−Ｒ間隔に対して低減することを示す。図１２に、ＱＲＳ群の下の面積の平均値が、あり得る呼吸イベントの過程にわたり、イベントの前後のそれぞれの値に対して変動することを示す。図１３に、複数のＥＣＧ信号の対に対する、Ｒ波振幅の平均値の、あり得る呼吸イベントの過程にわたる、イベントの前後のそれぞれの値に対する変動の相違を示す。図１４に、Ｒ−Ｒ間隔の振動の平均のピーク間の大きさが、あり得る呼吸イベントの過程にわたり、イベントの前後のそれぞれの値に対して低減することを示す。図１５に、２つの直交するＥＣＧ誘導から検出されたＱＲＳ群の下の面積の比の平均のピーク間の大きさが、あり得る呼吸イベントの過程にわたり、イベントの前後のそれぞれの値に対して低減することを示す。図１６に、あり得る呼吸イベントの過程にわたり最小二乗誤差線形近似を当てはめることによる、Ｒ波振幅の線形傾向の推定を示す。
【００３６】
呼吸イベント候補の各々に対し、呼吸イベントのＥＣＧから導出された相関の各々を生成する（図２、ボックス１７、１９、２１〜２７）。これら相関を用いて、呼吸イベント候補の各々を確定するかまたは拒否する。呼吸イベントのＥＣＧから導出された相関が予測値を生成することを示す呼吸イベント候補が許容される。
【００３７】
呼吸ピークが検証されると、それらを用いて、呼吸努力の推定を生成することができる（図２、ボックス２８）。呼吸努力は、Ｒ波振幅およびＱＲＳ群面積の値のシーケンスを、単一または複数のＥＣＧ信号からのそれぞれのチャネルに組み込むことによって推定する。これら組み込まれたチャネルは、呼吸努力のおよそのマーカとしての役割を果たす。呼吸努力の個々のおよそのマーカからの呼吸ピークのシーケンスを、それらマーカにおいて、間隔、信号変動および勾配が類似するピークを検出することによって確定する。個々のマーカから確定されるピークシーケンスを整合することによって、呼吸ピークのシーケンスを微調整して、誤検知の呼吸ピークを除去し、失われた呼吸ピークを回復する。
【００３８】
呼吸によってもたらされる、電極に対する心臓の回転は、Ｒ波の振幅とＲ波の面積（これは、ＱＲＳセグメント全体の面積に含まれる）との両方の振動において明白である。イベント中に呼吸が低減するかまたは完全に停止すると、信号のピーク間振動の低減が予測される。図９および図１０に示すように、イベント中のＲ波振幅／ＱＲＳ面積の平均ピーク間振動の値から、イベント前後のＲ波振幅／ＱＲＳ面積の平均ピーク間振動の値を引くことにより、Ｒｐｐ相関およびＡＲＥＡｐｐ相関を計算する。図９は、呼吸イベント中、Ｒ波の振幅の振動が低減することが分かることを示している。線Ａ、Ｂはそれぞれ、呼吸イベントの開始および終了をマークしている。図１０は、呼吸イベント中、ＱＲＳセグメントの面積の振動が低減することが分かることを示している。線Ａ、Ｂはそれぞれ、呼吸イベントの開始および終了をマークしている。
【００３９】
呼吸イベントの後、ＲＲ−間隔が低減することは周知である。これは、覚醒中の交感神経活動の亢進に関連する。覚醒中、交感神経活動の亢進が伴う、心拍数、収縮期血圧および拡張期血圧が上昇する短い期間がある。心拍数変動もまた、交感神経活動の亢進によって低減する。イベント中のＲＲＩの平均値とイベント後のＲＲＩの平均値との差を用いて、本方法で用いることができるＲＲＩ＿ｐｏｓｔ相関を計算した。図１１は、呼吸イベントの直後のＲ−Ｒ間隔の低減を観察することができることを示す。線Ａ、Ｂはそれぞれ、呼吸イベントの開始および終了をマークしている。
【００４０】
心臓の向きがずれることにより、ＱＲＳ面積およびＲ波振幅が変化する。イベント中の平均ＱＲＳ面積とイベント前後の平均ＱＲＳとの差を用いて、ＡＲＥＡ相関を計算する。図１２は、呼吸イベント中の平均ＱＲＳ面積の変化を示す。線Ａ、Ｂはそれぞれ、呼吸イベントの開始および終了をマークしている。
【００４１】
イベント中の平均Ｒ波振幅とイベント直前および直後の平均Ｒ波振幅との差を、２つのチャネルに対して計算し、その後それらを掛け合わせてＡｎｔｉｐｈａｓｅ相関を計算した。ＥＣＧ誘導が直交している場合、ずれは反対方向であるはずである。
【００４２】
図１３は、Ｒ波振幅の急な変化が明白であることを示し、各誘導の変化は逆方向である。線Ａ、Ｂはそれぞれ、呼吸イベントの開始および終了をマークしている。
【００４３】
呼吸性洞性不整脈（ＲＳＡ）は、ＥＣＧの心拍数（ＨＲ）またはＲ−Ｒ間隔（ＲＲＩ）において観察される自然に発生するリズムである。このリズムは、呼吸と心血管系との相互作用の直接の結果である。ＲＳＡは、呼吸数と同様のＲＲＩの最大および最小を示す周期的信号であることを特徴とする。通常、ＨＲは、吸気によって増大し、呼気によって低減する。示すように、イベント中のＲＲＩの平均ピーク間振動の値と、イベント前後の平均ピーク間振動の値との差を用いて、ＲＳＡ相関を計算した。
【００４４】
図１４は、呼吸イベント中のＲ−Ｒ間隔の振動の低減を観察することができることを示している。線Ａ、Ｂはそれぞれ、呼吸イベントの開始および終了をマークしている。
【００４５】
ＥＣＧから導出された呼吸信号（ＥＤＲ）を、直交するＥＣＧ誘導からのＱＲＳセグメントの面積の比として計算することにより、呼吸信号を確定する。イベント中のＥＤＲ信号の平均ピーク間振動と、イベント前後のＥＤＲ信号の平均ピーク間振動の値との差を用いて、ＥＤＲｐｐ相関を計算する。図１５は、呼吸イベント中のＥＤＲ信号の振動の低減を観察することができることを示している。線Ａ、Ｂはそれぞれ、呼吸イベントの開始および終了をマークしている。
【００４６】
図１６は、イベント中、Ｒ波の振幅が漸進的に増大した、ＥＣＧ信号データ例を示す。イベント中のＲ波振幅に対する、破線で示す線形最良適合の勾配を用いて、Ｒｓｌｏｐｅ相関を計算することができる。
【００４７】
ＱＲＳ面積およびＲ波振幅からの呼吸努力により、イベント分類が可能になる。これら信号を用いて呼吸努力を推定することができるため、それらを用いて、中枢型イベントと閉塞型イベントとを識別するのに役立てることができる（図２、ボックス２９）。呼吸努力は、閉塞型イベント中は、明白でありかつ恐らくは増大するはずであり、中枢型イベント中は、非常に低いままであるかまたは存在しないはずである。この例を図１７ａおよび図１７ｂに提示する。それらの図は、ＯＳＡイベントにおけるＥＣＧパラメータとＣＳＡイベントにおけるＥＣＧパラメータとの相違を示している。図１７ａは、閉塞型イベント中にＱＲＳ面積およびＲ波がいかに変化するかを示している。イベント中、呼吸努力が続いていることは明らかである。図１７ｂは、ＣＳＡイベント中にＱＲＳ面積およびＲ波が変化することを示している。呼吸努力は大幅に低減するかまたは存在しないことが明らかである。また、呼吸イベントのＥＣＧから導出された相関により、検査されているイベントのタイプに応じて異なる結果がもたらされる。たとえば、ＡＲＥＡｐｐ相関は、閉塞型イベントより中枢型イベントの方が、大きさが大きいことが予測される。それは、呼吸と相関するＱＲＳ面積の振動が大幅に低減するはずであるためである。
【００４８】
低呼吸イベントに比較して無呼吸イベント中の呼吸の方が大幅に低減するため、低呼吸イベントより無呼吸イベントの方が、呼吸イベントのＥＣＧから導出された相関のより大きい変化が予測される。呼吸イベントのＥＣＧから導出された相関のこれら変化を定量化することにより、本方法は、無呼吸イベントと低呼吸イベントとを識別することができる。本方法は、呼吸努力のＥＣＧから導出された測度に対して、呼吸イベント内の平均の振動の大きさの、呼吸イベント前後のそれぞれの値に対する比を計算し、これらの比が所定閾値を下回るイベントを無呼吸として分類し、他のイベントを低呼吸として分類することにより、無呼吸イベントと低呼吸イベントとを識別することができる。
【００４９】
呼吸イベントの総数を、ＥＣＧ信号記録期間の時間数で割ることにより、ＡＨＩを計算する（図２、ボックス３１）。睡眠評価時間間隔での無呼吸／低呼吸指数を計算することにより、閉塞型、中枢型ならびに閉塞型および混合型呼吸イベントの組合せに対して計算される別個の指数を用いて、睡眠時無呼吸の存在、深刻さおよびタイプの臨床推定が可能になる。指数は、本方法の使用者が見るために好適なモニタに表示されることが好ましい。
【００５０】
【表１】
2010540124.tif 000003
【００５１】
参考文献
[1] "Sleep-related
breathing disorders in adults: recommendations for syndrome definition and
measurement techniques in clinical research. The Report of an American Academy
of Sleep Medicine Task Force," Sleep, vol. 22, pp. 667- 89, Aug 1 1999.
[2] M. Schrader, C. Zywietz,
V. von Einem, B. Widiger, and G. Joseph, "Detection of sleep apnea in
single channel ECGs from the PhysioNet data base," 2000, pp. 263-266.
[3] C. Maier, H. Dickhaus, M.
Bauch, and T. Penzel, "Comparison of heart rhythm and morphological ECG
features in recognition of sleep apnea from the ECG," in Computers in
Cardiology, 2003, 2003, pp. 311-314.
[4] S. Vijendra, K. Behbehani,
E. A. Lucas, J. R. Burk, D. N. Burli, and D. H. Dao, "The use of R-wave
morphology in the detection of sleep-disordered breathing using the
electrocardiogram - a comparison between leads," 2004, pp. 3881-3884
Vol.6.
[5] F. Roche, V. Pichot, E.
Sforza, I. Court-Fortune, D. Duveｍey, F. Costes, M. Garet, and
J. C. Barthelemy, "Predicting sleep apnoea syndrome from heart period: a
time-frequency wavelet analysis," Eur Respir J, vol. 22, pp. 937-42, Dec
2003.
[6] J. D. Lattimore, D. S.
Celermajer, and I. Wilcox, Obstructive sleep apnea and cardiovascular
disease," J Am Coll Cardiol, vol. 41 , pp. 1429-37, May 7 2003.
[7] J. Wolf, J. Lewicka, and
K. Narkiewicz, "Obstructive sleep apnea: An update on mechanisms and
cardiovascular consequences," Nutr Metab Cardiovasc Dis, Feb 19 2007.
[8] M. T. Naughton,
"Heart failure and obstructive apnoea," Sleep Medicine Reviews, vol.
2, pp. 93-103, 1998.
[9] P. Grossman and E. W.
Taylor, "Toward understanding respiratory sinus arrhythmia: relations to
cardiac vagal tone, evolution and biobehavioral functions," Biol Psychol,
vol. 74, pp. 263-85, Feb 2007. 
[10] J. W. Denver, S. F. Reed,
and S. W. Porges, "Methodological issues in the quantification of
respiratory sinus arrhythmia," Biol Psychol, vol. 74, pp. 286-94, Feb
2007.
[11] F. Yasuma and J. Hayano,
"Respiratory sinus arrhythmia: why does the heartbeat synchronize with
respiratory rhythm?," Chest, vol. 125, pp. 683-90, Feb 2004.
[12] H. J. Burgess, J.
Kleiman, and J. Trinder, "Cardiac activity during sleep onset,"
Psychophysiology, vol. 36, pp. 298-306, May 1999.
[13] M. J. Carrington, R.
Barbieri, I. M. Colrain, K. E. Crowley, Y. Kim, and J. Trinder, "Changes
in cardiovascular function during the sleep onset period in young adults,"
J Appl Physiol, vol. 98, pp. 468-76, Feb 2005.
[14] R. J. Davies, P. J. Belt,
S. J. Roberts, N. J. AIi, and J. R. Stradling, "Arterial blood pressure
responses to graded transient arousal from sleep in normal humans," J Appl
Physiol, vol. 74, pp. 1123-30, Mar 1993.
[15] K. Narkiewicz, N.
Montano, C. Cogliati, P. J. van de Borne, M. E. Dyken, and V. K. Somers,
"Altered cardiovascular variability in obstructive sleep apnea,"
Circulation, vol. 98, pp. 1071-7, Sep 15 1998.
[16]
G. B. Moody, R. G. Mark, A. Zoccola, and S. Mantero, "Derivation of
Respiratory Signals from Multi-lead ECGs," Computers in Cardiology, vol.
12, pp. 113-116, 1985.
【図１】
2010540124.tif 000004
【図２】
2010540124.tif 000005
【図３】
2010540124.tif 000006
【図４】
2010540124.tif 000007
【図５】
2010540124.tif 000008
【図６】
2010540124.tif 000009
【図７】
2010540124.tif 000010
【図８】
2010540124.tif 000011
【図９】
2010540124.tif 000012
【図１０】
2010540124.tif 000013
【図１１】
2010540124.tif 000014
【図１２】
2010540124.tif 000015
【図１３】
2010540124.tif 000016
【図１４】
2010540124.tif 000017
【図１５】
2010540124.tif 000018
【図１６】
2010540124.tif 000019
【図１７ａ】
2010540124.tif 000020
【図１７ｂ】
2010540124.tif 000021
【国際調査報告】
2010540124.tif 000022
2010540124.tif 000023
2010540124.tif 000024
2010540124.tif 000025
2010540124.tif 000026