HRVの非線形指標とは?自分に必要?

DFA α1やサンプルエントロピーなどの非線形の心拍変動(HRV)指標は、心拍の間隔のばらつきではなく、そのパターンを表します。拍が前の拍を覚えているかどうか、次の拍がどれだけ予測しやすいかです。SDNNやRMSSDではわからないことが見え、研究用ソフトウェアにも登場します。ウェアラブルデバイスで自分の回復を追うだけなら必要ありません。2015年の欧州のポジションステートメントは、臨床での利用はこれまでのところかなり限られていると結論づけています。

このページの内容 5セクション

要点

  • 非線形指標は、拍のばらつきではなく、拍のパターンを表します。記憶(DFA)、予測しやすさ(エントロピー)、期外収縮への反応(タービュランス)です。
  • DFA α1が1に近いことは、拍がさまざまなスケールで互いを覚えていることを意味し、0.5に近いと、拍はノイズのようになります。エントロピーが0に近いのは繰り返しのパターン、2以上はほぼランダムです。
  • RMSSDの上位版ではありません。異なる種類の情報を加えるもので、設定と記録時間に大きく左右されます。
  • Sassiらの2015年のポジションステートメントは、これらの方法は信号の理解を進めたものの、臨床の道具としての成功はかなり限られていた、と述べています。
  • ウェアラブルデバイスによる毎日の測定には、自分のベースラインと比べたRMSSDかlnRMSSDで十分です。

非線形の心拍変動(HRV)指標は、心拍の間隔のばらつきではなく、そのパターンを表します。SDNNは間隔がどれだけ広く動くか、RMSSDは拍ごとにどれだけ跳ねるかを示します。非線形指標が問うのは、拍が前の拍を覚えているか( トレンド除去変動解析 DFA α1 トレンド除去変動解析の短期スケーリング指数。1に近いと、拍はあらゆる尺度で記憶をもち、0.5に近いとノイズのようにふるまいます。α2は、より長い区間での同じ傾きです。 DFA α1についてさらに詳しく (DFA))、次の拍がどれだけ予測しやすいか(サンプルエントロピー( SampEn サンプルエントロピー(SampEn) 心拍の系列がどれだけ予測しにくいか。一致した拍の連なりが、一拍長くしても一致する確率の対数に、マイナスを付けた値です。0付近は繰り返しのパターン、2以上はほぼランダムです。 サンプルエントロピー(SampEn)についてさらに詳しく )と近似エントロピー( ApEn 近似エントロピー サンプルエントロピーより古い、単位のない規則性の指標。連続する拍の各並びを自分自身とも一致するものと数えるため、短い記録では値が低いほうへ偏ります。 近似エントロピーについてさらに詳しく ))、そしてリズムが期外収縮にどう応えるか( 心拍タービュランス(HRT) 心拍タービュランス(HRT) 心室期外収縮の後に現れる、心拍リズムの短い加速と、それに続くゆっくりした回復。24時間心電図(ECG)から、タービュランス・オンセット(%)とタービュランス・スロープ(ms/拍)の二つの数値で測ります。 心拍タービュランス(HRT)についてさらに詳しく 、 減速能 減速能 心臓が遅くなるときの、その遅くなり方の平均的な大きさ(ミリ秒)。直前の拍より遅いすべての拍のまわりのパターンを平均して求めます。心筋梗塞の後に使われる、24時間の迷走神経の指標です。 減速能についてさらに詳しく (DC))です。ウェアラブルデバイスで自分の回復を追うだけなら、これらは必要ありません。研究用ソフトウェアには登場します。これらに関する2015年の欧州のポジションステートメントは、信号の理解は進んだものの、臨床の道具としての利用はこれまでのところかなり限られていると結論づけました。

なぜ「非線形」なのか

記録を一つ用意し、拍の順序をでたらめに入れ替えて、SDNNを計算してみてください。値は変わりません。ヒストグラム、 三角指数 HRV三角指数 記録内のNN間隔の数を、そのヒストグラムで最も高い棒の高さで割った値。全体の変動を表す単位のない幾何学的指標で、乱れた拍が少しあっても、ほとんど動きません。 HRV三角指数についてさらに詳しく 、バエフスキー指数も変わりません。これらの指標が見ているのは、間隔のばらつきだけです。シャッフルすると順序が壊れますが、生理学が宿っているのは順序のほうです。呼吸の波、もっとゆっくりした血圧の波、遅い拍の次にまた遅い拍が来やすいという性質などです。順序を読み取り、結果が値の単なる合計ではない方法を、非線形と呼びます。Peng、Havlin、Stanley、Goldbergerは1995年に、健康な心臓には長距離の相関があり、それが重い心不全では崩れることを示して、この一群の方法を有名にしました。

よく目にする二つの指標

DFA α1:拍は覚えているか?

リズムは過去を覚えているか

教材用の記録から5分間

α1(4〜16拍)
0.49
α2(16〜64拍)
0.65

休息十分の記録は合成データで、4拍ごとのきれいな呼吸の波に、拍ごとのランダムなノイズが重なっています。そのため4〜16拍ではほとんどノイズのように見え(α1 0.49)、ゆっくりした波によって16〜64拍にはいくらか記憶が現れます(α2 0.65)。安静時の実際の記録では、通常はα1が約1か、それ以上になります。

JavaScriptをオンにすると操作できます。

教材用の記録5分間のトレンド除去変動解析(DFA)です。各点は、一定の拍数のボックスの中での累積和の典型的な揺れで、両対数の軸では、点が直線に沿って並び、その傾きがα1(4〜16拍)とα2(16〜64拍)です。ホワイトノイズは0.5、1/fノイズは1になります。教材用の記録は合成データで、約0.5になります。安静時の実際の記録は、通常は約1か、それ以上です。系列を切り替えたり、ノイズを混ぜたりして、傾きの変化を見てください。

4〜64拍のボックスサイズに対する変動を示す両対数のグラフで、ボックスサイズごとの点と、2本の近似直線があります。4〜16拍の直線がα1、16〜64拍の直線がα2です。表示欄にはα1とα2が表示され、その意味を説明する文が続きます。4つのボタンで、休息十分の記録、ホワイトノイズ、1/fノイズ、異所性拍動ありの記録を選べます。スライダーでは、選んだ系列にホワイトノイズを混ぜられます。

トレンド除去変動解析は、時系列のふらつきが、より長い区間を見るにつれてどう大きくなるかを、両対数軸で調べます。記憶のないホワイトノイズの傾きは0.5、すべてのスケールで記憶のある1/fノイズは1、ランダムウォークは1.5になります。教材用の記録の五分間では、4〜16拍でα1 = 0.49、16〜64拍でα2 = 0.65になります。ホワイトノイズに切り替えると、どちらの傾きも0.5付近になり(α1 0.48)、1/fノイズに切り替えると1付近になります(α1 1.04)。安静時の実際の記録は、通常は約1かそれ以上になります。教材用の記録の値が低く出るのは、合成データだからです。四拍に一回のきれいな呼吸の波に、拍と拍の間のランダムなノイズが重なっており、短いスケールではノイズのように見えます。強く規則的な呼吸の揺れは、実際のα1の値も同じ方向に引き下げることがあります。運動では、強度が上がるにつれて約1から0.5に向かって下がり、そのため一部のトレーニングアプリはこれをリアルタイムで表示しています。

サンプルエントロピー(SampEn):次の拍はどれだけ予測しやすいか?

次の拍はどれだけ予測できるか

教材用の記録から5分間

SampEn
1.67
ApEn
1.02
許容範囲 r
5.7 ms
m拍の一致数
13
m+1拍の一致数
1

休息十分、m = 2、r = 0.20 × SD = 5.7 ms。 強調した2拍の連なりは、5.7 ms以内で他の13個の連なりと一致します。そのうち1個は、両方の連なりを1拍長くしても一致したままです。 記録内のすべての連なりで数えると、サンプルエントロピーは1.67です。呼吸のリズムによって、次の拍がある程度予測できます。

JavaScriptをオンにすると操作できます。

教材用の記録の最初の64拍です。m拍の連なりが1つ強調され、各拍の前後に±rの帯が付いています。帯の中に収まる他の連なりが「一致」で、色の濃い箱は、さらに1拍長くしても一致する連なりです。サンプルエントロピー(SampEn)は、記録全体で数えた一致のうち、1拍長くしても残る割合の対数にマイナスを付けたものです。近似エントロピー(ApEn)は、それより古い類似の指標です。mとrを変えたり、周期的な系列やランダムな系列に切り替えたりしてみてください。

心拍間隔64個を順に並べた折れ線グラフです。強調された箱は、m拍のテンプレートの連なりを示し、各拍の前後に±rミリ秒の許容の帯が付いています。薄い箱は、帯の範囲内でテンプレートに一致する、他のm拍の連なりです。濃い箱は、両方の連なりを1拍長くしても一致するものです。表示欄には、サンプルエントロピー、近似エントロピー、許容範囲(ミリ秒)、2つの一致数が表示され、その意味を説明する文が続きます。2本のスライダーで、m(1〜4拍)とr(標準偏差の0.1〜0.6倍)を設定できます。3つのボタンで、周期的な系列、休息十分の記録、ランダムな拍を選べます。

サンプルエントロピーは、二拍からなるすべての並びについて、許容範囲内で一致する並びを探し、その並びのうち一拍後にもまだ一致するものがどれだけあるかを調べます。その割合の対数にマイナスをつけたものがエントロピーです。教材用の記録のスコアは1.67、完全に周期的な時系列は、一致したものがずっと一致し続けるため0.01、ランダムな拍は2.35です。RichmanとMoormanは2000年に、旧来の近似エントロピーを改良するためにこれを考案しました。近似エントロピーは、各並びを自分自身とも一致するものと数え、記録の長さによってずれていました。どちらも、並びの長さmm、許容範囲rr、拍の数に大きく左右されるため、異なるプログラムや異なる記録時間の値は比べられません。

臨床で使われるもの

心拍タービュランスと減速能は、リスクの指標として開発され、おもに心筋梗塞を経験した人の24時間心電図で検証されました。タービュランスは、心室期外収縮に対するリズムの反応から圧受容器反射を読み取ります。減速能は、一日のうちの心拍の減速を平均したものです。Bauerら(2006)の、心筋梗塞から生還した1,455人の研究は、さらに二つのコホートでも検証され、減速能は駆出率やSDNNより死亡をよく予測することを示しました。これらは循環器の医師のためのリスク指標であり、体調管理の数値ではありません。ウェアラブルデバイスも表示しません。

総説の結論

Sassiらは2015年に、欧州心臓病学会のe-Cardiologyワーキンググループと欧州不整脈学会のために、この一群の方法全体を総説しました。フラクタル指標、エントロピー、シンボリックダイナミクス、カオスです。結論は、これらの方法は信号の理解を深めたこと、専門的な研究室が選ばれた患者で得た結果は興味深いものの、臨床での利用を確立するものではないこと、リスクの特定のような臨床の道具をつくる試みの成功はかなり限られていること、そして新たな前向き研究が必要であることです。ShafferとGinsberg(2017)の概説は、これらの指標を、欧州心臓病学会と北米ペーシング・電気生理学会のタスクフォース(Task Force)の指標と並べて挙げ、たとえばDFAが数時間にわたる記録向けに設計されたことを指摘しています。

これらは必要か?

ウェアラブルデバイスによる朝や夜間の測定には、必要ありません。自分のベースラインと比べたRMSSDやlnRMSSDが、それらの測定が得意とする問いに答えてくれます。一方、非線形の数値には、手首のセンサーが出すより、きれいで長い記録と、固定しておかなければならない設定が必要です。研究論文を読む人、Kubios HRVのような同様のプログラムを使う人、トレーニングアプリでDFA α1を目にする人は、知っておくと役に立ちます。自分の胸部ストラップの記録で試すには、RRコンバーターを使うと、書き出したデータをKubiosで開けるファイルに変換できます。用語集のDFA α1、サンプルエントロピー、近似エントロピー、心拍タービュランス、減速能のページには、それぞれの計算式と計算例があります。通常これらの隣に並ぶ周波数の数値については、HRVアプリのLF、HF、VLF、n.u.とは?で説明しています。

よくある質問

アプリが運動中にDFA α1を表示します。何のためのものですか?

一部のトレーニングアプリは、α1を運動強度の指標として使っています。楽な運動では1付近にあり、強度が上がるにつれて0.5に向かって下がります。Rogersらは2021年に、15人のランナーのデータから、0.75を横切る点が第1換気閾値(VT1)に当たると提案しました。この方法にはクリーンな胸部ストラップの信号が必要で、検証が続いています。

サンプルエントロピーは高いほどよいのですか?

原則としてそうではありません。完全に規則的なリズムは0付近、ランダムな拍は高い値になり、どちらの極端も健康的なパターンではありません。研究では、エントロピーの指標は年齢とともに低下しますが、値は設定に大きく左右されるため、比べられるのは、同じプログラムと同じ記録時間の値だけです。

時計でこれらを計算できますか?

計算自体はどのスマートフォンでも動きますが、ウェアラブルデバイスはこれらを表示しません。多くの指標は、手首のセンサーが出すより、きれいで長い間隔の時系列を必要とします。Kubios HRVのような研究用プログラムなら、胸部ストラップや心電図(ECG)の書き出しデータから計算できます。

出典

  1. Advances in heart rate variability signal analysis: joint position statement by the e-Cardiology ESC Working Group and the European Heart Rhythm Association co-endorsed by the Asia Pacific Heart Rhythm Society Europace, 2015
  2. Quantification of scaling exponents and crossover phenomena in nonstationary heartbeat time series Chaos, 1995
  3. Physiological time-series analysis using approximate entropy and sample entropy American Journal of Physiology: Heart and Circulatory Physiology, 2000
  4. Deceleration capacity of heart rate as a predictor of mortality after myocardial infarction: cohort study The Lancet, 2006
  5. An overview of heart rate variability metrics and norms Frontiers in Public Health, 2017
  6. Heart rate turbulence: standards of measurement, physiological interpretation, and clinical use: International Society for Holter and Noninvasive Electrophysiology consensus Journal of the American College of Cardiology, 2008
  7. A new detection method defining the aerobic threshold for endurance exercise and training prescription based on fractal correlation properties of heart rate variability Frontiers in Physiology, 2021