モデルを構造化し、ベイズモデルのベータをRコードで見つける方法は？

Question

私はJagsモデルを構造化し、ベイズのデータ​​分析のためのベータとプリオリを見つけようとしています。 私のモデルは3つの予測変数x1, x2, x3を持ち、その結果はベルヌーイ分布変数Yです。 事前確率の定義方法P(Y=1|X1), P(Y=1|X2), P(Y=1|X3)結果に影響を及ぼす3つの予測子は？ 私のデータはnXr, n=1920 r=4 columnsX1, X2, X3 and Yの行列に基づいています。

Answer 1

JAGSでGLMを構造化する方法を聞いているようです。 NXRは、データフレームではなく、マトリックスであることが予想される

library('runjags') 
template.jags(Y ~ X1 + X2 + X3, data=nXr, family='binomial', write.data=FALSE) 
## Inspect and edit the JAGSmodel.txt file ## 
results <- run.jags("JAGSmodel.txt", data=nXr) 

注：もしそうなら、最も簡単な方法は、あなたのようなものを使用するため、これを行うためにrunjagsパッケージにtemplate.jags機能を使用することが考えられますあなたが現在持っていると言っているように。これは、関連するパラメータに対して（比較的標準的な）情報提供の最小プライオリティーを使用しますが、必要に応じてプリオーダーを調整するために作成されたモデルファイルを編集することを強くお勧めします。私はあなたが固定効果の切片と係数を意味しない限り、あなたが「ベータを見つける」ことによってあなたが意味するものを知りませんか？

これで問題が解決しない場合は、質問にさらに情報を追加する必要があると思います。たとえば、プレディクタX1-3の詳細、正確に達成しようとしていること、最初の数行あなたのデータのどちらも傷つけません。