メインコンテンツへスキップ
教材一覧に戻る
数検1級

確率統計 — 確率分布・期待値分散・回帰相関

最終確認日:

この章で学ぶこと

不確かさを数理的に扱う確率統計を学びます。期待値・分散と代表的な分布、データから関係を読みとる回帰・相関が中心です。

  • 確率変数と確率分布(離散・連続)
  • 期待値と分散
  • 二項分布・正規分布
  • 相関係数と回帰直線(最小二乗法)

ポイント: 期待値は「平均的にどうなるか」、分散は「ばらつきの大きさ」を表す最も基本的な量です。連続分布では「確率密度関数を積分して確率を出す」のがポイント。回帰・相関は2変数データの関係を測ります。

1. 確率変数と期待値・分散

確率変数 XX期待値(平均) E[X]E[X]分散 V[X]V[X] は、離散・連続でそれぞれ次のように定義します。

離散連続(密度 ff)
期待値E[X]=ixipiE[X] = \sum_i x_i p_iE[X]=xf(x)dxE[X] = \displaystyle\int_{-\infty}^{\infty} x f(x)\, dx
分散V[X]=i(xiμ)2piV[X] = \sum_i (x_i - \mu)^2 p_iV[X]=(xμ)2f(x)dxV[X] = \displaystyle\int (x - \mu)^2 f(x)\, dx

分散は V[X]=E[X2](E[X])2V[X] = E[X^2] - (E[X])^2 でも計算できます。標準偏差は σ=V[X]\sigma = \sqrt{V[X]}

例題: さいころの目 XX(161\sim6 が等確率 16\tfrac16)の期待値と分散を求めよ。

E[X]=1+2+3+4+5+66=216=72E[X] = \frac{1+2+3+4+5+6}{6} = \frac{21}{6} = \frac{7}{2} E[X2]=1+4+9+16+25+366=916E[X^2] = \frac{1+4+9+16+25+36}{6} = \frac{91}{6} V[X]=E[X2](E[X])2=916494=18214712=3512V[X] = E[X^2] - (E[X])^2 = \frac{91}{6} - \frac{49}{4} = \frac{182 - 147}{12} = \frac{35}{12}

検算: 916=18212\dfrac{91}{6} = \dfrac{182}{12}494=14712\dfrac{49}{4} = \dfrac{147}{12}、差 35122.92\dfrac{35}{12} \approx 2.92。分散は正で妥当な大きさ。正しい。

2. 連続分布の例

例題: 確率密度 f(x)=2x (0x1)f(x) = 2x\ (0 \le x \le 1)、それ以外 00 の確率変数 XX の期待値を求めよ(全確率が 11 になることも確認せよ)。

まず全確率: 012xdx=[x2]01=1\displaystyle\int_0^1 2x\, dx = [x^2]_0^1 = 1(\circ、密度として正当)。期待値は

E[X]=01x2xdx=012x2dx=[2x33]01=23E[X] = \int_0^1 x \cdot 2x\, dx = \int_0^1 2x^2\, dx = \left[\frac{2x^3}{3}\right]_0^1 = \frac{2}{3}

検算: 密度が右上がり(0x10 \le x \le 1 で大きい xx ほど起こりやすい)なので、平均は区間の中点 0.50.5 より大きい 230.67\tfrac23 \approx 0.67 になるのは妥当。正しい。

3. 二項分布と正規分布

二項分布 B(n,p)B(n, p): 成功確率 pp の試行を nn 回行ったときの成功回数 XX の分布。

P(X=k)=(nk)pk(1p)nk,E[X]=np,V[X]=np(1p)P(X = k) = \binom{n}{k} p^k (1-p)^{n-k}, \qquad E[X] = np, \quad V[X] = np(1-p)

例題: コインを 1010 回投げたときの表の回数 XX(p=12p = \tfrac12)の期待値と分散を求めよ。

E[X]=1012=5,V[X]=101212=104=2.5E[X] = 10 \cdot \frac{1}{2} = 5, \qquad V[X] = 10 \cdot \frac{1}{2} \cdot \frac{1}{2} = \frac{10}{4} = 2.5

検算: 1010 回中だいたい半分の 55 回が表になるのは直感に合う。分散 2.52.5、標準偏差 2.51.58\sqrt{2.5} \approx 1.58。正しい。

正規分布 N(μ,σ2)N(\mu, \sigma^2): 平均 μ\mu・分散 σ2\sigma^2 の連続分布で、密度は

f(x)=12πσexp ⁣((xμ)22σ2)f(x) = \frac{1}{\sqrt{2\pi}\,\sigma} \exp\!\left(-\frac{(x-\mu)^2}{2\sigma^2}\right)

平均を中心とした左右対称の釣鐘形で、XXZ=XμσZ = \dfrac{X - \mu}{\sigma}標準化すると標準正規分布 N(0,1)N(0,1) に従います。二項分布は nn が大きいとき正規分布で近似できます。

大事: 正規分布では「平均 ±\pm 標準偏差」の範囲に約 68%68\%、「平均 ±2σ\pm 2\sigma」に約 95%95\% が入ります(6868-9595-99.799.7 則)。標準化して標準正規分布表を引くのが確率計算の基本手順です。

4. 相関係数と回帰直線

2変数データ (xi,yi)(x_i, y_i) の関係の強さは相関係数

r=(xixˉ)(yiyˉ)(xixˉ)2(yiyˉ)2=sxysxsyr = \frac{\sum (x_i - \bar{x})(y_i - \bar{y})}{\sqrt{\sum (x_i - \bar{x})^2}\,\sqrt{\sum (y_i - \bar{y})^2}} = \frac{s_{xy}}{s_x s_y}

で測ります(1r1-1 \le r \le 1)。rr11 に近いほど右上がりの直線関係が強い。最小二乗法による回帰直線 y=ax+by = a x + b の傾きは

a=sxysx2=(xixˉ)(yiyˉ)(xixˉ)2,b=yˉaxˉa = \frac{s_{xy}}{s_x^2} = \frac{\sum(x_i-\bar x)(y_i-\bar y)}{\sum(x_i-\bar x)^2}, \qquad b = \bar{y} - a\bar{x}

例題: データ (1,2),(2,2),(3,4),(4,4)(1,2),(2,2),(3,4),(4,4) の回帰直線を最小二乗法で求めよ。

平均は xˉ=1+2+3+44=2.5\bar{x} = \dfrac{1+2+3+4}{4} = 2.5yˉ=2+2+4+44=3\bar{y} = \dfrac{2+2+4+4}{4} = 3。偏差の積和と平方和を計算する。

(xixˉ)(yiyˉ)=(1.5)(1)+(0.5)(1)+(0.5)(1)+(1.5)(1)=1.5+0.5+0.5+1.5=4\sum(x_i-\bar x)(y_i-\bar y) = (-1.5)(-1) + (-0.5)(-1) + (0.5)(1) + (1.5)(1) = 1.5 + 0.5 + 0.5 + 1.5 = 4 (xixˉ)2=(1.5)2+(0.5)2+(0.5)2+(1.5)2=2.25+0.25+0.25+2.25=5\sum(x_i-\bar x)^2 = (-1.5)^2 + (-0.5)^2 + (0.5)^2 + (1.5)^2 = 2.25 + 0.25 + 0.25 + 2.25 = 5

傾き a=45=0.8a = \dfrac{4}{5} = 0.8、切片 b=30.82.5=32=1b = 3 - 0.8 \cdot 2.5 = 3 - 2 = 1。よって回帰直線は

y=0.8x+1y = 0.8\,x + 1

検算: 予測値は x=1,2,3,4x=1,2,3,41.8,2.6,3.4,4.21.8, 2.6, 3.4, 4.2。実測 2,2,4,42,2,4,4 との残差は 0.2,0.6,0.6,0.20.2, -0.6, 0.6, -0.2 で和が 00(最小二乗の性質)。また直線は点 (xˉ,yˉ)=(2.5,3)(\bar x, \bar y) = (2.5, 3) を通る(0.82.5+1=30.8\cdot2.5+1 = 3)。正しい。

どう問われるか

  • 一次では「期待値・分散の計算」「二項分布・正規分布の確率」が頻出。標準化と確率表の利用が鍵。
  • 二次では「連続分布の期待値・分散を積分で求める」「相関係数・回帰直線を計算する」「仮説検定・推定の基礎」が問われます。

よくまちがえるところ

確率統計は、言葉の定義の混同がそのまま計算の誤りになる分野です。

まちがい直し方
独立と排反を混同する同時に起こらないから独立だとする排反は同時に起こらないこと、独立は同時に起こる確率が各確率の積になること。別の概念
分散の加法を無条件に使うV[X + Y] = V[X] + V[Y] をいつでも使う独立(または無相関)のときに限る。一般には共分散の2倍が加わる
定数項を分散に足すV[2X + 3] を 4V[X] + 3 とするV[aX + b] = a² V[X]。定数 b は分散を変えない
相関を因果と読む相関係数が 0.9 だから x が y の原因だとする相関は直線的な関係の強さで、因果を意味しない

もっとも根が深いのは独立と排反の混同です。2つの事象が排反(同時に起こらない)であれば、同時に起こる確率は 0 です。一方、独立であれば同時に起こる確率は各確率の積になります。どちらの確率も正であるとき、この2つは両立しません。積が正なのに 0 になることはないからです。つまり排反な2事象は(確率が正なら)独立ではないのです。「関係がない」という日常語で両者を一括りにしないでください。

もう1つ、相関係数が 0 に近いことは「関係がない」ことを意味しません。相関係数が測るのは直線的な関係の強さだけです。放物線のうえに並んだデータのように、明確な関係があっても相関係数がほぼ 0 になる場合があります。散布図を描いて形を見るのが基本です。

たしかめ方

計算結果は次の観点で検算します。どれも短時間で確認できます。

見るところ成り立つべきこと破れていたら
確率の合計離散なら合計が 1、連続なら密度の全積分が 1場合の数え落とし、または規格化定数の誤り
分散の符号分散は必ず 0 以上E[X²] と (E[X])² の順序が逆
期待値の位置期待値は取りうる値の範囲の中に入る重みづけの取りちがえ
相関係数の範囲-1 以上 1 以下分母の標準偏差の計算ミス
回帰直線x と y それぞれの平均を座標とする点を通る切片の計算ミス
二項分布分散 np(1-p) は期待値 np より小さい1 - p を掛け忘れている

密度関数の例で確かめます。f(x) = 2x(0 ≦ x ≦ 1)なら全積分は 1 で密度として正当、期待値は 2/3 です。密度が右上がりなので平均が区間の中点 0.5 より大きくなるのは妥当で、値の位置からも整合が取れます。

回帰直線の例では、データ (1,2)、(2,2)、(3,4)、(4,4) から y = 0.8x + 1 が得られました。x の平均 2.5、y の平均 3 を代入すると 0.8 × 2.5 + 1 = 3 となり、確かに平均の点を通ります。残差の合計が 0 になることも、最小二乗法の性質からくる検算になります。

大事: 正規分布の確率を求めるときは、標準化 Z = (X - 平均)/標準偏差 を経由して標準正規分布表を引きます。標準偏差で割るのであって分散で割るのではありません。分散のまま割ると単位が合わず、値が大きくずれます。

自分でチェック

  • 独立と排反のちがいを、確率の式で説明できる
  • 分散の加法が使える条件を言える
  • V[aX + b] = a² V[X] と書ける
  • 相関係数が 0 でも関係がある場合があると説明できる
  • 標準化で割るのが標準偏差だと言える

まとめ

  • 期待値は平均、分散 V[X]=E[X2](E[X])2V[X] = E[X^2] - (E[X])^2、標準偏差 σ=V\sigma = \sqrt{V}
  • 二項分布 E=np, V=np(1p)E = np,\ V = np(1-p)、正規分布は標準化して表を引く
  • 相関係数 r=sxysxsyr = \dfrac{s_{xy}}{s_x s_y}1r1-1 \le r \le 1
  • 回帰直線は最小二乗法、傾き a=sxysx2a = \dfrac{s_{xy}}{s_x^2}、点 (xˉ,yˉ)(\bar x, \bar y) を通る

次章では、数値解析・初等整数論と二次対策を学び、本教材を仕上げます。

この教材きょうざいやくちましたか?

このしょうれんしゅう

ほかのぶんれんしゅうもチェック