第6・7章の確率変数は、人数や回数のように「とびとびの値」をとるものでした。この章では、身長や重さのように切れ目なく値が変わる量を確率変数として扱い、確率を曲線の下の**面積**(定積分)で表します。その代表が、左右対称な釣り鐘形の**正規分布**です。標準化と正規分布表を使って確率を求める方法を身につけ、テストの得点や製品の重さのようなデータに当てはめます。最後に、第7章の二項分布が、回数が多いと正規分布でよく近似できることを確かめます。定積分(微分積分の分野 第8章)を使います。
とびとびの値から、切れ目のない値へ
第6・7章で扱った確率変数は、さいころの目、芽が出た粒の数、表が出た回数のように、0 , 0, 0 , 1 , \ 1, 1 , 2 , \ 2, 2 , … \ \ldots … ととびとびの値 をとるものでした。このような確率変数を離散型 といいます。確率分布は表で書け、確率は表の値を足して求めました。
一方、身長・重さ・時間のような量は、165.3 165.3 165.3 cm、165.31 165.31 165.31 cm、…と、いくらでも細かい値をとりえます。この章では、こうした切れ目なく変わる値 をとる確率変数を扱います。表の値を「足す」代わりに、曲線の下の面積、つまり定積分 で確率を表すのが新しい点です(定積分は微分積分の分野 第8章・第9章)。
連続型確率変数と確率密度関数
第1章のヒストグラムを思い出してください。縦軸を「相対度数 ÷ 階級の幅」にしておくと、柱の面積がちょうど相対度数になり、柱の面積の合計は 1 1 1 になります。データを増やしながら階級の幅を細かくしていくと、柱の上の端はなめらかな曲線に近づきます。
階級の幅 0.5
0
1
2
0.25
0.5
0.75
階級の幅 0.25
0
1
2
0.25
0.5
0.75
赤い曲線:y = f(x) 柱の高さ=区間の確率 ÷ 幅
連続型では、ちょうど 1 1 1 つの値をとる確率は 0 0 0 です。P ( X = a ) = ∫ a a f ( x ) d x = 0 P(X = a) = \displaystyle\int_a^a f(x)\,dx = 0 P ( X = a ) = ∫ a a f ( x ) d x = 0 だからです。そのため P ( a ≦ X ≦ b ) P(a \leqq X \leqq b) P ( a ≦ X ≦ b ) と P ( a < X < b ) P(a < X < b) P ( a < X < b ) は同じ値になり、「以上」と「より大きい」を区別しなくてかまいません。f ( a ) f(a) f ( a ) の値そのものは確率ではなく、「a a a のあたりに値がどれだけ集まっているか」を表す密度 です。
人口密度の地図 にたとえると分かりやすくなります。地図の上の 1 点を指さして「ここに何人住んでいるか」と聞かれても、点には広さがないので答えは 0 0 0 人です。けれども「この区域には何人か」なら、区域の広さに人口密度をかけて足し合わせれば求まります。確率密度関数も同じで、1 点の確率は 0 0 0 でも、区間を決めればその上の面積として確率が決まります。密度が高い場所ほど、同じ幅の区間に多くの確率が集まっています。
連続型確率変数では、確率は確率密度関数 f ( x ) f(x) f ( x ) のグラフの下の面積 ∫ a b f ( x ) d x \displaystyle\int_a^b f(x)\,dx ∫ a b f ( x ) d x で表され、全体の面積は 1 1 1 、 1 点の確率は 0 0 0 になるということです。
例題1:確率密度関数と確率
確率変数 X X X のとる値の範囲が 0 0 0 ≦ X {}\leqq X ≦ X ≦ 2 {}\leqq 2 ≦ 2 で、確率密度関数が f ( x ) f(x) f ( x ) = a x ( 2 − x ) {}= ax(2 - x) = a x ( 2 − x ) (a a a は正の定数)です。定数 a a a の値と、確率 P ( 0 ≦ X ≦ 1 2 ) P\left(0 \leqq X \leqq \dfrac{1}{2}\right) P ( 0 ≦ X ≦ 2 1 ) を求めなさい。
【解答】
全体の面積が 1 1 1 なので
∫ 0 2 a x ( 2 − x ) d x \displaystyle \int_0^2 ax(2 - x)\,dx ∫ 0 2 a x ( 2 − x ) d x = a [ x 2 − x 3 3 ] 0 2 \displaystyle {}= a\Big[\,x^2 - \frac{x^3}{3}\,\Big]_0^2 = a [ x 2 − 3 x 3 ] 0 2 = a ( 4 − 8 3 ) \displaystyle {}= a\left(4 - \frac{8}{3}\right) = a ( 4 − 3 8 ) = 4 3 a \displaystyle {}= \frac{4}{3}a = 3 4 a = 1 \displaystyle {}= 1 = 1 より a = 3 4 ‾ \underline{a = \dfrac{3}{4}} a = 4 3 です。0 0 0 ≦ x {}\leqq x ≦ x ≦ 2 {}\leqq 2 ≦ 2 で f ( x ) f(x) f ( x ) ≧ 0 {}\geqq 0 ≧ 0 も満たしています。
P ( 0 ≦ X ≦ 1 2 ) \displaystyle P\left(0 \leqq X \leqq \frac{1}{2}\right) P ( 0 ≦ X ≦ 2 1 ) = 3 4 ∫ 0 1 2 ( 2 x − x 2 ) d x \displaystyle {}= \frac{3}{4}\int_0^{\frac{1}{2}} (2x - x^2)\,dx = 4 3 ∫ 0 2 1 ( 2 x − x 2 ) d x = 3 4 [ x 2 − x 3 3 ] 0 1 2 \displaystyle {}= \frac{3}{4}\Big[\,x^2 - \frac{x^3}{3}\,\Big]_0^{\frac{1}{2}} = 4 3 [ x 2 − 3 x 3 ] 0 2 1 = 3 4 ( 1 4 − 1 24 ) \displaystyle {}= \frac{3}{4}\left(\frac{1}{4} - \frac{1}{24}\right) = 4 3 ( 4 1 − 24 1 ) = 3 4 × 5 24 \displaystyle {}= \frac{3}{4} \times \frac{5}{24} = 4 3 × 24 5 = 5 32 ‾ \displaystyle {}= \underline{\frac{5}{32}} = 32 5
連続型確率変数の期待値と分散
離散型の期待値は「値 × 確率」の和 ∑ x k p k \sum x_k p_k ∑ x k p k でした。連続型では、幅 d x dx d x の小さな区間の確率が f ( x ) d x f(x)\,dx f ( x ) d x なので、和を積分に置きかえます。
V ( X ) V(X) V ( X ) = E ( X 2 ) {}= E(X^2) = E ( X 2 ) − m 2 {}- m^2 − m 2 は、( x − m ) 2 (x - m)^2 ( x − m ) 2 = x 2 {}= x^2 = x 2 − 2 m x {}- 2mx − 2 m x + m 2 {}+ m^2 + m 2 を展開して積分すると
∫ α β x 2 f ( x ) d x \displaystyle \int_\alpha^\beta x^2 f(x)\,dx ∫ α β x 2 f ( x ) d x − 2 m ∫ α β x f ( x ) d x \displaystyle {}- 2m \int_\alpha^\beta x f(x)\,dx − 2 m ∫ α β x f ( x ) d x + m 2 ∫ α β f ( x ) d x \displaystyle {}+ m^2 \int_\alpha^\beta f(x)\,dx + m 2 ∫ α β f ( x ) d x = E ( X 2 ) \displaystyle {}= E(X^2) = E ( X 2 ) − 2 m ⋅ m \displaystyle {}- 2m \cdot m − 2 m ⋅ m + m 2 ⋅ 1 \displaystyle {}+ m^2 \cdot 1 + m 2 ⋅ 1
となることから分かります。離散型のときの式変形の ∑ \sum ∑ が ∫ \int ∫ に替わっただけです。
離散型と連続型の関係は、次のように対応させると覚えやすくなります。
連続型の期待値・分散は、離散型の ∑ ( 値 ) × ( 確率 ) \sum (\text{値}) \times (\text{確率}) ∑ ( 値 ) × ( 確率 ) を ∫ x f ( x ) d x \int x f(x)\,dx ∫ x f ( x ) d x のように積分に置きかえて求め、V V V = E ( X 2 ) {}= E(X^2) = E ( X 2 ) − m 2 {}- m^2 − m 2 などの性質はそのまま使えるということです。
例題2:期待値と分散
例題1 の確率変数 X X X (確率密度関数 f ( x ) f(x) f ( x ) = 3 4 x ( 2 − x ) {}= \dfrac{3}{4}x(2 - x) = 4 3 x ( 2 − x ) ,0 0 0 ≦ X {}\leqq X ≦ X ≦ 2 {}\leqq 2 ≦ 2 ) の期待値、分散、標準偏差をそれぞれ求めなさい。
【解答】
E ( X ) \displaystyle E(X) E ( X ) = 3 4 ∫ 0 2 ( 2 x 2 − x 3 ) d x \displaystyle {}= \frac{3}{4}\int_0^2 (2x^2 - x^3)\,dx = 4 3 ∫ 0 2 ( 2 x 2 − x 3 ) d x = 3 4 [ 2 x 3 3 − x 4 4 ] 0 2 \displaystyle {}= \frac{3}{4}\Big[\,\frac{2x^3}{3} - \frac{x^4}{4}\,\Big]_0^2 = 4 3 [ 3 2 x 3 − 4 x 4 ] 0 2 = 3 4 ( 16 3 − 4 ) \displaystyle {}= \frac{3}{4}\left(\frac{16}{3} - 4\right) = 4 3 ( 3 16 − 4 ) = 1 ‾ \displaystyle {}= \underline{1} = 1 分布曲線は直線 x x x = 1 {}= 1 = 1 について対称なので、期待値が 1 1 1 になるのは納得できます。
E ( X 2 ) \displaystyle E(X^2) E ( X 2 ) = 3 4 ∫ 0 2 ( 2 x 3 − x 4 ) d x \displaystyle {}= \frac{3}{4}\int_0^2 (2x^3 - x^4)\,dx = 4 3 ∫ 0 2 ( 2 x 3 − x 4 ) d x = 3 4 [ x 4 2 − x 5 5 ] 0 2 \displaystyle {}= \frac{3}{4}\Big[\,\frac{x^4}{2} - \frac{x^5}{5}\,\Big]_0^2 = 4 3 [ 2 x 4 − 5 x 5 ] 0 2 = 3 4 ( 8 − 32 5 ) \displaystyle {}= \frac{3}{4}\left(8 - \frac{32}{5}\right) = 4 3 ( 8 − 5 32 ) = 6 5 \displaystyle {}= \frac{6}{5} = 5 6 より
V ( X ) \displaystyle V(X) V ( X ) = 6 5 \displaystyle {}= \frac{6}{5} = 5 6 − 1 2 \displaystyle {}- 1^2 − 1 2 = 1 5 ‾ , \displaystyle {}= \underline{\frac{1}{5}}, = 5 1 , σ ( X ) \displaystyle \sigma(X) σ ( X ) = 1 5 \displaystyle {}= \sqrt{\frac{1}{5}} = 5 1 = 5 5 ‾ \displaystyle {}= \underline{\frac{\sqrt{5}}{5}} = 5 5
正規分布
連続型の確率分布の中で最も大切なのが、次の正規分布 です。第6章の小話で、さいころの個数を増やすと目の和の分布が釣り鐘形になることを紹介しました。その釣り鐘形を式で表したものです。
N N N は正規分布を表す英語 normal distribution の頭文字で、N ( m , σ 2 ) N(m,\ \sigma^2) N ( m , σ 2 ) の 2 つ目には標準偏差の 2 乗、つまり分散を書きます。N ( 50 , 10 2 ) N(50,\ 10^2) N ( 50 , 1 0 2 ) なら平均 50 50 50 、 標準偏差 10 10 10 です。E ( X ) E(X) E ( X ) = m {}= m = m ,V ( X ) V(X) V ( X ) = σ 2 {}= \sigma^2 = σ 2 となることの証明には高校の範囲を少し超える積分が必要なので、厳密定義のページで扱います。e e e は微分積分の分野 第11章で学ぶ数で、ここでは式の形よりもグラフの性質を使えれば十分です。
-5
-4
-3
-2
-1
0
1
2
3
4
5
6
N(0, 1)
N(3, 1)
N(0, 2²)
服の型紙 にたとえると、正規分布の曲線はたった 1 枚の型紙です。平均 m m m は型紙を布のどこに置くかを、標準偏差 σ \sigma σ は型紙を何倍に拡大するかを決めます。拡大すると横に広がる分だけ高さは低くなり、面積はいつも 1 1 1 に保たれます。どの正規分布も同じ型紙から作られているので、「m m m から σ \sigma σ の何倍以内か」だけで確率が決まる、というのが性質 4 の意味です。
m−3σ
m−2σ
m−σ
m
m+σ
m+2σ
m+3σ
約 68.3%
約 95.4%
約 99.7%
正規分布 N ( m , σ 2 ) N(m,\ \sigma^2) N ( m , σ 2 ) は平均 m m m を中心に左右対称な釣り鐘形で、m m m が位置を、σ \sigma σ が広がりを決め、m m m ± σ {}\pm \sigma ± σ ,m m m ± 2 σ {}\pm 2\sigma ± 2 σ ,m m m ± 3 σ {}\pm 3\sigma ± 3 σ に入る確率は約 0.683 0.683 0.683 ,0.954 0.954 0.954 ,0.997 0.997 0.997 だということです。
例題3:σ の何倍かで確率を読む
ある地域の高校 3 年生の男子の身長 X X X cm は、正規分布 N ( 170 , 6 2 ) N(170,\ 6^2) N ( 170 , 6 2 ) に従うものとします。公式3 の表の値を使って、次の確率を求めなさい。
(1) P ( 158 ≦ X ≦ 182 ) P(158 \leqq X \leqq 182) P ( 158 ≦ X ≦ 182 )
(2) P ( X ≧ 176 ) P(X \geqq 176) P ( X ≧ 176 )
上の (1)、(2) の値をそれぞれ求めなさい。
【解答】
(1) 158 158 158 = 170 {}= 170 = 170 − 2 × 6 {}- 2 \times 6 − 2 × 6 ,182 182 182 = 170 {}= 170 = 170 + 2 × 6 {}+ 2 \times 6 + 2 × 6 なので、m m m ± 2 σ {}\pm 2\sigma ± 2 σ の範囲です。P ( 158 ≦ X ≦ 182 ) ≒ 0.954 ‾ \underline{P(158 \leqq X \leqq 182) \fallingdotseq 0.954} P ( 158 ≦ X ≦ 182 ) ≒ 0.954
(2) 176 176 176 = 170 {}= 170 = 170 + 6 {}+ 6 + 6 は m m m + σ {}+ \sigma + σ です。m m m ± σ {}\pm \sigma ± σ の外側の確率は 1 1 1 − 0.683 {}- 0.683 − 0.683 = 0.317 {}= 0.317 = 0.317 で、分布は左右対称なので、その半分が 176 176 176 以上の側にあります。
P ( X ≧ 176 ) \displaystyle P(X \geqq 176) P ( X ≧ 176 ) ≒ 0.317 2 \displaystyle {}\fallingdotseq \frac{0.317}{2} ≒ 2 0.317 ≒ 0.16 ‾ \displaystyle {}\fallingdotseq \underline{0.16} ≒ 0.16 第3章の小話で「正規分布なら偏差値 60 60 60 以上は上位約 16 16 16 %」と紹介したのは、この計算です(偏差値 60 60 60 は平均 + 1 + 1 + 1 標準偏差)。
標準化と正規分布表
m m m ± σ {}\pm \sigma ± σ のようなきりのよい範囲以外の確率は、積分を手で計算することができません(e − x 2 e^{-x^2} e − x 2 の不定積分は、知っている関数では書けないことが分かっています)。そこで、1 つの正規分布についてだけ数値の表を作っておき、どの正規分布もその表に直して読みます。
Z Z Z は第3章 公式6、第6章で学んだ標準化 そのものです。E ( Z ) E(Z) E ( Z ) = m − m σ {}= \dfrac{m - m}{\sigma} = σ m − m = 0 {}= 0 = 0 ,σ ( Z ) \sigma(Z) σ ( Z ) = σ σ {}= \dfrac{\sigma}{\sigma} = σ σ = 1 {}= 1 = 1 となり、さらに正規分布の形も保たれます(最後の一文の a a a = 1 σ {}= \dfrac{1}{\sigma} = σ 1 ,b b b = − m σ {}= -\dfrac{m}{\sigma} = − σ m の場合)。
正規分布表は 0 0 0 から u u u までの面積だけを載せています。それ以外の範囲は、分布が z z z = 0 {}= 0 = 0 について対称であることと、右半分の面積が 0.5 0.5 0.5 であることを組み合わせて求めます。
P ( Z ≧ u ) P(Z \geqq u) P ( Z ≧ u ) = 0.5 {}= 0.5 = 0.5 − p ( u ) {}- p(u) − p ( u )
P ( − u ≦ Z ≦ 0 ) P(-u \leqq Z \leqq 0) P ( − u ≦ Z ≦ 0 ) = p ( u ) {}= p(u) = p ( u )
P ( − a ≦ Z ≦ b ) P(-a \leqq Z \leqq b) P ( − a ≦ Z ≦ b ) = p ( a ) {}= p(a) = p ( a ) + p ( b ) {}+ p(b) + p ( b ) (a , a, a , b \ b b > 0 {}> 0 > 0 )
P ( a ≦ Z ≦ b ) P(a \leqq Z \leqq b) P ( a ≦ Z ≦ b ) = p ( b ) {}= p(b) = p ( b ) − p ( a ) {}- p(a) − p ( a ) (0 0 0 ≦ a {}\leqq a ≦ a < b {}< b < b )
-3
-2
-1
1
2
3
0
u
p(u)
p(u) = P(0 ≦ Z ≦ u)
右半分の面積は 0.5
**世界共通の時刻(協定世界時)**を思い浮かべてください。東京・ロンドン・ニューヨークの時刻表をそれぞれ用意しなくても、どの都市の時刻も協定世界時に直してしまえば、1 枚の表で比べられます。標準化も同じで、平均も標準偏差もばらばらな正規分布を、すべて「平均 0 0 0 ・ 標準偏差 1 1 1 」 の Z Z Z に直すことで、たった 1 枚の正規分布表で確率が読めるようになります。
正規分布の確率は、Z Z Z = X − m σ {}= \dfrac{X - m}{\sigma} = σ X − m で標準正規分布 N ( 0 , 1 ) N(0,\ 1) N ( 0 , 1 ) に直し、対称性と右半分の面積 0.5 0.5 0.5 を使って正規分布表 p ( u ) p(u) p ( u ) を読めば求まるということです。
例題4:正規分布表を使う
確率変数 X X X が正規分布 N ( 50 , 10 2 ) N(50,\ 10^2) N ( 50 , 1 0 2 ) に従うとき、公式4 の正規分布表を使って、次の確率を求めなさい。
(1) P ( X ≧ 62 ) P(X \geqq 62) P ( X ≧ 62 )
(2) P ( 38 ≦ X ≦ 55 ) P(38 \leqq X \leqq 55) P ( 38 ≦ X ≦ 55 )
上の (1)、(2) の値をそれぞれ求めなさい。
【解答】
Z Z Z = X − 50 10 {}= \dfrac{X - 50}{10} = 10 X − 50 とおくと、Z Z Z は N ( 0 , 1 ) N(0,\ 1) N ( 0 , 1 ) に従います。
(1) X X X = 62 {}= 62 = 62 のとき Z Z Z = 62 − 50 10 {}= \dfrac{62 - 50}{10} = 10 62 − 50 = 1.2 {}= 1.2 = 1.2 なので
P ( X ≧ 62 ) \displaystyle P(X \geqq 62) P ( X ≧ 62 ) = P ( Z ≧ 1.2 ) \displaystyle {}= P(Z \geqq 1.2) = P ( Z ≧ 1.2 ) = 0.5 \displaystyle {}= 0.5 = 0.5 − p ( 1.2 ) \displaystyle {}- p(1.2) − p ( 1.2 ) = 0.5 \displaystyle {}= 0.5 = 0.5 − 0.3849 \displaystyle {}- 0.3849 − 0.3849 = 0.1151 ‾ \displaystyle {}= \underline{0.1151} = 0.1151 (2) X X X = 38 {}= 38 = 38 のとき Z Z Z = − 1.2 {}= -1.2 = − 1.2 、X X X = 55 {}= 55 = 55 のとき Z Z Z = 0.5 {}= 0.5 = 0.5 なので
P ( 38 ≦ X ≦ 55 ) \displaystyle P(38 \leqq X \leqq 55) P ( 38 ≦ X ≦ 55 ) = P ( − 1.2 ≦ Z ≦ 0.5 ) \displaystyle {}= P(-1.2 \leqq Z \leqq 0.5) = P ( − 1.2 ≦ Z ≦ 0.5 ) = p ( 1.2 ) \displaystyle {}= p(1.2) = p ( 1.2 ) + p ( 0.5 ) \displaystyle {}+ p(0.5) + p ( 0.5 ) = 0.3849 \displaystyle {}= 0.3849 = 0.3849 + 0.1915 \displaystyle {}+ 0.1915 + 0.1915 = 0.5764 ‾ \displaystyle {}= \underline{0.5764} = 0.5764
データの分布に正規分布を当てはめる
テストの得点、同じ機械で作った製品の重さ、測定の誤差などのヒストグラムは、釣り鐘形に近くなることがよくあります。そのようなデータでは、平均値 x ‾ \overline{x} x と標準偏差 s s s をもつ正規分布 N ( x ‾ , s 2 ) N(\overline{x},\ s^2) N ( x , s 2 ) で分布を近似し、ある範囲に入る人数や、上位何%の境目を見積もることができます。
2 は、1 と逆の向きに表を使う読み方です。確率から u u u を探し、X X X = m {}= m = m + σ Z {}+ \sigma Z + σ Z で元の値に戻します。
実際のデータがぴったり正規分布になることはありません。得点には満点があり、身長がマイナスになることもありません。それでも、分布が左右対称な山形に近ければ、正規分布は手軽で見積もりのよい模型 になります。縮尺つきの地図 で、道の長さを定規で測れば実際の距離がおよそ分かるのと同じです。地図は道の細かい曲がりまでは描いていませんが、「だいたい何 km か」を知るには十分に役立ちます。
釣り鐘形のデータは正規分布 N ( x ‾ , s 2 ) N(\overline{x},\ s^2) N ( x , s 2 ) で近似でき、「範囲に入る人数=全体 × 確率」「上位 r r r の境目=表で p ( u ) p(u) p ( u ) = 0.5 {}= 0.5 = 0.5 − r {}- r − r となる u u u を読んで m m m + u σ {}+ u\sigma + u σ 」 と見積もれるということです。
例題5:テストの得点の見積もり
ある試験を 500 500 500 人が受け、得点の分布は正規分布 N ( 58 , 12 2 ) N(58,\ 12^2) N ( 58 , 1 2 2 ) で近似できるとします。公式4 の正規分布表を使って、次の値を見積もりなさい。
(1) 70 70 70 点以上の人数(小数第1位を四捨五入して整数で)
(2) 上位 2.5 2.5 2.5 % に入るための得点の境目(小数第1位を四捨五入して整数で)
上の (1)、(2) の値をそれぞれ求めなさい。
【解答】
得点を X X X 、Z Z Z = X − 58 12 {}= \dfrac{X - 58}{12} = 12 X − 58 とします。
(1) X X X = 70 {}= 70 = 70 のとき Z Z Z = 1 {}= 1 = 1 なので
P ( X ≧ 70 ) \displaystyle P(X \geqq 70) P ( X ≧ 70 ) = 0.5 \displaystyle {}= 0.5 = 0.5 − p ( 1.0 ) \displaystyle {}- p(1.0) − p ( 1.0 ) = 0.5 \displaystyle {}= 0.5 = 0.5 − 0.3413 \displaystyle {}- 0.3413 − 0.3413 = 0.1587 \displaystyle {}= 0.1587 = 0.1587 500 × 0.1587 500 \times 0.1587 500 × 0.1587 = 79.35 {}= 79.35 = 79.35 より、およそ 79 人 ‾ \underline{79 \text{ 人}} 79 人 です。
(2) P ( Z ≧ u ) P(Z \geqq u) P ( Z ≧ u ) = 0.025 {}= 0.025 = 0.025 、 つまり p ( u ) p(u) p ( u ) = 0.5 {}= 0.5 = 0.5 − 0.025 {}- 0.025 − 0.025 = 0.475 {}= 0.475 = 0.475 となる u u u は、表から u u u = 1.96 {}= 1.96 = 1.96 です。
c \displaystyle c c = 58 \displaystyle {}= 58 = 58 + 1.96 × 12 \displaystyle {}+ 1.96 \times 12 + 1.96 × 12 = 58 \displaystyle {}= 58 = 58 + 23.52 \displaystyle {}+ 23.52 + 23.52 = 81.52 \displaystyle {}= 81.52 = 81.52 より、上位 2.5 2.5 2.5 % の境目はおよそ 82 点 ‾ \underline{82 \text{ 点}} 82 点 です。
二項分布の正規分布による近似
第7章の二項分布 B ( n , p ) B(n,\ p) B ( n , p ) は、n n n が大きくなると確率の計算が大変になります。ところが n n n が大きいとき、二項分布のグラフは釣り鐘形になり、同じ平均・標準偏差をもつ正規分布の曲線にぴったり重なってきます。
n p np n p ,n p q npq n pq は第7章 公式2 の期待値と分散です。「同じ期待値・同じ分散の正規分布で置きかえる」と覚えれば十分です。
0
2
4
6
8
10
12
14
16
18
20
B(20, 1/2)
N(10, 5)
点描画 の絵を思い出してください。スーラの絵のように小さな色の点を並べた絵は、近くで見ると点の集まりですが、離れて見るとなめらかな色の広がりに見えます。二項分布の確率も、回数ごとの「とびとびの点」ですが、n n n が大きいと点が細かく密に並び、全体としてなめらかな正規分布の曲線に見えてくるのです。
第7章の厳密定義のページで、B ( 20 , 1 2 ) B\left(20,\ \dfrac{1}{2}\right) B ( 20 , 2 1 ) で期待値から標準偏差の 2 2 2 倍以内(6 6 6 ≦ X {}\leqq X ≦ X ≦ 14 {}\leqq 14 ≦ 14 ) に入る確率は約 0.959 0.959 0.959 だと計算しました。正規分布なら公式3 の約 0.954 0.954 0.954 で、n n n = 20 {}= 20 = 20 でもかなり近い値です。
回数 n n n が大きい二項分布 B ( n , p ) B(n,\ p) B ( n , p ) は、期待値 n p np n p ・ 分散 n p q npq n pq が同じ正規分布 N ( n p , n p q ) N(np,\ npq) N ( n p , n pq ) で近似でき、Z Z Z = X − n p n p q {}= \dfrac{X - np}{\sqrt{npq}} = n pq X − n p として正規分布表で確率が求まるということです。
例題6:硬貨 400 回の表の回数
硬貨を 400 400 400 回投げるとき、表が出る回数を X X X とします。正規分布で近似して、公式4 の正規分布表を使い、次の確率を求めなさい。
(1) P ( X ≧ 220 ) P(X \geqq 220) P ( X ≧ 220 )
(2) P ( 190 ≦ X ≦ 215 ) P(190 \leqq X \leqq 215) P ( 190 ≦ X ≦ 215 )
上の (1)、(2) の値をそれぞれ求めなさい。
【解答】
X X X は B ( 400 , 1 2 ) B\left(400,\ \dfrac{1}{2}\right) B ( 400 , 2 1 ) に従うので
n p \displaystyle np n p = 400 × 1 2 \displaystyle {}= 400 \times \frac{1}{2} = 400 × 2 1 = 200 , \displaystyle {}= 200, = 200 , n p q \displaystyle \sqrt{npq} n pq = 400 × 1 2 × 1 2 \displaystyle {}= \sqrt{400 \times \frac{1}{2} \times \frac{1}{2}} = 400 × 2 1 × 2 1 = 10 \displaystyle {}= 10 = 10 n n n = 400 {}= 400 = 400 は大きいので、Z Z Z = X − 200 10 {}= \dfrac{X - 200}{10} = 10 X − 200 は近似的に N ( 0 , 1 ) N(0,\ 1) N ( 0 , 1 ) に従います。
(1) X X X = 220 {}= 220 = 220 のとき Z Z Z = 2 {}= 2 = 2 なので
P ( X ≧ 220 ) \displaystyle P(X \geqq 220) P ( X ≧ 220 ) ≒ P ( Z ≧ 2 ) \displaystyle {}\fallingdotseq P(Z \geqq 2) ≒ P ( Z ≧ 2 ) = 0.5 \displaystyle {}= 0.5 = 0.5 − 0.4772 \displaystyle {}- 0.4772 − 0.4772 = 0.0228 ‾ \displaystyle {}= \underline{0.0228} = 0.0228 (2) X X X = 190 {}= 190 = 190 のとき Z Z Z = − 1 {}= -1 = − 1 、X X X = 215 {}= 215 = 215 のとき Z Z Z = 1.5 {}= 1.5 = 1.5 なので
P ( 190 ≦ X ≦ 215 ) \displaystyle P(190 \leqq X \leqq 215) P ( 190 ≦ X ≦ 215 ) ≒ p ( 1.0 ) \displaystyle {}\fallingdotseq p(1.0) ≒ p ( 1.0 ) + p ( 1.5 ) \displaystyle {}+ p(1.5) + p ( 1.5 ) = 0.3413 \displaystyle {}= 0.3413 = 0.3413 + 0.4332 \displaystyle {}+ 0.4332 + 0.4332 = 0.7745 ‾ \displaystyle {}= \underline{0.7745} = 0.7745 二項分布のまま計算した正確な値は、(1) が約 0.0255 0.0255 0.0255 、 (2) が約 0.793 0.793 0.793 です。とびとびの値を曲線で置きかえた分の小さなずれがあり、それを小さくする工夫(半整数補正)を厳密定義のページで紹介します。
※ここは発展ページ です。離散型の確率変数の期待値・分散・独立は第6章の厳密定義、二項分布は第7章の厳密定義で扱いました。ここでは、値が切れ目なく変わる確率変数を定積分で定め、正規分布の確率密度関数の全体の面積が 1 1 1 になること(ガウス積分)、期待値が m m m ・ 分散が σ 2 \sigma^2 σ 2 になることを証明します。最後に、二項分布を正規分布で近似する定理と、その近似をよくする工夫を紹介します。積分の計算は微分積分の分野 第8・14章の範囲を使います。
広義積分
正規分布はすべての実数値をとるので、無限に広い範囲での積分が必要です。
定義1:広義積分
関数 f ( x ) f(x) f ( x ) が x x x ≧ a {}\geqq a ≧ a で連続であり、極限 lim t → ∞ ∫ a t f ( x ) d x \displaystyle\lim_{t \to \infty} \int_a^t f(x)\,dx t → ∞ lim ∫ a t f ( x ) d x が有限の値として存在するとき、その値を
∫ a ∞ f ( x ) d x \int_a^{\infty} f(x)\,dx ∫ a ∞ f ( x ) d x と書く。∫ − ∞ a f ( x ) d x \displaystyle\int_{-\infty}^{a} f(x)\,dx ∫ − ∞ a f ( x ) d x も同様に定める。両方が存在するとき ∫ − ∞ ∞ f ( x ) d x \displaystyle \int_{-\infty}^{\infty} f(x)\,dx ∫ − ∞ ∞ f ( x ) d x = ∫ − ∞ a f ( x ) d x \displaystyle {}= \int_{-\infty}^{a} f(x)\,dx = ∫ − ∞ a f ( x ) d x + ∫ a ∞ f ( x ) d x \displaystyle {}+ \int_{a}^{\infty} f(x)\,dx + ∫ a ∞ f ( x ) d x と定める(a a a の選び方によらない)。
実践 問19 の lim t → ∞ ∫ 0 t z φ ( z ) d z \displaystyle\lim_{t \to \infty}\int_0^t z\,\varphi(z)\,dz t → ∞ lim ∫ 0 t z φ ( z ) d z は、この記号で ∫ 0 ∞ z φ ( z ) d z \displaystyle\int_0^{\infty} z\,\varphi(z)\,dz ∫ 0 ∞ z φ ( z ) d z と書けます。
連続型確率変数
定義2:確率密度関数と分布関数
区間 I I I (実数全体でもよい)で連続な関数 f ( x ) f(x) f ( x ) が
f ( x ) \displaystyle f(x) f ( x ) ≧ 0 \displaystyle {}\geqq 0 ≧ 0 ( x ∈ I ) , \displaystyle (x \in I), ( x ∈ I ) , ∫ I f ( x ) d x \displaystyle \int_I f(x)\,dx ∫ I f ( x ) d x = 1 \displaystyle {}= 1 = 1 を満たすとき、f f f を I I I 上の確率密度関数 という。確率変数 X X X が、I I I に含まれるすべての a a a ≦ b {}\leqq b ≦ b について
P ( a ≦ X ≦ b ) \displaystyle P(a \leqq X \leqq b) P ( a ≦ X ≦ b ) = ∫ a b f ( x ) d x \displaystyle {}= \int_a^b f(x)\,dx = ∫ a b f ( x ) d x を満たすとき、X X X を確率密度関数 f f f をもつ連続型確率変数 という。また、F ( x ) F(x) F ( x ) = P ( X ≦ x ) {}= P(X \leqq x) = P ( X ≦ x ) を X X X の分布関数 (累積分布関数)という。
確率変数を「標本空間の上の関数」として厳密に定めるには、無限個の結果を扱う確率論(測度論)が必要になるので、ここでは上の性質をもつものとして話を進めます。
定理1:分布関数の性質
X X X が I I I = [ α , β ] {}= [\alpha,\ \beta] = [ α , β ] 上の確率密度関数 f f f をもつとき、次が成り立つ。
(1) すべての a a a ∈ I {}\in I ∈ I で P ( X = a ) P(X = a) P ( X = a ) = 0 {}= 0 = 0 。 したがって P ( a ≦ X ≦ b ) P(a \leqq X \leqq b) P ( a ≦ X ≦ b ) = P ( a < X < b ) {}= P(a < X < b) = P ( a < X < b )
(2) I I I 上で F ( x ) = ∫ α x f ( t ) d t F(x) = \displaystyle\int_\alpha^x f(t)\,dt F ( x ) = ∫ α x f ( t ) d t であり、F F F は増加関数(広い意味)で、F ′ ( x ) F'(x) F ′ ( x ) = f ( x ) {}= f(x) = f ( x )
証明 (1) 幅が 0 0 0 の区間を考えると
P ( X = a ) \displaystyle P(X = a) P ( X = a ) = P ( a ≦ X ≦ a ) \displaystyle {}= P(a \leqq X \leqq a) = P ( a ≦ X ≦ a ) = ∫ a a f ( x ) d x \displaystyle {}= \int_a^a f(x)\,dx = ∫ a a f ( x ) d x = 0 \displaystyle {}= 0 = 0 である。P ( a ≦ X ≦ b ) P(a \leqq X \leqq b) P ( a ≦ X ≦ b ) = P ( X = a ) {}= P(X = a) = P ( X = a ) + P ( a < X < b ) {}+ P(a < X < b) + P ( a < X < b ) + P ( X = b ) {}+ P(X = b) + P ( X = b ) から後半が従う。
(2) X X X は I I I の外の値をとらないので
F ( x ) \displaystyle F(x) F ( x ) = P ( α ≦ X ≦ x ) \displaystyle {}= P(\alpha \leqq X \leqq x) = P ( α ≦ X ≦ x ) = ∫ α x f ( t ) d t \displaystyle {}= \int_\alpha^x f(t)\,dt = ∫ α x f ( t ) d t である。x x x < y {}< y < y なら F ( y ) − F ( x ) = ∫ x y f ( t ) d t ≧ 0 F(y) - F(x) = \displaystyle\int_x^y f(t)\,dt \geqq 0 F ( y ) − F ( x ) = ∫ x y f ( t ) d t ≧ 0 。F ′ ( x ) F'(x) F ′ ( x ) = f ( x ) {}= f(x) = f ( x ) は、微分積分の分野 第8章の「定積分で表された関数の微分」そのものである。(証明終)
(2) から、確率密度関数は分布関数の導関数、つまり「確率が x x x のまわりでどれだけの速さで増えるか」だと分かります。本文の「密度」という言葉は、この意味です。
期待値と分散
定義3:連続型確率変数の期待値と分散
X X X が I I I 上の確率密度関数 f f f をもつとき
E ( X ) \displaystyle E(X) E ( X ) = ∫ I x f ( x ) d x , \displaystyle {}= \int_I x f(x)\,dx, = ∫ I x f ( x ) d x , V ( X ) \displaystyle V(X) V ( X ) = ∫ I ( x − m ) 2 f ( x ) d x \displaystyle {}= \int_I (x - m)^2 f(x)\,dx = ∫ I ( x − m ) 2 f ( x ) d x ( m = E ( X ) ) \displaystyle (m = E(X)) ( m = E ( X )) と定める(I I I が無限に広いときは、広義積分が存在する場合に限る)。一般に、連続関数 g g g について E ( g ( X ) ) = ∫ I g ( x ) f ( x ) d x E(g(X)) = \displaystyle\int_I g(x) f(x)\,dx E ( g ( X )) = ∫ I g ( x ) f ( x ) d x とする。
定理2:1 次式の変換
X X X が I I I 上の確率密度関数 f f f をもち、a a a ≠ 0 {}\neq 0 = 0 ,b b b を定数、Y Y Y = a X {}= aX = a X + b {}+ b + b とする。
(1) E ( Y ) E(Y) E ( Y ) = a E ( X ) {}= aE(X) = a E ( X ) + b {}+ b + b ,V ( Y ) V(Y) V ( Y ) = a 2 V ( X ) {}= a^2 V(X) = a 2 V ( X ) ,V ( X ) V(X) V ( X ) = E ( X 2 ) {}= E(X^2) = E ( X 2 ) − { E ( X ) } 2 {}- \{E(X)\}^2 − { E ( X ) } 2
(2) Y Y Y は連続型確率変数で、その確率密度関数は g ( y ) g(y) g ( y ) = 1 ∣ a ∣ f ( y − b a ) {}= \dfrac{1}{|a|}\, f\!\left(\dfrac{y - b}{a}\right) = ∣ a ∣ 1 f ( a y − b ) である
証明 (1) 積分の線形性と ∫ I f ( x ) d x \displaystyle \int_I f(x)\,dx ∫ I f ( x ) d x = 1 \displaystyle {}= 1 = 1 から、第6章 定理1・定理2 の証明の ∑ \sum ∑ を ∫ \int ∫ に置きかえて同じ計算で得られる。
(2) a a a > 0 {}> 0 > 0 のとき、c c c ≦ d {}\leqq d ≦ d に対して、y y y = a x {}= ax = a x + b {}+ b + b と置換すると
P ( c ≦ Y ≦ d ) \displaystyle P(c \leqq Y \leqq d) P ( c ≦ Y ≦ d ) = P ( c − b a ≦ X ≦ d − b a ) \displaystyle {}= P\left(\frac{c - b}{a} \leqq X \leqq \frac{d - b}{a}\right) = P ( a c − b ≦ X ≦ a d − b ) = ∫ c − b a d − b a f ( x ) d x \displaystyle {}= \int_{\frac{c-b}{a}}^{\frac{d-b}{a}} f(x)\,dx = ∫ a c − b a d − b f ( x ) d x = ∫ c d f ( y − b a ) 1 a d y \displaystyle {}= \int_c^d f\!\left(\frac{y - b}{a}\right)\frac{1}{a}\,dy = ∫ c d f ( a y − b ) a 1 d y a a a < 0 {}< 0 < 0 のときは不等号の向きが入れかわり、置換で積分の上端と下端が逆になるので、1 a \dfrac{1}{a} a 1 が 1 ∣ a ∣ \dfrac{1}{|a|} ∣ a ∣ 1 に替わる。(証明終)
ガウス積分
正規分布の確率密度関数の全体の面積が 1 1 1 になることは、次の定理にもとづきます。
定理3:ガウス積分
∫ 0 ∞ e − x 2 d x \displaystyle \int_0^{\infty} e^{-x^2}\,dx ∫ 0 ∞ e − x 2 d x = π 2 , \displaystyle {}= \frac{\sqrt{\pi}}{2}, = 2 π , ∫ − ∞ ∞ e − x 2 2 d x \displaystyle \int_{-\infty}^{\infty} e^{-\frac{x^2}{2}}\,dx ∫ − ∞ ∞ e − 2 x 2 d x = 2 π \displaystyle {}= \sqrt{2\pi} = 2 π 証明 e − x 2 e^{-x^2} e − x 2 の不定積分は知っている関数で書けないので、不等式ではさみうちにする。
① 2 つの不等式 すべての実数 t t t で e t e^{t} e t ≧ 1 {}\geqq 1 ≧ 1 + t {}+ t + t である(h ( t ) h(t) h ( t ) = e t {}= e^t = e t − 1 {}- 1 − 1 − t {}- t − t は t t t = 0 {}= 0 = 0 で最小値 0 0 0 )。t t t = − s {}= -s = − s として e − s e^{-s} e − s ≧ 1 {}\geqq 1 ≧ 1 − s {}- s − s 、 また t t t = s {}= s = s として e − s e^{-s} e − s ≦ 1 1 + s {}\leqq \dfrac{1}{1 + s} ≦ 1 + s 1 (s s s > − 1 {}> -1 > − 1 )。 n n n を 2 2 2 以上の整数とし、s s s = x 2 n {}= \dfrac{x^2}{n} = n x 2 を代入して n n n 乗すると
( 1 − x 2 n ) n \displaystyle \left(1 - \frac{x^2}{n}\right)^n ( 1 − n x 2 ) n ≦ e − x 2 \displaystyle {}\leqq e^{-x^2} ≦ e − x 2 ( 0 ≦ x ≦ n ) , \displaystyle (0 \leqq x \leqq \sqrt{n}), ( 0 ≦ x ≦ n ) , e − x 2 \displaystyle e^{-x^2} e − x 2 ≦ ( 1 + x 2 n ) − n \displaystyle {}\leqq \left(1 + \frac{x^2}{n}\right)^{-n} ≦ ( 1 + n x 2 ) − n ( x ≧ 0 ) \displaystyle (x \geqq 0) ( x ≧ 0 ) (左の不等式は両辺が 0 0 0 以上の範囲で n n n 乗している)。e − x 2 e^{-x^2} e − x 2 > 0 {}> 0 > 0 なので
∫ 0 n ( 1 − x 2 n ) n d x \displaystyle \int_0^{\sqrt{n}} \left(1 - \frac{x^2}{n}\right)^n dx ∫ 0 n ( 1 − n x 2 ) n d x ≦ ∫ 0 ∞ e − x 2 d x \displaystyle {}\leqq \int_0^{\infty} e^{-x^2}\,dx ≦ ∫ 0 ∞ e − x 2 d x ≦ ∫ 0 ∞ ( 1 + x 2 n ) − n d x \displaystyle {}\leqq \int_0^{\infty} \left(1 + \frac{x^2}{n}\right)^{-n} dx ≦ ∫ 0 ∞ ( 1 + n x 2 ) − n d x (右の 2 つの広義積分が存在することは、0 0 0 < e − x 2 {}< e^{-x^2} < e − x 2 ≦ ( 1 + x 2 n ) − n {}\leqq \left(1 + \frac{x^2}{n}\right)^{-n} ≦ ( 1 + n x 2 ) − n ≦ 1 1 + x 2 {}\leqq \frac{1}{1 + x^2} ≦ 1 + x 2 1 と、∫ 0 t \int_0^t ∫ 0 t が t t t について増加し上に有界であることから分かる)。
② 置換 J k = ∫ 0 π 2 cos k θ d θ J_k = \displaystyle\int_0^{\frac{\pi}{2}} \cos^k\theta\,d\theta J k = ∫ 0 2 π cos k θ d θ とおく。左の積分で x x x = n sin θ {}= \sqrt{n}\sin\theta = n sin θ とおくと 1 1 1 − x 2 n {}- \frac{x^2}{n} − n x 2 = cos 2 θ {}= \cos^2\theta = cos 2 θ ,d x dx d x = n cos θ d θ {}= \sqrt{n}\cos\theta\,d\theta = n cos θ d θ なので、左辺は n J 2 n + 1 \sqrt{n}\,J_{2n+1} n J 2 n + 1 。 右の積分で x x x = n tan θ {}= \sqrt{n}\tan\theta = n tan θ とおくと 1 1 1 + x 2 n {}+ \frac{x^2}{n} + n x 2 = 1 cos 2 θ {}= \frac{1}{\cos^2\theta} = c o s 2 θ 1 ,d x dx d x = n cos 2 θ d θ {}= \frac{\sqrt{n}}{\cos^2\theta}\,d\theta = c o s 2 θ n d θ なので、右辺は n J 2 n − 2 \sqrt{n}\,J_{2n-2} n J 2 n − 2 (t t t → ∞ {}\to \infty → ∞ は θ \theta θ → π 2 {}\to \frac{\pi}{2} → 2 π に対応)。
③ ワリスの積分 θ \theta θ を π 2 \frac{\pi}{2} 2 π − θ {}- \theta − θ に置換すると、J k J_k J k は微分積分の分野 第14章 実践 問20 の I k = ∫ 0 π 2 sin k x d x I_k = \displaystyle\int_0^{\frac{\pi}{2}} \sin^k x\,dx I k = ∫ 0 2 π sin k x d x と等しく、J k J_k J k = k − 1 k J k − 2 {}= \dfrac{k-1}{k}J_{k-2} = k k − 1 J k − 2 を満たす。すると k J k J k − 1 kJ_kJ_{k-1} k J k J k − 1 = ( k − 1 ) J k − 1 J k − 2 {}= (k-1)J_{k-1}J_{k-2} = ( k − 1 ) J k − 1 J k − 2 なので、k J k J k − 1 kJ_kJ_{k-1} k J k J k − 1 は k k k によらず 1 ⋅ J 1 J 0 1 \cdot J_1 J_0 1 ⋅ J 1 J 0 = π 2 {}= \dfrac{\pi}{2} = 2 π に等しい。0 0 0 ≦ cos θ {}\leqq \cos\theta ≦ cos θ ≦ 1 {}\leqq 1 ≦ 1 より J k J_k J k は k k k について減少するので
π 2 \displaystyle \frac{\pi}{2} 2 π = k J k J k − 1 \displaystyle {}= kJ_kJ_{k-1} = k J k J k − 1 ≧ k J k 2 , \displaystyle {}\geqq kJ_k^2, ≧ k J k 2 , π 2 \displaystyle \frac{\pi}{2} 2 π = ( k + 1 ) J k + 1 J k \displaystyle {}= (k+1)J_{k+1}J_k = ( k + 1 ) J k + 1 J k ≦ ( k + 1 ) J k 2 \displaystyle {}\leqq (k+1)J_k^2 ≦ ( k + 1 ) J k 2 よって k k + 1 ⋅ π 2 \dfrac{k}{k+1}\cdot\dfrac{\pi}{2} k + 1 k ⋅ 2 π ≦ k J k 2 {}\leqq kJ_k^2 ≦ k J k 2 ≦ π 2 {}\leqq \dfrac{\pi}{2} ≦ 2 π となり、はさみうちの原理より k J k 2 kJ_k^2 k J k 2 → π 2 {}\to \dfrac{\pi}{2} → 2 π (k k k → ∞ {}\to \infty → ∞ )。 したがって
( n J 2 n + 1 ) 2 \displaystyle \left(\sqrt{n}\,J_{2n+1}\right)^2 ( n J 2 n + 1 ) 2 = n 2 n + 1 ( 2 n + 1 ) J 2 n + 1 2 \displaystyle {}= \frac{n}{2n+1}(2n+1)J_{2n+1}^2 = 2 n + 1 n ( 2 n + 1 ) J 2 n + 1 2 → 1 2 ⋅ π 2 , \displaystyle {}\to \frac{1}{2}\cdot\frac{\pi}{2}, → 2 1 ⋅ 2 π , ( n J 2 n − 2 ) 2 \displaystyle \left(\sqrt{n}\,J_{2n-2}\right)^2 ( n J 2 n − 2 ) 2 = n 2 n − 2 ( 2 n − 2 ) J 2 n − 2 2 \displaystyle {}= \frac{n}{2n-2}(2n-2)J_{2n-2}^2 = 2 n − 2 n ( 2 n − 2 ) J 2 n − 2 2 → 1 2 ⋅ π 2 \displaystyle {}\to \frac{1}{2}\cdot\frac{\pi}{2} → 2 1 ⋅ 2 π ①の両端はどちらも π 2 \dfrac{\sqrt{\pi}}{2} 2 π に収束するので、はさみうちの原理より前半の式を得る。
後半は、e − x 2 2 e^{-\frac{x^2}{2}} e − 2 x 2 が偶関数であることと、x x x = 2 u {}= \sqrt{2}\,u = 2 u の置換により 2 × 2 × π 2 2 \times \sqrt{2} \times \dfrac{\sqrt{\pi}}{2} 2 × 2 × 2 π = 2 π {}= \sqrt{2\pi} = 2 π 。 (証明終)
大学では、2 重積分と極座標を使って、( ∫ e − x 2 d x ) 2 \left(\int e^{-x^2}dx\right)^2 ( ∫ e − x 2 d x ) 2 を平面全体の積分に直す数行の証明を学びます。正規分布の式に π \pi π が現れるのは、この「円」が隠れているからです。
正規分布の性質
定義4:正規分布
m m m ∈ R {}\in \mathbb{R} ∈ R ,σ \sigma σ > 0 {}> 0 > 0 とする。実数全体の上の確率密度関数
f ( x ) \displaystyle f(x) f ( x ) = 1 2 π σ e − ( x − m ) 2 2 σ 2 \displaystyle {}= \frac{1}{\sqrt{2\pi}\,\sigma}\, e^{-\frac{(x - m)^2}{2\sigma^2}} = 2 π σ 1 e − 2 σ 2 ( x − m ) 2 をもつ連続型確率変数の分布を正規分布 N ( m , σ 2 ) N(m,\ \sigma^2) N ( m , σ 2 ) という。N ( 0 , 1 ) N(0,\ 1) N ( 0 , 1 ) を標準正規分布 、その確率密度関数を φ ( z ) \varphi(z) φ ( z ) = 1 2 π e − z 2 2 {}= \dfrac{1}{\sqrt{2\pi}}e^{-\frac{z^2}{2}} = 2 π 1 e − 2 z 2 と書く。
f ( x ) f(x) f ( x ) > 0 {}> 0 > 0 は明らかで、x x x = m {}= m = m + σ z {}+ \sigma z + σ z と置換すると ∫ − ∞ ∞ f ( x ) d x \displaystyle \int_{-\infty}^{\infty} f(x)\,dx ∫ − ∞ ∞ f ( x ) d x = ∫ − ∞ ∞ φ ( z ) d z \displaystyle {}= \int_{-\infty}^{\infty} \varphi(z)\,dz = ∫ − ∞ ∞ φ ( z ) d z = 1 \displaystyle {}= 1 = 1 となります(定理3)。したがって定義4 の f f f は確かに確率密度関数です。
定理4:正規分布の期待値と分散
X X X が N ( m , σ 2 ) N(m,\ \sigma^2) N ( m , σ 2 ) に従うとき、E ( X ) E(X) E ( X ) = m {}= m = m ,V ( X ) V(X) V ( X ) = σ 2 {}= \sigma^2 = σ 2 である。また、Z Z Z = X − m σ {}= \dfrac{X - m}{\sigma} = σ X − m は N ( 0 , 1 ) N(0,\ 1) N ( 0 , 1 ) に従い、一般に a X aX a X + b {}+ b + b (a a a ≠ 0 {}\neq 0 = 0 ) は N ( a m + b , a 2 σ 2 ) N(am + b,\ a^2\sigma^2) N ( am + b , a 2 σ 2 ) に従う。
証明 まず Z Z Z について。定理2 (2) で a a a = 1 σ {}= \frac{1}{\sigma} = σ 1 ,b b b = − m σ {}= -\frac{m}{\sigma} = − σ m とすると、Z Z Z の確率密度関数は σ f ( m + σ z ) \sigma f(m + \sigma z) σ f ( m + σ z ) = φ ( z ) {}= \varphi(z) = φ ( z ) となる。a X aX a X + b {}+ b + b についても同様に、確率密度関数は
1 ∣ a ∣ ⋅ 1 2 π σ e − ( y − b a − m ) 2 2 σ 2 \displaystyle \frac{1}{|a|}\cdot\frac{1}{\sqrt{2\pi}\,\sigma}\, e^{-\frac{\left(\frac{y-b}{a} - m\right)^2}{2\sigma^2}} ∣ a ∣ 1 ⋅ 2 π σ 1 e − 2 σ 2 ( a y − b − m ) 2 = 1 2 π ∣ a ∣ σ e − ( y − ( a m + b ) ) 2 2 a 2 σ 2 \displaystyle {}= \frac{1}{\sqrt{2\pi}\,|a|\sigma}\, e^{-\frac{(y - (am + b))^2}{2a^2\sigma^2}} = 2 π ∣ a ∣ σ 1 e − 2 a 2 σ 2 ( y − ( am + b ) ) 2 で、これは N ( a m + b , a 2 σ 2 ) N(am + b,\ a^2\sigma^2) N ( am + b , a 2 σ 2 ) の確率密度関数である。
次に E ( Z ) E(Z) E ( Z ) = 0 {}= 0 = 0 ,V ( Z ) V(Z) V ( Z ) = 1 {}= 1 = 1 を示す。実践 問19 と同じ計算で ∫ 0 ∞ z φ ( z ) d z \displaystyle \int_0^{\infty} z\varphi(z)\,dz ∫ 0 ∞ z φ ( z ) d z = 1 2 π \displaystyle {}= \frac{1}{\sqrt{2\pi}} = 2 π 1 が存在し、z φ ( z ) z\varphi(z) z φ ( z ) は奇関数なので ∫ − ∞ 0 z φ ( z ) d z \displaystyle \int_{-\infty}^0 z\varphi(z)\,dz ∫ − ∞ 0 z φ ( z ) d z = − 1 2 π \displaystyle {}= -\frac{1}{\sqrt{2\pi}} = − 2 π 1 。 よって E ( Z ) E(Z) E ( Z ) = 0 {}= 0 = 0 。
E ( Z 2 ) E(Z^2) E ( Z 2 ) は、z 2 φ ( z ) z^2\varphi(z) z 2 φ ( z ) = z ⋅ z φ ( z ) {}= z \cdot z\varphi(z) = z ⋅ z φ ( z ) ,{ − φ ( z ) } ′ \{-\varphi(z)\}' { − φ ( z ) } ′ = z φ ( z ) {}= z\varphi(z) = z φ ( z ) として部分積分すると
∫ 0 t z 2 φ ( z ) d z \displaystyle \int_0^t z^2\varphi(z)\,dz ∫ 0 t z 2 φ ( z ) d z = [ − z φ ( z ) ] 0 t \displaystyle {}= \Big[\,-z\varphi(z)\,\Big]_0^t = [ − z φ ( z ) ] 0 t + ∫ 0 t φ ( z ) d z \displaystyle {}+ \int_0^t \varphi(z)\,dz + ∫ 0 t φ ( z ) d z = − t φ ( t ) \displaystyle {}= -t\varphi(t) = − tφ ( t ) + ∫ 0 t φ ( z ) d z \displaystyle {}+ \int_0^t \varphi(z)\,dz + ∫ 0 t φ ( z ) d z t t t → ∞ {}\to \infty → ∞ で t φ ( t ) t\varphi(t) tφ ( t ) = 1 2 π ⋅ t e t 2 / 2 {}= \dfrac{1}{\sqrt{2\pi}}\cdot\dfrac{t}{e^{t^2/2}} = 2 π 1 ⋅ e t 2 /2 t → 0 {}\to 0 → 0 (微分積分の分野 第12章の大切な極限)、∫ 0 t φ \displaystyle \int_0^t \varphi ∫ 0 t φ → 1 2 \displaystyle {}\to \frac{1}{2} → 2 1 なので、∫ 0 ∞ z 2 φ ( z ) d z \displaystyle \int_0^\infty z^2\varphi(z)\,dz ∫ 0 ∞ z 2 φ ( z ) d z = 1 2 \displaystyle {}= \frac{1}{2} = 2 1 。 偶関数なので E ( Z 2 ) E(Z^2) E ( Z 2 ) = 1 {}= 1 = 1 、V ( Z ) V(Z) V ( Z ) = 1 {}= 1 = 1 − 0 2 {}- 0^2 − 0 2 = 1 {}= 1 = 1 。
最後に X X X = m {}= m = m + σ Z {}+ \sigma Z + σ Z なので、定理2 (1) より E ( X ) E(X) E ( X ) = m {}= m = m ,V ( X ) V(X) V ( X ) = σ 2 {}= \sigma^2 = σ 2 。 (証明終)
定理5:分布曲線の形
N ( m , σ 2 ) N(m,\ \sigma^2) N ( m , σ 2 ) の確率密度関数 f ( x ) f(x) f ( x ) について
(1) f ( m + h ) f(m + h) f ( m + h ) = f ( m − h ) {}= f(m - h) = f ( m − h ) (直線 x x x = m {}= m = m について対称)
(2) f f f は x x x < m {}< m < m で増加、x x x > m {}> m > m で減少し、x x x = m {}= m = m で最大値 1 2 π σ \dfrac{1}{\sqrt{2\pi}\,\sigma} 2 π σ 1 をとる
(3) 変曲点は x x x = m {}= m = m ± σ {}\pm \sigma ± σ の 2 点
(4) lim x → ± ∞ f ( x ) \displaystyle \lim_{x \to \pm\infty} f(x) x → ± ∞ lim f ( x ) = 0 \displaystyle {}= 0 = 0
証明 u u u = x − m σ {}= \dfrac{x - m}{\sigma} = σ x − m ,C C C = 1 2 π σ {}= \dfrac{1}{\sqrt{2\pi}\,\sigma} = 2 π σ 1 とおくと f f f = C e − u 2 / 2 {}= Ce^{-u^2/2} = C e − u 2 /2 。 (1)(4) は明らか。
f ′ ( x ) \displaystyle f'(x) f ′ ( x ) = − C σ u e − u 2 2 , \displaystyle {}= -\frac{C}{\sigma}\,u\,e^{-\frac{u^2}{2}}, = − σ C u e − 2 u 2 , f ′ ′ ( x ) \displaystyle f''(x) f ′′ ( x ) = C σ 2 ( u 2 − 1 ) e − u 2 2 \displaystyle {}= \frac{C}{\sigma^2}\,(u^2 - 1)\,e^{-\frac{u^2}{2}} = σ 2 C ( u 2 − 1 ) e − 2 u 2 f ′ f' f ′ の符号は − u -u − u の符号なので (2) が成り立つ。f ′ ′ f'' f ′′ の符号は u 2 u^2 u 2 − 1 {}- 1 − 1 の符号で、u u u = ± 1 {}= \pm 1 = ± 1 の前後で変わるので (3) が成り立つ。(証明終)
(3) から、標準偏差 σ \sigma σ は「山の中心から、カーブの曲がり方が切りかわる点までの距離」として、グラフの上で見ることができます。本文の fig02 で σ \sigma σ が大きいほど山がなだらかに見えるのはこのためです。
本文の「m m m ± σ {}\pm \sigma ± σ で約 0.683 0.683 0.683 」 などの値は、2 ∫ 0 k φ ( z ) d z 2\displaystyle\int_0^k \varphi(z)\,dz 2 ∫ 0 k φ ( z ) d z (k k k = 1 , {}= 1, = 1 , 2 , \ 2, 2 , 3 \ 3 3 ) を数値計算した 0.6827 0.6827 0.6827 ,0.9545 0.9545 0.9545 ,0.9973 0.9973 0.9973 を丸めたものです。第6章 定理6(チェビシェフの不等式)は、どんな分布でも「2 σ 2\sigma 2 σ 以上離れる確率は 1 4 \frac{1}{4} 4 1 以下」としか言えませんでしたが、正規分布と分かっていれば約 0.046 0.046 0.046 と、ずっと鋭く言えます。
二項分布の近似と半整数補正
定理6:ド・モアブル=ラプラスの定理
0 0 0 < p {}< p < p < 1 {}< 1 < 1 を固定し、X n X_n X n が二項分布 B ( n , p ) B(n,\ p) B ( n , p ) に従うとする。q q q = 1 {}= 1 = 1 − p {}- p − p とすると、すべての実数 a a a < b {}< b < b について
lim n → ∞ P ( a ≦ X n − n p n p q ≦ b ) \displaystyle \lim_{n \to \infty} P\left(a \leqq \frac{X_n - np}{\sqrt{npq}} \leqq b\right) n → ∞ lim P ( a ≦ n pq X n − n p ≦ b ) = ∫ a b φ ( z ) d z \displaystyle {}= \int_a^b \varphi(z)\,dz = ∫ a b φ ( z ) d z が成り立つ。
証明には、n C k {}_n\mathrm{C}_k n C k の大きさを n ! n! n ! の近似式(スターリングの公式。微分積分の分野 第15章の小話)で評価する計算が必要で、この本では省略します。本文の公式6 は、この定理を「n n n が大きいとき」の近似として使ったものです。
X n X_n X n は整数値しかとらないので、P ( X = k ) P(X = k) P ( X = k ) を「k k k − 1 2 {}- \frac{1}{2} − 2 1 から k k k + 1 2 {}+ \frac{1}{2} + 2 1 までの幅 1 1 1 の柱」の面積と見ると(本文 fig05)、曲線で置きかえるときは柱の端までを積分するほうが自然です。a a a ,b b b を整数とすると
P ( a ≦ X ≦ b ) \displaystyle P(a \leqq X \leqq b) P ( a ≦ X ≦ b ) ≒ P ( a − 1 2 − n p n p q ≦ Z \displaystyle {}\fallingdotseq P\biggl(\frac{a - \frac{1}{2} - np}{\sqrt{npq}} \leqq Z ≒ P ( n pq a − 2 1 − n p ≦ Z ≦ b + 1 2 − n p n p q ) \displaystyle {}\qquad\quad \leqq \frac{b + \frac{1}{2} - np}{\sqrt{npq}}\biggr) ≦ n pq b + 2 1 − n p )
これを半整数補正 (連続修正)といいます。本文の例で比べると、次のようになります。
(「補正なし」は端の値 a a a ,b b b をそのまま標準化したもの。値は小数第5位を四捨五入。)
n n n = 20 {}= 20 = 20 では補正なしのずれが目立ちますが、補正を入れると n n n = 20 {}= 20 = 20 でも小数第2位まで合います。n n n が大きくなるほど柱の幅 1 1 1 は n p q \sqrt{npq} n pq にくらべて小さくなるので、補正の有無の差は小さくなります。高校の教科書では補正なしで計算するのがふつうで、この章の問題も補正なしで答えています。
中心極限定理と次の章へ
定理6 は、X n X_n X n が「0 0 0 か 1 1 1 をとる独立な確率変数 n n n 個の和」(第7章 定義1)であることから、次の大きな定理の特別な場合になっています。
定理7:中心極限定理(証明は省略)
X 1 , X_1, X 1 , X 2 , \ X_2, X 2 , … \ \ldots … が互いに独立で、どれも同じ分布に従い、期待値 m m m , 分散 σ 2 \sigma^2 σ 2 (> 0 > 0 > 0 ) をもつとする。S n S_n S n = X 1 {}= X_1 = X 1 + ⋯ {}+ \cdots + ⋯ + X n {}+ X_n + X n とすると、すべての a a a < b {}< b < b について
lim n → ∞ P ( a ≦ S n − n m σ n ≦ b ) \displaystyle \lim_{n \to \infty} P\left(a \leqq \frac{S_n - nm}{\sigma\sqrt{n}} \leqq b\right) n → ∞ lim P ( a ≦ σ n S n − nm ≦ b ) = ∫ a b φ ( z ) d z \displaystyle {}= \int_a^b \varphi(z)\,dz = ∫ a b φ ( z ) d z が成り立つ。
もとの分布がさいころの目でも、第6章のバスの所要時間でも、独立にたくさん足し合わせると、標準化した和は標準正規分布に近づきます。第6章の小話で見た「さいころの個数を増やすと和の分布が釣り鐘形になる」現象の正体で、小話 kb1 のラプラスが示したのもこの定理の原型です。
S n S_n S n を n n n で割った X ‾ \overline{X} X = S n n {}= \dfrac{S_n}{n} = n S n で書き直すと、X ‾ − m σ / n \dfrac{\overline{X} - m}{\sigma/\sqrt{n}} σ / n X − m が近似的に N ( 0 , 1 ) N(0,\ 1) N ( 0 , 1 ) に従う、ということになります(E ( X ‾ ) E(\overline{X}) E ( X ) = m {}= m = m , σ ( X ‾ ) \sigma(\overline{X}) σ ( X ) = σ n {}= \frac{\sigma}{\sqrt{n}} = n σ は第6章 定理5)。また、もとの分布が正規分布なら、n n n によらず X ‾ \overline{X} X はちょうど正規分布 N ( m , σ 2 n ) N\left(m,\ \frac{\sigma^2}{n}\right) N ( m , n σ 2 ) に従うことが知られています(独立な正規分布の和は正規分布になる、という再生性 。証明は大学で学びます)。
第9章では、この X ‾ \overline{X} X を「標本平均」と呼び、実践 問15 の「1.96 1.96 1.96 倍以内に 95 95 95 %」を使って、一部のデータから全体の平均を推定する方法を学びます。