今後の課題

本方法の問題点は，解析に用いるデータのサンプル数がある一定数以上必要であることと，

解析の際に各信号をバンドパスフィルタに通し狭帯域信号にしていることである．データのサ ンプル数が少なすぎると，確率密度関数の推定精度が低くなり，また，バンドパスフィルタの 通過域を広くしすぎると様々な情報を含んでしまい，通過域を狭くしすぎると確定信号に近づ いてしまうので精度よく解析できなくなると考えられる．これらの問題点については，データ のサンプル数やバンドパスフィルタの通過域を変化させて同様の解析を行い，より適切な条件 の検討を行うべきであると考えられる．

本研究では，血圧調整システム内における3つの信号間における因果関係を調べたが，先に 述べたように血圧調整システムは非常に多くの要素が複雑に絡み合っている．したがって，ITE は4変数以上にも拡張できるので，変数を拡張したITEを用いてより多くの信号間における因 果関係を調べることで，血圧調整システムの機序をより詳しく解明することができると考えら れる．また，日常生活中において非拘束で測定できる心電図や脈波などの生体信号を長時間モ ニタリングし，得られた複数の信号をITEを用いて解析することで，血圧調整システム内の長 期的な変化を評価できる可能性があると考えられる．

参考文献

[1] 二宮 石雄：循環系の制御機構． 計測と制御．38(10)：33–38，1985．

[2] 吉澤 誠，田中 明，阿部 健一，竹田 宏，山家 智之，仁田 新一，阿部 裕輔，井街 宏：人 工心臓の制御．計測と制御．38(5)：328–333，1999．

[3] 檜垣 實男：超高齢者高血圧の血圧管理の意義と方策．日本老年医学会雑誌．43(6)：697–699， 2006．

[4] Takagi S, Saitoh S, Nakano M, Hayashi Y, Obara F, Onishi H, Shimamoto K : Rela-tionship between blood pressure level and mortality rate: an 18-year study conducted in two rural communities in Japan. Hypertension.18(2): 139–144, 2000.

[5] 大西 浩文，齋藤 重幸，島本 和明：端野・壮瞥町研究レビュー2007．Therapeutic Research 2007． 28(4)：513–525，2007．

[6] 本郷 利憲，廣重 力，豊田 順一：標準生理学，第6版．医学書院，p. 617，2006． [7] Guo GB, Abboud FM : Impaired central mediation of the arterial baroreflex in chronic

renal hypertension. Am J Physiol Heart Circ Physiol. 246: 720–727, 1984.

[8] Grassi G, Cattaneo BM, Seravalle G, Lanfranchi A, Mancia G : Baroreflex control of sympathetic nerve activity in essential and secondary hypertension. Hypertension. 31:

68–72, 1998.

[9] 本郷 利憲，廣重 力，豊田 順一：標準生理学，第6版．医学書院，pp. 407–413，p. 519， p. 559，2006．

[10] Eppinger H, Hess L : Zur Pathologie des vegetativen Nervensystems. Zeitschrift fur Klinische Medizin.67: 345–351, 1909.

[11] Guyton AC, Hall JE : Textbook of Medical Physiology 11 Edition, http://www.studentconsult.com/content/default.cfm?ISBN=0721602401

[12] 小澤 瀞司，福田 康一郎：標準生理学，第7版．医学書院，pp. 624–626，2009． [13] Mulder G : The heart of mental effort.,Thesis. University of Groningen. 1980.

[14] Korner PI, Tonkin AM, Uther JB : Valsalva constrictor and heart rate reflexes in subjects with essential hypertension and with normal blood pressure. Clinical and Experimental Pharmacology & Physiology.6(1): 97–110, 1979.

[15] 高田 重男：若年本態性高血圧患者の圧受容体反射を介する心拍数調節に関する研究．金 沢大学十全医学会雑誌．93(4)：655–667，1984．

[16] Parati G, Mancia G : The neck chamber technique. Giornale Italiano di Cardiologia.

22(4): 511–516, 1992.

[17] Pagani M, Somers V, Furlan R, Dell’Orto S, Conway J, Baselli G, Cerutti S, Sleight P, Malliani A : Changes in autonomic regulation induced by physical training in mild hypertension. Hypertension.12(6): 600-610, 1988.

[18] Persson PB, DiRienzo M, Castiglioni P, Cerutti C, Pagani M, Honzikova N, Akselrod S, Parati G : Time versus frequency domain techniques for assessing baroreflex sensitivity.

Hypertension.19(10): 1699-1705, 2001.

[19] 杉田 典大，吉澤 誠，田中 明，阿部 健一，山家 智之，仁田新一：乗り物酔いしやすさ と3次元映像感受性の関連性の定量的評価．日本バーチャルリアリティ学会第６回大会論 文集．2001，pp. 121–124．

[20] Yoshizawa M, Sugita N, Tanaka A, Abe K, Yambe T, Nitta S : Quantitative and physio-logical evaluation of three dimensional images. Proceedings of the Seventh International Conference on Virtual Systems and Multimedia, Berkeley, CA, USA, 2001, pp. 864–871.

[21] 杉田 典大，吉澤 誠，田中 明，阿部 健一，山家 智之，仁田 新一：血圧-心拍数間の最大相 互相関係数を用いた映像刺激の生態影響評価．ヒューマンインタフェース学論誌．4(4)： 227–234，2002．

[22] 阪田 治，椎名 毅，斎藤 陽一：多次元有向情報量とその脳波因果性解析への応用．電子情 報通信学会論文誌A．J83–A(5)：466–475，2000．

[23] Tanaka N, Okamoto H, Naito M : Detecting and evaluating intrinsic nonlinearity present in the mutual dependence between two variables. Physica D. Nonlinear Phenomena.147:

1–11, 2000.

[24] Schreiber T : Measuring information transfer. Physical Review Letters.85(2): 623–639, 2000.

[25] 是永 智宏，田中 尚樹，阿部 正紀：相互情報量および移動エントロピーに基づく情報移動 解析．電子情報通信学会総合大会講演論文集，2004–03–08．2004，p. 53．

[26] Katura T, Tanaka N, Obata A, Sato H, Maki A : Quantitative evaluation of interrelations between spontaneous low-frequency oscillations in cerebral hemodynamics and systemic cardiovascular dynamics. NeuroImage. 31(4): 1592–1600, 2006.

[27] Silverman BW : Density Estimation for Statistics and Data Analysis. Chapman and Hall, London, 1986.

[28] Jones MC, Marron JS, Sheather SJ : A brief survey of bandwidth selection for density estimation. Journal of the American Statistical Association.91(433): 401–407, 1996.

[29] Terrell GR : The maximal smoothing principle in density estimation. Journal of the American Statistical Association.85(410): 470–477, 1990.

[30] 谷崎 久志：密度関数のカーネル推定量におけるバンド幅の選択について：モンテカルロ 実験による小標本特性．国民経済雑誌．191(1)：59–70，2005．

[31] Seidel H, Herzel H : Bifurcations in a nonlinear model of the baroreceptor-cardiac reflex.

Physica D.115: 145–160, 1998.

[32] Seidel H, Herzel H : Modelling heart rate variability due to respiration and baroreflex.

in: Mosekilde E, Mouritsen OG. ed, Modelling the Dynamics of Biological Systems.

Springer, Berlin, 1995, pp. 205–229.

[33] Reiner VS, Antzelevitch C : Phase resetting and annihilation in a mathemetical model of sinus node. American Journal of Physiology.249(6): 1143–1153, 1985.

[34] Cavalcanti S, Belardinelli E : Modeling of Cordiovascular Variability Using a Differential Delay Equation. IEEE Transactions on Biomedical Engineering.43(10): 982–989, 1996.

[35] Press WH, Teukolsky SA, Vetterling WT, Flannery BP : Numerical Recipes in C. Cam-bridge University Press, CamCam-bridge, UK, 1992.

[36] Parmer RJ, Cervenka JH, Stone RA : Baroreflex sensitivity and heredity in essential hypertension. Circulation.85(2): 497–503, 1992.

[37] 浅井 宏祐: 自律神経機能検査，第4版．文光堂，pp. 159–195，2007．

[38] 小澤 利男：脈圧測定の臨床，Arterial stiffness動脈壁の硬化と老化．臨床動脈波研究会，

no. 8，2005．

[39] 大橋 俊夫：血液循環，標準生理学，第6版．医学書院，pp. 561-595，2006．

[40] Gribbin B, Steptoe A, Sleight P : Pulse wave velocity as a measure of blood pressure change. The Society for Psychophysiological Research.13(1): 86–90, 1976.

[41] Kazanavicius E, Gircys R, Vrubliauskas A, Lugin S : Mathematical methods for deter-mining the foot point of the arterial pulse wave and evaluation of proposed methods.

Information Technology And Control.34(1): 29–36, 2005.

[42] Chen W, Kobayashi T, Ichikawa S, Takeuchi Y, Togawa T : Continuous estimation of systolic blood pressure using the pulse arrival time and intermittent calibration. Medical and Biological Engineering and Computing.38(5): 569–574, 2000.

謝辞

本論文を結ぶにあたり，終始適切な御指導，御助言，御鞭撻を賜るとともに，学生生活にお きましても一方ならぬご配慮を頂きました東北大学サイバーサイエンスセンター 先端情報技 術研究部 吉澤 誠 教授に深く御礼申し上げます．

本論文の審査に際し，本研究に関して種々の観点から大変貴重な御意見を賜りました東北大 学加齢医学研究所 心臓病電子医学分野 山家 智之 教授，東北大学大学院 医工学研究科医工学 専攻 西條 芳文 教授に深く御礼申し上げます．

研究活動に対する適切な御指導，御助言，御鞭撻を賜り，また研究室におきましても御配慮 を賜りました東北大学サイバーサイエンスセンター 先端情報技術研究部 本間 経康 准教授に 心より御礼申し上げます．

本研究を遂行するにあたり，適切な御指導，貴重な御意見を賜り，日頃から御討論を頂き，

また実験等にも協力して頂きました東北大学大学院 工学研究科技術社会システム専攻 杉田 典 大 准教授に心より御礼申し上げます．

本研究に関して，日頃から貴重な御意見，御討論を頂き，実験等にも協力して頂きました東 北大学大学院 工学研究科電気・通信専攻 阿部 誠 助教に心から感謝致します．

事務処理などをはじめとして，研究生活上の様々な御世話をして頂きました曽根 知佳子 秘 書に心から感謝致します．

同じ分野の研究に携わり，日頃から貴重な御意見，御討論を頂き，また実験等に協力して頂 きました狩野 佑介氏（東北大学大学院 工学研究科電気・通信工学専攻博士課程前期2年），新 沼 大樹氏（東北大学大学院 工学研究科電気・通信工学専攻博士課程前期1年），松岡 成己氏

（東北大学大学院 工学研究科電気・通信工学専攻博士課程前期1年），荻原 健氏（東北大学工 学部 情報知能システム総合学科4年）に心から感謝致します．

また，本研究を行なうにあたり，貴重な時間を割いて実験に快く参加して下さった方々に心 より御礼申し上げます．

このほか吉澤・本間研究室の皆様には，学生生活全般にわたり日頃から大変お世話になりま した．心から感謝の意を表します．