Rで時系列のパワースペクトルを推定したいとき，もっともお気楽なのは，"spectrum"というコマンドを使うことだと思います．今回は，心拍変動 (RR間隔時系列)のスペクトル解析で"spectrum"コマンドを使った実例を紹介します．ちなみに，例として使用したRR間…

心拍変動を研究や商品開発で使いたいとき，まず，参照すべき論文があります．それが，以下のリンクにあります．https://www.ahajournals.org/doi/10.1161/01.CIR.93.5.1043 この論文のタイトルは， Heart Rate Variability Standards of Measurement, Physio…

心拍変動解析では，パワースペクトルの推定を使った心拍変動指標がよく使われます．HFパワーとか，LFパワーとか呼ばれるやつです．今回は，これら指標について説明しません．心拍変動指標を計算する前に，絶対に理解しておかなければならないことがあるから…

今回は実際の時系列データを読み込んでグラフを描いでみます． 濃厚接触者として自宅待機中（現在）の私の心拍変動 (RRI間隔)時系列をサンプルとして使ってください．以下のリンクからダウンロードできます． https://drive.google.com/file/d/1i0VIooI1-IKA…

今回は，何で「1/f（えふぶんのいち）ゆらぎ」とか，1/fノイズに注目するのか，何が不思議なのか，何がおもしろいのか，というお話です．1/fゆらぎについて，歴史と最新の成果を正確に記述することは，私には難しいので，まずは，現状の私の知識，印象，思い…

代表的な心拍変動指標（特に時間領域指標）の名前には，NNが入っていることが多いです．例えば，meanNN，SDNN，pNN50などです． meanNN すべてのNN間隔の平均値 SDNN すべてのNN間隔の標準偏差 pNN50 隣り合うNN間隔の差が50 msより長い拍の割合（％） 心臓…