密度比)、よくわからんかったので探して読んでみた。密度比を使ってはいけない場面があるよ、特に回帰係数の検定のときには気をつけなさい、というもの。以下BF)は直感的かつ筋の通った指標を提供してくれる。ベイズ・ファクターを求める方法としてよく使われているのがSavage-Dickey密度比(以下SDDR)である。=0\))のもとでの事前密度と事後密度の比になっている、ということである。積分しなくてもいいわけで、これは助かる。問題はこの「ある条件」が成り立たない場合である。実は、検定したいパラメータと局外パラメータが独立でない場合には、この「ある条件」は成り立たない。いう事態は、たとえば、Jeffreys-Zellner-Siow(JZS)事前分布を使った線形重回帰で、いくつかの共変量がネストされたモデルと一般的モデルの両方に含まれているときによく生じる。ベイジアンANOVAも、一般化線形モデルも、潜在過程と行動データを回帰構造でリンクさせる認知過程モデルも、みなこの問題と無縁でいられない。\Psi)\)とする。\theta_0\)について検定したい。BFは、\(H_1\)に対する\(H_0\)の証拠の強さを表す値で、周辺尤度の比として定義される。は\(H_0\)のもとでの局外パラメータの事前分布である。[分母はメモの誤記ではなくて、二重積分のはずなのにインテグラルはひとつだ。数学からっきしだめなのでわかんないんだけど、こういう略記ってありなんすか]これは積分があるので計算が難しい。そこでSDDRの登場である。(包含するほうのモデルのなかにある)検定したいパラメータの、事後の周辺密度と、制約の下での事前の周辺密度の比である。ふつう事前分布はよく知られた分布族に属しているから(t検定とか回帰の場合は単変量コーシー分布を使うことが多い)、分母を求めるのはかんたんである。\theta_0\)の密度を求める。\(H_0\)の下での\(\Psi\)の事前分布が、\(H_1\)の下での\(\Psi\)の条件付き事前分布と同一であるという仮定が必要である。\)のフルモデル\(H_1\)は事前分布\(\pi_0(\Psi)\)の帰無仮説\(H_0\)に縮約されなければならない。いいかえると、局外パラメータは\(H_0\)でも\(H_1\)でも全く同じ役割を果たさなければならない。この仮定は維持されることも多い。ふつう\(\theta\)は効果パラメータで(平均の差とか)、\(\Psi\)は残差分散のような局外パラメータで、後者は\(H_0\)でも\(H_1\)でも同じ役割を持つからだ。\pi_1(\Psi)$$のときには自動的に満たされる。次節では、心理学によくあるシナリオで、かつこの仮定が満たされない場合を示しましょう。\theta_0)\)を通じた期待値である。この修正式は元のBFと違って\(H_0\)にも依存するので、2つのモデルからのサンプリングが必要になり、ちょっと遅くなる。次の重回帰モデルについて考えよう。\(X_i\)は計画行列\(X\)の\(i\)行目。共変量は\(P\)個ですべて中心化済(\(X\)の列は平均0)としよう。\(\epsilon_i\)は\(N(0,\sigma^2)\)にiidに従うとする。傾きパラメータのうち一部だけについて、それがゼロでないかどうかに関心があるとしよう。関心がある\(K\)個のパラメータを\(\beta_K\)、残り\(R\)個のパラメータを\(\beta_R\)と書く。Zベイジアンモデル選択では、すべてのモデルに含まれているパラメータを共通パラメータと呼ぶことが多い。共通パラメータの事前分布はBFにはあまり影響しないことが知られている。\sigma^2\)が共通パラメータである。1/\sigma^2\)が用いられることが多い。ご存じJeffrey事前分布である。0\)について検定したい場合、\(\beta\)の正則事前分布が必要になる。\(H_R\)と\(H_P\)を比べたいだけなら\(\beta_K\)の事前分布さえ正則であればよいが、複数のネストされたモデルを比較したり、\(H_0\)とも比較したい場合には、\(\beta\)のすべてについて正則でなければならない。MVCとは多変量コーシー分布。\(\gamma\)はユーザがアプリオリに決める定数。N\)は中心化ずみ計画行列の共分散行列。上のは\(H_P\)の事前分布だけど、\(H_R\)だったら\(R\)次元のJZS事前分布をつかうわけである。JZS事前分布は中心が0で対称。尺度不変(従属変数や共変量のスケールを変えてもBFが変わらない)。\sqrt{2}/4\)になっている)。他にもいろいろ美点がある。ここで\(B_RP\)を求めるとき、検定したいパラメータ\(\beta_K\)と残りのパラメータ\(\beta_R\)は独立でないから、SDDRは使えない。細かく見てみてみよう。\pi_1(\Psi)$$が満たされていればいいわけだ。\)は多変量コーシー分布にならない。自由度\(1+K\)の多変量t分布になってしまうのだ。\(K\)が大きくなるほど、多変量コーシー分布に比べてゼロ近辺が盛り上がることになる。[面倒なの読み飛ばした]考察。このように、JZS事前分布で一部の傾きパラメータについて検定する際にSDDRは使えない。すべての傾きパラメータが0かどうかを検定する際にはSDDRでよい。g事前分布の場合もSDDRでよい。勉強にはなったけど、正直なところ偏回帰係数の検定にはいま興味がなくて、ほんとはSavage-Dickey密度比の導出について知りたかったのであった。CogPsy)をみるといいらしい。