医療従事者・セラピストの方へ

共鳴周波数呼吸法とは？

共鳴周波数呼吸法（RFB）は、ゆっくりとしたペースで行う呼吸法の一種で、患者は個人の共鳴周波数で呼吸します。共鳴周波数とは、心血管系の共鳴特性を利用して心拍変動を最大にする呼吸数（通常1分間に4.5〜6.5回）のことです（Vaschillo et al., 2006; Lehrer & Gevirtz, 2014）。

共鳴時には呼吸性洞性不整脈（RSA）が最大になります。心拍数は吸気時に上がり呼気時に下がりますが、その振幅が最も大きくなり、心拍数の振動は安静時のベースラインの4〜10倍に達します。繰り返し練習することで、圧受容器反射感受性が高まり、自律神経の調節が改善すると考えられています（Lehrer & Gevirtz, 2014; Shaffer & Meehan, 2020）。

RFBは、一般的なゆっくり呼吸のアプリとは異なります。1分間に1回ずれるだけでも、生理的な効果は大きく低下します（Steffen et al., 2017）。Precise Breathは、集団の標準値を使うのではなく、HRV分析によって患者ごとの呼吸ペースを特定します。

エビデンス

RFBとHRVバイオフィードバックについては、さまざまな臨床像を対象に研究が発表されています。

これらは、共鳴周波数呼吸法に関する既発表の研究結果です。Precise Breathはウェルネスと呼吸トレーニングのためのツールであり、医療機器ではありません。効果には個人差があります。

Precise Breathにおけるプロトコルの実装

患者の共鳴周波数を特定するための確立された臨床プロトコルでは、1回のセッションで複数の呼吸ペースを、心拍変動の分析を用いてテストします（Lehrer, 2000; Shaffer & Meehan, 2020）。Precise Breathは、このプロトコルを複数セッションにわたって継続的に精緻化できるよう応用しています：

• センサー：RR間隔に対応したBluetooth LE胸ベルトであれば使用できます。Polar H10（推奨）は、研究で検証された1拍ごとの精度を備え、HRV測定の参照基準となっています（Task Force, 1996; Gilgen-Ammann et al., 2019）。動作確認済み：Polar H10、Garmin HRM Dual。HR Service 0x180Dを使用するその他の標準的なBLE胸ベルトも動作するはずです。
• ブロック構成：各呼吸ブロックは、30秒の移行期間、15秒の安定化期間、120秒の測定ウィンドウで構成されます。これはLFパワー推定に必要な最低2分というTask Force（1996）の基準を満たしています。
• スコアリング：呼吸周波数におけるFFTのスペクトル振幅（重み85%）と、ヒルベルト変換による位相コヒーレンス（重み15%）を組み合わせます。信号品質によるゲーティング（アーチファクト率＋変動係数）で、信頼できないブロックを除外します。
• 調整モード：あらかじめ設定された5つのペース（4.5〜6.5 BPM）を、約14分間の1回のセッションでランダムな順にテストし、ベースラインの推定値を確立します。
• 探索モード：セッションごとに前後のペースを探索する、複数セッションにわたる適応型アルゴリズムです。確度を高めながら、患者個人の共鳴周波数へと収束していきます。
• IBIによる適応：十分なデータがあれば、アプリは安静時の心拍間隔と相関する共鳴周波数のセッション間の変動を考慮します。これはCapdevila et al.（2021）の知見と一致しています。
• プライバシー：すべてのデータは患者の端末内に保存されます。データが送信されることはなく、アカウントも不要です。

HRV分析とスコアリングの詳細

測定方法を理解したい臨床家・研究者向けの情報です。

信号処理パイプライン

各呼吸ブロックでは、胸ベルト型センサーから一連のRR間隔が得られます。これを5段階のパイプラインで処理します：

1. アーチファクト除去：2段階の適応型パイプラインです。(1) 範囲フィルターで300〜1500 msの範囲外の間隔にフラグを付けます。(2) 適応型の逐次差分フィルターで、拍間差の移動中央値と両側基準を用い、前後の両方から外れている一時的なスパイクだけにフラグを付けます。フラグが付いた間隔は削除せずに線形補間で置き換え、タイミングを保ちます。使用可能な間隔が30未満のブロックは除外されます。
2. 心拍タコグラム：アーチファクトを除去したRR間隔を瞬時心拍数（bpm）に変換し、連続した時間軸上に配置します。
3. 等間隔リサンプリング：タコグラムを自然3次スプラインで補間し、4 Hzの等間隔グリッドに変換します。これは対象となる最も高い呼吸周波数を大きく上回るもので、FFTによるスペクトル解析に適した形になります。
4. トレンド除去と窓処理：周波数領域の解析の前に、直流成分を除去し、スペクトル漏れを抑えるためにハニング窓を適用します。
5. スペクトル解析と位相解析：窓処理した信号にFFTを適用します。これとは別に、ヒルベルト変換によって、呼吸の基準に対する心拍の瞬時位相を抽出します。

総合スコア

120秒の測定ウィンドウごとに、総合スコア（0〜1）を算出します：

• スペクトル振幅（85%）：目標呼吸周波数（±0.015 Hzの帯域）におけるスペクトルの大きさを、スペクトル全体の大きさに対して正規化した値です。目標ペースでの心血管系の共鳴の強さ、つまりRSAが最大化されているかを示す主要な信号をとらえます。絶対値ではなく比率で自己正規化されるため、安静時心拍数やHRVの大きさがセッション間で異なっても影響を受けにくくなっています。
• 位相コヒーレンス（15%）：呼吸の基準に対する心拍の瞬時位相の平均合成ベクトル長（MRL）で、ヒルベルト変換によって算出します。MRL = 1.0は完全な位相同期、MRL = 0は位相関係がランダムであることを示します。120秒のウィンドウ（約10呼吸周期）は、安定した位相推定に理想的とされる20周期以上を下回るため、重みは15%とし、主要な指標ではなく参考的な指標としています。注：位相コヒーレンスには、吸気と呼気の比が対称（1:1）である必要があります。アプリではユーザーの選択肢として呼気を長くする比率（1:1.5、1:2）にも対応していますが、そのセッションでは位相の測定を無効にしています。非対称な呼吸では心拍の位相に系統的なずれが生じ、MRLを解釈できなくなるためです。呼気を長くしたセッションは、スペクトル振幅のみでスコア化されます。

信号品質によるゲーティング

信号品質の指標（0〜1）は、アーチファクト率と心拍数の変動係数から算出します：

信号品質が0.50未満のブロックは、共鳴周波数の推定から除外されます。この基準値により、体動によるアーチファクト、電極の接触不良、異所性の拍動が周波数の推定をゆがめるのを防ぎます。除外されたブロックはセッション概要に印が付くため、患者はセンサーの装着位置の問題に早く気づくことができます。

共鳴周波数の推定

各探索セッションの後、アプリは過去のすべてのセッションの有効なブロックをすべて使って、患者の共鳴周波数を推定し直します：

• 重み付きローレンツ関数フィット：過去の有効なすべてのブロック（最低5ブロック、信号品質 ≥ 0.50）の（ペース、スコア）の組に、品質で重み付けしたローレンツ関数（減衰振動子）の曲線を当てはめます。ローレンツ関数は共鳴系にとって物理的に正しいスペクトル形状で、中心周波数と幅についてグリッドサーチでフィットし、各グリッド点で振幅とベースラインを解析的に求めます。品質の2乗による重み付けで、過去の記録をすべて保持しつつ、品質の高いブロックの影響を大きくします。ローレンツ関数のグリッドサーチで有効なピークが見つからない場合は、代わりに重み付き二次関数を使用します。
• IBI共変量モデル：十分なデータがある場合は、拡張モデルでローレンツ関数の中心周波数をIBIに依存させます。x₀(IBI) = x₀ + β · IBI_cとすることで、安静時の心拍間隔と最適な呼吸ペースの相関を考慮します。これはCapdevila et al.（2021）の知見と一致しています。これにより、患者の安静時心拍数が過去のベースラインと異なる場合にセッションの途中で適応でき、ブロック2以降で探索するペースがその日の患者の状態に合わせて調整されます。
• 代替手法：ローレンツ関数と二次関数のどちらのフィットも十分に制約されない場合は、品質で重み付けしたビン法を代替の推定手法として使用します。
• 信頼度：各推定値には、データ量、当てはまりの良さ、ピークのカバー状況に基づく信頼度スコア（0〜1）が付きます。これはセッション概要に表示されるため、患者は推定値が時間とともにどれだけ安定してきたかを確認できます。

患者向けセットアップガイド

次の4つのステップで、患者は最初のダウンロードから実際の練習まで進められます。