HOME PAGE

Power Point ed.

Spike Density Estimation

- Histogram Method

- Kernel Method

Share it with your friend

Tweet

定常リニューアル過程

前章では記憶なしの特徴をもつスパイク時系列，すなわちポアソン過程について調べた．ポアソン過程のスパイク間隔は指数分布になり，一定時間内のスパイク数はポアソン分布に従うことをみた．また瞬間スパイク生成率を導入して定常ポアソン過程を定義することもできた．本章ではスパイク生成に記憶がある場合のモデルを導入する．特に大事なのはスパイク生成率が直前のスパイク生成時刻だけに依存する場合でリニューアル過程と呼ばれる．瞬間生成率が２つ以上前のスパイク生成時刻に依存する一般的な点過程モデルに対して，直前のスパイクのみに依存することからこの名称で呼ばれる．

神経細胞のスパイク時系列は強いリニューアル性があることが知られている．これは神経細胞のスパイク生成機構の生理学上の制約による．特にスパイク生成後$5$ [ms]程度までの間は次のスパイクを生成できない状態になる．これを絶対不応期（Absolute refractory period）と呼ぶ．このためリニューアル過程を用いて神経スパイク発火を高い精度で記述することができる．

ここで本章で使う統計量をまとめておく．用語は基本的に生存時間解析(Survival analysis)に従うが，生存時間解析ではイベントが個体の死亡を意味するため関連語句がスパイクの解析にはなじまない．そこでいくつかは違う名称で呼ぶこともある．まずスパイク発生時刻の確率密度関数(probability density function, p.d.f.)を

$$f\left( x\right) =\lim_{\Delta\rightarrow0+}\frac{P\left( t<X\leq t+\Delta\right) }{\Delta}$$

とする．これをスパイク密度分布(spike density funtion)，もしくはスパイク間隔分布と呼ぶ Inter-spike Interval (ISI) distribution)と呼ぶ．スパイクの発生時刻が$t$ 以上である確率を生存関数(survival function)と呼び

$$\overline{F}\left( t\right) =P\left( X>t\right) =\int_{t}^{\infty}f\left( u\right) du%$$ (1.15)

と表す．生存関数と相補関係にある累積分布関数(cumulative distribution function, c.d.f.)

$$F\left( t\right) =P\left( X\leq t\right) =\int_{0}^{t}f\left( u\right) du$$ (1.16)

は生存時間分布(lifetime distribution function)と呼ばれる．スパイク密度分布と分布関数は

$$f\left( t\right) =F^{\prime}\left( t\right)$$ (1.17)

$$=-\overline{F}^{,}\left( t\right)$$ (1.18)

の関係にある．ただし記号プライムは微分の意味．

ハザード関数 (瞬間スパイク生成率)

リニューアル過程では瞬間スパイク生成率は，最後のスパイクの時刻を基準にして時間的に変動する．少なくとも$t$ 秒間スパイクが生成されず，$t$ においてスパイクを生成する確率として，瞬間スパイク生成率

$$r\left( t\right) =\lim_{\Delta\rightarrow0+}\frac{P\left( t<X\leq t+\Delta\vert t<X\right) }{\Delta}%$$

を導入しよう．点過程理論では条件付強度関数(conditional intensity function)という．生存時間解析(Survival analysis)ではハザード関数(hazard function)，年齢別故障率(age-specific failure rate)，回復関数(recovery function)などの名称でも呼ばれる[Cox, 1962]．

従ってハザード関数はスパイク密度分布(ISI分布)を用いて次のようにかける．

生存関数で表記すると

$$r\left( t\right) =-\frac{\overline{F}^{,}\left( t\right) }{\overl... ...ft( t\right) }=-\frac{d}{dt}\left\{ \log\overline{F}\left( t\right) \right\} .$$ (1.20)

両辺を0 から$t$ まで積分して，生存関数とハザード関数の関係式

$$\overline{F}\left( t\right) =\exp\left\{ -\int_{0}^{t}r\left(u\right) du\right\}$$ (1.21)

を得る．これから分布関数は

微分すると，スパイク密度分布(ISI分布)は

この式は確率密度分布とハザード関数の関係式1.20の逆関数になっている．1.20と1.24は併せて非常に重要な公式なので覚えておこう．定常ポアソン過程のスパイク密度分布が指数分布であることから，定常ポアソン過程のハザード関数は $r\left(t\right)=\lambda$ であることがわかる．

スパイク密度分布

それでは，代表的なスパイク密度分布(スパイク間隔分布，ISI分布)を見てみよう．

ガンマ分布

ガンマ分布はアーラン分布の拡張であることは既に述べた．リニューアル過程の中でも数学的な取り扱いが比較的容易で使い勝手が良い分布である．ガンマ分布のハザード関数を求めてみよう．

$$f\left( t;\lambda,\kappa\right) =\frac{\left( \kappa\lambda\right) ^{\kappa}}{\Gamma\left( \kappa\right) }t^{\kappa-1}e^{-\kappa\lambda t}%$$ (1.24)

ガンマ分布の平均は$1/\lambda$ ，分散は $1/\kappa\lambda$ ，また $CV=1/\sqrt{\kappa}$ ．ただし，平均が$1/\lambda$ になるように1.5を規格化してある．図にいくつかの$\kappa$ について分布を表示した．．

ゾ 1.3: 繧ャ繝ウ繝槫ッ?コヲ蛻?ク?
ゾ 1.4: 繧ャ繝ウ繝槫?蟶??繝上じ繝シ繝蛾未謨ー縲!-- MATH $\lambda=1，\kappa=0.5，1，2$ -->

ワイブル分布

ガンマ分布のハザード関数が定数に収束するのに対し，工業製品の故障率のようなイベントの場合は経年劣化するからハザード関数は時間とともに増加すると考えられる．時間とともにハザード関数がべき的に増加する場合を考えよう．

$$r\left( t;\lambda,\kappa\right) =\kappa\lambda\Gamma\left( 1+1/\kappa\right) \left\{ \Gamma\left( 1+1/\kappa\right) \lambda t\right\} ^{\kappa-1}$$

$$=\kappa\left\{ \Gamma\left( 1+1/\kappa\right) \lambda\right\} ^{\kappa }t^{\kappa-1}.$$

ただし$\kappa<1$ のときはハザード関数が減少関数になる．これは時間が経つにつれてスパイク発生率は低くなることを表す．べき乗のハザード関数を持つ分布はワイブル分布(Weibull distribution)とよばれ

$$f\left( t;\lambda,\kappa\right) =\kappa\lambda\Gamma\left( 1+1/\k... ...[ -\left\{ \Gamma\left( 1+1/\kappa\right) \lambda t\right\} ^{\kappa}\right]$$

で与えられる．この表記では，期待値は$1/\lambda$ ，分散は $1/\lambda^{2}\left[ \Gamma\left( 1+2/\kappa\right) /\Gamma\left( 1+1/\kappa\right) ^{2}-1\right]$ である． $\kappa=1$ のときにワイブル分布は指数分布になる．分布関数は

$$F\left( t;\lambda,\kappa\right) =1-\exp\left[ -\left\{ \Gamma\left( 1+1/\kappa\right) \lambda t\right\} ^{\kappa}\right].$$

ゾ 1.5: 繝ッ繧、繝悶Ν蟇?コヲ蛻?ク?$\lambda =1$ ， $\kappa =0.5,1,1,3$
ゾ 1.6: 繝ッ繧、繝悶Ν蛻?ク?未謨ー縲IMG WIDTH="43" HEIGHT="15" ALIGN="BOTTOM" BORDER="0" SRC="img3.png" ALT="$ \lambda =1$"> ， $\kappa =0.5,1,1,3$
ゾ 1.7: 繝ッ繧、繝悶Ν蛻?ク??繝上じ繝シ繝蛾未謨ー縲IMG WIDTH="43" HEIGHT="15" ALIGN="BOTTOM" BORDER="0" SRC="img3.png" ALT="$ \lambda =1$"> ， $\kappa =0.5,1,1,3$

逆ガウス分布

神経細胞スパイクの場合は死亡率と全く逆のハザード関数が描ける．神経細胞のスパイク活動のISIを良く表すとされる逆ガウス分布(Inverse Gaussian distribution)[]を見てみよう．

$$f\left( t;\lambda,\kappa\right) =\sqrt{\frac{\kappa}{2\pi t^{3}}}\exp\left[ -\frac{\kappa}{2t}\left( \lambda t-1\right) ^{2}\right]$$    for $$x>0$$，0 for $$x<0$$

逆ガウス分布の期待値は$1/\lambda$ ，分散は $1/\kappa \lambda^{3}$ ，また また累積分布関数は

$$F\left( t;\lambda,\kappa\right) =\frac{1}{2}\left\{ 1+\operatorna... ...}\left[ \sqrt {\frac{\kappa}{2t}}\left( \lambda t+1\right) \right] \right\}.$$

逆ガウス分布のハザード関数を見てみると増加から減少に転じている．スパイク発生率はゼロから始まり，時間とともに増加する．これはガンマ分布のハザード関数と同じ性質である．ガンマ分布と違うのはハザード関数が減少に転じる点にある．つまり時間が経つにつれてスパイク発生頻度は減少していく．

ゾ 1.8: 騾?ぎ繧ヲ繧ケ蟇?コヲ蛻?ク?平均$\lambda =1$ ．
ゾ 1.9: 騾?ぎ繧ヲ繧ケ蛻?ク??繝上じ繝シ繝蛾未謨ー縲均$\lambda =1$ ．

神経細胞のモデルとしてどのISI分布を用いたモデルが適しているかについて一般的なコンセンサスはないと思われる．神経細胞の種類やその細胞の存在する部位によって発火活動の特徴は大きく異なることが知られている．逆ガウス分布はブラウン運動のファーストパッセージタイムとして導出されるという背景があり，不応期をよく捉えるとする研究がある．ガンマ分布はバースト活動から規則発火まで表現することができ，数学的な取り扱いも用意であるという利点がある．