JMP® Student Editionについてご質問があります。
このソフトで、RRI(心拍間隔)のエクセルのデータからストレス指標のLF/HF(低周波成分/高周波成分)を知ることは可能かどうということと、もし可能であれば、その手順もご教示頂きたいです。
JMP® Student Edition について質問があります。
このソフトを使って、RRI(心拍間隔)のExcelデータからストレス指標LF/HF(低周波成分/高周波成分)を判定することは可能でしょうか? もし可能であれば、手順を教えてください。
SAS Institue Japan株式会社 JMPジャパン事業部でシニアテスターをしています小野と申します.
私自身,「HF/LF比」というものをご質問で初めて知った状況ですので,かなり曖昧な回答になります.
基本的には,HF/LF比の計算はJMPの標準機能で提供されておらず,かつ,プログラミングの関数でも用意されていない処理をしなければいけないので,JMPでやろうとするとかなり面倒だと思います.
まず,JMPもしくはJMP Pro(以降,まとめて"JMP"と述べます)でマウス操作もしくは1つの関数だけによって,RR間隔データからストレス指標であるHF/LF比を計算するようなものはありません.なお,JMP Student Editionは,JMP Proの機能を持ち合わせています.他にも例えばBMIの計算などもJMPではGUIで提供されていないのですが,そのような標準提供されていない計算をJMPで行おうとする場合,自分でプログラミングする必要があります.
(今回の「HF/LF比」計算を行うアドインはないのですが,JMP Marketplaceというサイト(https://marketplace.jmp.com/ja)に,「アドイン」という機能を用いた,現状有姿での提供の無料プログラムがいくつか公開されています.)
HF/LF比をJMPで計算するとしたら,次のような方法でプログラミングする(自分でコードを書く)必要があります.
- JSL (JMPスクリプト言語; JMP Scripting Language)というJMP専用の言語でプログラミングする.
- JMP内にて,(JMP 18以降の)Pythonのコードでプログラミングする.
RR間隔データからHF/LF比を計算する際には,いくつかの処理方法や計算方法があります.すでにどれがいいのか定番の方法もあるのかもしれませんが,ChatGptにて,以下のプロンプトで尋ねたところ,いくつかの論文が紹介されました.
「RR間隔からLF/HF比を求めるのに,(a) 等間隔データにするために線形補間やスプライン補完などをする,(b) 周波数領域でのパワーを求めるのに,そのままの値を用いるか,滑らかにするか,また,(c)最終的にパワースペクトルをどう数値積分するか,といったところでいろいろな手法が考えられます.これらを比較して検討した論文はありますか?」
以上のプロンプトへの論文の中で簡単な計算方法(例えば,線形補間により等間隔データを作り,生の高速Fourier変換を行い,台形法則で積分布する)は,JSLでもいくつか関数が用意されています(線形補間はInterpolate関数, 高速Fourier変換はFFTで行えます).しかし,その他の処理はJMPでは特に関数が用意されておらず,その部分を自分でスクラッチで記述する必要があります.どの計算方法が良いのかが分からない場合には,いずれの方法も試して結果が変わらないことを示した方がいいと思いますので,JMPのJSLで記述するのは手間がかかると思います.
以上のような計算する上で考慮すべき細かいこと以外にも,そもそもLF/HF比そのものが妥当な指標か推定対象にマッチしたものかどうかをまずは考慮すべきかもしれません.(時間領域だけでなく)周波数領域でも現実の心電図データを眺めてきた専門家の意見を聞くのが早いのではないかと思います(周波数領域でのデータの解釈がわからず,いずれかの何かしらの方法でLF/HF比を求めただけだと,計算が正しいかなども含めて不安が残るかもしれません.)
(なお,ARIMAモデルをあてはめてそのパワースペクトル密度を見るのであれば,「時系列分析」プラットフォームで行えます.ただし,このプラットフォームはLFパワー,HFパワー,LF/HF比などはまったく計算しません.)
前向きな良い答えじゃなくて申し訳ないのですが,何か不明な点などがありましたらご返信ください.
これは私の専門分野ではないので、データからLF/HF比を適切に測定するにはかなりの背景知識が必要だと思います。多くの尺度について説明している次のサイトを見つけました: https://pmc.ncbi.nlm.nih.gov/articles/PMC5624990/ 。また、この比率の使用に関する問題点を挙げている参考文献もいくつか見つけました。しかし、JMPの機能に関しては、これらの尺度はどれも簡単に構築でき、比較的単純なJMPの式ですべて実行できるようです。データが時系列データであることを考えると、単純なスプラインやより複雑な時系列モデルなど、いくつかの組み込み手法を使用でき、これらの結果はデータセット内の変数として簡単に保存できます。そのため、適切な計算とその解釈には多くの専門知識が必要であるように思われますが、JMPの実装は私にとって複雑ではありません。
したがって、ご質問への私の回答は、「はい、JMPを使えば心拍数データからLF/HFを簡単に測定できます」となります。ただし、手順を説明するには、具体的にどのような測定方法を用いるのか、もう少し背景を説明する必要があります。yusuke_onoさんの回答はJMPの計算が複雑に思えるかもしれませんが、それは私がLF/HFについて彼らよりも理解が浅いからかもしれません。
返信が遅くなり申し訳ございません。ご返信ありがとうございます。
マウス操作だけでHF/LF比を計算する方法はないと承知しております。プログラミングが必要で、「アドイン」という機能があることを教えていただき、ありがとうございます。
JMP では少し難しそうなので、他の方法も検討してみようと思います。