Power Point ed.

Spike Density Estimation

- Histogram Method

- Kernel Method

Share it with your friend

Tweet

Next: スパイク時系列モデルの推定（暫定版） Up: スパイク時系列モデルの実現 Previous: 一般点過程の時間伸縮理論と尤度関数目次索引

Subsections

非定常リニューアル過程

変動レートの時間伸縮

非定常リニューアル過程は変動ポアソン過程の場合を拡張して考えるとわかりやすい．非定常ポアソン過程と同様に時間伸縮した軸上で定常リニューアル過程を考える．例として $\Lambda$ 上の定常ガンマ過程を考え，そのイベント系列を $\left\{ \Lambda_{1},\Lambda_{2},\cdots\right\}$ とする．イベント間隔 $\Lambda_{i}-\Lambda_{i-1}$ が平均が，分散が $1/\kappa$ である標準ガンマ分布

$\displaystyle g\left( z;\kappa\right) =\frac{1}{\Gamma\left( \kappa\right) }% \kappa\left( \kappa z\right) ^{\kappa-1}e^{-\kappa z}%$

に従うとする[Berman, 1981,Barbieri et al., 2001]． $\kappa$ はshape parameterと呼ばれる．実時間との対応は $\Lambda_{i}=\Lambda\left( t_{i}\right)$ だから，ガンマ分布に従う乱数を発生させ，変数変換の公式

$\displaystyle z\left( t_{i}\vert t_{i-1},\lambda_{t_{i-1}:t_{i}}\right) =\Lambd... ...Lambda\left( t_{i-1}\right) =\int_{t_{i-1}}^{t_{i}} \lambda\left( u\right) du$

を満たす $t_{i}$ を逐次的に求めれば，実時間でのスパイク時系列を求めることができる．

条件付きスパイク密度分布

ガンマ分布

条件付きスパイク密度分布 $p\left( t_{i}\vert t_{i-1},\left\{ \lambda _{t}\right\} \right)$ を求めよう．分布 $p\left( t_{i}\vert t_{i-1},\left\{ \lambda _{t}\right\} \right)$ は次のような $t_{i}$ からへの変数変換（タイムリスケーリング）を施すことで，標準ガンマ分布になるとする．すなわち，変数変換の公式から

$\displaystyle g\left( z\vert t_{i-1},\lambda_{t_{i-1}:t_{i}}\right) =\left\vert... ...dt}\right\vert ^{-1}p\left( t_{i}\vert t_{i-1},\lambda_{t_{i-1}:t_{i}}\right)$

$dz/dt_{i}=\lambda\left( t_{i}\right)$ に注意して

$\displaystyle p\left( t_{i}\vert t_{i-1},\lambda_{t_{i-1}:t_{i}}\right)$	$\displaystyle =\left\vert \frac{dz}{dt_{i}}\right\vert g\left( z\vert t_{i-1},\lambda_{t_{i-1}:t_{i}}\right)$
	$\displaystyle =\lambda\left( t_{i}\right) \frac{1}{\Gamma\left( \kappa\right) }... ...ppa-1}\exp\left[ -\kappa\int_{t_{i-1}}^{t_{i}}\lambda\left( u\right) du\right]$

となる．

ワイブル分布

平均・分散 $\Gamma\left( 1+2/\kappa\right) /\Gamma\left( 1+1/\kappa\right) ^{2}-1$ の標準ワイブル分布は

$\displaystyle g\left( z;\kappa\right) =\kappa\Gamma\left( 1+1/\kappa\right) \le... ...exp\left[ -\left\{ \Gamma\left( 1+1/\kappa\right) t\right\} ^{\kappa}\right]$

$dz/dt_{i}=\lambda\left( t_{i}\right)$ に注意して

$\displaystyle p\left( t_{i}\vert t_{i-1},\lambda_{t_{i-1}:t_{i}}\right)$	$\displaystyle =\left\vert \frac{dz}{dt_{i}}\right\vert g\left( z\vert t_{i-1},\lambda_{t_{i-1}:t_{i}}\right)$
	$\displaystyle =\lambda\left( t_{i}\right) \kappa\Gamma\left( 1+1/\kappa\right) ... ...ppa\right) \int_{t_{i-1}}^{t_{i}}\lambda\left( u\right) du\right\} ^{\kappa-1}$
	$\displaystyle \cdot\exp\left[ -\left\{ \Gamma\left( 1+1/\kappa\right) \int_{t_{i-1}% }^{t_{i}}\lambda\left( u\right) du\right\} ^{\kappa}\right]$

となる．

逆ガウス分布

同様にして，逆ガウス分布に従うリニューアル過程は，平均・分散 $1/\kappa$ の標準逆ガウス分布

$\displaystyle g\left( z;\kappa\right) =\sqrt{\frac{\kappa}{2\pi z^{3}}}\exp\left[ -\frac{\kappa}{2z}\left( z-1\right) ^{2}\right]$

を用いてISI分布は

$\displaystyle p\left( t_{i}\vert t_{i-1},\lambda_{t_{i-1}:t_{i}}\right)$	$\displaystyle =\left\vert \frac{dz}{dt_{i}}\right\vert g\left( z\vert t_{i-1},\lambda_{t_{i-1}:t_{i}}\right)$
	$\displaystyle =\lambda\left( t_{i}\right) \sqrt{\frac{\kappa}{2\pi\left\{ \int... ...ht) du-1\right\} ^{2}}{\int_{t_{i-1}}^{t_{i}}\lambda\left( u\right) du}\right]$

となる．

時間伸縮理論

条件付きスパイク密度分布が求まれば，これから瞬間スパイク生成率(ハザード関数)が求まる．

$\displaystyle r\left( t_{i}\vert t_{i-1},\lambda_{t_{i-1}:t_{i}}\right) =\frac{... ...t_{t_{i-1}}^{t_{i}}p\left( u\vert t_{i-1},\lambda_{t_{i-1}:t_{i}}\right) du}%$

瞬間スパイク生成率が求まれば，時間を十分小さな幅 $\Delta$ のビンに区切って確率 $r\left( t_{i}\vert t_{i-1},\lambda_{t_{i-1}:t_{i}}\right) \Delta$ でスパイクを発生させればよい．もちろんこの方法は効率が悪い．

瞬間性スパイク生成率 $r\left( t_{i}\vert t_{i-1},\lambda_{t_{i-1}:t_{i}}\right)$ が求まれば，非定常リニューアル過程であっても非定常ポアソン過程としてあつかえる．非定常ポアソン過程の時と同様に $r\left( t_{i}\vert t_{i-1},\lambda_{t_{i-1}:t_{i}}\right) \Delta$ を積分すると

$\displaystyle \zeta=\int_{t_{i-1}}^{t_{i}}r\left( u\vert t_{i-1},\lambda_{t_{i-1}:t_{i}}\right) du$

$t_{i}$ が $p\left( t_{i}\vert t_{i-1},\lambda_{t_{i-1}:t_{i}}\right)$ に従うとき， $\zeta$ は指数分布に従う．したがって指数分布に従う乱数 $\zeta$ を生成し，数値計算で上式を満たす $t_{i}$ を求めればよい．時間伸縮理論を用いたこの方法は非常に汎用性のある方法で，より一般の点過程について成り立つから証明はそこで示す．ここでは時間伸縮理論に基づいた数値計算アルゴリズムをまとめておく．

アルゴリズム 17 時間伸縮理論を用いた非定常リニューアル過程の実現

$i=1,t_{1}=0.$
レートの指数分布に従う乱数 $\eta$ を生成する．
時間伸縮 $\Lambda_{t}=\int_{0}^{t}\lambda\left( u\right) du$ により，伸縮したISIを求める． $z_{t}=\Lambda_{t}-\Lambda_{t_{i-1}}$
ISI分布を求める $p\left( t\vert t_{i-1},\lambda_{t_{i-1}:t_{i}}\right) =\lambda\left( t\right) g\left( z_{t}\vert t_{i-1},\lambda_{t_{i-1}:t_{i}}\right)$
ハザード関数を求める
$r\left( t\vert t_{i-1},\lambda_{t_{i-1}:t}\right) =\frac{p\left( t\vert t_{i-... ...t) }{1-\int_{t_{i-1}}^{t}p\left( u\vert t_{i-1},\lambda_{t_{i-1}:u}\right) du}$
$\eta=\int_{t_{i-1}}^{t}r\left( u\vert t_{i-1}\right) du$ を満たす $t_{i}$ を求める．
$t_{i}\geq T$ ならば終了．
$t_{i}<T$ ならば $i\leftarrow i+1.\ 2$ へ．