3. 統計的推定 保健統計 2010年度 Ⅰ 母集団と標本 a) 標本調査の利点 b) 標本調査における誤差 Ⅱ 確率と確率分布 a) 確率の公理 b) 確率の計算定理 1) 加法定理 2) 条件つき確率と乗法定理 c) 確率分布 1) 2) 3) 4) 確率変数 期待値と分散 2項分布 正規分布 Ⅲ 統計的推定 a) 標本平均の標本分布 b) 点推定 c) 区間推定 1) 母平均の区間推定 ⅰ) 中心極限定理 ⅱ) 母分散が既知の場合の区間推定 ⅲ) 母分散が未知の場合の区間推定 2) 母比率の区間推定 ⅰ) 標本比率の標本分布 ⅱ) 母比率の区間推定 Ⅰ 母集団と標本 母集団(個体数N) 標本(個体数n) × × × × × × × × × × • ある集団についての調査をおこなうとき、調査対象となる集 団(母集団)からその一部を標本として選び、調査する方法 がある。これを標本調査という。 • 標本調査の例として次のようなものが挙げられる。 労働力調査(完全失業率はこの調査の結果求められ る) ⇒ 全国の15歳以上(約1億1千万人)の母集団から、 約10万人を標本として選ぶ調査 内閣支持率調査などの世論調査 ⇒ 全国の有権者(20歳以上の日本国民)(約1億人)の 母集団から、約1000人(新聞社のおこなう内閣支持率 調査の場合) その他、視聴率調査、街頭でのアンケート、製品の品質 管理のための抜き取り調査など、数多くの標本調査がお こなわれている。 a) 標本調査の利点 標本調査をおこなうメリットとして、次のようなことが挙げられる。 • 費用・時間の削減 → 調査票を配布回収する調査では、調査票の印刷費、集計にか かる機械処理費用、人件費などと全部を集計しおえるまでの時 間がだいぶ削減できる。 • 得られる情報の増加、精度の向上 → 調査には調査員が使われることが多いが、ベテランの調査員 は調査の内容をきちんと説明できるので、答えづらい内容を聞い たり、正しい結果を導いたりすることができる。 • 全数調査が不可能な場合にも調査可能 → ガラスの耐久性についての品質管理を調査するなどの場合、 全数調査をおこなうことは不可能である。 b) 標本調査における誤差 標本調査の結果と、真の状態との間にはズレがある。こ のズレのことを誤差というが、標本調査における誤差に は次の2つの種類のものが組み合わさったものである。 1. 非標本誤差 - 調査もれ、無回答、記入ミスなど ⇒ 全数調査でも起こりうる 統計理論によりコントロール不可能 2. 標本誤差 - 標本の偏りによるもの ⇒ 標本調査に固有のもの 統計理論によりコントロール可能 • 標本の偏りによる誤差がどの程度の範囲に収まるかを、 統計理論によって知ることができる。⇒確率の問題 <菅内閣発足直後の支持率の例> 母集団(有権者1億人) × × × × 標本1(朝日1088人) 60% × × × × × × × × 2010年6月10日付の朝刊各 紙に掲載された菅内閣支持 率を見ると、異なった結果に なっている。 同じ対象に同じ調査をおこ なっても、標本によってその結 果が異なる。 これが、標本の偏りである。 × × × × × × 標本2(読売1057人) 64% 標本3(毎日1018人) 66% 標本4(日経888人) 68% 標本5(共同1021人) 61% Ⅱ 確率と確率分布 a) 確率の公理 1. どのような事象Aに対しても、確率の値は常に0と1の間の値 をとる。すなわち、 0 P( A) 1 2. おこりうる事象全体の集合をSとすれば、Sの確率は1である。 P( S ) 1 3. A,B,… が同時に起こらない事象(このとき、A,B,… を排反 事象という)のとき、A,B,… のいずれかが起こる確率はそれ ぞれの事象が起こる確率の和に等しい。すなわち P( A B ) P( A) P( B) b) 確率の計算定理 松中がホームランを 松中がホームランを 打つ(A1) 打たない(A2) ホークスが勝つ(B1) 0.1 0.495 引き分け(B2) 0.01 0.05 ホークスが負ける(B3) 0.04 0.305 計 0.15 0.85 計 0.595 0.06 0.345 1 • 松中がホームランを打ち、ホークスが勝つ確率 → A1とB1が ともに起きる確率である。これをA1とB1の同時確率といい、 P(A1∩B1)とあらわす。(∩は「かつ」(and)を表す記号。capとよぶ。) • 松中がホームランを打つかどうかに関わらず、ホークスが勝 つ確率 → A1が起こるかどうかに関わらず、B1が起きる確率 である。これをB1の周辺確率といい、P(B1)とあらわす。 1) 加法定理 松中がホームランを 松中がホームランを 打つ(A 1) 打たない(A 2) ホークスが勝つ(B 1) 0.1 0.495 引き分け(B 2) 0.01 0.05 ホークスが負ける(B 3) 0.04 0.305 計 0.15 0.85 (例) 松中がホームランを打つか、ホークスが勝つ確率 P( A1 B1 ) P( A1 ) P( B1 ) P( A1 B1 ) 0.15 0.595 0.1 0.645 計 0.595 0.06 0.345 1 加法定理 (∪は「または」(or)を表す記号。cupとよぶ。) <排反事象の場合> (例) ホークスが勝つか、引き分ける確率 P( B1 B2 ) P( B1 ) P( B2 ) 0.595 0.06 0.655 排反事象の場合の 加法定理 2) 条件つき確率と乗法定理 • P(E)>0のとき、事象Eの起こることを条件として、事象Fが起こることを、 (Eを条件とする)Fの条件つき確率といい、P(F|E)であらわす。 (例) 袋の中に、赤球3個、白球2個の計5個の球が入っている。この袋から 球を続けて2個取り出すとき、2個とも赤球となる確率を考えてみよう。 1個目が赤球となる確率は、 P(赤1 ) 3 5 1個目が赤球であったという条件のもとで、 2個目が赤球となる確率は、 2 P (赤 2 | 赤1 ) 4 よって、2個とも赤球となる確率は、 3 2 3 P(赤1 赤2 ) P(赤1 ) P(赤 2 | 赤1 ) 5 4 10 1個目 2個目 乗法定理 (例) 松中がホームランを打ったときに、ホークスが勝つ確率 ⇒ A1を条件とするB1の条件つき確率P(B1|A1)である。 この条件つき確率を用いて、松中がホームランを打ち、ホークスが勝 つ確率を考えると、乗法定理により P( A1 B1 ) P( A1 ) P( B1 | A1 ) となる。よって条件つき確率P(B1|A1)は同時確率を周辺確率で割ること によって求めることができ、 P( B1 | A1 ) となる。 P( A1 B1 ) 0.1 0.67 P( A1 ) 0.15 松中がホームランを 松中がホームランを 打つ(A 1) 打たない(A 2) ホークスが勝つ(B 1) 0.1 0.495 引き分け(B 2) 0.01 0.05 ホークスが負ける(B 3) 0.04 0.305 計 0.15 0.85 計 0.595 0.06 0.345 1 <独立事象の乗法定理> • 事象Eが起こっても起こらなくても事象Fの確率に変化がないとき、すな わちP(F|E) = P(F|Ec) = P(F)のとき、事象Eと事象Fは独立であるという。 ( Ec はEが起こらないという状況をあらわす) 雨が降る(A1) 白鵬が勝つ(B1) 白鵬が負ける(B3) 計 雨が降らない(A2) 0.2 0.6 0.05 0.15 0.25 0.75 計 0.8 0.2 1 この例で雨が降った場合の白鵬が勝つ条件つき確率は P( B1 | A1 ) P( A1 B1 ) 0.2 0.8 P( A1 ) 0.25 雨が降らない場合の白鵬が勝つ条件つき確率は P( B1 | A2 ) P( A2 B1 ) 0.6 0.8 P( A2 ) 0.75 となり、 P(B1|A1) = P(B1|A2) = P(B1)であることから、雨が降るか降らな いかと、白鵬が勝つか負けるかは独立である。 • 事象Eと事象Fが独立である場合、乗法定理は P( E F ) P( E ) P( F ) となる。 c) 確率分布 1) 確率変数 • サイコロを3回振る実験を考える。 • 1の目が出た場合を○、1の目以外が出た場合を×とあらわ すと、起こりうる結果は ○○○, ○○×, ○×○, ×○○, ○××, ×○×, ××○, ××× の8通りである。 • ここで、1の目が何回出たかによって分類するなら • 2回目に振ったサイコロの目は1回目に振ったさいころの目と は独立であるので、独立事象の乗法定理が用いられる。 • 1の目が出た回数を x 回とし、それに対応する確率を P(x) とあらわすと、次のように整理できる。 サイコロを3回振った時の1の目の出る回数 0.8 確率 0.6 0.4 0.2 0 0 1 2 1の目の出る回数 3 • このようにとりうる値†のそれぞれにある確率が対応してい る変数を確率変数といい、その対応関係を確率分布という。 † 連続変数の場合はその値を含む微小な区間を考える。 2) 期待値と分散 等 • 右のような確率で賞金がもらえるくじ 1等 があったとする。 • このくじを1枚購入した時点で、いくら 2等 の賞金がもらえるかはわからない。 3等 • しかし、大体いくらぐらいもらえるか を知りたい。 はずれ • そのとき、 もらえる金額×当たる確率 の総和がもらえると期待できる金額 となる。 1000000 もらえる金額 1000000円 20000円 100円 0円 当たる確率 1 50000 1 1000 1 10 44949 50000 1 1 1 44949 20000 100 0 20 20 10 0 50(円) 50000 1000 10 50000 このくじの期待値は50(円)であるという • このことは、次のように考えることができる。 • 主催者が、全部で5万本のくじを作成したとする。当たる確率 を考えると、このときくじの中に、1等を1本、2等を50本、3等 を5000本入れる必要がある。このくじが、全部で5万本あっ たとすると、下のような度数分布表であらわすことができる。 1等 2等 3等 はずれ 計 xi 1000000 20000 100 0 fi fixi 1 1000000 50 1000000 5000 500000 44949 0 2500000 • ある人がこのくじを5万本全部買い占めたとする。くじの当選 番号が発表された後で当選金の払い戻しを受ける場合、そ の合計金額は確実に2500000(円)であり、1枚あたりの当選 金(すなわち算術平均)を考えると、2500000÷50000=50 (円)であり、期待値に一致する。 期待値=確率変数の算術平均 † このことから、期待値のことを、「平均」「平均値」などと呼ぶこともある。 • サイコロを3回振る実験で1の目が出た回数をxとするなら、x の期待値は 0 125 75 15 1 75 30 3 108 1 1 2 3 0 216 216 216 216 216 216 216 216 2 となり、1の目が出る回数の期待値は0.5回である。 • またサイコロを6回振る実験をおこなうと x P(x) 0 1 2 3 4 5 6 15625 18750 9375 2500 375 30 1 46656 46656 46656 46656 46656 46656 46656 となるので、 1の目が出る回数の期待値は 15625 18750 9375 2500 375 30 1 1 2 3 4 5 6 46656 46656 46656 46656 46656 46656 46656 18750 18750 7500 1500 150 6 46656 0 1 46656 46656 46656 46656 46656 46656 46656 0 となり、6回ふれば1の目が1回ぐらい出るという直感に一致 する。 • 期待値は E(x) x P(x) とあらわすことができる。 • 分散は V( x) ( x E(x))2 P(x) となる。 • 連続型確率変数の場合は E( x) x f( x)dx V( x) ( x E( x)) 2 f( x)dx となる。 • 確率分布は、いくつかの種類に分類することができる。 – 離散型確率分布 2項分布、ポアソン分布、負の2項分布、超幾何分布、・・・ – 連続型確率分布 正規分布、t分布、カイ2乗分布、・・・ 3) 2項分布 [定義] 起こりうる結果がAかBかという2つの結果しか起こらな い試行† をn回繰り返したとき、Aという結果がx回おこったと する。このxの確率分布を2項分布という。 † このような試行をベルヌーイ試行という [分布関数] Aが起こる確率をp、Bが起こる確率をq(=1-p)とす ると、2項分布は p(x)=nCxpxqn-x という式であらわすことができる。この式を2項分布の分布関 数という。 (例) サイコロを3回振る実験では、A(1の目が出る)かB(1の 目が出ない)かという2つの結果しか起こらない試行をn(=3) 回繰り返したとき、A (1の目が出る)という結果がx回おこっ た。このxの確率分布は2項分布(にしたがう)といわれる。 • 1 5 p ,q ,n 3 この例では、 6 であるので、分布関数にあては 6 めると、 p(x) C 1 x 5 3 x となる。 3 x 6 6 • xのとりうる値は0,1,2,3の4つであるので、この分布関数は次 のような関係を表している。 ◎数学補足 nCxについて • nCxはn個の中からx個を選ぶ組み合わせの数であり、次の ように定義される。 n Cx n! x!(n x)! • ここで、!は階乗を表す記号であり、次のようなものである。 n! = n ×(n-1)×・・・×2×1 よって、nCxは次のように計算できる。 n (n 1) (n x 1) (n x) 2 1 n Cx x ( x 1) 2 1 (n x) 2 1 x個 n (n 1) (n x 1) x ( x 1) 2 1 x個 たとえば、5人の班の中から2人の委員を選ぶ組み合わせは 5 C2 5 4 20 10(通り ) 2 1 2 となる。 • サイコロを3回振る実験において、1の目が1回出るパターン は、 ○××, ×○×, ××○ の3通りあるが、これはサイコロを 振る3回のうち、何回目に1の目が出るかを考えたものであり、 3 C1 3 3(通り ) 1 である。 • また、nC0は定義のように計算できないので、 nC0=1と特別 に定義する。 [期待値と分散] 2項分布の期待値(平均)は E(x)=np 分散は V(x)=npq となる。 • 離散型確率変数の期待値は、一般に E(x) x P(x) によっ て求めることができるので、 E(x) 0 0.579 1 0.347 2 0.069 3 0.005 0 0.347 0.138 0.015 0.5 となる。 • 確率変数が2項分布にしたがう場合、期待値は E(x) np と して求めることができる。すなわち、すべてのとりうる値と対 応する確率が得られなくても、期待値が計算できるのである。 • この例の場合 E( x) 3 1 1 0.5 となる。 6 2 1 5 5 • また分散は、 V( x) 3 となる。 6 6 12 2項分布(n=10) 2項分布(n=5) 0.2 0.1 0 0 1 2 3 4 5 0 2項分布(n=100) 1 2 3 4 5 6 7 8 9 10 20 0.3 18 0.4 16 0.35 0.3 0.25 0.2 0.15 0.1 0.05 0 0.5 0.12 0.1 0.08 14 12 10 8 30 27 24 21 33 6 20 18 16 14 12 10 8 6 4 2 0 0 4 0.05 2 0.1 0.16 0.14 0.12 0.1 0.08 0.06 0.04 0.02 0 0 0.15 18 15 0.2 12 9 2項分布(n=50) 6 0 0.25 3 0.06 0.04 0.02 2項分布(n=20) 0 4) 正規分布 • 2項分布において、nを大きくしていくと、左右対称のつりがね 型の正規分布といわれる分布に近づく。 • 2項分布は離散型確率変数の分布であるが、nを無限に大き くしたとき、xのとりうる値は無限に大きくなる。すなわちxは連 続型確率変数として扱われる。 n=500のとき P (x) 0.06 0.05 0.04 0.03 0.02 0.01 144 136 128 120 112 104 96 88 80 72 64 56 48 40 32 24 16 8 0 0 x • 正規分布は数学的に望ましい性質を持った分布 • 身長や知能指数などがこの分布にしたがうといわれている。 • 密度関数 f ( x) 1 e 1 x 2 2 2 e 2.718(自然対数の底) 2 正規分布の平均は、分散は 2 • 正規分布は平均μ、分散σ2の値によって、中心の位置や山 の高さが変わってくる。 <平均の異なる正規分布> σ=1の正規分布 0.5 μ=0 μ=3 μ=-4 0.4 0.3 0.2 0.1 0 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 <分散の異なる正規分布> μ=0の正規分布 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 σ=1 σ=2 σ=1/2 -7 -6 -5 -4 -3 -2 -1 0 1 2 3 4 5 6 7 • これらの正規分布は、中心の位置を移動させたり、目盛りの 幅を変える(横に伸ばしたり、縮めたりする)ことによって、全 て同じ正規分布となる。 <標準化と標準正規分布> • A君は、あるテストで英語が90点、数学が65点であった。 ⇒ 英語の方が数学より成績が良かった?? • 英語の平均点が80点、数学の平均点が50点だった。⇒ 英 語は平均点より10点高い、数学は平均点より15点高い。数 学の方が良い?? • 英語と数学のどちらが成績が良かったのだろうか?⇒ 標準 化の必要性(これを応用したものが偏差値) • 英語が平均80、標準偏差10の正規分布、数学が平均50、 標準偏差20の正規分布にそれぞれしたがうとする。 英語と数学の成績の分布 f(x) 0.05 数学 英語 0.04 0.03 0.02 0.01 0 0 10 20 30 40 50 60 70 80 90 100 点数 • 平均や分散の異なるものを比較するとき、平均や分散をそろ え、その相対的な位置によって比較しようというのが標準化 の考えである。 • 標準化は次のような変換である。 z x • この例で、英語は(90-80)/10=1 数学は(65-50)/20=0.75 となり英語の方が成績が良いことになる。 • 偏差値は、このzを用いて 50+10×z で求められる。この人 の英語の偏差値は60、数学の偏差値は57.5である。 <標準正規分布> • 正規分布にしたがう変数について、このような変換をおこなう と、標準正規分布(平均0、分散1の正規分布)になる。 • 標準正規分布では±1の範囲に68.3%、±2の範囲に95.4%、 ±3の範囲に99.7%が含まれる。 標準正規分布 0.45 0.4 0.35 0.3 0.25 0.2 0.15 0.1 0.05 0 -3.5 -3 -2.5 -2 -1.5 -1 -0.5 0 0.5 1 1.5 2 2.5 3 3.5 Ⅲ 統計的推定 a) 標本平均の標本分布 母集団(大きさ N) 標本(大きさ n) × × × 標本平均 x × × × × × × × × 標本平均 x × × 標本平均 x × × × 母平均 μ • • 標本調査をおこなう場合、通常は1つの標本についての特性値(標本平 均や標本平均など)がわかり、それから母集団の特性値についての推論 をおこなう。母集団全体の情報はわからない。 しかし母集団全体の情報が分かり、とりうるすべての標本について考え ることができたなら、標本の特性値についての分布を考えることができる。 これを標本分布という。 • 500人受講している科目の採点に、25人だけ採点して全体 の平均点を推定しようとするとき、25人の組み合わせ全てか ら標本平均が計算でき、その分布を考えることができる。 • 一般にN個の母集団からn個の標本を選ぶ組み合わせの数 はNCnとあらわすことができる。 N Cn N! n!( N n)! N ( N 1) 1 n (n 1) 1 ( N n) ( N n 1) 1 N ( N 1) ( N n 1) ( N n) ( N n 1) 1 n (n 1) 1 ( N n) ( N n 1) 1 N ( N 1) ( N n 1) 分母も分子も n個ずつ n (n 1) 1 <簡単な例> 中国地方5県の乗用車登録台数(2010年4月末現在、軽自動車 除く)は次のようになっている。(単位: 台) 鳥取 島根 鳥取 島根 197247 216744 232149 184958 岡山 岡山 684788 651448 広島 広島 910219 863016 山口 山口 504297 474582 出典: 中国運輸局『管内保有車両数』 これを10万台単位で四捨五入し、各都道府県の頭文字をア ルファベットで表すと T 2 となる。 母平均、母分散は S 2 O 7 H 9 Y 5 22795 5 5 (2 5) 2 (2 5) 2 (7 5) 2 (9 5) 2 (5 5) 2 2 5 9 9 4 16 0 7.6 5 • この5県を母集団とし、その中から2県を選んで標本とする と、選び方は5C2=10通りとなる。それぞれの標本につい て、標本平均を求め、その分布をあらわすと次のようにな る。 x 2 4.5 5.5 3.5 4.5 5.5 3.5 8 6 7 標本平均の標本分布 2.5 2 度数 パターン T,S 2,2 T,O 2,7 T,H 2,9 T,Y 2,5 S ,O 2,7 S ,H 2,9 S ,Y 2,5 O ,H 7,9 O ,Y 7,5 H ,Y 9,5 1.5 1 0.5 0 0 0.5 1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 6.5 7 7.5 8 標本平均 • 次に標本平均の平均と分散について考えよう。 標本平均の度数分布表から、次のように計算できる。 x fi(度数) 2 3.5 4.5 5.5 6 7 8 計 E( x ) 1 2 2 2 1 1 1 10 f i x i2 f i xi 2 7 9 11 6 7 8 50 4 24.5 40.5 60.5 36 49 64 278.5 f i xi 50 5 f i 10 f i xi2 278.5 2 V (x) ( E ( x ))2 5 2.85 f i 10 ※ 度数分布表からの平均の計算は、(度数×階級値)の総和を度数 の合計で割れば良い なお、この分散の式は計算式であり、次のようにして求 めたものである。 V (x) f i ( xi E ( x ))2 m f i xi2 2 E ( x )f i xi m( E ( x ))2 m f i xi2 2m E( x ) m( E ( x ))2 m f i xi2 ( E ( x ))2 m ※ 分散については、{度数×(階級値-平均)2}の総和を度数の 合計で割ったものとなる • 標本平均の平均、分散と、母平均、母分散の関係として E( x) V (x) 2 が成り立つ。分散に関しては N n 2 V (x ) N 1 n である。この例では、 V (x) 5 2 7.6 3 7.6 2.85 5 1 2 4 2 ※全国規模の統計調査などを考えた場合、母集団の大きさNは非常に 2 大きいので、N n は1に近くなり、V ( x ) とみなせる。 N 1 n 視聴率調査の場合、関東地区1580万世帯から600世帯を選ぶので N n 15800000 600 0.999962≒1 N 1 15800000 1 <補足> 母分散σ2について、ここでは個々の値から平均を引いた ものを2乗して加え、個数で割った。すなわち、 (2 5) 2 (2 5) 2 (7 5) 2 (9 5) 2 (5 5) 2 5 2 とした。 教科書の分散の定義によると、この分母は5-1=4になる はずである。 この教科書の定義は標本不偏分散といわれることもあり、 あとで説明する「母分散がわからない場合の区間推定」 をおこなうときに、計算が簡単になる。 b) 点推定 母集団(個体数 N) 標本(個体数 n) × × × × × × × × × × × × 母平均μ 母分散σ2 母数θ 推論 標本平均x 標本分散s2 標本統計量t 標本から計算された1つ の数値によって、母集団 の数値を推定することを 点推定という。 たとえば、標本平均を母 平均の推定値と考えるこ とや、標本メディアンを母 集団のメディアンの推定 値と考えることである。 ただし、一般に t≠θであ る。 c) 区間推定 • 点推定で母数θをピタリと推定することは難しい。そのため、 標本統計量tの近くの区間を設定し、その区間内に母数θが 含まれることを推定する。これを区間推定という。 ※ たとえば、日本全国全世帯の家計の平均年収を知りたいとき、1万世帯を標 本として調査し、500万円という標本平均を得たとする。この500万円±10万 円という区間をとればよいのか、±30万円という区間を取ればよいのかを考 えていく。(区間が広がれば母平均が含まれる可能性は高くなるが、実用性 に劣る) 1) 母平均の区間推定 ⅰ) 標本平均の標本分布の形状 x の標本分布について、 E (x ) N n 2 V (x ) N 1 n が成り立っていた。母集団の個体数が十分大きいとき、 V (x) が成り立つ 2 n 次に、標本平均 x の分布がどのような形になるのか考えてみよう。 ア) 母集団の分布が正規分布の場合 母集団が平均μ、分散σ2の正規分布にしたがっているとする。 標本平均 x は n x x xn x 1 2 n x i 1 i n であり、正規分布にしたがう変数の和(をnという定数で割ったもの)と なっている。 したがって、正規分布の再生性†より、 x は正規分布にしたがう。 † 確率変数XとYがそれぞれN(μx,σ2x), N(μy,σ2y) にしたがうとき、その1次結 合α X+βY はN(αμx+βμy,α2σ2x+β2σ2y )にしたがう。これを正規分布の再生性と いう。 イ) 母集団の分布が正規分布ではない場合 母集団の分布が正規分布でない場合でも、標本の個体数 n が大きいと き、次のような定理によって標本平均 x の分布は正規分布となる。 <中心極限定理> 算術平均μ, 分散σ2をもつ母集団からとられた大きさ n の標 本の平均 x の分布は、母集団の分布がどのようなもので あっても、 n が大きくなるとき、正規分布 N(μ, )に近づく。 2 n ※ 以上 ア), イ) より、nが大きい時には母集団の分布にかかわらず、標 本平均 x の分布は正規分布となり、標準化された変数 x n の分布は、標準正規分布 N(0, 1) に近づく。 z ⅱ) 母分散が既知の場合の区間推定 標準正規分布にしたがう変数が、-1.96と1.96の間の値をとる確率は 95%である。よって、z x はnが大きいときには、中心極限定理によ n り標準正規分布にしたがうので、 P(1.96 x 1.96) 0.95 n となる。この式のカッコ内を変形すると 1.96 x 1.96 n となり、標本平均 x は 1.96 n n の区間内に95%の確率で含まれる。 x の分布 標準化 z 1.96 また P(1.96 n μ 1.96 zの分布 x n -1.96 0 1.96 n x 1.96) 0.95 のカッコ内は次のようにも変形できる。 n 1.96 x 1.96 1.96 x 1.96 n n n 1.96 x 1.96 x 1.96 n n x 1.96 n n 1.96 n x 1.96 n と x 1.96 n x 1.96 なことを意味している。 1.96 n x 1.96 × μ 1.96 × x n n x 1.96 × n n は次のよう x を中心に、 x 1.96 n という区間を考えると、とりうる標本のうち95% がこの区間内に母平均μを含む。 • このように母数が含まれると考えられる区間を信頼区間とい い、その区間に母数が入ると信頼できる程度を信頼係数と いう。 • この場合、 ( x 1.96 区間である。 n , x 1.96 n ) はμの信頼係数95%の信頼 (例) 20歳男性の身長を調べるために、100人を標本として選 んだところ、標本平均 x =170であった。σ=8であるとき、母 平均μの95%信頼区間を求めよ。 (解) μの95%信頼区間は ( x 1.96 n , x 1.96 n ) 8 8 ,170 1.96 ) 100 100 (170 1.568,170 1.568) (170 1.96 (169.43,171.57) となる。 ⅲ) 母分散が未知の場合の区間推定 母集団(大きさ N) 標本(大きさ n) 信頼区間を求める場合、 z × × × × × × × × × × × × 母平均 μ 母分散 σ2 標本平均 x 標本分散 s2 x が標準正規分布 n にしたがうという性質を用いる。 しかし、母平均の推定をおこ なう場合に、母分散σ2が分 かっているということは、あま り多くない。 (過去の調査に おいて母分散のおおよその 値が分かり、それを用いるな どの例外はあるが) 母分散σ2がわからないとき、代わりに標本分散s2を用いる。 このとき、 t x が自由度n-1のt分布にしたがう。 s/ n 正規分布とt分布 0.45 0.40 0.35 0.30 normal t1 t5 t10 0.25 0.20 0.15 0.10 0.05 0.00 -3.00 -2.00 -1.00 0.00 1.00 2.00 3.00 ※ t分布は標準正規分布を上からつぶしたような、左右対称の形をしている。 自由度が小さいほどつぶれ具合が大きく、自由度が大きいほど標準正規 分布に近くなっている。 <自由度について> 自由度とは、自由に値を取ることのできる個体数のことである。 この場合は、t統計量の自由度は標本分散 s2 の分子に含まれる xi のうち、自由 に値を取ることのできる個数である。 n s2 ( x1 x ) ( x2 x ) ( xn x ) n 1 2 2 とができるが、xn は 2 x i n x (x x) i 1 i n 1 2 なので、x1, …, xn-1 は自由に値をとるこ を満たすように決められ、自由度はn-1となる。 • 母集団の分散が分からないとき、母平均μの95%信頼区 間は、t分布の95%点をt0.95とあらわすと、 ( x t0.95 s s , x t0.95 ) n n となる。 t0.95はt分布表からその値を求める。 x ※ より正確には、母集団の分布が正規分布にしたがうとき、t s/ n が自由度n-1のt分布にしたがう。 しかし、母集団の分布が正規分布にしたがわない場合でも、標本の 大きさがある程度大きければ、 t x は近似的に自由度n-1のt分 s/ n 布にしたがうとみなせる。 x また、nが十分大きい場合、t分布は正規分布に近づくので、 t s/ n が正規分布にしたがうと考えることもある。 x の分布 zの分布 標準化 z 1.96 n μ 1.96 x n -1.96 0 1.96 n tの分布 母分散が分からない場合、 t x s n 変換 が自由度n-1の t分布にしたがう。 t統計量の95%が含まれる区 間の境界値であるt0.95の値を、 t分布表から探し出す。 t (自由度n-1のt分布) x s n -t0.95 t0.95 (例) 20歳女性の身長を調べるために、9人を標本として選ん だところ、標本平均 x =160であった。s=9であるとき、母平 均μの95%信頼区間を求めよ。 (解) 自由度9-1=8のt分布のt0.95=2.306なので、 μの95%信 頼区間は s s , x t0.95 ) n n 9 9 (160 2.306 ,160 2.306 ) 9 9 (160 6.92,160 6.92) ( x t0.95 (153.08,166.92) となる。 2) 母比率の区間推定 ⅰ) 標本比率の標本分布 母集団(大きさ N) × 標本(大きさ n) × × × × × × × × × × × 母比率 p 標本比率 pˆ まず、標本比率 pˆ の標本 分布を考えよう。 内閣支持率を例にとると、 標本比率 pˆ とは、標本n 人のうちのx人が「内閣を 支持する」と答えた割合 であり、 pˆ x である。 n よって pˆ の標本分布を考えるためには、まずxの標本分布を 考えればよい。 • 標本として選ばれた人の答えは、それぞれ「内閣を支持す る」か「内閣を支持しない」かのいずれである。 また選ばれた人が 「内閣を支持する」人である確率は、母 比率pに等しい。 よって、n人の標本を選ぶことは、AかBかという2つの結果し か起こらない試行 をn回繰り返すこととみなすことができ、 「内閣を支持する」人の人数xは2項分布にしたがう。 • 2項分布の期待値は E(x) = np、分散は V(x) = npq である ので、これを用いて、 pˆ の平均、分散を考えてみると、 x E ( x) np E ( pˆ ) E ( ) p n n n x V ( x) npq pq ˆ V ( p) V ( ) 2 2 n n n n となる。 • また、「内閣を支持する」人を1、「内閣を支持しない」人を0と 表すことを考える。n人の標本の中に「内閣を支持する」人は をx人含まれるので、このようにあらわした場合、 pˆ x n は大 きさnの標本の平均とみなすことができ、中心極限定理が適 用できる。 よって、 pˆ の分布は、平均p、分散 pq の正規分布にしたがう。 n 標準化された変数 z pˆ p は標準正規分布にしたがう。 pq n ⅱ) 母比率の区間推定 pˆ p pq が標準正規分布にしたがうことから、母比率pの n pq pq ( pˆ 1.96 , pˆ 1.となる。 96 ) 95%信頼区間は n n z (例) NHK連続テレビ小説「ゲゲゲの女房」最終回(2010.9.25放送)の視 聴率は23.6%であった。この数値は関東地区の約1600万世帯から600 世帯をサンプルとして選んだ結果である。このデータから、関東地区全世 帯の視聴率の95%信頼区間を求めよう。 (解) pの代わりに pˆ を用いてpの95%信頼区間を計算すると ( pˆ 1.96 (0.236 1.96 pq pq , pˆ 1.96 ) n n 0.236 0.764 0.236 0.764 ,0.236 1.96 ) 600 600 (0.236 0.034,0.236 0.034) (0.202,0.270) となる。
© Copyright 2024 ExpyDoc