Poisson過程(Poisson process)とは、時刻 $t$ までに起こった出来事の回数 $N(t)$ を表す計数過程で、重ならない時間区間での回数が独立であり、長さ $h$ の区間での回数が平均 $\lambda h$ の Poisson 分布に従うものである。独立に同じ指数分布に従う待ち時間を順に並べると Poisson 過程が得られ、逆にどの Poisson 過程でも出来事の間隔は独立に同じ指数分布に従う。増分が独立で定常であり、短い時間 $h$ に 1 回起こる確率が $\lambda h+o(h)$、2 回以上起こる確率が $o(h)$ であるという局所的な条件からも特徴づけられ、ばらばらに起こる出来事の回数の基本の模型になる。
Poisson 過程は、時間の中でばらばらに起こる出来事(電話の着信、放射性原子の崩壊、店への来客など)の回数を、時刻 $t$ までの回数 $N(t)$ として表す確率過程の基本の模型である。重ならない時間区間での回数が独立で、長さ $h$ の区間での回数が平均 $\lambda h$ の Poisson分布 に従う、という 2 つの条件で定める。
この記事では、Poisson 過程を定義し、次の 3 つを証明する。独立に同じ指数分布に従う待ち時間を並べると Poisson 過程ができること(主定理 1)、逆にどの Poisson 過程でも出来事の間隔は独立に同じ指数分布に従うこと(主定理 2)、そして「短い時間に 1 回起こる確率は長さに比例し、2 回以上起こる確率は無視できる」という局所的な条件だけから Poisson 分布が出てくること(主定理 3)である。最後に、点に独立に色を塗ると色ごとに独立な Poisson 過程に分かれること(間引き)を示す。
以下、確率空間 $(\Omega,\mathcal{F},P)$ を固定し、$\mathbb{Z}_{\ge0}=\{0,1,2,\dots\}$ とする。$\mathrm{Po}(\mu)$ は平均 $\mu$ の Poisson 分布である(Poisson分布 の定義「Poisson 分布」)。
$\mathbb{Z}_{\ge0}$ に値をとる確率変数の族 $(N(t))_{t\ge0}$ で、すべての $\omega\in\Omega$ について次を満たすものを計数過程という。
$N(t)$ は「時刻 $t$ までに起こった出来事の回数」、増分 $N(t)-N(s)$ は「区間 $(s,t]$ で起こった回数」と読む。
$\lambda>0$ とする。計数過程 $(N(t))_{t\ge0}$ が次の 2 条件を満たすとき、強度(rate)$\lambda$ の Poisson 過程(Poisson process)という。
これは Dur19 §3.7(p. 175)の定義に、計数過程であること($N(0)=0$、単調、右連続)を加えたものである。増分の分布が区間の長さだけで決まるので、Poisson 過程の増分は定常である。$N(t)$ 自身は $t=0$ からの増分なので $N(t)\sim\mathrm{Po}(\lambda t)$、平均は $E[N(t)]=\lambda t$ であり、$\lambda$ は単位時間あたりの平均の回数を表す。
出来事の時刻と間隔を、計数過程から次のように取り出す。
計数過程 $N$ について、$k\ge1$ に対し
$$
S_k:=\inf\{t\ge0\mid N(t)\ge k\}\in[0,\infty]
$$
(集合が空なら $S_k=\infty$)を第 $k$ の到着時刻、$S_0:=0$ とし、$T_k:=S_k-S_{k-1}$($S_{k-1}< \infty$ のとき)を第 $k$ の到着間隔という。
計数過程 $N$ と $k\ge1$、$t\ge0$ について、$\{S_k\le t\}=\{N(t)\ge k\}$ である。とくに $S_k$ は $[0,\infty]$ に値をとる確率変数であり、$S_1\le S_2\le\cdots$ である。
要点:$N$ の右連続性により、$N(u)\ge k$ となる $u$ の集合は $S_k$ 自身を含む。
$\omega$ を固定し、$A:=\{u\ge0\mid N(u,\omega)\ge k\}$ とおく。$N(t)\ge k$ なら $t\in A$ なので $S_k=\inf A\le t$ である。逆に $S_k\le t$ とする。このとき $A$ は空でない。$N$ は単調なので、$u\in A$ かつ $u'\ge u$ なら $u'\in A$ である。$u>S_k$ なら、下限の定義により $u$ より小さい $A$ の元があるので $u\in A$ である。右連続性により $N(S_k)=\lim_{u\downarrow S_k}N(u)\ge k$ なので $S_k\in A$、よって $t\ge S_k$ も $A$ に入り、$N(t)\ge k$ である。
$\{N(t)\ge k\}\in\mathcal{F}$ なので、すべての $t$ で $\{S_k\le t\}\in\mathcal{F}$ であり、$S_k$ は確率変数である。$N(t)\ge k+1$ なら $N(t)\ge k$ なので $\{S_{k+1}\le t\}\subset\{S_k\le t\}$ がすべての $t$ で成り立ち、$S_k\le S_{k+1}$ である。$\square$
最後に、待ち時間の分布を定める。$\lambda>0$ について、$P(T>t)=e^{-\lambda t}$($t\ge0$)を満たす確率変数 $T$ の分布を、強度 $\lambda$ の指数分布 $\mathrm{Exp}(\lambda)$ という。これは $(0,\infty)$ 上の密度 $\lambda e^{-\lambda t}$ をもち、平均は $1/\lambda$ である(高校数学の範囲での導入と無記憶性は 一様分布と指数分布 にある)。
構成の鍵は、到着時刻の組の密度が、最後の時刻だけで決まる簡単な形をしていることである。
$T_1,T_2,\dots$ を独立で、どれも $\mathrm{Exp}(\lambda)$ に従う確率変数とし、$S_n:=T_1+\cdots+T_n$ とおく。$n\ge0$ について $(S_1,\dots,S_{n+1})$ は $\mathbb{R}^{n+1}$ 上の密度
$$
g(s_1,\dots,s_{n+1})=\begin{cases}\lambda^{n+1}e^{-\lambda s_{n+1}} & (0< s_1< s_2<\cdots< s_{n+1})\\ 0 & (\text{それ以外})\end{cases}
$$
をもつ。
$T_1,\dots,T_{n+1}$ は独立なので、$(T_1,\dots,T_{n+1})$ の分布はそれぞれの分布の積測度であり、Fubiniの定理 により密度 $\prod_{i=1}^{n+1}\lambda e^{-\lambda t_i}$($t_i>0$)をもつ。線形写像 $\Phi(t_1,\dots,t_{n+1}):=(t_1,\ t_1+t_2,\ \dots,\ t_1+\cdots+t_{n+1})$ は下三角で対角成分がすべて $1$ の行列で表されるので行列式は $1$ であり、$(0,\infty)^{n+1}$ を $\{0< s_1<\cdots< s_{n+1}\}$ の上へ全単射に写す。行列式 $1$ の線形写像は Lebesgue 測度を保つので、$(S_1,\dots,S_{n+1})=\Phi(T_1,\dots,T_{n+1})$ の密度は $\Phi^{-1}(s)$ での元の密度の値、すなわち $\lambda^{n+1}e^{-\lambda(t_1+\cdots+t_{n+1})}=\lambda^{n+1}e^{-\lambda s_{n+1}}$ である。$\square$
$\lambda>0$ とし、$T_1,T_2,\dots$ を独立で、どれも $\mathrm{Exp}(\lambda)$ に従う確率変数とする。$S_n:=T_1+\cdots+T_n$ とおき、
$$
N(t):=\#\{n\ge1\mid S_n\le t\}\qquad(t\ge0)
$$
と定める(ただし、ある $T_k$ が $0$ 以下であるか、$S_n$ が $\infty$ に発散しないという確率 $0$ の事象の上では $N\equiv0$ とおく)。このとき $N$ は強度 $\lambda$ の Poisson 過程であり、その到着時刻は $S_n$、到着間隔は $T_n$ である(確率 $0$ の事象を除く)。さらに、$0=t_0< t_1<\cdots< t_k$ と $n_1,\dots,n_k\in\mathbb{Z}_{\ge0}$ について
$$
P\bigl(N(t_j)-N(t_{j-1})=n_j\ (j=1,\dots,k)\bigr)=\prod_{j=1}^ke^{-\lambda(t_j-t_{j-1})}\frac{(\lambda(t_j-t_{j-1}))^{n_j}}{n_j!}
$$
である。
段 1(計数過程であること)。 $P(T_k\le0)=0$ なので、確率 1 ですべての $T_k$ は正である。また $E[e^{-S_n}]=\prod_{i=1}^nE[e^{-T_i}]=\bigl(\tfrac{\lambda}{\lambda+1}\bigr)^n$ であり、Markovの不等式 の定理「Markov の不等式」を $e^{-S_n}$ に使うと $P(S_n\le t)=P(e^{-S_n}\ge e^{-t})\le e^{t}\bigl(\tfrac{\lambda}{\lambda+1}\bigr)^n\to0$ である。$S_n$ は $n$ について増加するので $P(\lim_nS_n\le t)\le P(S_n\le t)$ がすべての $n$ で成り立ち、$P(\lim_nS_n\le t)=0$ である。$t=1,2,\dots$ について合わせると、確率 1 で $S_n\to\infty$ である。この事象の上では $0< S_1< S_2<\cdots\to\infty$ であり、$N(t)=\max\{n\mid S_n\le t\}$($S_1>t$ なら $0$)は有限で、$N(0)=0$、単調非減少であり、$N$ は区間 $[S_n,S_{n+1})$ の上で一定の値 $n$ をとるので右連続である。また $N(t)\ge k\iff S_k\le t$ なので、到着時刻は $S_k$ である。
段 2(増分の同時分布)。 $0=t_0<\cdots< t_k=:t$ と $n_1,\dots,n_k$ をとり、$n:=n_1+\cdots+n_k$ とおく。事象「すべての $j$ で $N(t_j)-N(t_{j-1})=n_j$」は、「$S_n\le t< S_{n+1}$ であり、$S_1,\dots,S_n$ のうち区間 $I_j:=(t_{j-1},t_j]$ に入るものがちょうど $n_j$ 個」という事象である($n=0$ なら「$S_1>t$」)。$S_1<\cdots< S_n$ は小さい順に並んでいるので、これは「最初の $n_1$ 個が $I_1$、次の $n_2$ 個が $I_2$、…に入り、$S_{n+1}>t$」と同じである。この $(s_1,\dots,s_n)$ の集合を $D$ とすると、lem-pp-joint-density により
$$
P(\cdots)=\int_D\Bigl(\int_t^\infty\lambda^{n+1}e^{-\lambda s_{n+1}}\,ds_{n+1}\Bigr)ds_1\cdots ds_n=\lambda^ne^{-\lambda t}\,\mathrm{vol}(D)
$$
である。$D$ は、各 $j$ について「$I_j$ の中の $n_j$ 個の小さい順の点」の集合の直積であり、長さ $\ell$ の区間の中の $m$ 個の小さい順の点の集合の体積は $\ell^m/m!$ である(立方体 $I^m$ は、座標の大小の順序 $m!$ 通りで、体積の等しい $m!$ 個の部分に分かれ、座標の一致する部分は体積 $0$)。したがって $\mathrm{vol}(D)=\prod_j(t_j-t_{j-1})^{n_j}/n_j!$ であり、$\lambda^n=\prod_j\lambda^{n_j}$、$e^{-\lambda t}=\prod_je^{-\lambda(t_j-t_{j-1})}$ と合わせて主張の式を得る。
段 3(Poisson 過程であること)。 $0\le s_0< s_1<\cdots< s_r$ が与えられたとき、$s_0>0$ なら $t_0=0$ を先頭に加えた分点で段 2 を使い、最初の区間 $(0,s_0]$ の回数 $n_1$ について和をとる(その因子の和は $1$)。すると増分 $N(s_i)-N(s_{i-1})$ の同時分布は $\prod_i\mathrm{Po}(\lambda(s_i-s_{i-1}))$ の積の形になる。同時分布が周辺分布の積なので増分は独立であり、それぞれ $\mathrm{Po}(\lambda(s_i-s_{i-1}))$ に従う。$\square$
段 2 の計算は Dur19 pp. 175–176 の計算($N(t)$ が Poisson 分布に従うことを $S_n$ の密度で確かめ、増分の独立性は指数分布の無記憶性から導く)とは別の、同時密度による証明である。
主定理 1 は Poisson 過程を 1 つ作った。次の定理は、定義を満たす過程はどれもこの構成と同じ形をしていることを示す。
$N$ を強度 $\lambda$ の Poisson 過程とする。このとき確率 1 で $0< S_1< S_2<\cdots$、$S_k\to\infty$ であり、到着間隔 $T_1,T_2,\dots$ は独立で、どれも $\mathrm{Exp}(\lambda)$ に従う。
$N'$ を thm-pp-construction で作った強度 $\lambda$ の Poisson 過程とし、その到着時刻を $S'_k$ とする。
段 1(有限個の時刻での値の分布は同じ)。 $0< a_1<\cdots< a_m$ について、$(N(a_1),\dots,N(a_m))$ は増分 $N(a_1)-N(0),\dots,N(a_m)-N(a_{m-1})$ の和で書け、増分は独立で分布は $\mathrm{Po}(\lambda(a_i-a_{i-1}))$($a_0:=0$)と決まっている。よって $(N(a_1),\dots,N(a_m))$ の同時分布は $\lambda$ と $a_i$ だけで決まり、$N'$ についての同時分布と等しい。
段 2(到着時刻の分布は同じ)。 $n\ge1$ と実数 $a_1,\dots,a_n\ge0$ について、lem-pp-switch により
$$
P(S_1\le a_1,\dots,S_n\le a_n)=P\bigl(N(a_1)\ge1,\dots,N(a_n)\ge n\bigr)
$$
であり、右辺は段 1 により $N'$ で計算したものと等しい($a_i$ に同じ値があっても、異なる値だけを並べて段 1 を使えばよい)。したがって $[0,\infty]^n$ に値をとる $(S_1,\dots,S_n)$ と $(S'_1,\dots,S'_n)$ の同時分布関数は、有限の $a_i$ で一致する。$a_i=\infty$ を含む値は、いつも成り立つ条件 $S_i\le\infty$ を除いた残りの $i$ についての同じ形の値なので、同じ理由で一致する。$\prod_i[0,a_i]$($a_i\in[0,\infty]$)の形の集合は積について閉じ、$[0,\infty]^n$ の Borel 集合族を生成するので、2 つの同時分布は等しい。
段 3(結論)。 $(S'_1,\dots,S'_n)$ は確率 1 で $0< S'_1<\cdots< S'_n< \infty$ を満たすので、$(S_1,\dots,S_n)$ も同じである。到着間隔 $(T_1,\dots,T_n)$ は $(S_1,\dots,S_n)$ の同じ線形変換なので、$(T'_1,\dots,T'_n)$ と同じ同時分布、すなわち独立な $\mathrm{Exp}(\lambda)$ の組の分布をもつ。$n$ は任意なので、$T_1,T_2,\dots$ は独立で、どれも $\mathrm{Exp}(\lambda)$ に従う。最後に $P(S_n\le t)=P(N(t)\ge n)\to0$($n\to\infty$)なので、確率 1 で $S_n\to\infty$ である。$\square$
主定理 1 と 2 を合わせると、強度 $\lambda$ の Poisson 過程と、「独立な $\mathrm{Exp}(\lambda)$ の待ち時間を順に並べた出来事の列」は、同じものの 2 つの見方である。とくに Poisson 過程の分布(有限個の時刻での値の同時分布)は $\lambda$ だけで決まる。
強度 $\lambda$ の Poisson 過程と $t>0$、$n\ge1$ について、$N(t)=n$ という条件のもとで $(S_1,\dots,S_n)$ は集合 $\Delta:=\{0< s_1<\cdots< s_n\le t\}$ の上の一様分布に従う。すなわち、$[0,t]$ の一様分布に従う独立な $n$ 個の点を小さい順に並べたものと同じ分布である。
要点:thm-pp-construction の証明の段 2 と同じ計算で、条件付きの密度が定数になる。
thm-pp-interarrival により、どの Poisson 過程の $(S_1,S_2,\dots)$ も thm-pp-construction の構成と同じ分布をもち、$N(t)=n\iff S_n\le t< S_{n+1}$ なので、構成した過程で示せばよい。Borel 集合 $B\subset\Delta$ について、lem-pp-joint-density により$$P\bigl((S_1,\dots,S_n)\in B,\ N(t)=n\bigr)=P\bigl((S_1,\dots,S_n)\in B,\ S_{n+1}>t\bigr)=\lambda^ne^{-\lambda t}\,\mathrm{vol}(B)$$である。$B=\Delta$ とすると $P(N(t)=n)=\lambda^ne^{-\lambda t}t^n/n!$ なので、条件付き確率は $\mathrm{vol}(B)/\mathrm{vol}(\Delta)$ である。独立な一様な $n$ 点 $U_1,\dots,U_n$ を小さい順に並べた組が $B$ に入るのは、ある置換 $\sigma$ で $(U_{\sigma(1)},\dots,U_{\sigma(n)})\in B$ となるときで、$B\subset\Delta$ なのでそのような $\sigma$ は高々 1 つである。各 $\sigma$ でその確率は $\mathrm{vol}(B)/t^n$ なので、合わせて $n!\,\mathrm{vol}(B)/t^n=\mathrm{vol}(B)/\mathrm{vol}(\Delta)$ となり、一致する。$\square$
Poisson 分布という具体的な形を仮定せずに、次の 4 つの性質だけから Poisson 過程が出てくる。これが Poisson 過程を「ばらばらに起こる出来事」の模型として使う根拠である。
$\lambda>0$ とし、計数過程 $N$ が次を満たすとする。$o(h)$ は $h\downarrow0$ で $o(h)/h\to0$ となる量を表す。
逆。 1 と 2 は定義から出る。$N(h)\sim\mathrm{Po}(\lambda h)$ なので $P(N(h)=1)=\lambda he^{-\lambda h}=\lambda h+\lambda h(e^{-\lambda h}-1)$ で、第 2 項を $h$ で割ると $0$ に近づく。また $P(N(h)\ge2)=1-e^{-\lambda h}-\lambda he^{-\lambda h}$ であり、$e^{-\lambda h}=1-\lambda h+O(h^2)$ を代入すると $O(h^2)=o(h)$ である。
本題。 $0\le s< t$、$h:=t-s$、$k\in\mathbb{Z}_{\ge0}$ を固定し、$P(N(t)-N(s)=k)=\mathrm{Po}(\lambda h)(\{k\})$ を示せばよい(独立増分は仮定されている)。$n\ge1$ について区間 $(s,t]$ を $n$ 等分し、$m$ 番目の小区間での増分を $X_{n,m}$($m=1,\dots,n$)とする。1 により $X_{n,1},\dots,X_{n,n}$ は独立であり、2 によりどれも $N(h/n)$ と同じ分布をもつ。$B_{n,m}:=1$($X_{n,m}\ge1$ のとき)、$0$(それ以外)とおき、$S_n:=\sum_mB_{n,m}$ とする。$B_{n,1},\dots,B_{n,n}$ は独立で(独立な確率変数の関数だから)、どれも
$$
p_n:=P(N(h/n)\ge1)=\lambda\frac hn+o\Bigl(\frac hn\Bigr)
$$
を成功確率とする。よって $np_n\to\lambda h$、$p_n\to0$ であり、Poisson分布 の系「少数の法則」により $P(S_n=k)\to\mathrm{Po}(\lambda h)(\{k\})$ である。
一方、$N(t)-N(s)=\sum_mX_{n,m}$ と $S_n$ が異なるのは、ある $m$ で $X_{n,m}\ge2$ のときに限るので、4 により
$$
P\bigl(N(t)-N(s)\ne S_n\bigr)\le\sum_{m=1}^nP(X_{n,m}\ge2)=n\cdot o\Bigl(\frac hn\Bigr)=h\cdot\frac{o(h/n)}{h/n}\to0
$$
である。したがって $\lvert P(N(t)-N(s)=k)-P(S_n=k)\rvert\le P(N(t)-N(s)\ne S_n)\to0$ であり、左辺の第 1 項は $n$ によらないので $P(N(t)-N(s)=k)=\mathrm{Po}(\lambda h)(\{k\})$ である。$\square$
これは Dur19 Theorem 3.7.2(p. 175)であり、そこでも区間を $n$ 等分して Poisson 近似の定理を使っている。仮定 4 を外すと結論は成り立たない(反例の表の 3 行目)。
$N$ を強度 $\lambda$ の Poisson 過程、$m\ge1$、$p_1,\dots,p_m\ge0$、$p_1+\cdots+p_m=1$ とする。$\{1,\dots,m\}$ に値をとる確率変数 $C_1,C_2,\dots$ が独立で、どれも $P(C_k=i)=p_i$ を満たし、族 $(C_k)_k$ と族 $(N(t))_t$ は独立であるとする。第 $k$ の出来事に色 $C_k$ を塗り、
$$
N_i(t):=\#\{k\mid 1\le k\le N(t),\ C_k=i\}
$$
とおく。このとき $N_1,\dots,N_m$ は、それぞれ強度 $\lambda p_i$ の Poisson 過程であり($p_i=0$ なら確率 1 で $N_i\equiv0$)、互いに独立である(有限個の時刻での値の同時分布が、過程ごとの分布の積になる)。
要点:分点ごとの回数と色の数え方を掛け合わせると、Poisson分布 の定理「Poisson 分布の色分け」と同じ計算になる。
各 $N_i$ が計数過程であることは定義から明らかである。$0=t_0< t_1<\cdots< t_r$ をとり、$\Delta_j:=t_j-t_{j-1}$、区間 $(t_{j-1},t_j]$ での $N_i$ の増分を $Y_{ij}$ とする。$k_{ij}\in\mathbb{Z}_{\ge0}$ を与え、$n_j:=\sum_ik_{ij}$ とおく。事象「すべての $i,j$ で $Y_{ij}=k_{ij}$」は、「すべての $j$ で $N(t_j)-N(t_{j-1})=n_j$」と「各 $j$ について、番号 $n_1+\cdots+n_{j-1}+1$ から $n_1+\cdots+n_j$ までの $C_k$ のうち値 $i$ のものがちょうど $k_{ij}$ 個」の共通部分である。後者は $(C_k)$ の重ならない番号の組で決まる事象の共通部分で、$(C_k)$ と $N$ は独立なので、確率は$$\prod_{j=1}^r\Bigl(e^{-\lambda\Delta_j}\frac{(\lambda\Delta_j)^{n_j}}{n_j!}\cdot\frac{n_j!}{\prod_ik_{ij}!}\prod_ip_i^{k_{ij}}\Bigr)=\prod_{j=1}^r\prod_{i=1}^me^{-\lambda p_i\Delta_j}\frac{(\lambda p_i\Delta_j)^{k_{ij}}}{k_{ij}!}$$である(多項係数 の命題「結果の回数の分布」と $\sum_ip_i=1$ による)。右辺は $(i,j)$ ごとの $\mathrm{Po}(\lambda p_i\Delta_j)$ の確率の積なので、$Y_{ij}$ はすべて独立で $Y_{ij}\sim\mathrm{Po}(\lambda p_i\Delta_j)$ である。$t_0=0$ でない分点は、thm-pp-construction の証明の段 3 と同じく $0$ を加えて足し上げればよい。したがって各 $N_i$ は強度 $\lambda p_i$ の Poisson 過程である。どの有限個の時刻も 1 組の分点に含められ、$N_i$ の値はその増分の和なので、有限個の時刻での $(N_1,\dots,N_m)$ の値の同時分布は過程ごとの分布の積になる。$\square$
これは Dur19 Theorem 3.7.4(p. 178)である。驚くべき点は独立性で、来客が強度 $10$(1 時間あたり)の Poisson 過程で、客が独立に確率 $1/2$ ずつ男女に分かれるなら、ある 1 時間に男性客が多かったことは、同じ時間の女性客の数について何の情報も与えない。逆向きの操作が重ね合わせで、独立な強度 $\lambda$ と $\mu$ の Poisson 過程の和は強度 $\lambda+\mu$ の Poisson 過程である。実際、和の増分は 2 つの独立な Poisson 分布に従う増分の和なので Poisson分布 の定理「Poisson 分布の再生性」により $\mathrm{Po}((\lambda+\mu)(t-s))$ に従い、重ならない区間の増分の独立性も 2 つの過程の独立性と独立増分から出る。
ある窓口への電話が強度 $\lambda=3$(1 分あたり)の Poisson 過程に従うとする。2 分間に 1 本も来ない確率は $P(N(2)=0)=e^{-6}\approx0.00248$ である。最初の 1 分に 2 本、次の 2 分に 4 本来る確率は、独立増分により
$$
e^{-3}\frac{3^2}{2!}\cdot e^{-6}\frac{6^4}{4!}=\frac92\cdot54\cdot e^{-9}\approx0.0300
$$
である。最初の着信までの待ち時間は $\mathrm{Exp}(3)$ に従い、平均 $1/3$ 分である(thm-pp-interarrival)。
バスが強度 $\lambda$ の Poisson 過程に従って来るとき、バスの間隔の平均は $1/\lambda$ である。ところが、時刻 $t>0$ に停留所に着いた人が「自分が着いた時刻を含む間隔」$S_{N(t)+1}-S_{N(t)}$ の長さを測ると、その平均は $1/\lambda$ より長い。
次のバスまでの待ち時間 $W:=S_{N(t)+1}-t$ は $\mathrm{Exp}(\lambda)$ に従う。実際、thm-pp-interarrival により thm-pp-construction の構成で計算してよく、$u\ge0$ について、lem-pp-joint-density の密度を $n$ ごとに積分すると
$$
P(W>u)=\sum_{n\ge0}P(S_n\le t,\ S_{n+1}>t+u)=\sum_{n\ge0}\lambda^ne^{-\lambda(t+u)}\frac{t^n}{n!}=e^{-\lambda u}
$$
である(各項は thm-pp-construction の証明の段 2 で $S_{n+1}>t$ を $S_{n+1}>t+u$ に替えたもの)。前のバスからの経過時間 $A:=t-S_{N(t)}$ は、$0\le a< t$ について $P(A>a)=P(N(t)-N(t-a)=0)=e^{-\lambda a}$、$a\ge t$ で $P(A>a)=0$ なので、$E[A]=\int_0^te^{-\lambda a}\,da=(1-e^{-\lambda t})/\lambda$ である。よって着いた時刻を含む間隔の平均は
$$
E[W]+E[A]=\frac{2-e^{-\lambda t}}{\lambda}
$$
で、$t$ が大きいとほぼ $2/\lambda$ になる。長い間隔ほど「着いた時刻を含む」確率が高いので、平均より長い間隔を測りやすい。待ち時間の平均が $1/\lambda$ のままで変わらないのは指数分布の無記憶性(一様分布と指数分布 の定理「無記憶性と指数分布」)の表れであり、時刻表どおりに間隔 $1/\lambda$ で来るバスなら待ち時間の平均はその半分の $1/(2\lambda)$ になる。
強度 $\lambda$ の Poisson 過程で、区間 $[0,1]$ にちょうど 1 回だけ出来事が起こったと分かっているとき、その時刻 $S_1$ は $[0,1]$ の一様分布に従う(cor-pp-uniform で $n=1$)。強度 $\lambda$ の値によらないことに注意する。たとえば $P(S_1\le 1/4\mid N(1)=1)=1/4$ である。
| 外す条件 | 反例 | 成り立たなくなること |
|---|---|---|
| thm-pp-construction の、間隔が指数分布であること | 間隔がすべて定数 $1$($T_k\equiv1$) | $N(t)=\lfloor t\rfloor$ は確率 $1$ で決まった値をとり、Poisson 分布に従わない |
| thm-pp-construction の、間隔が独立であること | $T\sim\mathrm{Exp}(\lambda)$ を 1 つとり、$T_1=T_2=\cdots=T$ | $N(t)=\lfloor t/T\rfloor$ で、$P(N(t)=0)=e^{-\lambda t}$ は合うが、$P(N(t)\ge2)=1-e^{-\lambda t/2}$ が Poisson 分布の値 $1-e^{-\lambda t}-\lambda te^{-\lambda t}$ と異なる |
| thm-pp-local の 4(2 回以上が $o(h)$) | 独立な強度 $\lambda$ と $\mu>0$ の Poisson 過程 $\tilde N_1,\tilde N_2$ について $N:=\tilde N_1+2\tilde N_2$ | 1・2・3 は満たすが、$P(N(h)\ge2)\ge1-e^{-\mu h}$ は $o(h)$ でなく、$E[N(h)]=(\lambda+2\mu)h$ なので $N$ は強度 $\lambda$ の Poisson 過程でない |
| thm-pp-local の 2(定常増分) | 強度 $1$ の Poisson 過程 $\tilde N$ と $\Lambda(t):=\lambda t+\lambda\max\{t-1,0\}$ について $N(t):=\tilde N(\Lambda(t))$ | 1・3・4 は満たすが($h\le1$ で $N(h)\sim\mathrm{Po}(\lambda h)$)、$N(2)-N(1)\sim\mathrm{Po}(2\lambda)$ なので $N$ は強度 $\lambda$ の Poisson 過程でない |
| thm-pp-thinning の、色が独立に同じ分布に従うこと | 色を順に $C_1=1,C_2=2,C_3=1,\dots$ と交互に塗る | $N_1(t)-N_2(t)\in\{0,1\}$ がいつも成り立ち、$N_1$ と $N_2$ は独立でない |
要点:2〜5 行目について、外した条件以外の条件を満たすことと、結論が崩れることを直接の計算で確かめる。
2 行目:$N(t)\ge2\iff 2T\le t$ なので $P(N(t)\ge2)=1-e^{-\lambda t/2}$ である。$\lambda t=1$ のとき、これは $1-e^{-1/2}\approx0.393$、Poisson 分布の値は $1-2e^{-1}\approx0.264$ で異なる。
3 行目:$N$ は計数過程で、増分は独立な 2 つの過程の同じ区間の増分で決まるので独立で定常である。$P(N(h)=1)=P(\tilde N_1(h)=1)P(\tilde N_2(h)=0)=\lambda he^{-(\lambda+\mu)h}=\lambda h+o(h)$ であり、$P(N(h)\ge2)\ge P(\tilde N_2(h)\ge1)=1-e^{-\mu h}$ を $h$ で割ると $\mu$ に近づく。強度 $\lambda$ の Poisson 過程なら $E[N(h)]=\lambda h$ である。
4 行目:$\Lambda$ は増加で連続なので $N$ は計数過程であり、重ならない区間は重ならない区間に写るので増分は独立で、$N(t)-N(s)\sim\mathrm{Po}(\Lambda(t)-\Lambda(s))$ である。$h\le1$ では $N(h)\sim\mathrm{Po}(\lambda h)$ なので、3・4 は thm-pp-local の証明の「逆」と同じ計算で成り立つ。一方 $N(2)-N(1)\sim\mathrm{Po}(2\lambda)$ は $N(1)\sim\mathrm{Po}(\lambda)$ と分布が異なる。
5 行目:$N(t)=n$ のとき、色 $1$ の個数は $\lceil n/2\rceil$、色 $2$ は $\lfloor n/2\rfloor$ である。$P(N_2(t)=0)>0$ だが、$N_1(t)\ge2$ なら $N(t)\ge3$ で $N_2(t)\ge1$ なので $P(N_1(t)\ge2,N_2(t)=0)=0\ne P(N_1(t)\ge2)P(N_2(t)=0)$ である。
Mathpediaは寄付と、参考文献の書籍リンク(Amazonアソシエイト)の紹介料で運営されています。 支援について / 寄付する