戻る

統計学入門 第5章の問題 5.1 一様分布

一様分布の問題

問題

いくつかの一様分布に対して、次の計算を行え。
(i) [0,6][0, 6]上の一様分布の密度関数、期待値、分散。
(ii) 同じく、チェビチェフの不等式の成立
(iii) [0,1][0, 1]上の一様分布の歪度、尖度

解答

(i)

密度関数f(x)f(x)は、ある定数をCCとして

f(x)={C(0x6)0(otherwise)\begin{align} f(x) = \left\{ \begin{aligned} C &\quad (0 \leq x \leq 6)\\ 0 &\quad (otherwise) \end{aligned} \right. \end{align}

です。確率の規格化条件により、

1=f(x)dx=C06dx=6C\begin{align*} 1= \int_{-\infty}^{\infty} f(x) dx = C \int_{0}^{6} dx = 6C \end{align*}

CCが満たさなければならないので、C=1/6C = 1/6となることが分かります。よって、求める密度関数は

f(x)={16(0x6)0(otherwise)\begin{align} f(x) = \left\{ \begin{aligned} \frac{1}{6} &\quad (0 \leq x \leq 6)\\ 0 &\quad (otherwise) \end{aligned} \right. \end{align}

です。

期待値μ\mu

μ=06xf(x)dx=1606xdx=3\begin{align*} \mu &= \int_{0}^{6} x f(x) {\rm d}x = \frac{1}{6} \int_{0}^{6} x {\rm d}x = 3 \end{align*}

です。

分散σ2\sigma^2を求めるために2次のモーメントは

E(X2)=06x2f(x)dx=1606x2dx=16[x33]06=12\begin{align*} E(X^2) &= \int_{0}^{6} x^2 f(x) {\rm d}x = \frac{1}{6} \int_{0}^{6} x^2 {\rm d}x = \frac{1}{6} \left[ \frac{x^3}{3} \right]_{0}^{6} = 12 \end{align*}

です。よって、分散は

σ2=E(X2)μ2=1232=129=3\begin{align*} \sigma^2 &= E(X^2) - \mu^2 = 12 - 3^2 = 12 - 9 = 3 \end{align*}

となります。

(ii)

チェビチェフの不等式が

P(Xμkσ)1k2\begin{align} P(|X - \mu| \geq k\sigma) &\leq \frac{1}{k^2} \end{align}

であることを思い出すと、確率変数xxの範囲は

 x33kx33kまたはx3+3k\begin{align*} &~\quad|x - 3| \geq \sqrt{3} k \\ &\Leftrightarrow x \leq 3 - \sqrt{3} k \quad \text{または} \quad x \geq 3 + \sqrt{3} k \end{align*}

となります。この窓に非ゼロの確率含まれるのかどうかで条件を分けることができ、(a) 0<k30 < k \leq \sqrt{3}、(b) k>3k > \sqrt{3}とします。

(a) 0<k30 < k \leq \sqrt{3}で、

P(X33k)=033k16dx+3+3616dx=16(33k)=133k\begin{align} P(|X - 3| \geq \sqrt{3}k) &= \int_0^{3-\sqrt{3}k} \frac{1}{6}\, {\rm d}x + \int_{3 + \sqrt{3}}^6 \frac{1}{6}\, {\rm d}x \notag \\ &= \frac{1}{6}(3 - \sqrt{3}k) \notag \\ &= 1 - \frac{\sqrt{3}}{3} k \end{align}

となります。これをチェビチェフの不等式(3)(3)と比較すると、

133k1k2 1 - \frac{\sqrt{3}}{3} k \leq \frac{1}{k^2}

となるので、うまく式変形して

g(k)=k233k31 g(k) = k^2 - \frac{\sqrt{3}}{3} k^3 \leq 1

となります。g(k)g(k)の最大値を求めます。これをkkで微分すると

g(k)=2k(23k) g^\prime(k) = 2k(2 - \sqrt{3} k)

なので、極大値(兼最大値)はk=2/3k = 2/\sqrt{3}です。この値をg(k)g(k)に代入すればg(2/3)=4/9<1g(2/\sqrt{3}) = 4/9 < 1であり、チェビチェフの不等式を満たします。

一方、(b) k>3k > \sqrt{3}のとき明らかにP(X33k)=0P(|X - 3| \geq \sqrt{3}k) = 0です。なので、0<1/k2(k>3)0 < 1/k^2\, (k > \sqrt{3})なので、自明にチェビチェフの不等式を満たします。

(iii)

歪度・尖度はそれぞれ3次のモーメントα3\alpha_3、4次のモーメントα4\alpha_4を用いてそれぞれ

S=E[(Xμ)3]σ3,K=E[(Xμ)4]σ43, S = \frac{E[(X - \mu)^3]}{\sigma^3}\,, \quad K = \frac{E[(X - \mu)^4]}{\sigma^4} - 3\,,

で定義されます。

まずは、3次のモーメントを計算します。期待値がμ=1/2\mu = 1/2なので、

E[(Xμ)3]=01(x12)3dx=1/21/2t3dt=0 E[(X - \mu)^3] = \int_{0}^{1} \left(x - \frac{1}{2}\right)^3 \, {\rm d}x = \int_{-1/2}^{1/2} t^3 \, {\rm d}t = 0

となります。よって、歪度S=0S = 0です。

次に4次のモーメントを計算します:

E[(Xμ)4]=01(x12)4dx=201/2t4dt=180\begin{align} E[(X - \mu)^4] &= \int_{0}^{1} \left(x - \frac{1}{2}\right)^4 \, {\rm d}x = 2 \int_{0}^{1/2} t^4 \, {\rm d}t = \frac{1}{80} \end{align}

2次のモーメントが1/31/3であることから分散はσ2=1/12\sigma^2 = 1/12です。よって4次のモーメントは

α4=180122=95 \alpha_4 = \frac{1}{80} \cdot 12^2 = \frac{9}{5}

となります。よって尖度は

K=α43=953=65 K = \alpha_4 - 3 = \frac{9}{5} - 3 = -\frac{6}{5}

となります。

※尖度の定義にある3-3はガウス分布の4次のモーメントの値です。なので、K<0K < 0であるので、[0,1][0, 1]の一様分布はガウス分布よりも「鈍く丸っこい」分布であることが分かります。

参考資料

戻る