JP4832289B2 - 呼吸事象検出システムおよび方法 - Google Patents

Description

本発明は、生理的パラメータの移動性記録、特に心肺パラメータの誘導性体積変動記録からのデータ処理を改良するシステムおよび方法を提供し、更に、この改良した処理によって可能となる、咳や嘆息のような間欠的な生理的事象の検出システムおよび方法も提供する。

２００３年４月１０日に出願された、Sackner InmanおよびWilhelmによる"System and Methods for Respiratory Event Detection"（呼吸事象検出システムおよび方法）と題する米国仮特許出願第６０／４６１，７３８号、および２００３年９月２６日に出願された、Sackner Coyle、Keenan Inman、およびWilhelmによる"Systems and Methods for Respiratory Event Detection"（呼吸事象検出システムおよび方法）と題する米国仮特許出願第６０／５０６，９０４４号の優先権を主張する。これら２件の特許出願の内容は、ここで引用したことにより、全ての目的に対してその全体が本願にも含まれるものとする。

生理的パラメータの連続的移動性監視は、日常生活の臨床関連徴候の基準や、いかにその経験が被験者の併発活動および挙動によって形作られるかということについて、我々の理解を広げることができる。有望な新技術が開発されているが、移動性監視データは解釈が難しいことが多い。主な問題の１つは、肉体的活動が、研究所または診療所にいるときのように制御されないことであり、気が付かないと、通常の運動によって誘発される肉体的変化と病気の指示との混乱を招く可能性がある。したがって、臨床関連移動性監視システムは、運動性活動を登録し、評価できる内容を提供することによって、医師がいずれの異常値（例えば、ＥＣＧまたは呼吸パターン）でも肉体的運動に帰することができるのか、または生理的不規則性(dysregulation)を表すのか判断できる有益性を備えてしかるべきである。日常的に、監視対象被験者の活動、姿勢、および居所を紙と鉛筆で、または電子的に記録すれば、より多くの関連情報を提供することができ、踏み込んだ医学的解釈に役立てることができる。しかしながら、日記だけを用いるのでは、居所や活動レベルの変化を報告する際に記録が不正確なことが多いため、不十分である。

発話(speech)は、特に呼吸を移動性監視と混乱する可能性のある重要な活動である。発話状態の出現は、最も頻繁に見られる肉体的活動および挙動の一種であり、肺機能に加えて、種々の生理的系統を変化させる可能性がある。例えば、脈拍は、会話的発話によって、毎分５回から１０回増加するのが通例である。脈拍の増加は、社会的に要請がある状況では一段と高くなる可能性がある。マイクロフォンによる聴覚的記録が、発話活動を定量化するために用いられてきたが、監視対象の人の声以外にも、周囲の音や他の人の発話というような音も拾い上げるので、限界がある。咽喉マイクロフォン(throat microphone)の方が選択度が高いが、長時間装着することは不便であり、監視対象被験者が周囲の注目を引くので望ましくない。

誘導性体積変動記録法（ＩＰ：inductive plethysmography）は、心肺機能の非影響呼吸監視のための科学的および医療的に受け入れられている業界標準(gold standard)であり、医療および研究用設定において広く用いられている。呼吸については、この技法は、胸郭および腹区画双方の拡大および縮小を、正弦波状に配した電気配線から成りこれを可撓性帯に埋め込んだＩＰセンサを用いて、測定することにより、呼吸器系によって移動させられる空気量を近似する。高周波低電圧の発振電流が配線を通過して、コイルの自己インダクタンスを測定するために必要な磁界を発生する。これは、帯が包囲する断面積に比例する。胸郭および腹帯(abdominal band)の較正後、２つの信号の加重和は、一回換気量に対応または比例する。

したがって、当技術分野では、肉体的活動、特に発話を登録または検出し、活動情報を利用して、移動性監視データの解釈を改良し一層の信頼性向上を図るためにシステムおよび方法の改良が必要となっている。このようなシステムは、誘導性体積変動記録法による移動性監視に直接適用可能であるとよい。

この章または本願のいずれの章においても、いずれの参考文献の引用または特定も、このような参考文献が本発明の従来技術として使用可能であるとは解釈しないものとする。

本発明の目的は、呼吸パラメータ、例えば、肺容積（Vt)を求めることができる肺モデル、好ましくは二区画肺モデルを特徴付けるのに十分なパラメータの移動性記録からのデータ処理を改良するシステムおよび方法を提供することによって、従来技術における欠陥を克服することである。つまり、好ましい入力パラメータは、胸郭（ＲＣ）および腹部（ＡＢ）の断面積、円周、または半径（あるいは同様の幾何学的変数）の指示を得るというようなことにより、ＲＣ体積およびＡＢ体積を特徴づける。このような好ましいパラメータは、外皮内に混み込むというようなことにより、検査対象被験者の胴体に対して固定したセンサを用いた、電気的、光学的、またはその他の技術に基づく呼吸体積変動記録法によって決定することが好ましい。

呼吸誘発体積変動法（ＲＩＰ）は好ましい測定技術であるが、本発明のシステムおよび方法は、他のセンサ技術への適用も容易である。このようなセンサ技術には、例えば、身体インピーダンス・センサ、水銀含有シラスティック(silastic)歪みゲージ、ベローズ呼吸運動描写器、体積呼吸運動描写器、差動線形変換器、身体周囲の誘導性変換器、身体の直径を検知する磁気記録計、局部的な移動を測定する圧電変換器、光の反射による移動分析等が含まれる。更に、本発明および方法は、身体周囲、直径、距離、および移動の複数のパラメータを反映する複数の信号を発生するセンサ技術への適合化も容易であり、これらを合わせると、少なくとも、二区画呼吸モデル、および／または胸郭および腹部の形状変化のような、その改善の証印(indicia)が得られる。加えて、これらの方法は、熱電対、サーミスタ、終末呼吸二酸化酸素センサ、鼻圧力(nasal pressure)および流套管(flow cannulas)呼吸マスク、空気流差圧センサ等を含む、種々の空気流モニタが発生する信号のような、空気流を反映する信号にも適用することができる。

本発明の別の目的は、好ましくは、肺容積の時間経過、連続呼吸の一回換気量等のような、呼吸を反映する生理的データの処理を改良し、呼吸停止、減呼吸、咳、嘆息等のような、生理的事象を検出するシステムおよび方法を提供することによって、従来技術における欠陥を克服することである。これら別のシステムおよび方法に対して入力するデータは、多くのセンサおよび多くのセンサ技術によって求めることができる。好ましくは、この入力データは、二区画肺モデル、即ち、呼吸モデルを特徴付けるのに十分な移動性記録から求める。

ここでは、主に移動性記録に応用した場合について本発明を説明するが、そのシステムおよび方法の一部または全ては、研究所、診療所、または病院におけるような、その他の設定にも適用可能であることは当然認められよう。

ある種の実施形態は、添付した特許請求の範囲に要約されている。
特許および印刷刊行物を含む種々の引例が本願全体を通じて引用されている。このような引用された引例の内容は全て、ここで引用したことにより、あらゆる目的に対して、その全体が本願にも含まれるものとする。

本発明は、以下の本発明の好適な実施形態の説明、本発明の具体的な実施形態の図示例、および添付図面を参照することにより、一層深く理解することができよう。

以下に、好ましい呼吸誘発体積変動記録（ＲＩＰ）技術から得られる入力データを主に参照しながら本発明について説明するが、記載するシステムおよび方法は、他のセンサ技術から得られるデータとでも使用できることは理解されよう。例えば、どのようにして求めたかには係わらず、二区画呼吸モデルを特徴付けるデータに、暫定的信号処理を適用することができる。また、不可欠な入力データをどのようにして得ようとも、肺容積の指標の時間過程（連続または離散）からの呼吸事象の決定は有用である。

更に、本発明には様々な特徴や形態があり、以下では主に単一のシステムおよび方法の構成要素として説明する。また、本発明は、そのような単一の実施形態には限定されるのではなく、これらの特徴および形態のその他の組み合わせや部分的組み合わせを有するその他の実施形態も含み、これらには独自の有用性があることも理解されよう。また、本発明のシステムおよび方法の記載および図は、ある順序になっているが、この順序は唯一でも限定でもないことは理解されよう。ステップ、特徴、および形態は、記載するデータの従属性が一致するのであれば、他の時間的順序でも実施することができる。
５．１ 好ましい体積変動記録システム
端的な背景という目的で、呼吸体積変動記録法（ＲＰ）および呼吸誘発体積変動記録法（ＲＩＰ）について次に説明し、続いて、このセンサ技術の好ましい実施形態について端的に説明する。「体積変動記録法」は、ここでは、例えば、肺、血管中の血液等のような、含まれる流体の体積変化による器官の容積判定を意味するものとする。
Ａ．ＲＩＰ技術の概要
ＲＰは、呼吸機能、例えば、被験者の胸郭（ＲＣ）を通る１つ以上の断面または被験者の腹部（ＡＢ）を通る１つ以上の断面の面積（直径、半径等）を特徴付けるために十分な、時々刻々の容積を反映する指標を判定する。このような判定から、呼吸の二区画モデルを参照して、肺の容積、心臓の容積等を判定することができる。

ＲＩＰは、このような指標を決定する際に、被験者の身体の対象の断面周囲に配した配線ループの自己インダクタンスを測定する。自己インダクタンスは、ループによって囲まれる断面積に大きく依存する。好ましくは、ループ自己インダクタンスを測定する際に、ループを発振器に含ませ、発振周波数を測定する。例えば、２００１年４月１７日出願の米国特許出願第０９／８３６，３８４号（移動性誘導性体積変動記録システムの改良）、２０００年４月４日発行の米国特許第６，０４７，２０３号（センサ外皮を含む移動性誘導性体積変動記録システム）、２００２年１月２９日発行の米国特許第６，３４１，５０４号（誘導性体積変動記録センサのための伸張可能な導電性繊維）、１９８９年２月２８日発行の米国特許第４，８０７，６４０号（伸張可能な誘導性体積変動記録変換器）、１９９４年７月２６日発行の米国特許第５，３３１，９６８号（誘導性体積変動記録センサおよび回路）、および１９９４年４月１２日発行の米国特許第５，３０１，６７８号（伸張可能な誘導性体積変動記録変換器）を参照のこと。

具体的には、呼吸の測定では、ＲＣおよびＡＢ ＲＩＰループから得た信号をフィルタ処理し、平滑化し、較正して、時々刻々の肺容積を反映する指標を求めると有益である。例えば、１９９２年１１月３日発行の米国特許第５，１５９，９３５号（個々の肺機能の測定）、１９８９年３月２８日発行の米国特許第４，８１５，４７３号（呼吸容積の監視方法）、１９８３年１月５日発行の米国特許第４，３０８，８７２号（呼吸容積監視方法）、２００２年７月２日発行の米国特許第６，４１３，２２５号（誘導性体積変動記録呼吸モニタの較正方法）、１９８９年５月３０日発行の米国特許第４，８３４，１０９号（誘導性体積変動記録呼吸モニタの較正方法）、および１９８３年２月１５日発行の米国特許第４，３７３，５３４号（誘導性体積変動記録呼吸モニタの較正方法）を参照のこと。

加えて、誘導性体積変動記録法を用いると、他の生理的指標を判定することもできる。例えば、時々刻々の心臓容積、心臓出力、および心室壁の運動を反映するデータを、中胸部センサ・ループから抽出することができ、センサ・ループから、先端部または首に関するデータを抽出することができ、動脈および静脈拍を反映するデータを抽出することができる。例えば、２００２年３月２６日出願の米国特許出願第１０／１０７，０７８号（心室容積信号の抽出のための信号処理技術）、１９９３年１月１２日発行の米国特許第５，１７８，１５１号（心臓出力の誘導性体積変動記録測定方法）、１９９１年８月２０日発行の米国特許第５，０４０，５４０号（中央静脈圧の誘導性体積変動記録測定）、１９９１年１月２２日発行の米国特許第４，９８６，２７７号（中央静脈圧の誘導性体積変動記録測定）、１９８４年６月２６日発行の米国特許第４，４５６，０１５号（首部容積変化の測定）、および１９８４年６月５日発行の米国特許第４，４５２，２５２号（首および口の容積測定による心臓パラメータの判定）を参照のこと。

Ｂ．好ましいＲＩＰシステム
図１は、本発明のシステムおよび方法によって分析する信号を測定するのに適した、好ましい（しかし限定的ではない）ＲＩＰシステムを示す。この好ましいシステムは、外皮１、データ記録部３、および分析用ソフトウェア・パッケージと共に構成され、オフラインおよび対話式信号分析、表示、および報告の作成を行う分析用コンピュータ４を含む。図示した好ましい外皮（ここでは、「シャツ」とも呼ぶ）は、選択可能で再利用可能な伸縮性材料で作られ、ぴったり装着し被験者の身体と共に伸縮することができ、１つ以上の埋め込みＲＩＰセンサ帯を含む。センサは、導電性配線ループで構成され、身体の断面周囲に配置され、例えば、正弦波状配列を有することにより、伸縮自在に構成されている。ＲＩＰセンサのインダクタンスは、電気および電子技術分やでは周知の方法を用いて、電気信号、通例では周波数信号に変換される。

センサをシャツの中に埋め込むのは、被験者の胴体に対するそれらの正確で持続性のある配置を確保するためである。呼吸を測定するためには、２つのセンサが好ましい。胸郭の高さにあるセンサ５、および腹部の高さにあるセンサ６である。心臓パラメータを測定するためには、剣状突起の高さの真下に位置する中央胸ＲＩＰセンサ（図示せず）が好ましい。シャツは、他の種類のセンサを含んでもよく、即ち、取り付けてもよく、電磁センサ７、身体の姿勢および運動を検知するための１つ以上の加速度センサ（シャツ内にあるように図示されている一例の加速度計８、他は図示されていない）等、マイクロフォン、パルス酸素計（図示せず）、カプノグラフ、１つ以上のＥＥＧ電極等を含む。病院、診療所、または研究所では、広範囲の生理的センサから他の信号を得ることもできる。図示したＲＩＰ技術は現時点における好適な実施形態であり、今後の外皮構成およびセンサのパッケージングによっては、例えば、追加のＲＩＰセンサおよび他の技術の体積変動記録センサを内蔵することもあり得る（例えば、光センサを採用したり、被験者との相互作用の後における電磁または音響エネルギにおけるドプラ・シフトを検知する）。そして、以下に記載する方法は、追加のセンサを含むこのようなデータにも適用可能である。

図示のシャツに装着したセンサは、ケーブル２によってデータ記録部に取り付けてもよく、他の周辺診断デバイス（図示せず）を補助ケーブル１０によってデータ記録部に取り付けてもよい。あるいは、ローカル・ワイヤレス・リンク、光リンク、あるいはシャツまたはその繊維内に埋め込んだリンクを、離散ケーブルの代わりに用いてもよい。データ記憶ユニット３は、帯上に装着したり、ポケットに入れたり、またはシャツ内に埋め込むことができるように、小型で軽量であることが好ましい。このユニットは、完全公開格納(full disclosure storage)およびオフライン分析のために、十分な精度および確度でセンサ波形を格納すること、および／または内部である種の処理および分析を行うこともできる。データは、着脱自在の標準化されたメモリ・カード１１（例えば、フラッシュ・メモリ・カード）によって分析コンピュータ４に転送されるか、あるいはセル電話技術を用いた伝送のようなワイヤレス・リンクによって転送される。データ記録ユニットは、ディジタル日記帳を実現するために、タッチ・スクリーンを含んでもよく、そのデータも分析データに転送し、センサ読み取り値と相関付けるとよい。

本発明のシステムおよび方法は、分析コンピュータ４（および、任意に一部データ・ユニット３において）実行する分析ソフトウェアによって実現する。このソフトウェアは、メモリ・カードからデータを読み取り、あるいはそれ以外の方法でセンサ・データを受け取り、本発明の方法にしたがって、毎分換気量、一回換気量、呼吸速度(respiratory rate)、呼吸流量(respiratory flow rate)等のような、被験者の呼吸パターンを特徴付ける種々の生理的パラメータを抽出し表示する。分析ソフトウェア・パッケージは、分析コンピュータのユーザが対話的に制御する利点がある。また、本発明は、本発明の方法を実現する分析ソフトウェアの選択部分全てを、プログラム製品として入手可能とし、これをＣＤ−ＲＯＭのようなコンピュータ読み取り可能媒体上で分析コンピュータに転送することも想定している。

加えて、分析ソフトウェアは、ＲＩＰ信号と共に、またはこれらと組み合わせて、他のセンサからの信号も処理するという利点がある。例えば、パルス酸素計の信号は、同時発生の動脈酸素飽和情報を提供することができ、加速度計信号（および同様のセンサからの信号）を処理して、被験者の姿勢および活動の指標を提供することができる。姿勢および活動は、同時に測定する心肺測定値について挙動内容を定めるので重要である。

５．２ 呼吸信号処理
図２を参照して、呼吸信号処理方法の好適な実施形態についてこれより説明する。簡略化のため、そして限定ではなく、この記載を組織化するために見出しを用いる。

５．２．１ 入力信号
必須の入力信号は、非侵襲的に、好ましくは、二要素呼吸モデルにしたがって、呼吸容積を決定するためのパラメータを定める。例えば、Stern et al. Phychophysiology Recording（心理生理学記録）、 Second Edition, Oxford University Press, Oxford, U.K. の第１０章（呼吸器系）を参照のこと。端的に言えば、二要素モデルは、胸郭（大きさが可変の円筒と同様）の直径、および横隔膜（円筒内のピストンと同様）の位置から、肺の容積を概念化する。横隔膜の位置は、腹部の容積によって反映され、その比較的圧縮できない内蔵が横隔膜によって圧縮される。したがって、好適な実施形態では、１対の呼吸信号２１がＲＣおよびＡＢの大きさを反映し、前述のＲＩＰ技術センサからというように、呼吸体積変動記録センサ（ＲＰ）から入力する。あるいは、入力信号は、出力Vtおよび／またはVol信号（後に説明する）、または同等物を出力する別の呼吸モデルに対するデータを定めることもできる。このような出力信号は、二要素モデルから発生したか否かには係わらず、呼吸事象認識方法、即ち、本発明によって処理することができる。

他の実施形態では、複数の入力信号が、胴体の寸法、動き、または形状を反映することができ、ＩＰまたはＲＩＰ技術以外によって発生することもできる。例えば、ＲＩＰと同様に、光または光ファイバ系ＲＰも使用することができる。また、胴体形状の数個のパラメータを反映するセンサも、本発明に入力することができる。このような形状反映信号を処理して、ＲＣおよびＡＢ信号を発生し、後に、例えば、周知の再帰法によって決定される係数との線形結合によって用い、再帰法によりこのような線形結合と同時ＲＰまたはＲＩＰ信号との間の差を最少にすることができる。あるいは、このような形状反映信号は、例えば、ＱＤＣ法の拡張版によって処理して直接Vt信号（以下で説明する）を発生したり、例えば、同時実空気流情報に適用する線形再帰法によってVol信号（以下で説明する）を発生することもできる。具体的には、このような形状反映信号が被験者の胴体の表面上で規定される距離の適当な複数の指標を含む場合、このような信号は、数学の技術分野では既知の方法（具体的には、微分幾何法）によって処理し、胴体の形状の表現を得ることができる。二次元表面（メトリック）上に十分な距離を定義すれば、その表面の形状を三次元空間内に構築できることは公知である。比較的小さな形状変化のみを伴う呼吸監視では、この構造は、静止している胴体の形状を中心に線形化すればよい。したがって、この幾何学的構造は、線形統計モデル構築方法によって、距離反映信号を較正空気体積に線形的に当てはめることによって、回避することができれば好ましい。

１つ以上の加速度計、パルス酸素計、およびマイクロフォンからの追加の信号は、呼吸信号と共に処理すると有利である。被験者上に位置し、例えば、胴体上、または脚部上、または足部上（例えば、靴の中）に装着した１つ以上の加速度計（図１における８）からの信号は、被験者の姿勢および運動即ち活動の指標を与える。このような情報は、同時生理的データの挙動内容を定め、（以下で説明するように）ＲＰ信号からのモーション・アーチファクト(motion artifact)をフィルタ処理または適応的フィルタ処理するために用いることができる。周囲の音響からの妨害を最少に抑えて被験者の発声を捕獲するために装着する頚部マイクロフォン（図１における９）のようなマイクロフォンを用いると、発話の有無を判定することができ、咳、嘆息、およびその他の呼吸事象を判別する際に役立てることができる。その他の任意のセンサも、呼吸の解釈または処理に関連する追加情報を提供することができる。

５．２．２ 入力信号のディジタル化
信号処理における最初のステップは、信号のディジタル化２５およびフィルタ処理２７である。以下の説明は、好ましい一連のステップに従っており、２５においてアナログ・ドメイン信号を最初にディジタル化し、次いで２７において公知のディジタル・フィルタリング技法にしたがってディジタル化した信号をフィルタ処理する。余り好ましくはないが、信号を全体的にアナログ・ドメインでフィルタ処理し、次いでディジタル化してもよい。しかしながら、好ましい一連のステップにおいてでさえ、ディジタル化するアナログ信号がサンプリング・レートの半分よりも大きい周波数の成分を有さない（ナイキスト条件）ように、予備的な明示または暗示アナログ・フィルタ処理ステップを実行することは理解されよう。

好ましくは、生のアナログ信号をディジタル化するには、シャツに装着した電子回路、またはディジタル形態での記録および転送のためのデータ記憶ユニットにおける電子回路のいずれかを用いる。信号は、少なくとも８ビットで、更に好ましくは１６ビットで量子化（Ａ／Ｄ変換）する。必要であれば、呼吸信号のみを１６ビットで量子化し、他の信号（加速度計、マイクロフォン等）を８ビットで量子化してもよい。好ましいサンプリング・レートは、信号の種類によって異なる。

ＲＩＰセンサから発生したアナログ信号は、通常、約５０Ｈｚサンプリング・レート（２０ミリ秒のサンプル時間）での処理に合わせてディジタル化するので、実際の信号は、２５Ｈｚまで、および４０ミリ秒期間まで表される。更に周波数が高い呼吸現象の処理のため、そして更に精度の高いディジタル化信号の時間微分計算のためには、４００Ｈｚまでのサンプリング・レート（２．５ミリ秒のサンプリング時間）を用いることもできる。あるいは、３０〜４０Ｈｚという低いサンプリング・レートを用いると、携帯型データ記憶ユニットにおいて記憶の必要性を制限することができる。これらのサンプリング・レートは、例えば、２００１年４月１７日に出願した同時係属中の特許出願第０９／８３６，３８４号に記載されている。

加速度計信号は、好ましくは、約１０Ｈｚ（１００ミリ秒のサンプル時間）（８〜１６Ｈｚの範囲）でディジタル化することが好ましく、加速度「Ｇ」の単位で定量化することができる。音質（強度）の判定のために、マイクロフォン信号を、好ましくは、約１０Ｈｚでディジタル化し、デシベル単位で定量化するが、発話、いびき、咳等を区別するのに十分な音響特性を得るためには、１ｋＨｚから１０ｋＨｚまでのサンプリング・レートが好ましい。

５．２．３ 入力信号のフィルタ処理
信号のフィルタ処理２７は、信号を測定しているときに、そしてそれらを移転可能なメモリ・カードに記録する前に、好ましくはディジタルで行うとよい。あるいは、ディジタル化したがフィルタ処理しない信号をメモリ・カードに記録しておき、後に分析コンピュータ上においてソフトウェアでフィルタ処理してもよい。

加速度計信号は、実施形態の中には、呼吸信号の適応フィルタ処理を行うために用いると有益な場合もあり、ディジタル化し、処理して、被験者の姿勢、および／または活動、および／または運動の指示を与える。好ましい処理では、入力信号をフィルタ処理して、姿勢および／または位置に対して敏感なローパス成分と、活動および／または運動に対して敏感なハイパス成分とを得る。ハイパス／ローパスの帯域分離は、姿勢変化（恐らく、被験者毎に特性がはっきり分かれる）に特有の周波数、例えば、０．１Ｈｚと１Ｈｚ（好ましくは０．４−０．６Ｈｚ）との間で選択することが好ましい。次いで、姿勢の指標を、ローパス・フィルタ処理した加速度計信号から判定することがでできる。例えば、胴体において加速度計が観察する重力加速度の方向は、着座／起立位置と仰臥位置との間で区別することができ、一方腿に装着した加速度計は、同様に、着座／仰臥位置と起立位置との間で区別することができる。したがって、外皮に装着した加速度計および腿に装着した加速度計からの信号を組み合わせると、仰臥、着座、または起立位置の指標を得ることができる。また、活動または運動の指標は、ハイパス・フィルタ処理した信号から判定することができる。例えば、活動強度は、ハイパス信号における全エネルギによって反映される。活動の速度または種類は、ハイパス信号内にあるモード周波数(modal frequency)において反映される。更に、足または靴に装着した加速度計は、歩行または走行中に発生する力の指標を与えることができる。

ＲＰまたはＲＩＰによって測定するＲＣおよびＡＢ信号は、勿論、呼吸を反映し、心臓活動、アーチファクトによって誘発した活動または運動、およびノイズを反映する成分も含有する場合もある。図３Ａは、ＲＣおよびＡＢが発生したことがわかっているスペクトルを示す（一定比率ではない）。最初に、呼吸信号５０は、全体的に中程度の振幅を有し、基本周波数が約０．２〜０．５Ｈｚ（１４〜３０呼吸／分）で、高調波が約１．０〜１．５Ｈｚまで達している。信号５２に含まれる心臓信号は、全体的に振幅は呼吸振幅の１０〜２０％以下であり、基本周波数が１〜２Ｈｚ（６０〜１２０回／分（しかし、場合によっては０．７４Ｈｚ、即ち、４５回／分と、低くなることもある）、高調波は３Ｈｚ以上である。特に移動性監視では、モーション・アーチファクトが大抵の場合見られる。これらのアーチファクトは、信号５２にも含まれており、振幅および周波数スペクトルが広範囲にわたって変動するが、通常の運動の多くでは、モード周波数が約１〜２Ｈｚであり、これよりも高い高調波がある。最後に、その他の発生源からのノイズは、センサに固有のノイズ、および信号処理によって混入したノイズを含むが、一般に振幅は小さく、対象範囲ではスペクトルは平坦である。

その結果、好適な実施形態では、ＲＰまたはＲＩＰ信号をフィルタ処理して、呼吸信号およびその高調波を今後の処理のために保持しつつ、心臓、運動、およびノイズ成分を抑える。一実施形態では、呼吸成分を非呼吸成分全般から最良に分離するように、単一のローパス・フィルタを選択し、ＲＣおよびＡＢ信号に適用する。あるいは、一団の被験者がいる場合には、個々の被験者毎に個別化する。図３Ａは、理想的なこのようなローパス線形フィルタ５４を、全体的な周波数応答と共に模式的に示す。通過帯域は０．０Ｈｚから約１．３Ｈｚであり、減衰が増加するロール・オフ帯域は約１．３Ｈｚから約１．５Ｈｚであり、拒絶帯域は約１．５Ｈｚよりも上である。本発明のこのディジタル・フィルタおよびその他のディジタル・フィルタは、当技術分野では、重み係数選択型有限インパルス応答（ＦＩＲ）フィルタとして知られている。例えば、Smith, 2nd ed. 1999, The Scientist and Engineer's Guide to Digital Signal Processing（科学者および技術者のためのディジタル信号処理の指針）, California Technical Publishing, San Diego, CAを参照のこと。このようなＦＩＲフィルタは、１００〜２５０個の係数を有し、必要に応じて１０００個以上の係数を有することもできる。

別の実施形態では、ＲＰおよびＲＩＰ信号のフィルタが非呼吸信号の測定可能な特性に調節するか、またはこれに適合すると、特に移動性内容の場合に、呼吸信号の非呼吸信号からの分離の改善に至ることがわかっている。このような別の「オープン・ループ」型実施形態の１つでは、異なる通過帯域を有する複数のフィルタの１つ（または通過帯域が調節可能な単一のフィルタ）を、加速度計信号から決定される運動または活動特性に応じて、場合毎に選択する。運動および活動が活発になるに連れて、モーション・アーチファクトは一般に低い方の周波数まで延びていくことが分かっているので、このような実施形態では、運動または活動強度が増加（減少）するときには、通過帯域が短い（長い）フィルタを選択することができる。運動活動は、例えば、フィルタ処理した最新の加速度計信号の根二乗平均として判定することができる。

図３Ｂは、このような実施形態の１つを模式的に示す。ここでは、フィルタ特性５４（図３Ａにおける特性５４と同様）が選択されるのは、被験者の運動が殆どまたは全くないときであり、約１．１〜１．２Ｈｚから１．３〜１．４Ｈｚまでにロール・オフ帯域を有するフィルタ特性５８が選択されるのは、被験者に平均的な運動がある場合であり、被験者に大きな運動がある場合には、約０．９Ｈｚから１．１Ｈｚまでにロール・オフ帯域を有するフィルタ特性６０を選択する。ある種の被験者では、約０．１〜０．２ｇ（重力加速度）未満の過渡被験者加速度が、被験者に殆どまたは全く運動がないことを示し、約０．１〜０．２ｇから０．３〜０．５ｇの過渡加速度が平均的な運動を示し、約０．３〜０．５ｇを超える過渡加速度が大きな運動を示す。この実施形態は、３組以上のフィルタ重み係数を予め計算し、１つ以上の加速度計またはその他の運動センサからの過渡入力の（ＲＭＳ）平均振幅にしたがって係数集合から選択することによって実現することができる。

更に適応的な実施形態も、適応信号処理または信号推定の既知の原理にしたがって構築することができる。非呼吸信号成分を抑えるために適応的に構築した実施形態の１つでは、フィルタ処理した加速度計データと、モーション・アーチファクトによる信号成分との間の相関、およびＥＣＧまたはパルス酸素計信号からの心拍データと、心臓成分(cardiac component)による信号成分との間の相関を利用する。更に、加速度計の信号は、ＲＰまたはＲＩＰの呼吸成分とは無関係である。同様に、心拍データも呼吸成分とは総じて無関係であるが、呼吸によって誘発される脈拍の変動(variability)および振幅による何らかの相関がある場合もある。モーション・アーチファクトは、時々刻々の適応処理、または加速度計信号をノイズ信号基準として用いてＲＰまたはＲＩＰ（あるいはその他の）信号をフィルタ処理する（このようなフィルタ処理は一般にディジタルである）ことによって抑えることができ、一方心臓信号は、一般に、心拍データをノイズ信号基準として用い、更に時々刻々の適応処理またはフィルタ処理を行えば抑えることができる。適応干渉抑制手順が、当技術分野では公知であり、信号および多数のノイズ源双方に適用可能である。適応処理は、好ましくは、ＬＭＳアルゴリズムに基づく。例えば、Widrow eta al. 1985, Adaptive Signal Processing（適応信号処理）, Pearson Education, Singapore、特に第１２章（適応的干渉相殺）を参照のこと。

個々の被験者毎にフィルタおよび適応パラメータを個別化することが好ましい場合が多い。例えば、競技者についてのこのようなパラメータは、心肺疾患を患っている者とは異なると思われる。個々のパラメータは、較正期間において被験者に種々の選択した強さの運動または活動を行うように命令し、この期間における測定から決定すればよい。

５．２．４ 呼吸信号の較正
ディジタル化したＲＰ信号は、非呼吸成分を抑えるためにフィルタ処理されていることが好ましく、次にこれらを較正する（２９）。この処理は、これらの信号、好ましくは、ＲＣおよびＡＢ ＲＩＰセンサの処理済み出力を組み合わせ、実際の時々刻々の肺容積（ここではＶｔで示す）に比例する別の信号を生成する。２つよりも多いＲＰセンサからのデータが使用可能な場合、先のディジタル化およびフィルタ処理を全てのセンサ信号に対して行うことが好ましく、全ての信号に対する線形当てはめにより、肺容積に比例する単一の信号を決定するために、容易に分かる方法でこの較正ステップを拡張する。このステップおよび後続のステップは、センサ信号を受け取る毎に実行することができ、あるいは信号を用いて後から実行してもよい。

二区画呼吸モデルによれば、以下の線形組み合わせにより、時々刻々の肺容積に比例するＶｔ信号が得られる。

（あるいは、ＲＣおよびＡＢの全体的に線形の組み合わせを用いてもよい）。実際の呼吸空気体積（または流量）の測定値が使用できるのであれば、Ｖｔを更に較正して、実際の肺容積（ここではVolで示す）を密接に反映する信号とすることができる。Volは、Ｖｔ信号の線形スケーリングによって得ることができる。

Ｍはこれら実際の空気体積の同時適正較正Ｖｔ信号に対する比率への最良の当てはめである。実際の空気体積の測定は、移動性内容には好ましいバッグ(rebreathing bag)、肺活量計、統合呼吸気流計等を用いて測定することによって行うことができる。被験者には、ある回数（例えば、３〜１０回）の呼吸の間、最大限息を吸入し、そして最大限吐き出し、バッグ（例えば、８００ｍｌの）を最大限へこませ、次いで最大限膨らませるように命令する。次いで、例えば、最少二乗当てはめによって、既知のVolおよび較正したＶｔ信号からＭを求めることができる。

図５Ａは、静止時における安定した呼吸の間に監視した実際の被験者から得たＲＣおよびＡＢ信号の一例を示す。以下で説明する方法で決定したＫを用いて、これらの信号を組み合わせて、図示のＶｔ信号を得る。更に、Ｖｔ信号を調整(scale)してVol信号（図示せず）にすることもできる。後続の処理の多くは、このＶｔ信号（および／またはVol信号）と、これから抽出した呼吸パラメータとに基づくことができる。

Ａ．基本的な較正
Ｋは、好ましくは、改良した定量的診断較正（ＱＤＣ：Quantitative Diagnostic Calibration）手順によって決定する。まず、基本的ＱＤＣについて説明し、次いで、改良した手順の説明が続く。例えば、米国特許第６，４１３，２５５号および第４，８３４，１０９号、ならびにSchaf (ed.), 1989, Heart Lung Interactions in Health and Disease（健康体および病患体における心肺相互作用）, (ISBN 082477986X), pp676-678を参照のこと（ここで引用したことにより、あらゆる目的に対して全てが本願にも含まれるものとする）。

基本ＱＤＣは、後続の監視において目立つ挙動状態を表す挙動状態において正常に呼吸する被験者によって測定したＲＣおよびＡＢ信号から、Ｋを決定または較正する。例えば、睡眠の研究では休止仰臥(rest supine)、移動の研究では安定した歩行中の較正である。大抵の場合、被験者は適正な呼吸および活動について、ならびに較正手順について命令を受け、ＲＣおよびＡＢ信号自体によって、または同時加速度計信号によって、遵守を確認することができる。較正測定は、好ましくは、少なくとも１００回の呼吸（約５分）を含み、更に好ましくは２５０回の呼吸（約１０〜１５分）、そして更に一層好ましくは、それよりも多い呼吸（１時間および１２００回の呼吸まで）を含む。較正測定期間が長い程、特に、移動性の設定が実際に近い程、精度向上には有益である。信号は、モーション・アーチファクトが比較的少ないか、あるいはフィルタ処理して抑えなければならない。

次に、ＲＣおよびＡＢ信号における全ての吸入および吐出差を判定する。これは、各最少信号とその直後の最大信号との間の差、および各最大信号とその直後の最少信号との間の差である。少なくとも２００から５００（または２０００まで）通りの差の集合が得られなければならない。最後に、ＡＢ信号において測定した全ての吸入および吐出差の標準偏差の、全てのＲＣ信号において測定した全ての吸入および吐出差のＳＤに対する比として、Ｋを推定する。式で表すと次の通りである。

（数３）
K = SD （ＡＢ吸入および吐出差）／SD（ＲＣ吸入および吐出差）

Ｂ．較正の改良
Ｋの決定および較正を改良することにより、較正の目的で用いられるＲＣおよびＡＢ信号において測定した吸入および吐出差の集合の中にある異常値の排除を更に進める方法が得られる。これらの較正方法は、振子空気現象、即ち、各肺または２つの肺間の内部呼吸に依存するので、差が、小さい標準偏差で一定する周期性呼吸を表すように、差の集合から異常な差の値を厳格に排除することは有益である。本発明者は、以下のマルチステップ異常排除方法を発見した。

改良した方法の一実施形態では、ＲＣおよびＡＢ信号において測定した全ての吸入および吐出差の集合の平均および標準偏差（ＳＤ）を別個に判定する。次いで、平均から±３ＳＤよりも大きな差の値（第１異常値）を破棄し、吸入および吐出差の第１縮小集合を得る。次に、これら吸入および吐出差の第１縮小集合の平均および標準偏差（ＳＤ）を判定し、±２ＳＤよりも大きな全ての差（第２異常値）を破棄して、ＲＣおよびＡＢ吸入および吐出差の第２縮小集合を得る。ＲＣおよびＡＢ吸入および吐出差の第２縮小集合の平均および標準偏差（ＳＤ）を判定し、平均から±１ＳＤよりも大きい差を全て（第３異常値）破棄し、ＲＣおよびＡＢ吸入および吐出差の最終集合を得る。次いで、これらの最終集合を用い、前述の式にしたがってＫを判定する。

更に一般的には、この改良した方法は、異常値排除のステップを削減または増加する方法や、差排除判断基準を有する方法も含む（排除判断基準は、後のステップ程、前のステップよりも増々厳格になる）。例えば、安定性が高いデータには二ステップ方法を適用可能であり、約±２．５ＳＤおよび±１．５ＳＤの排除判断基準を採用する。安定性が低いデータに適用可能な別の実施形態では、連続的な判断基準±４ＳＤ、±３ＳＤ、±２ＳＤ、および最後に±１ＳＤを採用した４つの排除ステップを有する。

更に較正を改良するには、複数の較正係数集合（Ｋおよび任意にＭ）を採用する。各集合は、被験者の特定の挙動状態に当てはまり、被験者がその挙動状態を取るときに用いられる。これが有益なのは、較正係数（複数の係数）が監視対象被験者の状態に依存することが観察されているからである。即ち、Ｋは、被験者の姿勢、例えば、仰向け、傾斜(reclining)、着座、および起立に影響され、更に被験者の活動の強さにも影響される（呼吸速度に影響を及ぼすためである）。この改良の一実施形態では、較正に適したＲＣおよびＡＢ信号を時々刻々測定する。被験者には、周期的に較正プロセスを、例えば、３０分毎、１時間毎等で実行するように命令し、記録されるＲＰデータが、較正係数（複数の係数）を判定し直すのに必要な較正情報を有するようにするとよい。また、恐らくはデータ記録ユニット（図１における３）における入力によって、または加速度計データの処理から、挙動状態の大きな変化が通知されたときに、被験者に較正プロセスを実行するように命令してもよい。

更に別の実施形態では、後続の監視期間中に予期されるそれぞれの挙動状態において、被験者の監視の開始時に、呼吸較正(breathing calibration)を行う。例えば、呼吸較正データは、被験者の仰向け、着座、静止起立、安定した歩行等によって得ることができる。このデータから較正係数（複数の係数）を計算し、被験者の状態とインデックスさせる。後に監視中に、加速度計（またはその他の運動検出器）の信号から、現在の挙動状態を判定し、この現在の挙動状態に対応する所定の係数を用いて、ディジタル化しフィルタ処理したＲＣおよびＡＢ ＲＰセンサ信を、Ｖｔ（および／またはVol）信号に組み合わせる。挙動状態−姿勢、および活動は、前述のように、フィルタ処理した加速度計信号から判定することができる。

更に別の実施形態では、最適線形推定技術、特に、カルマン・フィルタ技術を採用する。ここでは、挙動状態によりインデックスした較正係数（複数の係数）（Ｋおよび任意にＭ）の多数の集合が、線形呼吸モデルを形成する。このモデルに対する入力は、ＲＰ信号、好ましくは、ＲＣおよびＡＢ ＲＩＰセンサ信号、ならびに、例えば、同時加速度計データから得られるような挙動状態の指標である。これによって較正された出力は、高い分解能レベルで、時々刻々、監視対象被験者の現在の挙動状態に応答する。モデルの出力（複数の出力）は、Ｖｔ（および／またはVol）信号であり、これを実際の測定値または現在判定されているＶｔ（またはVol）値と比較すればよい。したがって、カルマン・フィルタまたは予測−補正技術は、当技術分野では公知のように、繰り返し適用し、相対的または実際の現肺容積の最適な推定値を得ることができる。例えば、Maybeck, 1979, Stochastic models, estimation, and control（統計モデル、推定、および制御）, Academic Press, New Yorkを参照のこと。更に、較正信号によって、適応技術を繰り返し適用し、以前に決定したモデル係数を改良即ち洗練化することもできる。例えば、Widrow et al., 1985, Adaptive Signal Processing（適応信号処理）, Pearson Education, Singapore、特に、第９章（適応モデリングおよびシステム同定）を参照のこと。

５．２．５ アーチファクトの認識および基準線の決定
次に、直前のＲＰまたはＲＩＰ信号の処理から得るのが好ましいが必ずしもそうとは限らないＶｔ信号、および任意にVol信号を更に処理し（図２における３１）、呼吸パラメータを決定し、（直前のフィルタ処理後に残っている）アーチファクトを排除し、呼吸事象（例えば、咳、嘆息等）を認識する（図２における３５）ために有用な呼吸パラメータに対して基準線値を確立する。Ｖｔ信号は、肺容積に比例し、任意の単位でよい。Vol信号は、実際の肺容積変化を反映し、リットル単位とするとよい。また、元のディジタル化しフィルタ処理したＲＰセンサ信号も、次の処理のために利用可能である。図４Ａは、４回の正常な呼吸の場合におけるＶｔ信号の一例を示す（時間は右に向かって進んで行く）。また、その時間微分信号dV/dt（図４ＡではdVtで示す）も示す。これは、公知の数値微分法によって計算することができる。Ｖｔは肺容積に比例するので、dV/dtは呼吸空気流に比例する。

これらの信号から種々の呼吸パラメータ、例えば、Ｖｔ信号において呼吸７４についてのパラメータ、およびdV/dt信号において呼吸７６についてのパラメータを得ることができる。例えば、呼吸７４の吸入体積７８は、この呼吸のＶｔ最大値（最終吸入体積）と、直前のＶｔ最小値（直前の最終吐出体積）との間の差であり、この呼吸の吐出体積８０は、この呼吸のＶｔ最大値（最終吸入体積）とその直後のＶｔ最小値（現在の最終吐出体積）との間の差である。この呼吸のピーク吸引流（ＰＩＦ）８２は、dV/dtの最大値であり、そのピーク吐出流（ＰＥＦ）８４はdV/dtの最小値である。更に別の呼吸パラメータ（図示せず）に吸入／吐出率（ＩＥ率）があり、これは、呼吸の吸入時間のその吐出時間に対する比率であり、吸入時間とは、呼吸のＶｔ最大値（最終吸入体積）とその直前のＶｔ最小値（直前の最終吐出体積）との間の時間であり、吐出は、呼吸のＶｔ最大値（最終吸入体積）とその直後のＶｔ最小値（現在の最終吐出体積）との間の時間である。更に別の呼吸パラメータも、これらのパラメータおよび入力信号から機械的に得ることができる。

連続呼吸に対して決定した呼吸パラメータは、時間順に並べた値の系列を形成する。これら呼吸パラメータの時間的系列は、Ｖｔ（および／またはVol ）信号および任意に主要センサ・データと共に、更に別の処理、特に、アーチファクトの除去、基準線の決定、および呼吸事象の認識（図２における３５）のために利用することができる。

Ａ．アーチファクトの除去
アーチファクト呼吸およびそれから得られるパラメータを更に認識し、多くの場合は破棄することが好ましい。このようなアーチファクトの認識は、呼吸信号解釈および生理学の専門家が作成した評価を表す１つ以上の規則を適用することによって行うことが好ましい。一般に、これらの規則は、既知の規則に基づく処理の枠組によって適用するとよく、これにしたがって、規則エンジンが、コード化した規則を適用する。これによって、更に別の規則を追加し、特定の規則集合を特定の用途に選択することもできる。以下の規則は、広範囲で有用であることが認められている。

第１の「異常値排除」規則は、極値パラメータによって呼吸を破棄する。図３Ｃは、特定の呼吸パラメータの値の振幅のヒストグラム例であり、横軸上にパラメータの振幅、縦軸上にパラメータ値の相対数を有し、対応するパラメータ振幅が測定シーケンスにおいて見出される。通常、このようなヒストグラムによって、多数の小振幅事象と大振幅事象の末端との間の広範な振幅範囲が明らかとなる。振幅閾値７０よりも左側にある小振幅事象は、大抵の場合、残留ノイズ、心臓、およびその他の非呼吸成分に起因する。振幅閾値７２よりも右側にある大振幅事象は、大抵の場合、残留する運動および同様のアーチファクトに起因する。したがって、第１規則は、低振幅閾値未満のパラメータ値、または大振幅閾値よりも大きいパラメータ値によって、呼吸を破棄する。嘆息、咳等のような有効な呼吸事象を破棄しないために、これらの閾値は、大抵の場合、平均または振幅ヒストグラムのモードから約±３ＳＤ以上に設定する。

ここでは、「２５％未満規則」および「実在呼吸規則」として知られる更に別の規則がある。「２５％未満規則は、吸入または吐出、あるいは双方の体積が所定の閾値未満であることを、呼吸にマークする。マークされた呼吸は、以後の処理では考慮されず、これらの期間から得られたパラメータ値は、破棄することが好ましい。好ましくは、所定の閾値は、初期Ｖｔ較正に所定の比率を乗算して求める。初期Ｖｔ較正（Vt-cal）は、較正期間中（または初期較正期間における５分間の窓内）に測定した吸入または吐出あるいは双方の体積の平均である。したがって、Vt-calを求めるために用いられる体積は、単一ステップまたはマルチステップの異常値排除を克服したそれらである。一般的に好ましい所定の比率（ｒ１）は、２０％と３０％との間にあり、更に好ましくは２４％と２６％との間にあることが分かっており、以下では２５％を用いる。しかしながら、この比率は、被験者の測定値および評価を用いて低振幅アーチファクトを最良に認識することによって、被験者毎に個別化するとよい。

「実在呼吸規則」は、直前の最終吐出体積（直前の最少Ｖｔ）および現在の最終吐出体積（直後の最少Ｖｔ）の差が固定比にVt-calを乗算した値未満である場合、呼吸がアーチファクトよりも「真の」呼吸に近いと思われるものとして、呼吸に印を付ける。一般的に好ましい固定比率（ｒ２）は、２００％と３００％との間にあり、更に好ましくは２４０％と２６０％との間にあることが分かっており、以下では２５０％を用いる。この場合も、被験者の測定値および評価を用いて低振幅アーチファクトを最良に認識することによって、この比率を被験者毎に個別化するとよい。

Ｂ．基準線の決定
次に、呼吸行為(performance)における偏差を評価する際に続いて用いるために、呼吸信号およびパラメータ系列の基準線値を決定する。基準線値は、好ましくは、移動メディアン・フィルタ(moving median filter)によって決定する。これは、観察群の統計的中心値を返すことが知られている。適用すると、メディアン・フィルタで処理された、呼吸パラメータの現時点における値は、現時点を含む時間窓内で発生したこの１組のパラメータの値（有効であり、アーチファクトではない）の統計的中心値となる。

窓は、変化に対する応答性（窓の大きさが小さい程高い）と、ノイズおよびアーチファクトに対する不感性（窓の大きさが大きい程高い）との間の折衷案として選択する。各パラメータは、異なるレートの変化、ノイズ特性、アーチファクト感度等を有する可能性があるので、パラメータ毎に別個に窓特性を選択することが有益である。一般に、好ましい窓の大きさは、１から１．５〜３、そして５分までで様々な値を取り、１〜２分が大抵の場合に相応しい。また、パラメータによっては、窓の位置が、現時点の時刻よりも進んでいる場合、ここを中心とする場合、これよりも遅れている場合がある。例えば、１分を中心とする窓は、一般に呼吸体積には有用であることがわかっているが、０．５分（約９回の呼吸）窓（約０．５秒）は、一般に呼吸時間（例えば、吸入時間）に有用であることがわかっている。また、個々の被験者に記録した信号の特性にばらつきがあるような場合、これらの被験者の窓特性を個別化すると有益であると考えられる。

５．３ 事象の認識
呼吸事象は、呼吸信号の時間的経過、および既に決められている呼吸パラメータの系列から認識する。本発明は、呼吸停止３７、減呼吸３９、嘆息４１、咳４３、および発話４５の、多くの場合１つよりも多い方法による、自動事象認識３５を含む。追加の事象も、同様に、信号およびパラメータ系列から認識することができる。信号および処理済み医療データの、担当介護士による手作業での検討が、患者に適正な方策を求めるために必要となる場合が多い。

５．３．１ 呼吸停止認識
１つ以上の呼吸体積パラメータ、吸入、吐出、または一回換気量が、所定の閾値に当該呼吸パラメータの現中心基準線(running median baseline)を乗算した値よりも小さい場合、呼吸がほぼない期間を呼吸停止として認識する。呼吸パラメータの現中央基準線は、現在の呼吸よりも前の１から５分の窓を用いて決定することが好ましい。好ましくは、所定の閾値は、１と５０％との間であり、更に好ましくは、１５と３０％との間である。更に、呼吸減少期間が、５および６０秒の間の所定期間、更に好ましくは８および２０秒の間続くはずである。多くの被験者では、２分の窓、２５％の閾値、および１０秒の期間とすると、適当な呼吸停止認識が得られる。これらの値は、過去の監視データを参考にして、個々の被験者毎に個別化するとよい。

呼吸事象の呼吸体積パラメータが基準線よりも２５％以下の場合、以前のアーチファクト除去ステップによって破棄されている可能性があるので、代わりに、２回の真の呼吸（直前の章を参照のこと）が５秒よりも長い時間と６０秒よりも長い時間との間、更に好ましくは合計時間が８秒よりも長い時間、更に一層好ましくは１０秒よりも長い時間である場合に、呼吸停止を認識してもよい。呼吸停止期間は、呼吸パラメータの現中心値から除外することが好ましい。

図５Ｂは、実際に睡眠中の被験者から観察した呼吸停止期間の一例を示す。時点１１０および１１２の間、呼吸体積が実質的に直前の基準線の２５％未満となったことが容易に分かる。時点１１０および１１２の間の時間は、約１９秒である。

Ａ．呼吸停止分類
呼吸停止は、概略的に、閉鎖性呼吸停止（気道の閉鎖）、中枢性呼吸停止(central apneas)（ＣＮＳ呼吸駆動の減少のため）、または混合呼吸停止（閉鎖および駆動減少の要素による）に分類することができる。この判断基準にしたがって認識された呼吸停止は、追加のパラメータを用いて、呼吸位相関連（「ePhRL」または「位相関連」として知られている）、および呼吸努力評定（「努力」）(respiratory effort rate)を生ずるとして、分類することができる。

位相関連は、呼吸毎に、ディジタル化しフィルタ処理したＲＣおよびＡＢ信号を参照することによって決定する。ＲＣおよびＡＢは、双方の信号が同じ方向の移動を反映し、双方が吸入の間増加するか、吐出の間減少すれば、同相（または協調されている）と考えられる。これらの信号は、逆方向の移動を反映する場合、位相外れ（または非協調）と考えられる。この位相関連の好ましい数値尺度に、ＲＣおよびＡＢ信号が位相外れになっている呼吸の時間の割合がある。ePhRLが約４０％未満である場合、正常な呼吸から軽い呼吸非協調を意味し、約４０％よりも大きい値は、中程度から激しい呼吸非協調を意味する。したがって、閉鎖性呼吸停止は、１に近い割合によって特徴付けられ、中枢性呼吸停止は正常に近い割合によって特徴付けられる。低および高割合を分離するために好ましい閾値は、２０％と６０％との間であり、一般に４０％が相応しい。閾値は、個々の被験者毎に個別化するとよい。

また、図５Ｂもこれらの呼吸位相関連を示す。図示した呼吸停止期間の外側、即ち、時点１１０以前および時点１１２以降に、この図の観察では、ＲＣおよびＡＢ信号が同相となり、双方とも、ePhRLが正常な範囲内で、吸入および吐出において一緒に移動している。呼吸停止期間、時点１１０および１１２の間、これらの信号は位相外れとなり、ＲＣは吸入しようとし、一方ＡＢは吐出しようとし、あるいはその逆であり、ePhRLは１．０に近く、閉鎖性プロセスを意味する。

努力評定とは、呼吸努力の評定であり、呼吸努力は、呼吸体積パラメータが、当該パラメータ対する現中央値(running median)の１％および２５％の間であることによって認識される。好ましくは、呼吸努力は、一回換気量が現中央値の一回換気量の５％および２５％の間であることによって認識され、中央値は、好ましくは、先行する２分窓から取る。呼吸努力評定に対する基準線は、中央呼吸評定であり、中央値は先行する１０回の呼吸の窓からであることが好ましい。

これら追加のパラメータに関して、呼吸停止期間を、ePhRLが約４０％よりも大きく、呼吸努力評定が基準線呼吸の約７５％よりも大きい場合に、閉鎖性に分類することが好ましい。呼吸停止期間は、ePhRLが約４０％未満である場合、あるいは呼吸努力評定が基準線の約０から２５％よりも低い場合に、中枢性に分類する。

５．３．２ 減呼吸認識
減呼吸期間は、呼吸停止認識に用いた呼吸パラメータによって、減呼吸と認識することができるが、閾値は緩和されている。好ましくは、減呼吸、即ち、減少呼吸は、現呼吸に先行する１から５分の期間を有する窓からの現中央値との比較によって認識する。好ましくは、所定の閾値は２０および８０％の間であり、更に好ましくは、現中央値基準線の２５および５０％の間である。呼吸停止認識におけると同様、減呼吸の期間は、５および６０秒、更に好ましくは８および２０秒の所定期間続くはずである。多くの被験者では、２分の窓、５０％未満で２５％よりも大きい呼吸体積閾値、および１０秒の期間によって、適切な減呼吸認識が得られる。これらの値は、過去の監視データに鑑み、個々の被験者毎に個別化するとよい。

減呼吸期間は、呼吸停止と同様の分類のように、閉鎖性、中枢性、または混合に分類する。減呼吸期間は、呼吸パラメータの現中央値に含まれることが好ましい。

５．３．３ 嘆息認識
本発明は、嘆息も認識する。嘆息の頻度が普段よりも多い場合、心配事に対する反動または鬱状態のような、心理学的窮状を意味することが多い。心配事は、肺疾患の間に起こる呼吸困難の期間に増加することが多い。

最初に、モーション・アーチファクトを除外するために、嘆息は、吐出期間が好ましくは０．２５から３秒の範囲を有する所定の閾値よりも大きい、真の呼吸（直前の章で定義した）であるとする。好ましい閾値は約１秒であり、個別化してもよい。これらの要件を満たす呼吸は、その吸入または吐出あるいは一回換気量が所定の閾値に現中央点基準線値（好ましくは、先行する２分窓から判断する）を乗算した値よりも大きい場合、嘆息であると認識する。Ｖｔの好ましい閾値は、例えば、１００および１０００％の間であるが、多くの被験者では、閾値を約２５０％とすると、適切な嘆息の認識が得られた。この閾値は、過去の監視データに鑑み、個々の被験者事に個別化してもよい。

図４Ｂは、実際の被験者から記録したデータにおける嘆息９０を示す。ここでは、Ｖｔ９０が現中央値Ｖｔ基準線の約６００％となっている。この嘆息は、空気流信号dV/dtにおいても、容易に視認できるピーク９２として表されており、このピークは、呼吸体積の増大を比例的に反映し、確認する。

５．４ 発話認識
監視対象の被験者が発話しているか否かの認識は重要である。何故なら、発話は、呼吸、心臓、およびその他の生理的信号の意味を限定する内容を規定することが多いからである。発話中における生理的活性化(physiological activation)の量および質は、音の大きさや迅速性のような、認知的、感情的、そして多くの場合関連する身体的要因に依存することは知られている。これらは、発話および誘出される感情を引き起こす交流(interaction)の社会的および相互人間的な質によって変動する可能性がある。したがって、発話の生理的効果は、秘匿対公開設定、つまらない会話対魅力のある会話、および軽い会話対論争のような討論では異なると思われる。個性、社会的不安、内気等は、これらの反応に影響を及ぼすが、移動性設定においては現在では利用されていない。日常生活における殆どの社会的交流は発話を伴い、発話の発生を監視することによって、社会的交流を定量化することができる。ある不規則性は、社会的行動の増加（例えば、熱中）または減少（鬱状態）によって特徴付けることができる。社会的隔絶が高齢者において頻繁に観察され、自律神経系の機能の望ましくない変化に関連することが示されている。

発話は、一層複雑で潜在的に個別の行動であり、対応して呼吸に影響を及ぼし得るので、この章では最初に発話認識に適したパラメータおよび閾値を決定できる概略的な方法を、続いて健康な被験者群への適用を紹介する。

５．４．１ パラメータおよび閾値の決定
本発明は、個々の被験者および同様の被験者群において発話を認識するためのパラメータおよびそれらの閾値を選択する系統的な方法も提供する。以下では発話に関して記載するが、これらの方法は、咳、嘆息などのような他の呼吸事象を認識するための代わりのパラメータを選択するためにも適用することができる。一般に、これらの方法では、種々の発話および非発話状態における測定、次いで測定データの処理による、発話（またはその他の呼吸事象）の認識に最も適したパラメータおよび閾値の決定を伴う。

したがって、第１ステップでは、意図した被験者、または意図した被験者群の代表に、発話活動および非発話活動に加わるように命令する。非発話活動は、着座、休憩、歩行、軽い活動等を含むことができる。発話活動は、これら非発話活動によって達成されても、されなくてもよい。更に別の発話活動は、種々の社会的設定、例えば、朗読または多人数の会話において発生することもある。呼吸信号およびパラメータの移動性ＲＩＰ記録は、これら様々な種類の活動の間に行われ（好ましくは前述のように）、測定値を分別し活動毎に分ける。

種々のパターン認識および分類技法を集合化データに適用し、データ群を適した精度で区別することが可能なパラメータおよび値のパターンを識別することができる。生理的内容を得るために発話を認識するのに適した精度は、２０％および１０％以上の間の偽陰性または陽性率とすることができる。他の用途では、別の精度レベルを要求することもある。パターン分類技法は、分散技術の再帰、判別、および分析のような統計的手法、あるいは、例えば、ｋ−平均クラスタリングまたはニューラル・ネットワークによる自動集合化または整列を含むことができる。例えば、Duda et al, 2000 2nd ed., Pattern Classification（パターン分類）, Wiley Interscience, New Yorkを参照のこと。好ましい出力は、単一のパラメータまたは数個のパラメータの関数、および発話または非発話の認識、そして任意に発話種類の認識に用いることができる閾値である。

５．４．２ 発話認識
この方法は、移動性被験者からのＲＩＰ記録における発話期間および非発話期間を認識するためのものである。一般に、線形判別分析から、数個の単一パラメータ、および２〜４個のパラメータの線形結合を識別し、１〜１５％レベルの偽陰性または陽性認識精度が得られる。これらのパラメータおよび組み合わせは多大な計算を必要とせず、数時間から半日以上の移動性記録に定期的に適用することができる。

好ましい認識パラメータは、吸入／吐出（ＩＥ）率、分別吸入時間(fractional inspiratory time)、吸入流量、および吐出時間、ならびにこれらの変種の線形結合および割合係数を含む。これらのパラメータは全て発話を検出する際に良く機能したが、ＩＥ率、即ち、吸入時間の吐出時間に対する比率（約０．５２の閾値による）は、他よりも一貫性があり、最も好ましい。この単一パラメータの利点は、以下を含む。１）例えば、吸入流量のような体積較正を必要としないので、測定が容易である。２）個々の呼吸毎に判定することができるが、変動性パラメータ(variability parameter)は少なくとも１分の測定期間を必要とする。このような発話認識については、第６章において更に説明する。

カットオフを０．５２に固定したＩＥ率は、健康な人々には最も好ましいが、ある種の他の母集団、例えば、慢性閉塞性肺疾患または喘息の患者には、異なるＩＥ率のカットオフ、あるいは異なるパラメータまたは結合さえも必要となる場合もある。更に、発話認識パラメータは、他の母集団に対して個別化することができ、同じ方法、例えば、線形統計分析を用いる個人でさえも、個別化することができる。

５．５ 呼吸困難およびＦＥＶ_１の監視
本発明は、医療上重要である患者の呼吸困難およびＦＥＶ_１／ＶＣ双方に対する指標および代用品(surrogate)を連続的に監視する方法を提供する。

Ａ．呼吸困難の監視
呼吸困難は、困難なまたは苦しい呼吸の間隔、呼吸停止の感じ、または呼吸努力が完全に満たされないという体感であり、多くの人によって、肺過膨張の直接または間接的な結果であると考えられている。過膨張は、換気効果と実際の吸気との正常な関係、または吸気の知覚を乱し、患者が呼吸しようとしても、意図した空気量が移動しない。例えば、過膨張は、横隔膜の静止時の長さを変化させ、肋間筋肉をそれらの最適範囲の外側に変化させ、発生する力を弱め、これによって所与の神経呼吸駆動に対する空気流も減少させる。

呼吸困難は、喘息の患者または慢性閉塞性肺疾患（ＣＯＰＤ）、例えば、肺気腫および慢性気管支炎の患者には一般的であり、これらおよびその他の疾患の評価および管理には重要な兆候である。喘息は過膨張に至る虞れがある。何故なら、関連する気管支収縮が気道の抵抗を著しく高め、ＣＯＰＤもそうする可能性があるからである。何故なら、関連する肺容積弾性係数の増大が、肺の反跳を減少させ、こうして空気流が制限されるからである。これらの状態が悪化すると、空気流は、呼吸駆動の増大には正常に応答しなくなり、心配になり、場合によっては恐怖に陥る。

呼吸困難は、主観的な評価であり、患者に質問することによって、現在は評価している。しかし、患者は質問の繰り返しにうんざりしたり、呼吸停止の感覚に慣れてしまうこともある。したがって、換気駆動の客観的尺度が呼吸困難の感覚に繋がることが示されれば、医療的に貴重であろう。所与の活動レベルにおいて、このような換気駆動尺度は、比較的一定でなければならず、したがって換気駆動の増大は患者の困難を示すことになる。

本発明は、呼吸困難に関連する通気努力のこのような客観的指標を提供するが、重要なのは、これが、前述のＶｔおよび／またはVol信号から容易に決定可能な呼吸パラメータに依存するということである。好ましい呼吸困難指標は、微小換気体積(V_E)のピーク吸入流量(P₁F)に対する比率、即ち、比率V_E/P_lFである。この比率は、最終肺容積（ＥＥＬＶ）が増加するに連れて、即ち、進行する過膨張と共に高くなり、したがって、先に論じたように、呼吸困難と関連付けられる。図６は、肺過膨張に対するV_E/P_lFの依存性を実証する。ここでは、過膨張は、正の終端吐出圧力（ＰＥＥＰ）である１０ｃｍＨ_２Ｏに被験者を置くことによって誘発させた。ＰＥＥＰが開始した後、被験者の肺が正の終端吐出圧力によって膨張すると、ＥＥＬＶは劇的に約３リットル増加し、ＰＥＥＰが終了した後、ＥＥＬＶは正常に戻り始める。ＥＥＬＶの増加（過膨張）に伴う力学的変化によって、空気流がＶ_Ｅに減少する。次いで、被験者は、呼吸神経駆動を横隔膜および肋間筋肉に対して増大させることによって補おうとして、P_lF（空気流）の増加に至る。したがって、V_E/P_lF比は、ここでは、少なくとも１／３に減少する。このように、過膨張は、呼吸困難に伴い、V_E/P_lFの減少を招く。

要約すると、V_E /P_lFは、呼吸筋肉効率および呼吸停止のインデックスとして役立つ。これは、Pif/Vtによって測定した呼吸中枢からの駆動を、通気によって測定した呼吸系の出力に繋げる。動的肺過膨張のために非効率的な呼吸筋肉収縮を有する患者では、この比率が低下するように、通気のいずれの変化に対しても不釣り合いな程度まで呼吸駆動を増加させる。

被験者のV_E /P_lF比は、本発明のシステムおよび方法によって、連続的に、またはほぼ連続的に、または間欠的に、あるいは偶発呼吸困難や規定した診断閾値のように、必要に応じて、容易に監視することができる。V_Eは、１分のような時間期間に発生する吸入または吐出体積（図４Ａにおける７８および８０）の和として、または平均吸入または吐出体積の呼吸速度との積によって容易に決定することができる。P_lF（図４Ａにおける８２）は、呼吸体積曲線の導関数として決定することができる。V_E/P_lF比の値は、これらのパラメータを、肺容積に比例するＶｔ曲線から測定したのか、または実際の肺容積の変化を反映するVol曲線から測定したのかには実質的に依存してはならない。更に、V_E/P_lF比の閾値は、母集団または個人のいずれについても決定することができる。例えば、第１閾値を主観的な呼吸困難のいつもの開始時に選択することができ、第２閾値を危険な可能性のあるレベルの呼吸困難の開始時に選択すればよい。

Ｂ．ＦＥＶ_１の監視
ＦＥＶ_１として知られる１秒間の強制吐出体積は、気道の妥当性および開通性の尺度として容認されている。この標準的な測定を行うには、被験者に、できるだけ迅速に最大肺容積（最大吸入の後）から最大努力で吐出させ、残留容積までずっと続ける。肺疾患およびその進行の重要な指標の１つは、ＦＥＶ_１の肺活量（ＶＣ）に対する比率、即ち、ＦＥＶ_１/VCである。通例、健康な被験者の比率は少なくとも８０〜８５％であり、一方喘息またはＣＯＰＤの被験者のＦＥＶ_１/VC比率は７０％以下である。この値は、ＣＯＰＤが悪化するに連れて、または喘息の発作が始まり進行するに連れて低下する。

ＦＥＶ_１/VCは広く用いられ信頼されているが、この比の測定は、間欠的であっても疲れる。ＦＥＶ_１は、概念的には単純であるが、被験者からの最大の努力を要求する。同様に、単純であるが、ＶＣの測定では、被験者が完全に吐出することを要求するが、これは、肺容積が残留容積に近いと特に厳しい要求となる。したがって、ＦＥＶ_１/VCの測定を繰り返すことは、特に、そもそも被験者が肺に問題を抱えている場合には、重荷であり、したがって、精度が低い場合が多いと思われる。したがって、通常の呼吸中に容易に測定可能なＦＥＶ_１/VC（またはその切迫した変化）の代用物または指標があれば、広範囲にわたって有効であろう。

本発明の方法およびシステムは、組み合わせて発見した換気パラメータの連続的監視を行い、現在のまたは切迫した気道障害の、高感度であるが侵襲性の低い尺度とする。これらの組み合わせは、一般に、ＦＥＶ_１/VC比の代用または指標として用いることができる。このような組み合わせの１つに、ピーク吸入流P_lFの一回換気量VTに対する比率がある。この比率、即ち、P_lF/VTは、ＦＥＶ_１/VCを反映するように流量の変化に対して感度が高いことが発見されている。更に別の組み合わせは、ピーク吐出に達するまでの時間を吐出時間で除算したものである。ピーク吐出に達する時間とは、吐出の開始（Ｖｔ信号の最大値）からdVt/dtの最小値（最大外方空気流）までの時間であり、吐出時間とは、吐出の開始から終了までの全時間である。多数の線形回帰によって評価されている、２０％よりも大きなＦＥＶ_１.0の変化を追跡する変数の更に別の組み合わせ（感度および特定性は０．９０に等しく、ＲＯＣ分析は、０．８９の曲線（ＡＵＣ）の下のエリアを解明した）には、吐出周期性流量のピーク対平均の比率（P_EF/M_EF)、一回換気量に対する胸郭の寄与（％ＲＣ）、および胸腹収縮不同時性による吐出時間の分画(fraction)がある。これらのパラメータを、本発明のシステムおよび方法によって既に説明したように連続的にまたは間欠的に監視すれば、P_lF/V_T、P_EF/M_EF、％ＲＣ、および胸腹収縮不同時性による吐出時間の分画の呼吸毎の分析が可能となる。

これらのパラメータおよび指標は、既に説明したＶｔ（および／またはVol）信号の処理の検査から容易に判断できる。したがって、これらの指標は、患者の努力や、不適当な患者の努力に伴う要因が混同する可能性には無関係に、気道の障害を監視することを可能にする。その結果、この努力に無関係な気道開通性の評価は、実際の疾患や慢性疾患中の危機における気道変化の開始または初期段階を示す、早期警告を与えることができる。

５．６ 咳検出方法
これは、咳検出および分類方法（図２における４３）の代替実施形態を含む。咳の検出は重要である。何故なら、例えば、咳の頻度が高まることは、心不全を伴うことが多い、実際の肺水腫の早期兆候でもあるからである。一般に、これらの方法を進めると、ＡＢ、ＲＣ、およびＶ_Ｔ信号を含む入力呼吸パラメータからの呼吸事象候補、ならびに任意に音声入力からの音響事象候補を認識する。次いで、呼吸事象候補および関連する音響事象候補の特定の組み合わせから咳事象を検出する。咳の種類および重大度は、呼吸および音響事象パラメータの値によって判別することができる。

５．６．１ 第１の咳認識方法
第１の咳検出方法によれば、咳は、好ましくは、０．２５から３秒の範囲を有する所定の閾値よりも長い吐出期間によって、真の呼吸と認識しなければならない。有用で好ましい閾値は、約１秒であり、個別化してもよい。次いで、これらの判断基準を満たす呼吸を、咳として認識するのは、そのピーク吐出流（ＰＥＦ）が、先行する二分窓から判定した現中心基準線ＰＥＦ値の所定の閾値よりも大きいときである。好ましいＰＥＦ閾値は、現中心基準線ＰＥＦ値の１００％と１０００％以上との間である。多くの被験者について、ＰＥＦ閾値を約２５０％よりも大きくすると、適切な咳の認識が得られた。この値は、過去の監視データに鑑み、個々の患者毎に個別化してもよい。

図４Ｃは、実際の被験者から記録したデータにおける２回の咳９４および９８を示す。ＰＥＦは、dV/dt曲線から判定する。この曲線では、同じ２回の咳９６および１０２が、短く鋭い呼気として容易に視認可能である。ここで、咳９６に対するＰＥＦは現中心ＰＥＦ基準線の約４００％であり、一方咳１０２については、ＰＥＦは基準線の約３８０％である。

５．６．２ 代わりの咳認識方法
図７は、特定的に音響入力を咳検出の補助として組み込んだ、代わりの咳検出方法を示す。この副章におけるこの図以降では、入力データおよび派生データを、多くの場合以下の略語で示す。

要約すると、図７は、最初に、高周波数ノイズおよび低周波数移動アーチファクトを更に抑える（未だ十分に抑えられていない場合）ように設計された２つの帯域通過フィルタによって、帯域通過のフィルタ処理をＶｔおよびＡＢ信号に施すことを示す。フィルタ処理した信号のピーク・ピーク電力（または呼吸振幅、即ち、最大吐出および最大吸入の差）が既定の閾値−Ｔを超過している場合、呼吸および音響信号双方を更に詳しく検査して、咳事象の存在を検出する。

音声信号（例えば、咽喉マイクロフォンからの）を発話認識フロント・エンドで処理し、音声事象が有声または無声音声を含むか否か判定を行う。この判定に重要なのは、派生信号PITCHmである。これは有限期間におけるピッチ値の平均である。被験者が発話しているかまたは会話に加わっている場合、この平均レベルは大きく上昇するはずであり、咳の場合には上昇しないはずである。ピッチ値を計算するには、ケプストラムまたはメル周波数ケプストラム係数（ＭＦＣＣ：Mel Frequency Cepstral Coefficient）内にあるピーク・ピーク電力を測定する。別の重要な派生信号に、PITCH信号がある。音響信号処理から出力されるのは、図９におけるEVT形跡によって示されるようなパルスであり、タイミングおよび期間は、入力音響データにおいて検出される有意の音声事象のそれに等しい。

音響事象がない場合、咳は検出されない。音響事象がある場合、その期間によって、フィルタ処理した呼吸信号のどれを咳痕跡検出器に印加すべきか判断する。音響事象の期間が比較的長い場合（即ち、中心有意音響事象(median significant sound event)よりも長い）、例えば、≧６００ミリ秒の場合、低周波数帯域通過フィルタ処理した呼吸データLFBを、咳検出器によって分析する。音響の期間が比較的短い場合、即ち、中心有意音響事象よりも長い）、例えば、≦６００ミリ秒の場合、高周波数帯域通過呼吸データHFBを分析する。この信号選択は、移動およびモーション・アーチファクトの適切な濾過に繋がるので、咳の痕跡を一層明確に検出できることが分かっている。種々の咳の痕跡を、後に図９、図１２Ａ、図１２Ｂ、図１３Ａ、図１３Ｂ、図１４Ａおよび図１４Ｂに示す。

５．６．３ ステップの詳細−ディジタル・フィルタおよびピーク電力の判定
一回換気量形跡Ｖｔは、ＲＣおよびＡＢ帯域の線形加重和であり、並列な２つのＦＩＲ帯域通過フィルタを通され、（フィルタ処理した信号の最大値が反映する）ピーク電力を測定して、可能な事象の存在を判定する（ピーク電力が閾値Ｔを超過した場合）。入力呼吸信号のフィルタは、好ましくは、有限インパルス応答（ＦＩＲ）設計であるが、位相ずれや時間遅延を最少に抑えた無限インパルス応答（ＩＩＲ）フィルタを用いてもよい。ここでは、対応する音響信号と時間的に一致したままでいるように、呼吸信号の位相が十分に落ちついていなければならない。

１０２４というフィルタ長は、十分に尖鋭な特性を得るのに好ましい長さとして決定した。図８Ａおよび図８Ｂは、前述の低および高帯域通過フィルタの周波数および位相応答を示す。表２．１は、これら好ましいそれぞれのフィルタのパラメータを纏めた。

これらのフィルタは、Matlab^TMＦＩＲ最少二乗法を、モデル次数(model order)１０２４と共に用いて設計した。前述のフィルタのパラメータは、胸郭および腹部から取り込んだ十分な呼吸移動（ＲＣおよびＡＢ）を保持しつつ被験者の肉体の移動が可能な範囲で、フィルタ処理するように選択する。

Ａ．電力閾値
ピーク・ピーク電力は、ここでは正方向に向かう信号上の最大点から負に向かう信号上の最少点までを測定したものと定義することが好ましく、咳事象候補がフィルタ処理した呼吸信号の中にあるか否か判定を行うために用いられる。この閾値を超さない場合、有意な咳があるとは見なされない。信号ＬＦＢ、ＨＦＢ、およびＦＡＢを測定してこの判断を下す。ＦＡＢはＡＢをフィルタ処理した形跡からのフィルタ残余であり、ＲＣおよびＡＢが位相外れでＶｔに対して減算効果を有し、帯域内における真の成果(true effort)を減少させる場合に有益である。閾値−Ｔは、漠然とした全体的に約２００ｍｌの吐出体積であるが、特定の個人の個々の母集団毎に調節することができる。

５．６．４ 詳細なステップ−音声事象検出器
図９は、マイクロフォン入力の一例−形跡ＭＩＣから得た、音響包絡線の一例−形跡ＳＥを示す。音響包絡線は、全ての呼吸帯域と同じサンプリング周波数にダウン・サンプリングすることが好ましい。これは、５０Ｈｚであることが好ましい。これによって、フィルタ残余および呼吸信号の微分(derivation)の影響を最少に抑える。このダウン・サンプリングは、マイクロフォン・ストリームからの３０サンプル毎の平均化を伴い、１５００Ｈｚにおいてサンプリングして、５０Ｈｚの音響包絡線が得られる。

次に、音声事象検出および期間判定のために音響包絡線を処理する。音声事象の開始が認識されるのは、音響包絡線が、較正背景ノイズ閾値の倍数に決定してある閾値を超えたときである。好ましくは、ノイズ閾値は、長期間のマイクロフォン記録（２４０時間までを用いた）から較正し、音響包絡線信号の３０のレベルを表す＋１および−１の間の信号の変動を監視することによって決定する。ノイズ閾値の２倍、即ち、６０の値が、有効な事象閾値であることが分かっている。音声事象が終了するのは、音響包絡線がノイズ閾値未満（ここでは、３０の値）に再度低下したときである。咽喉マイクロフォンを使用することにより、背景ノイズの影響を最少に抑える。音声事象は、ＥＶＴ形跡において、振幅が１０で期間が音声事象の長さに等しいパルスで印される。図９は、音声事象と、関連するＨＦＢ信号も示す。

５．６．５ ステップの詳細−咳痕跡検出器
有意な音声事象が、呼吸事象の可能性がある事象と一致した場合、これらの信号から１つを、音声期間に応じて選択し、咳痕跡について更に分析する。有意な音声咳事象の期間を判定した後、ＬＦＢ信号またはＨＦＢ信号を更に分析して、咳痕跡の存在を調べる。分析する通過帯域の選択を補助するために、音声事象咳期間を測定する。音声事象の期間が短い場合、即ち、約６００ｍｓ未満の咳事象の場合、ＨＦＢ信号を分析する。何故なら、咳事象は、短い程、高い周波数成分を有すると思われるからである（咳をより短時間で終わらせるため）。逆に、咳の期間が長い程信号の周波数は低くなるので、咳痕跡検出のためにはＬＦＢ信号を選択する。

図９のＨＦＢ形跡に典型的な咳痕跡を示す。咳痕跡は、ＨＦＢまたはＬＦＢあるいは双方の痕跡において、終末が鋭く（高ピーク吐出流に対応する）、続いて鋭い吸入が生ずることが好ましい。この吸入は、咳事象に分類された音声事象に伴って発生する。ＨＦＢまたはＬＦＢが残す最も低いサンプル値は、関連する音声事象の中央領域近くに位置することが好ましい。中央位置は、音声事象の開始から音声事象期間の３３％よりも大きく、事象の終了から事象期間の３３％未満の時点として規定する。更に、この最小値は、−Ｔ値を超過し、当該事象の期間における中心サンプルのいずれの側にも一定の傾きを有していなければならない。この場合の−Ｔ値は、静かな呼吸または落ちついた呼吸と識別された領域において計算した、個々の被験者についての平均呼吸体積に基づいて計算すればよい。したがって、各サンプル間の差[x(n)-x(n-1)]は、痕跡の中心よりも前では負であり、中心の跡および終了の前では正のはずである。ノイズは信号から濾過されており、計算に影響を及ぼすことはない。

更に、最小値のいずれの側においても、ＨＦＢまたはＬＦＢの形跡の傾き（およびこれらの傾きの傾斜(gradient)）は、以下の制約に該当することが好ましい。第１に、痕跡は最小値の中心サンプルの各側において同様の傾きとほぼ対称であるとよい。中心サンプル即ち最小値の各側における各傾きの終点は、信号が最大振幅に達し減少し始める前の点である。これらの終点は、事象の終了の通過後または事象の終了前で、事象時間期間の５０％よりも長い時間期間を超過してはならない。これらの厳しい制約を適用することにより、咳状事象を誤って検出する可能性は大幅に低下する。あるいは、閾値を指定し、ピーク吐出流および連続するピーク吸入流がこの閾値を超えなければならないとしてもよい。

５．６．６ ステップの詳細−フロント・エンド処理
このステップは、音声波形を簡潔なパラメータ表現（好ましくは、周波数対置換の表現形態）に変換し、咳の音響を発話音響と区別できるようにする。一般的に、前者の方が周波数が低く、後者の方が周波数が高い。したがって、周波数関連閾値を簡潔な表現の中に規定し、閾値未満の信号が咳の音響である可能性が高くなるようにするとよい。

咳の呼吸痕跡を有する事象候補は、それに伴う音響事象が咳の音響を含むと判定されない場合、咳であると見なされない。逆に、咳の音響痕跡を有する事象候補は、それに伴う呼吸事象が十分な咳特性を有していない場合、咳であると見なされない。代わりの検査では、咳−発話閾値よりも高い信号電力があっても、咳−発話閾値よりも低い信号電力が増大すると、この音響事象を咳として受け入れる。事象候補は、PITCH値がある閾値よりも高い場合も、真の咳とは見なされない。PITCH値がこの閾値の直ぐ下である場合でも、PITCHm値がこの閾値よりも高ければ、事象候補は咳とは見なされない。PITCHmは既定の時間期間内における全てのPITCH値の平均である。これらPITCH値の平均がこの閾値よりも高い場合、この事象の前後に発話があり、したがってこの事象は恐らく発話であると推測する。

発話音声信号の特性は、約１０ミリ秒の時間増分の間固定であると見なされ、したがってこのようなセグメントの間に分析を行う。発話信号の固定部分の一例を図１１Ａに示す。長い期間では、発話信号特性は、異なる音声音響が生成されることを反映して、変化することが確実である。短時間のスペクトル分析は、音声信号をそのように特徴付ける公知の方法である。

パラメータを用いて音声信号のピッチ特性を抽出し表現する技法には、線形予測符号化（ＬＰＣ）、メル周波数ケプストラム係数（ＭＦＣＣ）等のように、数種類ある。ＭＦＣＣは、咳検出方法には好ましいことがわかっている。一般に、ＭＦＣＣは、人の耳の臨界帯域幅の既知の変形に基づいているので、これらの係数は、メル周波数目盛りで表され、１０００Ｈｚ未満の周波数では線形であり、１０００Ｈｚよりも高い周波数では対数表現である。これらのフィルタは、発話の音素的に重要な特性を取り込む。

Ａ．ＭＦＣＣの判定
図１０は、ＭＦＣＣを計算する好ましいプロセスのフローチャートである。これは、１５００Ｈｚでサンプリングした音声入力を処理する。このサンプリング周波数は、発話および咳の成分を解明するために選択したものである。このプロセスの最初のステップ、即ち、フレーム・ブロック化ステップは、連続する発話信号をＮサンプルずつのフレームにブロック化する。隣接するフレームは、Ｍ個のサンプル（Ｍ＜Ｎ）によって分離されている。最初のフレームは、最初のＮ個のサンプルから成る。２番目のフレームは、最初のフレーム後のＭ番目のサンプルから始まり、Ｎ−Ｍサンプルだけ重複する。同様に、３番目のフレームは、最初のフレームの後の２Ｍ番目のサンプル（即ち、２番目のフレーム後のＭ番目のサンプル）から始まり、Ｎ−２Ｍサンプルだけ重複する。このプロセスは、音声全体を１つ以上のフレームにブロック化するまで続ける。好ましいブロック化パラメータＮおよびＭは、Ｎ＝６４（〜４０ミリ秒の窓処理(windowing)と同等であり、高速radix-2ＦＦＴが容易に行われる）、およびＭ＝３２である。

この窓処理ステップは、個々の各フレームを窓処理し、フレーム境界における信号の不連続を最少に抑える。連続で円滑な窓を用い、各フレームの開始および終点において０になるように信号を徐々に低下させることによって、スペクトル歪みを最少に抑える。ウィンドウをw(n)、０≦ｎ≦Ｎ−１と定義すると、窓処理の結果は次の信号となる。ここで、Ｎは各フレーム内のサンプル数である。

好ましくは、本発明ではハミング・ウィンドウを用いる。これを次のように定義する。

テキスト処理ステップは、高速フーリエ変換であり、Ｎ個のサンプルから成る各フレームを時間ドメインから周波数ドメインに変換する。ＦＦＴは、離散フーリエ変換（ＤＦＴ）を実施するための周知のアルゴリズムであり、以下のように、Ｎ個のサンプルの集合{x_n}に対して定義する。

一般に、X_nは複素数である。得られる系列{X_n}は、以下のように解釈される。ゼロ周波数はｎ＝０に対応し、正の周波数０＜ｆ＜Ｆ_ｓ／２は値１≦ｎ≦Ｎ／２−１に対応し、負の周波数−Ｆ_ｓ／２＜ｆ＜０はＮ／２＋１≦ｎ≦Ｎ−１に対応する。ここで、Ｆ_ｓはサンプリング周波数を示す。このステップの結果は、多くの場合スペクトル・ピリオドグラムと呼ばれている。図１１Ｂは、図１１Ａの信号のスペクトルを示す。

次のステップは、メル周波数ラッピング (mel-frequency wrapping)である。精神物理学の研究によって、発話信号に対する音響の周波数内容の人による知覚は、線形の目盛りには従わないことが示されている。したがって、Ｈｚ単位で測定した実際の周波数ｆの各トーン毎に、「メル」目盛りと呼ばれる目盛り上で主観的ピッチ(subjective pitch)を測定する。「メル」目盛りは、１０００Ｈｚ未満では線形周波数間隔を有し、１０００Ｈｚを超えると対数的間隔を有する。基準点として、認知聴覚閾値(perceptual hearing threshold)よりも４０ｄＢ高い、１ｋＨｚトーンのピッチを、１０００メルと既定する。したがって、以下の近似式によって、Ｈｚ単位で所与の周波数ｆについてメルを計算する。

通常、主観的音声スペクトルのシミュレーションは、フィルタ・バンクによって行われ、図１１Ｃに示すように、メル目盛り上にフィルタを均一に離間して配置する。最初のバンクは、三角形状帯域通過周波数応答を有し、間隔および帯域幅は、一定のメル周波数インターバルによって決定される。入力信号S(ω)のメル・フィルタ処理したスペクトルは、したがって、S(ω)が入力であるときのこれらの周波数の出力電力から成る。メル・スペクトル係数Ｋの数は、通例、１８および２４の間で選択する。尚、このフィルタ・バンクは、周波数ドメインにおいて適用され、したがってこれは単に、図１１Ｃの三角形状窓係数に、図１１Ｂの時間周波数スペクトルを乗算したに過ぎないことに注意すること。この方法では、低サンプリング・レートのために、Ｋ＝１０メル・スケールのフィルタ・バンクを、入力信号周波数スペクトルに適用することが好ましいことがわかっている。

ケプストラム判定の最終ステップでは、対数メル・スペクトルを時間に逆変換し、メル周波数ケプストラム係数（ＭＦＣＣ）を求める。発話スペクトルのケプストラム表現では、所与のフレーム分析に対する信号の局部的スペクトル特性の表現が得られる。メル・スペクトル係数（およびそれらの対数）は実数であるので、これらは、離散余弦変換（ＤＣＴ）を用いて時間ドメインに変換することができる。したがって、最終ステップの結果であるメル電力スペクトル係数をＳ~_k, k=1,2,...,Kとすると、ＭＦＣＣのc~_nは次のように計算することができる。

尚、最初の成分c~₀をＤＣＴから除外すると好都合であることに注意すること。何故なら、これは話者特定情報を殆ど搬送しない入力信号の平均値を表すからである。
図１１Ｄは、図１１Ａ〜図１１Ｃに既に示した発話信号に対するケプストラム出力を示す。咳および無声発話音響は、一般に、メル周波数閾値の１．５〜２を下回ることが分かっている。有声発話が信号例にあることは明白である。何故なら、ピッチが高くなると、信号電力がこの閾値よりも高くなるからである。PITCHm信号は、メル周波数スペクトルの単純平均、または電力加重平均等として得ることができる。PITCH信号は、離散余弦変換から得られる最大メル周波数ケプストラム係数として得られる。

５．６．７ 咳の例
Ａ．慢性閉塞性肺疾患（ＣＯＰＤ）
慢性閉塞性肺疾患（ＣＯＰＤ）は、一般に、呼吸機能の漸進悪化に至る１群の肺障害のことである。肺への空気流を徐々に阻害するＣＯＰＤに共通する２つの原因は、気管支炎および気腫である。慢性気管支炎では、気道が遮断され膨張し、気管支内にある粘液生成腺が拡大し、過剰な量の粘液が肺内に分泌される。したがって、この形態のＣＯＰＤでは、この過剰な粘液を排出するために咳をする必要性が高くなる。

図１２Ａおよび図１２Ｂは、システムに入力した音声およびビデオ入力、ならびにソフトウェア・アプリケーションで実施した本発明の方法によって特定したＣＯＰＤ咳を示す。ＨＦＢおよびＬＦＢ形跡は、図１２Ａにおける真の咳が、鋭い（短期間および大きい空気流）吐出、およびそれに続く鋭い吸入によって特徴付けられていることを示す。更に、音声事象が、咽喉マイクロフォン入力から検出され、低いピッチを有することを特徴とし、咳の音響を含む可能性が最も高い。図１２Ｂは、ＣＯＰＤ患者からの数種類の咳でない事象および１つの真の咳事象を示す。咳でない事象は、低ピッチ音響事象として見られ、関連する呼吸上の咳指標（ＬＦＢまたはＨＦＢ信号における鋭い吸入および吐出）が欠如していた。一方、真の咳事象は、伴う音響および呼吸事象が適正な特性を有することを特徴とする。

Ｂ．膵胞性繊維症（ＣＦ）
膵胞性繊維症（ＣＦ）は、生命を脅かす多系統状態であり、主に肺および消化器系を冒す。ＣＦは、粘りのある粘液の分泌を招き、気道を閉塞し、気道から粘液を排出しようとするために、頻繁に咳をする必要性が生ずる。咳をすると、粘液が緩み、呼吸が楽にできるようになる。

図１３Ａおよび図１３Ｂは、二人のＣＦ患者からの咳を示す。関連する形跡、特にＨＦＢおよびＬＦＢ形跡の検査から、これらの咳は、ＣＯＰＤの咳よりも激しく、大きな振幅および／または高い空気流を有することが明らかである。更に、振幅が十分なので、咳の痕跡は、フィルタ処理しない一回換気量（Ｖｔ）の形跡において特定が容易である。

Ｃ．感染後の咳（ＰＩＣ）
感染後の咳（ＰＩＣ）は、上気道のビールス感染後に最もよく見られ、いかに粘液分泌が増大しようと、それに係わらずに永続する炎症により、咳を誘発する。図１４Ａおよび図１４Ｂは、ＰＩＣ咳を２例示す。

５．６．８ 咳の重大度および分類
検出した咳事象は、更に、バンドパス・フィルタで処理した肺容積データ、ＬＦＢおよびＨＦＢ信号の特定的な特性を抽出することによって分析することもできる。この特性には、咳痕跡の深さ即ち振幅、および咳痕跡終端における反射吸入駆動(reflex inspieratory drive)が含まれる。咳の病理的原因の判別を可能にする尺度には、咳の深さの、安静で緊張のないときの呼吸の特定した期間における、被験者毎に計算した平均吐出体積との比率が含まれる。これによって、個別の較正に基づいて重大度を判定することができ、したがって肺疾患の判定に役立つ。更に、このような尺度には、咳事象中における吐出および吸入体積双方の変化率が含まれる。更に別の処置は、咳のセグメントを分析し、咳事象の異なる間隔で体積の変化率を比較することである。

簡略化する場合には、これらの信号（咳の音量）の振幅およびそれらの傾き（空気流量）を組み合わせて、咳の種類を互いに判別するための診断用判断基準を得ることができる。これらのパラメータは、例えば、咳の異なる深さ、および咳事象の終端における反射的吸入活動を反映する。フィルタ処理していないＶｔにおける咳痕跡の外観は、特にひどい咳の別の指標となる。より単純なこれらの鋭さ判断基準を用いると、前述の例によって示したように、ＣＦ咳は、その重大度が高い可能性が高いために、ＣＯＰＤ咳は、その重大度が低い可能性が高いために、そしてＰＩＣ咳は、その重大度が中間である可能性が高いために、認識できることがわかっている。フィルタ処理していない一回換気量形跡Ｖｔに咳痕跡がある場合、最も高い重大度の咳が伴う。

６．例
この章では、第５．５章の方法による発話認識について説明する。
Ａ．方法
測定被験者には、９人の男と９人の女が含まれ、平均年齢（±ＳＤ）が２１．３（±１．２）歳であり、全員肉体的には非常に健康であり、現在喫煙しておらず、呼吸器系疾患の病歴もない。実験手順を全て説明した後、参加者全員が告知に基づく同意を承認し、その書式に署名した。次いで、被験者はＲＩＰ記録外皮（図１における１）−VivoMetrics, Inc. (Ventura, CA)からのLifeshirt（登録商標）を装着し、測定を開始するに当たり、８００ｍｌのバッグの吸入および吐出を７回行い、それを完全に充満させそして空にすることによって、デバイスの呼吸センサの較正を行った。較正は、着座および起立姿勢において行った。次いで、被験者は静かに着座し（安静着座、４分）、連続して発話（過去１週間の経験について）（発話、４分）、質問書に記入した（書き込み、＞４分）。次に、彼らは通常の１日を送り、翌朝モニタを取り外して帰宅した。

メモリ・カード上に記憶されているデータを、パーソナル分析コンピュータ（図１参照）にダウンロードし、本発明の分析および表示ソフトウェアによって処理した。記録上に較正期間を印し、自動的に分析してＫおよびＭ係数を求めた。これらの係数を用いて、ＶｔおよびVol信号をＲＩＰ信号から計算し、記録全体にわたって各呼吸毎に種々の呼吸パラメータを、ＶｔおよびVol信号から計算した。実験期間（安静着座、発話、書き込み）を印し、４つの１分セグメントに分割した。呼吸数、一回換気量、毎分換気量、吸入流速（平均）、吸入時間、吐出時間、ＩＥ率、端数吸入時間（全時間で除算した吸入時間）、および一回換気量に対する胸郭の寄与を計算し、それらの呼吸毎の変動性を、変動の係数（CV%)にインデックスさせた。つまり、CV%は、傾向線周囲の標準偏差を、あらゆる線形傾向を除去した後の平均で除算した比率であった。呼吸期間の最初の３回の１分セグメントの判別分析によって、各パラメータのＦ比率、期間の％補正分類、および最適カットオフ・スコアを計算した。パラメータおよびそれらの分類特性の線形結合も判定した。事前に省略した最後の１分セグメントのデータについて、判別関数の妥当性を判断した。

Ｂ．結果
表１および表２は、発話を他の２つの活動と区別するために、１〜３分における呼吸パラメータを、その有効サイズによってランク付けして示す。これらの表において、「Ｆ」は、Ｆ比率であり、ｄｆ＝１．１７、「cutoff」は条件を判別するための最適なカットオフ評価点であり、「%false」は、このカットオフを用いた偽分類の割合である。

分析によって、吸入／吐出（ＩＥ）率、分別吸入時間、吸入流量、および吐出時間（ならびにこれらのパラメータの変形の百分率係数）は全て、発話を認識するのに適することがわかった。発話対書き込みに対する判別分析では、ＩＥ率を最初に取り上げた。これは、発話を書き込みと区別し、カットオフ判断基準を０．５２としたときに、約９８．５％の正しい分類が行われ、発話を安静着座と区別する際にも同様に成功した。ＩＥ率は、発話を書き込みと対比するための最も高いＦ比率が得られた。これらの表から、タイミング、タイミングの可変性、および体積パラメータの複合線形関数も、発話を他の２つの状態から分離するのに有用であると言えることが示唆されている。呼吸パラメータの一部、例えば、呼吸速度または肋骨寄与は、発話をはっきりと示さなかった。発話を安静着座と対比するには、タイミング可変性の４つのパラメータが最も大きなＦ比率を有した。恐らく、これは、安静着座では、非常に規則的な呼吸パターンが呼吸毎に得られたからであろう。ヒストグラムから、ＩＥ率は、作業の範囲内ではほぼ正常に分布し、安静着座および書き込みについての分布は、発話中のそれらとは殆ど重複しなかったことが明らかとなった。

ここに記載し特許請求する発明は、ここに開示した好適な実施形態によってその範囲が限定される訳ではない。何故なら、これらの実施形態は、本発明のいくつかの形態の例示であることを意図するからである。同等の実施形態も、本発明の範囲に該当するものとする。実際、ここに示し説明した実施形態に加えて、本発明の種々の修正が、以上の説明から、当業者には明白となろう。このような修正も、添付した特許請求の範囲に該当するものとする。

この中で多数の引例を引用したが、その開示内容は、ここで引用したことにより、全ての目的に対してその全体が本願にも含まれるものとする。更に、これらの引例はいずれも、先にどのように特徴付けたかには係わらず、ここに特許請求する主題の本発明に対する従来技術としては認めないこととする。

図１は、好ましい体積変動記録および分析システムを示す。 図２は、好ましい信号処理エレメントを示す。 図３Ａは、信号分布およびフィルタ特性の一例を示す。 図３Ｂは、信号分布およびフィルタ特性の一例を示す。 図３Ｃは、信号分布およびフィルタ特性の一例を示す。 図４Ａは、呼吸信号の一例を示す。 図４Ｂは、呼吸信号の一例を示す。 図４Ｃは、呼吸信号の一例を示す。 図５Ａは、別の呼吸信号の例を示す。 図５Ｂは、別の呼吸信号の例を示す。 図６は、肺過呼吸の影響を示す。 図７は、咳検出の代替方法を示す。 図８Ａは、好ましいディジタル・フィルタ応答を示す。 図８Ｂは、好ましいディジタル・フィルタ応答を示す。 図９は、咳の間に記録したデータの一例を示す。 図１０は、音声信号における好ましいピッチ判定方法を示す。 図１１Ａは、音響信号の一例におけるピッチ判定を示す。 図１１Ｂは、音響信号の一例におけるピッチ判定を示す。 図１１Ｃは、音響信号の一例におけるピッチ判定を示す。 図１１Ｄは、音響信号の一例におけるピッチ判定を示す。 図１２Ａは、ＣＯＰＤ患者の咳の一例を示す。 図１２Ｂは、ＣＯＰＤ患者の咳の一例を示す。 図１３Ａは、ＣＦ患者の咳の一例を示す。 図１３Ｂは、ＣＦ患者の咳の一例を示す。 図１４Ａは、ＰＩＣ患者の咳の一例を示す。 図１４Ｂは、ＰＩＣ患者の咳の一例を示す。 図１４Ｃは、ＰＩＣ患者の咳の一例を示す。 図１４Ｄは、ＰＩＣ患者の咳の一例を示す。

Claims (10)

1. 監視対象被験者の呼吸を反映する信号を処理する方法であって、
   フィルタ特性に関連付けられた前記被験者の現在の挙動状態は、被験者の運動が殆どあるいは全くないか、被験者の平均的な運動があるか、または、被験者の大きな運動があるか、どうかを決定するステップと、
   決定された前記挙動状態に１対１に関連付けられたフィルタを選択するステップと、
   非呼吸起源の信号成分を低減するために、被験者の挙動状態に関連付けられ前記選択されたフィルタのフィルタ特性に基づいて前記呼吸信号をフィルタ処理するステップと、
   胸郭（ＲＣ）の大きさを示す少なくとも１つのフィルタ処理した信号と、腹部の大きさ（ＡＢ）を示す少なくとも１つの信号との組み合わせから肺の容積を示す信号（Ｖｔ）を求めるステップと、
   を備えている、方法。
2. 前記被験者の現在の挙動状態が仰臥、着座、起立、歩行、および走行の群から選択した少なくとも１つの状態であるかどうかを決定するステップを更に備えている、請求項１記載の方法。
3. 被験者の挙動状態を、前記被験者の位置および加速度に応答する少なくとも１つの加速度計からの信号によって示す、請求項１記載の方法。
4. ２つの異なるフィルタ特性を異なる時点において前記呼吸信号に適用し、前記適用するフィルタ特性を、被験者の挙動状態に応じて選択する、請求項１記載の方法。
5. 前記挙動状態が被験者の運動減少を示すときに適用する少なくとも１つのフィルタ特性の通過帯域は、前記挙動状態が被験者の運動増加を示すときに適用する少なくとも１つのフィルタ特性の少なくとも通過帯域を含む、請求項３記載の方法。
6. 前記被験者の運動は、該被験者の位置および加速度に応答する加速度計からの少なくとも１つの加速度計信号の全てまたは一部の電力によって示される、請求項５記載の方法。
7. モーション・アーチファクトを、前記被験者の運動を示す１つ以上の更に別の基準信号に時々刻々応じて低減するステップを更に備えている、請求項１記載の方法。
8. 被験者の運動を示す前記基準信号は、前記被験者の位置および加速度に応答する少なくとも１つの加速度計からのハイパス・フィルタ処理した信号から成る、請求項７記載の方法。
9. 前記被験者の心臓活動を示す１つ以上の更に別の基準信号に時々刻々応じて心臓信号成分を抑えるステップを更に備えている、請求項１記載の方法。
10. 被験者の心臓活動運動を示す前記基準信号は、ＥＣＧ入力および／またはパルス酸素計から得た信号から成る、請求項９記載の方法。

Images