確率分布

連続型確率分布

統計検定®ドリルの「連続型確率分布」に関する問題と解説です。全10問。

この分野をドリルで練習する →
  1. Q1 · 2級 · 計算 · 一問一答

    ある工場で生産される部品の重さ $X$ は、平均 $500\,\text{g}$、標準偏差 $20\,\text{g}$ の正規分布 $N(500, 20^2)$ に従うとする。 この部品の重さが $530\,\text{g}$ 以上 $540\,\text{g}$ 未満となる確率を、以下の標準正規分布表を用いて求めよ。 【標準正規分布表(必要な範囲のみ抜粋)】表中の数値は、標準正規分布に従う確率変数 $Z$ について上側確率 $P(Z \geq z)$ を表す。$z$ の値は「行(小数第1位まで)」と「列(小数第2位)」の交点から読み取る。 $$\begin{array}{c|cccccccccc} z & 0.00 & 0.01 & 0.02 & 0.03 & 0.04 & 0.05 & 0.06 & 0.07 & 0.08 & 0.09 \\ \hline 1.4 & 0.0808 & 0.0793 & 0.0778 & 0.0764 & 0.0749 & 0.0735 & 0.0721 & 0.0708 & 0.0694 & 0.0681 \\ 1.5 & 0.0668 & 0.0655 & 0.0643 & 0.0630 & 0.0618 & 0.0606 & 0.0594 & 0.0582 & 0.0571 & 0.0559 \\ 1.6 & 0.0548 & 0.0537 & 0.0526 & 0.0516 & 0.0505 & 0.0495 & 0.0485 & 0.0475 & 0.0465 & 0.0455 \\ 1.9 & 0.0287 & 0.0281 & 0.0274 & 0.0268 & 0.0262 & 0.0256 & 0.0250 & 0.0244 & 0.0239 & 0.0233 \\ 2.0 & 0.0228 & 0.0222 & 0.0217 & 0.0212 & 0.0207 & 0.0202 & 0.0197 & 0.0192 & 0.0188 & 0.0183 \\ 2.1 & 0.0179 & 0.0174 & 0.0170 & 0.0166 & 0.0162 & 0.0158 & 0.0154 & 0.0150 & 0.0146 & 0.0143 \end{array}$$

    答え

    $0.0440$(約$4.40\%$)

    標準化

    確率変数 $X$ が正規分布 $N(\mu, \sigma^2)$ に従うとき、$Z = \dfrac{X - \mu}{\sigma}$ とおくと $Z$ は標準正規分布 $N(0, 1)$ に従う。この変換を「標準化」と呼び、正規分布表を使うためには必ずこの変換を行う必要がある。 本問では $\mu = 500$、$\sigma = 20$ なので、 $X = 530$ のとき $Z = \dfrac{530 - 500}{20} = 1.50$ $X = 540$ のとき $Z = \dfrac{540 - 500}{20} = 2.00$ したがって、求める確率は $P(530 \leq X < 540) = P(1.50 \leq Z < 2.00)$ と書き換えられる。

    表の読み取りと引き算

    標準正規分布表は「$Z$ がある値 $z$ 以上になる確率(上側確率)」$P(Z \geq z)$ を表にしたものである。 表より、$P(Z \geq 1.50) = 0.0668$、$P(Z \geq 2.00) = 0.0228$ と読み取れる。 $P(1.50 \leq Z < 2.00)$ は「$Z \geq 1.50$ となる面積」から「$Z \geq 2.00$ となる面積」を引いた部分に対応する(裾の広い方から、より裾の狭い方を除く、というイメージ)。 $$P(1.50 \leq Z < 2.00) = P(Z \geq 1.50) - P(Z \geq 2.00) = 0.0668 - 0.0228 = 0.0440$$ よって、求める確率は約 $0.0440$(4.40%)である。
  2. Q2 · 準1級 · 暗記 · 4択

    確率変数 $X$ が平均 $\mu$、分散 $\sigma^2$ の正規分布 $N(\mu, \sigma^2)$ に従うとき、$X$ の確率密度関数 $f(x)$ として正しいものはどれか。
    • $f(x) = \dfrac{1}{\sigma\sqrt{2\pi}} \exp\left(-\dfrac{(x-\mu)^2}{2\sigma^2}\right)$
    • $f(x) = \dfrac{1}{\sigma\sqrt{2\pi}} \exp\left(-\dfrac{(x-\mu)^2}{\sigma^2}\right)$
    • $f(x) = \dfrac{1}{\sigma^2\sqrt{2\pi}} \exp\left(-\dfrac{(x-\mu)^2}{2\sigma^2}\right)$
    • $f(x) = \dfrac{1}{\sigma\sqrt{2\pi}} \exp\left(\dfrac{(x-\mu)^2}{2\sigma^2}\right)$

    答え

    $f(x) = \dfrac{1}{\sigma\sqrt{2\pi}} \exp\left(-\dfrac{(x-\mu)^2}{2\sigma^2}\right)$

    補足

    正規化定数 $\dfrac{1}{\sigma\sqrt{2\pi}}$ は、$f(x)$ を $-\infty$ から $\infty$ まで積分した値がちょうど $1$ になるように決められている。 指数部の $-\dfrac{(x-\mu)^2}{2\sigma^2}$ は、$x=\mu$ で最大値(頂点)をとり、$x$ が $\mu$ から離れるほど急速に $0$ に近づく釣鐘型の形状を作り出している。$\sigma$ が大きいほど山の広がりが大きく(なだらかに)なる。 特に $\mu=0,\ \sigma^2=1$ の場合が標準正規分布 $N(0,1)$ であり、標準正規分布表はこの特別な場合の面積(確率)をあらかじめ計算してまとめたものである。

    正規化係数の導出

    なぜ先頭の係数が $\dfrac{1}{\sigma\sqrt{2\pi}}$ になるのかを確認する。確率密度関数である以上、 $$\int_{-\infty}^{\infty} f(x)\,dx = 1$$ を満たす必要がある。そこで、まず係数を $C$ とおいた $$g(x) = C \exp\left(-\dfrac{(x-\mu)^2}{2\sigma^2}\right)$$ について、$\displaystyle\int_{-\infty}^{\infty} g(x)\,dx = 1$ となるように $C$ を定める。 置換 $t = \dfrac{x-\mu}{\sigma}$(すなわち $x = \mu + \sigma t$、$dx = \sigma\,dt$)を行うと、 $$\int_{-\infty}^{\infty} \exp\left(-\dfrac{(x-\mu)^2}{2\sigma^2}\right) dx = \sigma \int_{-\infty}^{\infty} e^{-t^2/2}\,dt$$ となるので、あとは基準となる積分 $\displaystyle I = \int_{-\infty}^{\infty} e^{-t^2/2}\,dt$ の値を求めればよい。
    $I$ の値は、$I$ を直接求める代わりに $I^2$ を考える、というのが定石の求め方である。$s$ という別の積分変数も使って $$I^2 = \left(\int_{-\infty}^{\infty} e^{-t^2/2}\,dt\right)\left(\int_{-\infty}^{\infty} e^{-s^2/2}\,ds\right) = \iint_{\mathbb{R}^2} e^{-(t^2+s^2)/2}\,dt\,ds$$ と、$1$変数の積分の積を$2$変数(平面)上の積分に書き換える。 これを極座標 $t = r\cos\theta,\ s = r\sin\theta$($dt\,ds = r\,dr\,d\theta$)に変換すると、$t^2+s^2 = r^2$ なので $$I^2 = \int_{0}^{2\pi}\int_{0}^{\infty} e^{-r^2/2}\, r\,dr\,d\theta = 2\pi \int_{0}^{\infty} r e^{-r^2/2}\,dr$$ ここで $u = r^2/2$ と置換すると $du = r\,dr$ なので、$\displaystyle\int_{0}^{\infty} r e^{-r^2/2}\,dr = \int_{0}^{\infty} e^{-u}\,du = 1$ となり、 $$I^2 = 2\pi \times 1 = 2\pi \quad \Rightarrow \quad I = \sqrt{2\pi}\ \ (I>0 \text{ より})$$ したがって $\displaystyle\int_{-\infty}^{\infty} \exp\left(-\dfrac{(x-\mu)^2}{2\sigma^2}\right) dx = \sigma\sqrt{2\pi}$ となり、$g(x)$ の積分を $1$ にするには $C = \dfrac{1}{\sigma\sqrt{2\pi}}$ とすればよいことがわかる。これが正規分布の確率密度関数の先頭に $\dfrac{1}{\sigma\sqrt{2\pi}}$ が現れる理由である。
  3. Q3 · 準1級 · 暗記 · 4択

    形状母数 $\alpha > 0$、尺度母数 $\beta > 0$ のガンマ分布 $\mathrm{Ga}(\alpha, \beta)$ の確率密度関数として正しいものはどれか。
    • $f(x) = \dfrac{1}{\Gamma(\alpha)\,\beta^{\alpha}}x^{\alpha-1}e^{-x/\beta}\ \ (x>0)$
    • $f(x) = \dfrac{1}{\Gamma(\alpha)\,\beta^{\alpha}}x^{\alpha}e^{-x/\beta}\ \ (x>0)$
    • $f(x) = \dfrac{1}{\Gamma(\alpha)\,\beta^{\alpha-1}}x^{\alpha-1}e^{-x/\beta}\ \ (x>0)$
    • $f(x) = \dfrac{\beta^{\alpha}}{\Gamma(\alpha)}x^{\alpha-1}e^{-x/\beta}\ \ (x>0)$

    答え

    $f(x) = \dfrac{1}{\Gamma(\alpha)\,\beta^{\alpha}}x^{\alpha-1}e^{-x/\beta}\ \ (x>0)$

    正規化定数の確認

    $y = x/\beta$ と置換すると $dx = \beta\,dy$ なので $$\int_0^\infty x^{\alpha-1}e^{-x/\beta}\,dx = \beta^{\alpha}\int_0^\infty y^{\alpha-1}e^{-y}\,dy = \Gamma(\alpha)\,\beta^{\alpha}$$ したがって、全体の積分を $1$ にするための係数は $\dfrac{1}{\Gamma(\alpha)\,\beta^{\alpha}}$ になる。$x$ の指数 $\alpha-1$ とガンマ関数の定義 $\Gamma(\alpha) = \int_0^\infty x^{\alpha-1}e^{-x}dx$ の指数が揃っていることがポイント。

    特別な場合

    $\alpha = 1$ のときは平均 $\beta$ の指数分布、$\alpha = k/2,\ \beta = 2$ のときは自由度 $k$ のカイ二乗分布 $\chi^2(k)$ になる。また、平均 $\beta$ の指数分布に独立に従う $n$ 個の和は $\mathrm{Ga}(n, \beta)$ に従う。
  4. Q4 · 準1級 · 暗記 · 4択

    確率変数 $X$ が、形状母数 $\alpha > 0$、尺度母数 $\beta > 0$ のガンマ分布に従い、確率密度関数 $f(x) = \dfrac{1}{\Gamma(\alpha)\,\beta^{\alpha}}x^{\alpha-1}e^{-x/\beta}\ \ (x>0)$ をもつとき、$X$ の期待値 $E[X]$ と分散 $V[X]$ の組として正しいものはどれか。
    • $E[X] = \alpha\beta,\ V[X] = \alpha\beta^2$
    • $E[X] = \dfrac{\alpha}{\beta},\ V[X] = \dfrac{\alpha}{\beta^2}$
    • $E[X] = \alpha\beta,\ V[X] = \alpha^2\beta^2$
    • $E[X] = \alpha\beta,\ V[X] = \alpha(\alpha+1)\beta^2$

    答え

    $E[X] = \alpha\beta,\ V[X] = \alpha\beta^2$

    モーメントの導出

    $y = x/\beta$ と置換すると $\displaystyle\int_0^\infty x^{s-1}e^{-x/\beta}\,dx = \beta^{s}\,\Gamma(s)$ となる。これを使うと $$E[X^k] = \frac{1}{\Gamma(\alpha)\,\beta^{\alpha}}\int_0^\infty x^{\alpha+k-1}e^{-x/\beta}\,dx = \frac{\beta^{\alpha+k}\,\Gamma(\alpha+k)}{\Gamma(\alpha)\,\beta^{\alpha}} = \frac{\Gamma(\alpha+k)}{\Gamma(\alpha)}\,\beta^{k}$$ $\Gamma(\alpha+1) = \alpha\Gamma(\alpha)$ より $E[X] = \alpha\beta$、$\Gamma(\alpha+2) = (\alpha+1)\alpha\Gamma(\alpha)$ より $E[X^2] = \alpha(\alpha+1)\beta^2$。 $$V[X] = \alpha(\alpha+1)\beta^2 - \alpha^2\beta^2 = \alpha\beta^2$$ $\alpha(\alpha+1)\beta^2$ は $E[X^2]$ であり、平均の2乗を引き忘れた値である。

    レート表記との対応

    尺度母数の逆数であるレート $\lambda = 1/\beta$ を用いて $f(x) = \dfrac{\lambda^{\alpha}}{\Gamma(\alpha)}x^{\alpha-1}e^{-\lambda x}$ と書く流儀もあり、その場合は $E[X] = \dfrac{\alpha}{\lambda}$、$V[X] = \dfrac{\alpha}{\lambda^2}$ となる。選択肢の $\dfrac{\alpha}{\beta},\ \dfrac{\alpha}{\beta^2}$ は、本問の尺度表記の密度に対してレート表記の公式を当てはめてしまった誤り。問題文の密度関数でどちらの表記かを必ず確認すること。
    $\alpha = 1$ は平均 $\beta$ の指数分布で $E = \beta,\ V = \beta^2$、$\alpha = k/2,\ \beta = 2$ は自由度 $k$ のカイ二乗分布で $E = k,\ V = 2k$ となり、いずれもこの公式に代入すればすぐに得られる。
  5. Q5 · 準1級 · 計算 · 4択

    確率変数 $X$ が確率密度関数 $f(x) = \dfrac{1}{\Gamma(\alpha)\,\beta^{\alpha}}x^{\alpha-1}e^{-x/\beta}\ \ (x>0)$ をもつガンマ分布に従うとする($\alpha > 1$)。$E\!\left[\dfrac{1}{X}\right]$ として正しいものはどれか。
    • $\dfrac{1}{(\alpha-1)\beta}$
    • $\dfrac{1}{\alpha\beta}$
    • $(\alpha-1)\beta$
    • $\dfrac{1}{(\alpha+1)\beta}$

    答え

    $\dfrac{1}{(\alpha-1)\beta}$

    ガンマ関数を使った計算

    $\displaystyle\int_0^\infty x^{s-1}e^{-x/\beta}\,dx = \Gamma(s)\,\beta^{s}$ を $s = \alpha-1$ として使うと $$E\!\left[\frac{1}{X}\right] = \frac{1}{\Gamma(\alpha)\,\beta^{\alpha}}\int_0^\infty x^{\alpha-2}e^{-x/\beta}\,dx = \frac{\Gamma(\alpha-1)\,\beta^{\alpha-1}}{\Gamma(\alpha)\,\beta^{\alpha}} = \frac{\Gamma(\alpha-1)}{\Gamma(\alpha)\,\beta} = \frac{1}{(\alpha-1)\beta}$$ 最後に $\Gamma(\alpha) = (\alpha-1)\Gamma(\alpha-1)$ を使った。積分 $\int_0^\infty x^{\alpha-2}e^{-x/\beta}dx$ が有限になるには $\alpha - 1 > 0$ が必要で、これが条件 $\alpha > 1$ の理由である。

    1/E[X] との違い

    $\dfrac{1}{\alpha\beta}$ は $\dfrac{1}{E[X]}$ であり、一般に $E[1/X] \neq 1/E[X]$ である。実際、$1/x$ は $x>0$ で下に凸なのでイェンセンの不等式より $E[1/X] \ge 1/E[X]$ となり、$\dfrac{1}{(\alpha-1)\beta} > \dfrac{1}{\alpha\beta}$ と整合している。
    この計算は、ベイズ統計で逆ガンマ分布(正規分布の分散の共役事前分布)の期待値を求めるときや、指数分布のレート(平均の逆数)の推定量 $1/\bar{X}$ の偏りを調べるときにそのまま使われる。
  6. Q6 · 準1級 · 暗記 · 4択

    確率変数 $X$ が確率密度関数 $f(x) = \dfrac{1}{B(\alpha,\beta)}x^{\alpha-1}(1-x)^{\beta-1}\ \ (0<x<1)$ をもつベータ分布 $\mathrm{Be}(\alpha, \beta)$ に従うとき($\alpha>0,\ \beta>0$)、$X$ の期待値 $E[X]$ と分散 $V[X]$ の組として正しいものはどれか。
    • $E[X] = \dfrac{\alpha}{\alpha+\beta},\ V[X] = \dfrac{\alpha\beta}{(\alpha+\beta)^2(\alpha+\beta+1)}$
    • $E[X] = \dfrac{\alpha}{\alpha+\beta},\ V[X] = \dfrac{\alpha\beta}{(\alpha+\beta)^2(\alpha+\beta-1)}$
    • $E[X] = \dfrac{\beta}{\alpha+\beta},\ V[X] = \dfrac{\alpha\beta}{(\alpha+\beta)^2(\alpha+\beta+1)}$
    • $E[X] = \dfrac{\alpha}{\alpha+\beta},\ V[X] = \dfrac{\alpha\beta}{(\alpha+\beta)^3}$

    答え

    $E[X] = \dfrac{\alpha}{\alpha+\beta},\ V[X] = \dfrac{\alpha\beta}{(\alpha+\beta)^2(\alpha+\beta+1)}$

    モーメントの導出

    ベータ関数の定義から $\displaystyle\int_0^1 x^{\alpha+k-1}(1-x)^{\beta-1}dx = B(\alpha+k, \beta)$ なので $$E[X^k] = \frac{B(\alpha+k,\beta)}{B(\alpha,\beta)} = \frac{\Gamma(\alpha+k)\,\Gamma(\alpha+\beta)}{\Gamma(\alpha)\,\Gamma(\alpha+\beta+k)}$$ $k=1$ で $E[X] = \dfrac{\alpha}{\alpha+\beta}$、$k=2$ で $E[X^2] = \dfrac{\alpha(\alpha+1)}{(\alpha+\beta)(\alpha+\beta+1)}$。よって $$V[X] = \frac{\alpha(\alpha+1)}{(\alpha+\beta)(\alpha+\beta+1)} - \frac{\alpha^2}{(\alpha+\beta)^2} = \frac{\alpha\beta}{(\alpha+\beta)^2(\alpha+\beta+1)}$$

    直感とベイズ統計での意味

    $x^{\alpha-1}$ が $1$ に近い側、$(1-x)^{\beta-1}$ が $0$ に近い側の重みを表すので、$\alpha$ が大きいほど分布は $1$ 寄りになり期待値は $\frac{\alpha}{\alpha+\beta}$ となる($\frac{\beta}{\alpha+\beta}$ は $1-X$ の期待値)。分母の $\alpha+\beta+1$ から、$\alpha+\beta$ が大きいほど分散が小さく、分布が期待値の周りに集中することがわかる。
    二項分布の成功確率 $p$ の事前分布を $\mathrm{Be}(\alpha,\beta)$ とし、$n$ 回中 $x$ 回成功を観測すると、事後分布は $\mathrm{Be}(\alpha+x,\ \beta+n-x)$ となる(共役事前分布)。事後平均 $\frac{\alpha+x}{\alpha+\beta+n}$ はこの期待値の公式から直ちに得られる。$\alpha=\beta=1$ は一様分布 $U(0,1)$ で、$E = \frac{1}{2},\ V = \frac{1}{12}$ と一致する。
  7. Q7 · 準1級 · 計算 · 4択

    確率変数 $X$ が確率密度関数 $f(x) = 12x(1-x)^2\ \ (0<x<1)$ をもつとき、$X$ の分散 $V[X]$ として正しいものはどれか。
    • $\dfrac{1}{25}$
    • $\dfrac{1}{5}$
    • $\dfrac{3}{50}$
    • $\dfrac{6}{125}$

    答え

    $\dfrac{1}{25}$

    ベータ分布として読む

    $x^{1}(1-x)^{2} = x^{2-1}(1-x)^{3-1}$ なので、これはベータ分布 $\mathrm{Be}(2, 3)$ である。実際 $B(2,3) = \dfrac{1!\,2!}{4!} = \dfrac{1}{12}$ で、係数 $12$ と一致する。 公式 $V[X] = \dfrac{\alpha\beta}{(\alpha+\beta)^2(\alpha+\beta+1)}$ に $\alpha=2,\ \beta=3$ を代入して $$V[X] = \frac{2\cdot 3}{5^2\cdot 6} = \frac{1}{25}$$

    モーメントからの確認

    $E[X] = \dfrac{2}{5}$、$E[X^2] = \dfrac{\alpha(\alpha+1)}{(\alpha+\beta)(\alpha+\beta+1)} = \dfrac{2\cdot 3}{5\cdot 6} = \dfrac{1}{5}$ なので $$V[X] = \frac{1}{5} - \left(\frac{2}{5}\right)^2 = \frac{5 - 4}{25} = \frac{1}{25}$$ $\frac{1}{5}$ は $E[X^2]$ のまま平均の2乗を引き忘れた値、$\frac{3}{50}$ は分母を $(\alpha+\beta-1)$ にした値、$\frac{6}{125}$ は分母を $(\alpha+\beta)^3$ にした値である。
  8. Q8 · 準1級 · 暗記 · 4択

    標準コーシー分布(位置母数 $0$、尺度母数 $1$)の確率密度関数として正しいものはどれか。
    • $f(x) = \dfrac{1}{\pi(1+x^2)}\ \ (-\infty<x<\infty)$
    • $f(x) = \dfrac{2}{\pi(1+x^2)}\ \ (-\infty<x<\infty)$
    • $f(x) = \dfrac{1}{2\pi(1+x^2)}\ \ (-\infty<x<\infty)$
    • $f(x) = \dfrac{1}{\sqrt{2\pi}(1+x^2)}\ \ (-\infty<x<\infty)$

    答え

    $f(x) = \dfrac{1}{\pi(1+x^2)}\ \ (-\infty<x<\infty)$

    正規化定数の確認

    $\dfrac{d}{dx}\arctan x = \dfrac{1}{1+x^2}$ なので $$\int_{-\infty}^{\infty}\frac{dx}{1+x^2} = \Big[\arctan x\Big]_{-\infty}^{\infty} = \frac{\pi}{2} - \left(-\frac{\pi}{2}\right) = \pi$$ よって係数は $\dfrac{1}{\pi}$ になる。係数 $\frac{2}{\pi}$ は $x>0$ の半分だけで積分した場合(半コーシー分布)の値。一般の位置 $\mu$・尺度 $\sigma$ では $f(x) = \dfrac{\sigma}{\pi\{\sigma^2 + (x-\mu)^2\}}$ となる。

    他の分布との関係

    標準コーシー分布は自由度 $1$ の $t$ 分布と一致する。また、独立な標準正規確率変数 $Z_1, Z_2$ の比 $Z_1/Z_2$ も標準コーシー分布に従う。
    形は正規分布と同じく $0$ を中心とする左右対称の山型だが、裾が $1/x^2$ のオーダーでしか減衰しない(正規分布は $e^{-x^2/2}$)ため、極端な値が出やすい「裾の重い」分布の代表例である。
  9. Q9 · 準1級 · 暗記 · 4択

    確率変数 $X$ が確率密度関数 $f(x) = \dfrac{1}{\pi(1+x^2)}\ \ (-\infty<x<\infty)$ をもつ標準コーシー分布に従うとき、$X$ の期待値と分散についての記述として正しいものはどれか。
    • 期待値も分散も存在しない
    • 期待値は $0$ で、分散は存在しない
    • 期待値は $0$ で、分散は $1$ である
    • 期待値は存在しないが、分散は $1$ である

    答え

    期待値も分散も存在しない

    期待値が存在しない理由

    期待値が存在するためには $E[|X|] < \infty$(絶対収束)が必要である。ところが $$E[|X|] = 2\int_0^\infty \frac{x}{\pi(1+x^2)}\,dx = \frac{1}{\pi}\Big[\log(1+x^2)\Big]_0^\infty = \infty$$ となり発散する。密度は $0$ について対称だが、$\int_{-\infty}^{0}$ と $\int_0^{\infty}$ がそれぞれ $-\infty,\ +\infty$ に発散するため「$\infty - \infty$」となって値が定まらない。対称性から「期待値は $0$」とするのは典型的な誤りである。

    分散も存在しない

    分散は期待値を前提に定義されるので、期待値が存在しなければ分散も存在しない。実際 $E[X^2] = \displaystyle\int_{-\infty}^\infty \frac{x^2}{\pi(1+x^2)}dx$ も、被積分関数が $x\to\pm\infty$ で $1/\pi$ に近づくため発散する。「分散は $1$」は、標準コーシー分布の尺度母数 $1$ を標準正規分布の分散 $1$ と同じものと取り違えた誤りである。 なお、中心の位置を表す指標としては中央値(= 最頻値 = $0$)が、散らばりの指標としては四分位範囲(標準コーシー分布では四分位点が $\pm 1$)が使える。
  10. Q10 · 準1級 · 暗記 · 4択

    $X_1, X_2, \ldots, X_n$ が互いに独立に標準コーシー分布に従うとき、標本平均 $\bar{X} = \dfrac{1}{n}\displaystyle\sum_{i=1}^n X_i$ の分布として正しいものはどれか。
    • $n$ によらず、標準コーシー分布に従う
    • $n$ が大きくなると、正規分布に近づく
    • $n$ が大きくなると、$0$ に確率収束する
    • 尺度母数 $1/\sqrt{n}$ のコーシー分布に従う

    答え

    $n$ によらず、標準コーシー分布に従う

    特性関数による証明

    標準コーシー分布の特性関数は $\varphi(t) = E[e^{itX}] = e^{-|t|}$ である。独立性より $$\varphi_{\bar{X}}(t) = \prod_{i=1}^n \varphi\!\left(\frac{t}{n}\right) = \left(e^{-|t|/n}\right)^n = e^{-|t|}$$ となり、$\bar{X}$ の特性関数は $X_1$ と全く同じである。つまり何個平均しても散らばりは小さくならない。(モーメント母関数は存在しないため、特性関数を用いる。)

    大数の法則・中心極限定理が成り立たない理由

    大数の法則は期待値の存在、中心極限定理は分散の存在(有限)を前提としている。コーシー分布はどちらも持たないため、$\bar{X}$ は $0$ に収束せず、正規分布にも近づかない。 $\bar{X}$ のばらつきが $1/\sqrt{n}$ のオーダーで小さくなるのは分散が有限な場合の話であり、コーシー分布にはあてはまらない。位置母数の推定には、標本平均ではなく標本中央値(一致推定量となる)などが使われる。

← 解説記事一覧に戻る