Tidak, perhitungan Anda tidak benar, karena:
Sebuah) b1 dan b3 mungkin berkorelasi dalam distribusi posterior, dan
b) bahkan jika tidak, itu bukan cara Anda menghitungnya (pikirkan hukum angka besar).
Tapi jangan takut, ada cara yang sangat mudah untuk melakukan ini di WinBUGS. Cukup tentukan variabel baru:
b1b3 <- b1 + b3
dan memonitor nilainya.
EDIT:
Untuk penjelasan yang lebih baik tentang poin pertama saya, misalkan posterior memiliki distribusi normal multivariat bersama (tidak dalam hal ini, tetapi berfungsi sebagai ilustrasi yang bermanfaat). Lalu parameternyabsaya memiliki distribusi N(μsaya,σ2saya), dan interval yang kredibel adalah 95% (μsaya- 1,96σsaya,μsaya+ 1,96σsaya) - perhatikan bahwa ini hanya tergantung pada mean dan varians.
Sekarang b1+b3 akan memiliki distribusi N(μ1+μ3,σ21+ 2ρ13σ1σ3+σ23). Perhatikan bahwa istilah varians (dan karenanya interval kredibel 95%) melibatkan istilah korelasiρ13 yang tidak dapat ditemukan dari interval untuk b1 atau b3.
(Poin saya tentang hukum angka besar adalah bahwa standar deviasi dari jumlah 2 variabel acak independen lebih kecil dari jumlah standar deviasi.)
Adapun cara menerapkannya di WinBUGS, sesuatu seperti ini adalah apa yang ada dalam pikiran saya:
model {
a ~ dXXXX
b1 ~ dXXXX
b2 ~ dXXXX
b3 ~ dXXXX
b1b3 <- b1 + b3
for (i in 1:N) {
logit(p[i]) <- a + b1*x[i] + b2*w[i] + b3*x[i]*w[i]
y[i] ~ dbern(p[i])
}
}
Pada setiap langkah sampler, simpul b1b3
akan diperbarui dari b1
dan b3
. Tidak memerlukan prior karena hanya fungsi deterministik dari dua node lainnya.
Beberapa pemikiran: 1) Saya tidak yakin apakah fakta bahwa ini adalah masalah Bayesian. 2) Saya pikir pendekatan Anda benar 3) Interaksi dalam regresi logistik rumit. Saya menulis tentang ini di sebuah makalah tentang SAS PROC LOGISTIC, tetapi gagasan umum berlaku. Makalah itu ada di blog saya dan tersedia di sini
sumber
Saat ini saya mengalami masalah yang sama. Saya juga percaya bahwa pendekatan untuk menghitung efek total dari w adalah benar. Saya percaya ini dapat diuji melalui
h0: b2 + b3 * rata-rata (x) = 0; ha: b2 + b3 * rata-rata (x)! = 0
Namun, saya menemukan sebuah makalah oleh Ai / Norton, yang mengklaim bahwa "besarnya efek interaksi dalam model nonlinier tidak sama dengan efek marginal dari istilah interaksi, dapat dari tanda yang berlawanan, dan signifikansi statistiknya tidak dihitung oleh perangkat lunak standar. " (2003, hal. 123)
Jadi mungkin Anda harus mencoba menerapkan formula mereka. (Dan jika Anda mengerti bagaimana melakukan itu, tolong beri tahu saya.)
PS. Ini tampaknya menyerupai chow-test untuk regresi logistik. Alfred DeMaris (2004, p. 283) menjelaskan tes untuk ini.
Referensi:
Ai, Chunrong / Norton, Edward (2003): Istilah interaksi dalam model logit dan probit, Economic Letters 80, hlm. 123-129
DeMaris, Alfred (2004): Regresi dengan data sosial: pemodelan variabel respon kontinu dan terbatas. John Wiley & Sons, Inc., Hoboken NJ
sumber