数理学第九回） 尤度比定とは？.pptVIP

数理統計学第９回 数理統計学(第九回）尤度比検定とは？ 浜田知久馬 胃がんのスクリーニング検査 正常 胃癌 計 検査の 陽性 1497 16 1513 判定 陰性 8483 4 8487 計 9980 20 10000 αエラー 1497/9980=0.15 βエラー 4/20=0.20 問題1 ２つの検査のα,βエラーどちらを選ぶ 検査Ａ 正常 胃癌 計 検査の 陽性 4990 20 5010 判定 陰性 4990 0 4990 計 9980 20 10000 検査Ｂ 正常 胃癌 計 検査の 陽性 50 10 60 判定 陰性 9930 10 9940 計 9980 20 10000 問題2 ２つの妊娠診断薬どちらを選ぶ 検査Ａ 非妊娠 妊娠 計 検査の 陽性 25 250 275 判定 陰性 475 250 725 計 500 500 1000 検査Ｂ 非妊娠 妊娠 計 検査の 陽性 250 475 725 判定 陰性 250 25 275 計 500 500 1000 妊娠診断薬 妊娠診断薬 この検査薬は、妊娠の早期判定の補助として用いるもので、確定診断は、他の所見とともに医師により総合的になされるものです。 妊娠すると妊婦の尿中に分泌されるヒト絨毛性性腺刺激ホルモン(hCG)を鋭敏な抗原?抗体反応で検出する診断薬 尿が希釈され過ぎているか反対に濃縮されている場合 ：妊娠してても陰性 病気の場合やホルモン値に影響を及ぼすお薬を飲んでいる場合：妊娠してなくても陽性 問題 ３　麻原彰晃を無罪と判定した場合についてαエラー，βエラーの観点から論ぜよ ４　美人は必ずしも幸せな人生をおくるとは限らない．この仮説について 　　背理法を用いて証明せよ. 検定の構成法 ?一様最強力検定は存在しないことが多い. ?推定問題の原理的構成法 　最尤法,最小二乗法,モーメント法 　最尤法では確率が最大になるように母数推定 ?ネイマン?ピアソンの基本定理 　確率の比(尤度比)に着目すればよい. 　最尤法に基づいた検定(尤度比検定） 尤度比検定の例 ダーウィンのデータ 仮説：自家受精群と他家受精群で 　　　母平均μが等しいか？ 帰無仮説：Ｈ0：μ1= μ2 対立仮説：Ｈ1：μ1≠μ2 σ2=32 (既知） 仮説の母数空間における表現 尤度比検定統計量 ｆ(Y;θ)＝ｆ(θ; Y)　 　Yをgivenとして確率をθ の関数と考えたとき 尤度(likelihood)とよぶ. maxH1ｆ(θ; Y)　　 maxθ∈ω1ｆ(θ; Y) ───────＝ ─────────＞ｃ maxH0ｆ(θ; Y)　　 maxθ∈ω0ｆ(θ; Y) log(maxH1ｆ(θ; Y))- log(maxH0ｆ(θ; Y))＞logｃ ならば,Ｈ0を棄却 正規分布の確率密度関数 σ２は既知 ｎ個Ｙ１ ，???，Ｙｎ のｎ個のデータの得られる確率ｆ f＝ｆ(y１) ?ｆ(y2) ???ｆ(yｎ) ＝Πｆ(yi) 尤度(Ｈ0) Ｈ0：μ1=μ2=μ 第ｈ群（ｈ＝１：自家受精群，ｈ＝２：他家受精群）の ｉ番目の観測値をｙｈｉ（i＝１,2,???15）で表すことにする． 対数尤度(Ｈ0) 対数尤度(Ｈ1) μ1≠μ2 対数尤度(Ｈ1) μ1≠μ2 対数尤度比 尤度比： maxH1f／maxH0f 対数尤度比：log(maxH1f)- log(maxH0f) 対数尤度比 対数尤度比 Z検定 自由度∞,σ既知の平均値の差の検定 　　　 Ｈ0の下でZは標準正規分布, Z ２ は自由度1のカイ2乗分布にしたがう. ZがZαを越えるときＨ0を棄却 2×対数尤度比がχ2αを越えるときＨ0を棄却 尤度の計算プロ