アーラン分布の生存関数を計算する方法を解説！！

アーラン分布の生存関数をポアソン分布から計算する方法を解説します。
最初にアーラン分布を思い出しておきましょう。

生存関数$S(x)=1–F(x)$を計算してみましょう。
素直に積分を計算する方法と、ポアソン分布との関係を使って計算する方法があります。

素直に積分で計算する方法

$$\begin{array}{rl}S(x)& ={\int }_x^{\mathrm{\infty }}\frac{{\left(\frac{1}{\mu }\right)}^k{\xi }^{k-1}{e}^{-\frac{1}{\mu }\xi }}{\mathrm{\Gamma }(k)}d\xi \\ & =\frac{{\left(\frac{1}{\mu }\right)}^k}{\mathrm{\Gamma }(k)}{\left[({\xi }^{k-1}\mu +{}_{k-1}{P}_1{\xi }^{k-2}{\mu }^2+{}_{k-1}{P}_2{\xi }^{k-3}{\mu }^3+\cdots (k-1)!{\mu }^k)e^{–\frac{1}{\mu }\xi }\right]}_x^{\mathrm{\infty }}\\ & =\frac{{\left(\frac{1}{\mu }\right)}^k}{\mathrm{\Gamma }(k)}(-x^{k-1}\mu -{}_{k-1}{P}_1{x}^{k-2}{\mu }^2-{}_{k-1}{P}_2{\xi }^{k-3}{\mu }^3-\cdots (k-1)!{\mu }^k)e^{–\frac{1}{\mu }x}\end{array}$$
と計算できます。指数分布に係る部分積分を効率的に計算する方法は下の記事で解説しています。

あわせて読みたい

指数分布の確率密度関数との積分をする時の部分積分公式 この記事では、指数分布の確率密度関数との積分をする時の部分積分公式を解説します。よくある状況を考えるために、$f(x)\in C^{\mathrm{\infty }}$とし、$latex ^\exists{N}…

ポアソン分布を利用して計算する方法

アーラン分布は、指数分布の$k$個の足し合わせでした。
つまり、$X\sim \mathrm{\Gamma }(k,\frac{1}{\mu })$は、単位時間あたり発生回数が平均$\frac{1}{\mu }$のポアソン分布に従うイベントが、$k$回発生するまでの時間を表現しています。
そこで、単位時間発生回数が平均$\frac{1}{\mu }$のポアソン分布は、$x$時間で平均$\frac{1}{\mu }x$回発生するので、$Y\sim Po(\frac{1}{\mu }x)$として
$P(X>x)=P(Y<k)=e^{-\frac{1}{\mu }x}(\frac{{\left(\frac{1}{\mu }\right)}^0}{0!}+\frac{{\left(\frac{1}{\mu }\right)}^1}{1!}+\cdots +\frac{{\left(\frac{1}{\mu }\right)}^{k-1}}{(k-1)!})$
と計算することができます。
これは普通に積分で計算したときの結果と一致しています(少し表現は違いますが、一致していることが確認できると思います)。

いかがでしたでしょうか。個人的には、素朴に計算した方が良いと思います。