本記事は「これなら分かる!はじめての数理統計学」シリーズに含まれます。目次へ不適切な内容があれば,記事下のコメント欄またはお問い合わせフォームよりご連絡下さい。\frac{\pi^2}{3}\end{align}ただし,$\Gamma(\cdot)$はガンマ関数を表す。\int_0^{\infty}t^{n-1}e^{-t}dt\end{align}標準ロジスティック分布は,累積分布関数が標準シグモイド関数となる確率分布として定式化されます。\bbR\end{align}であり,モーメント母関数の変数は$|t|<1$とします。ロジスティック分布は,累積分布関数が標準シグモイド関数となる確率分布として定式化されます。\frac{1}{1+e^{-x}}\end{align}したがって,確率密度関数は以下のように導出されます。\frac{e^{-x}}{\left(1+e^{-x}\right)^2}\end{align}モーメント母関数モーメント母関数の定義に従って計算していきます。\label{モーメント母関数}\end{align}ここで,標準ロジスティック分布の累積分布関数が標準シグモイド関数であったことに着目しましょう。以下のように置換します。\label{置換}\end{align}これは,標準シグモイド関数を$u$と置くことを意味しています。後のために,$e^x$を$u$で表しておきましょう。\label{シグモイドの微分}\end{align}となります。式($\ref{シグモイドの微分}$)はベータ関数の被積分関数の形をしていますので,式($\ref{置換}$)の置換により標準ロジスティック分布のモーメント母関数($\ref{モーメント母関数}$)がベータ関数を用いて表せることを期待します。\Gamma(1+t)\Gamma(1-t)\label{4}\end{align}と計算できます。平均と分散を求めるためには,モーメント母関数の性質を利用します。原点周りの一次モーメント,すなわち期待値を求めます。0\end{align}続いて,原点周りの二次モーメントを求めます。\label{二次モーメント}\end{align}ただし,$\Gamma(1)=1$を利用しました。ここで,やや天下り的なのですが,式($\ref{二次モーメント}$)はポリガンマ関数の級数表示を用いて計算することができます。順を追って説明していきます。ディガンマ関数$\psi(x)$は対数ガンマ関数の導関数として定義されます。\frac{\Gamma^{\prime}(x)}{\Gamma(x)}\end{align}ポリガンマ関数は,ディガンマ関数の$m$回微分として定義されます。\frac{d^{m}}{dx^{m}}\psi(x)\end{align}$m=1$のポリガンマ関数を計算してみましょう。\frac{\Gamma^{\prime\prime}(x)\Gamma(x)-\left\{\Gamma^{\prime}(x)\right\}^2}{\left\{\Gamma(x)\right\}^2}\end{align}分子に式($\ref{二次モーメント}$)の形が出現しました。$x=1$を代入してみましょう。\psi^{(1)}(1)\label{整理後}\end{align}ただし,$\Gamma(1)=1$を利用しました。\zeta(2)\end{align}となります。,$\zeta(s)$はリーマンゼータ関数を表します。$\zeta(2)$の値はバーゼル問題としても知られていて,オイラーによって値が求められました。証明は別ページで行うことにして,結果のみを用います。\label{オイラー}\end{align}式($\ref{オイラー}$)を式($\ref{整理後}$)に代入し,その結果を式($\ref{二次モーメント}$)に代入すると,以下が得られます。\frac{\pi^2}{3}\end{align}最後に,一次モーメントと二次モーメントから分散を求めます。再生性を示すためには,再生性を示したい分布に従う独立な二つの確率変数を考え,その和のモーメント母関数を計算したときに,パラメータが和の形になっていることを示します。ロジスティック分布のモーメント母関数の積をとっても同じモーメント母関数の形が現れないため,ロジスティック分布に再生性はありません。ロードマップ確率分布のロードマップロジスティック分布はロードマップには含まれていませんが,累積分布関数が標準シグモイド関数となる確率分布として導出されました。本稿の執筆にあたり参考にした文献は,以下でリストアップしております。※コメントに返信があったらメールで通知する。私はロボットではありません。