ベイズ統計の基礎.PDF

ベイズ統計の基礎 水産資源学におけるベイズ統計の 応用ワークショップ 2007年8月2-3 日，中央水研 遠洋水産研究所 外洋資源部 鯨類管理研究室 岡村 寛 なぜベイズ法を使うのか？ • 水産資源学・生態学では，必ずしも実験デザインを きちんとコントロールした上でデータが得られるわけ ではない • その場合，データとその扱いに様々な問題が起こる データ間の相関，個体差によるばらつき， コントロールされない自然変動， 小さな標本サイズ，多くの欠測値， 直接的なデータはないが間接的な情報ならある， 複数のデータセットの統合，複雑な意思決定 なぜベイズ法を使うのか？ • そのような場合，パラメータの数は多くなり， モデルは複雑な階層構造を持つものになり がち． • そんなときベイズ法は大変ありがたく，強力な 武器となる． • ベイズ法のひとつのネックは計算が大変であ るということだったが，理論と道具の発展によ り，敷居がどんどん低くなってきている． このトークの内容 • ベイズ法の基本的考え方 • 従来の統計学とベイズ統計学 • 水産資源管理に対するベイズ法 の応用例 • 階層ベイズモデル，モデル選択 ベイズ法の基本的考え方 釣り人の漁場選択 漁場A ：0.8 漁場B: 0.2 魚が釣れる確率0.1 魚が釣れる確率0.6 魚が釣れない確率0.9 魚が釣れない確率0.4 確率記号で • 漁場Aを選んだとき釣れる確率（条件付確率） P(釣れる|A) = 0.1 • 漁場Bを選んだとき釣れる確率（条件付確率） P(釣れる|B) = 0.6 • ある人が漁場Aを選んで，魚が釣れる確率 P(A, 釣れる) = P(A) ×P(釣れる|A) = 0.8 ×0.1 = 0.08 P(釣れる|A) = P(A, 釣れる)/P(A) 条件付確率 • Xが起こったという条件のもとで，Yが起 こる確率 P(Y|X) = P(X, Y)/P(X) • Yが起こったという条件のもとで，Xが起 こる確率 P(X|Y) = P(X, Y)/P(Y) • P(X, Y) = P(X|Y)P(Y) = P(Y|X)P(X) 結果から原因を推測する P(X, Y) = P(X|Y)P(Y) = P(Y|X)P(X) X ：原因，Y ：結果とせよ． P(X|Y) = P(原因|結果) = P(Y|X)P(X)/P(Y) = P(結果|原因)P(原因)/P(結果) 結果（データ）から原因（パラメータ，仮説）の推測へ サンプル 自然現象 ベイズの公式 P(X|Y) = P(Y|X)P(X)/P(Y) • P(X|Y) ：データYが与えられた時のパラメータXの確 率分布（事後分布） • P(Y|X) ：パラメータXが与えられた時のデータYが得 られる確率（尤度） • P(X) ：パラメータXのデータを得る前の確率分布（事 前分布） • P(Y) ：可能なすべてのパラメータの値のもとでデー タYが得られる確率分布 P(Y) = ΣX P(Y|X)P(X) ベイズの公式 P(X|Y) = P(X,Y)/Σ P(X,Y) X = P(Y|X)P(X)/ΣX P(Y|X)P(X) = P(Y|X)P(X)