問題タブ [r2jags]

JAGS でモデルを実行しようとしていますが、次のエラーが発生しました。

私が指定したモデルは

次に、データと初期値を指定します。最後に、モデルを実行します。

手伝って頂けますか？このエラーが発生する理由を知っていますか? 前もって感謝します。

私はJAGSとベイジアン統計に非常に慣れていません.Crawleyの第2版Rブックのベイジアン統計に関する第22章に従おうとしているだけです。単純な線形モデルの本に表示されているコードを正確にコピーします: growth = a + b *tannin で、成長とタンニンという 2 つの連続変数の 9 行があります。データとパッケージは次のとおりです。

ASCII ファイルは次のようになります。

しかし、このコードを使用すると:

次のようなエラー メッセージが表示されます。

私は入門ベイジアン分析コースのプロジェクトに取り組んでおり、R を定期的に使用するのもかなり初めてです。見つけた、またはまとめたデータセットを使用して、階層モデルを構築することになっています。私はデータセットをまとめて、ある国の社会的自由とその富がその幸福にどの程度影響するかという問題を分析しました。データセットは次のようになります。

ステータスはフリーダム ハウスを整数に変換したもので、1自由でない、2部分的に自由、3自由な社会を意味します。私の最初のモデルでは、GDP データを追加したくありません。

クラスで見た例と同じように見えるコードをまとめましたが、問題が発生しました。修正方法がわかりません。

リストの作成からモデリング、事前確率の設定、jagsコマンドの実行まで、すべてを設定する方法は次のとおりです。

次のようなエラーが発生します。

私のデータセットには 460 の観測値があります。私は作成したループを調べていますがi、1 から 460 に移動する必要があると思います。エラーはここにあると思います:mu[i] <- alpha[Country[i]] + beta1[Country[i]]*Year[i]間違った場所を見ているのでしょうか?

誰でもこれをデバッグするのを手伝ってもらえますか? どんな助けでも大歓迎です。

グループ レベルの予測変数を使用して、マルチレベルのロジスティック回帰モデルを適合させる作業を行っています。R 経由で JAGS を使用しています。モデルをパッケージrunjagsと比較すると、異なる動作が得られます。R2Jags

問題を示す再現可能な例を書き込もうとしました。以下では、二項モデルからデータをシミュレートし、データを 8 つのプロットと 2 つのブロックにインデックス付けしてから、マルチレベル ロジスティック回帰を適用して、以下のコードの成功確率 (b1および) を回復します。b2一番下までスクロールして、2 つの近似の概要を確認します。

私の質問は：

1. これら 2 つの適合の事後分布が異なるのはなぜですか? 同じデータ、単一のモデル仕様を使用し、それぞれの前に乱数ジェネレーターを設定しています。事後分布の平均が異なるのはなぜですか? また、Rhat 値がこれほど異なるのはなぜですか?

パッケージのバージョン情報:

rjags フィットからの出力は次のとおりです。

R2jags フィットからの出力は次のとおりです。

2 つの適合からの mu.alpha のトレース プロットを次に示します。まず、rjags フィットから:

次に、R2jags フィットから: