ロジスティック回帰は、線形結合として説明できます
η=β0+β1X1+...+βkXk
g
g(E(Y))=η
ここで、リンク関数はロジット関数です
E(Y|X,β)=p=logit−1(η)
Y{0,1}η
E(Y)=P(Y=1){0,1}E(Y|X,β)P(Y=1|X,β)P(Y=1|X,β)pY
yi∼Bernoulli(p)
βXηE(Y|X,β)=ηY{0,1}η[0,1]
βiηpY
tβiμiσ2i
model {
# setting up priors
a ~ dnorm(0, .0001)
b ~ dnorm(0, .0001)
for (i in 1:N) {
# passing the linear combination through logit function
logit(p[i]) <- a + b * x[i]
# likelihood function
y[i] ~ dbern(p[i])
}
}
ご覧のように、コードはモデル定義に直接変換されます。ソフトウェアが行うことは、a
およびの正規分布からいくつかの値を描画しb
、それらの値を使用して推定しp
、最後に尤度関数を使用して、これらのパラメーターが与えられたデータの可能性を評価します(これはベイズの定理を使用する場合です。ここを参照してください)詳細な説明)。
βiΣ
⎛⎝⎜⎜⎜⎜β0β1⋮βk⎞⎠⎟⎟⎟⎟∼MVN⎛⎝⎜⎜⎜⎜⎜⎡⎣⎢⎢⎢⎢μ0μ1⋮μk⎤⎦⎥⎥⎥⎥,⎡⎣⎢⎢⎢⎢⎢σ20σ1,0⋮σk,0σ0,1σ21⋮σk,1……⋱…σ0,kσ1,k⋮σ2k⎤⎦⎥⎥⎥⎥⎥⎞⎠⎟⎟⎟⎟⎟
...しかし、これは詳細に入るので、ここで終了しましょう。
ここでの「ベイズ」の部分は、事前確率を選択し、ベイズの定理を使用して、確率論的にモデルを定義しています。ここを参照してください「ベイズモデル」の定義は、ここでいくつかのためのベイズ的アプローチに関する一般的な直感。また、このアプローチを使用すると、モデルの定義が非常に簡単で柔軟性があることにも気づくでしょう。
Kruschke、JK、Aguinis、H.&Joo、H.(2012)。時が来ました:組織科学におけるデータ分析のためのベイズの方法。 組織研究方法、15(4)、722-752。
Gelman、A.、Jakulin、A.、Pittau、GM、およびSu、Y.-S。(2008)。ロジスティックモデルおよびその他の回帰モデルの、情報量の少ないデフォルトの事前分布。 応用統計の年報、2(4)、1360–1383。