「積分発火モデル」の版間の差分

14行目: 14行目:
 神経細胞の電気的特性については、HodgkinとHuxleyによって細胞膜上に発現しているイオンチャネルの膜電位依存性とそれらによる活動電位生成機構、およびその数理的な表現が明らかにされたが<ref name=Hodgkin1952><pubmed>12991237</pubmed></ref>[1]、それ以前にLapicqueによって、細胞膜のキャパシタとしての特性や神経興奮現象(活動電位生成)に対する閾値となる電位、および、閾値に到るまでの過程について詳細に調べられていた<ref name=Lapicque1907>'''Lapicque, L. (1907).'''<br>Recherches quantitatives sur l'excitation électrique des nerfs traitée comme une polarization. Journal de physiologie et de pathologie générale, 9, 620-635.</ref><ref name=Lapicque2007><pubmed>18046573</pubmed></ref> [2,3]。入力として与えられた電流はキャパシタとしての特性により膜電位に積算(integrate)される。膜電位が上昇して閾値に到達すると、活動電位を発生(fire)する。この膜電位の閾値に達するまでの積算過程をモデル化したのが、積分発火モデルである。このモデルでは活動電位生成中の膜電位変動は記述しない。それは非常に短時間(< 2msec)の過程であり、膜電位挙動のほとんどの時間が閾値に到るまでの積算過程であるとみなせるためである。神経細胞の状態を表す変数が膜電位のみの1変数であるため、計算量も多くない。このため、多くの研究において、採用されてきたモデルである。一方で、神経細胞応答の本質であるアクティブな膜伝導性を一切無視したモデルであるため、実際の神経細胞応答とは異なる特性を示すことから、その拡張モデルも多く提案されてきた。以下では、まず、基本となる積分発火モデルであるLeaky Integrate-and-Fireモデル(LIFモデル、と呼ばれる)を解説し、次にその拡張モデルとして代表的なものを紹介する。
 神経細胞の電気的特性については、HodgkinとHuxleyによって細胞膜上に発現しているイオンチャネルの膜電位依存性とそれらによる活動電位生成機構、およびその数理的な表現が明らかにされたが<ref name=Hodgkin1952><pubmed>12991237</pubmed></ref>[1]、それ以前にLapicqueによって、細胞膜のキャパシタとしての特性や神経興奮現象(活動電位生成)に対する閾値となる電位、および、閾値に到るまでの過程について詳細に調べられていた<ref name=Lapicque1907>'''Lapicque, L. (1907).'''<br>Recherches quantitatives sur l'excitation électrique des nerfs traitée comme une polarization. Journal de physiologie et de pathologie générale, 9, 620-635.</ref><ref name=Lapicque2007><pubmed>18046573</pubmed></ref> [2,3]。入力として与えられた電流はキャパシタとしての特性により膜電位に積算(integrate)される。膜電位が上昇して閾値に到達すると、活動電位を発生(fire)する。この膜電位の閾値に達するまでの積算過程をモデル化したのが、積分発火モデルである。このモデルでは活動電位生成中の膜電位変動は記述しない。それは非常に短時間(< 2msec)の過程であり、膜電位挙動のほとんどの時間が閾値に到るまでの積算過程であるとみなせるためである。神経細胞の状態を表す変数が膜電位のみの1変数であるため、計算量も多くない。このため、多くの研究において、採用されてきたモデルである。一方で、神経細胞応答の本質であるアクティブな膜伝導性を一切無視したモデルであるため、実際の神経細胞応答とは異なる特性を示すことから、その拡張モデルも多く提案されてきた。以下では、まず、基本となる積分発火モデルであるLeaky Integrate-and-Fireモデル(LIFモデル、と呼ばれる)を解説し、次にその拡張モデルとして代表的なものを紹介する。


[[ファイル:Kitano 積分発火モデルFig1.png|サムネイル|450px|'''図1. 積分発火モデルとMATモデル(拡張された積分発火モデルを参照)の模式図'''
[[ファイル:Kitano 積分発火モデルFig1.png|サムネイル|450px|'''図1. 積分発火モデルとMulti-timescale Adaptive Threshold (MAT)モデル(拡張された積分発火モデルを参照)の模式図'''
どちらのモデルも、膜電位 (黒) が閾値<math>V_{th}</math>(マゼンタ) に到達するとスパイクを生成するが、スパイク生成時のリセットルールが異なる。<br>
どちらのモデルも、膜電位 (黒) が閾値<math>V_{th}</math>(マゼンタ) に到達するとスパイクを生成するが、スパイク生成時のリセットルールが異なる。<br>
'''A.''' 積分発火モデル:スパイクを生成すると、膜電位をリセット電位<math>V_{reset}</math>(シアン) にリセットする。パラメータは<math>\tau_m=10\mbox{ }[ms]</math>, <math>R_m=10\mbox{ }[M\Omega]</math>, <math>E_L=-65\mbox{ }[mV]</math>, <math>V_{th}= -50\mbox{ }[mV]</math>, <math>V_{rest}= -60\mbox{ }[mV]</math>である。<br>
'''A.''' 積分発火モデル:スパイクを生成すると、膜電位をリセット電位<math>V_{reset}</math>(シアン) にリセットする。パラメータは<math>\tau_m=10\mbox{ }[ms]</math>, <math>R_m=10\mbox{ }[M\Omega]</math>, <math>E_L=-65\mbox{ }[mV]</math>, <math>V_{th}= -50\mbox{ }[mV]</math>, <math>V_{rest}= -60\mbox{ }[mV]</math>である。<br>
'''B.''' MATモデル:スパイクを生成すると、膜電位をリセットせずに閾値 (マゼンタ) を変動させる。閾値パラメータは<math>\omega=10\mbox{ }[mV]</math>, <math>\alpha_1=4\mbox{ }[mV]</math>, <math>\alpha_2=3\mbox{ }[mV]</math>である。その他 のパラメータは<math>\tau_m=5\mbox{ }[ms]</math>, <math>R_m=10\mbox{ }[M\Omega]</math>, <math>E_L=0\mbox{ }[mV]</math>である。]]
'''B.''' Multi-timescale Adaptive Thresholdモデル:スパイクを生成すると、膜電位をリセットせずに閾値 (マゼンタ) を変動させる。閾値パラメータは<math>\omega=10\mbox{ }[mV]</math>, <math>\alpha_1=4\mbox{ }[mV]</math>, <math>\alpha_2=3\mbox{ }[mV]</math>である。その他 のパラメータは<math>\tau_m=5\mbox{ }[ms]</math>, <math>R_m=10\mbox{ }[M\Omega]</math>, <math>E_L=0\mbox{ }[mV]</math>である。]]
 
== 積分発火モデル ==
== 積分発火モデル ==
Integrate and fire model
Integrate and fire model