問題タブ [rjags]

Rのベイジアンモデルで使用する次のJAGSモデルがあります。変数「R」の事後分布を推定しようとしています。R を除くすべての変数は、決定論的ノードであると想定されています。各変数 s_A、z_A、z_W、および d はベクトルです。tau_s は data.frame です。したがって、TTD_aquifer と O2s_all は、各 i のベクトルであると予想されます。

これを jags.model() で実行すると、次のエラーが表示されます: RUNTIME ERROR: Invalid vector argument to exp. そのため、R のようにベクトルを exp() に入力できないようです。決定論的な例として、TTD_aquifer と O2s_all の方程式は R で正常に実行されます。指数の問題を回避するには、JAGS で TTD_aquifer の式をどのように記述すればよいですか?

共変量およびyの関数としてモデル化できる従属変数に関するデータがあります。とは「プロット」レベルで観察され、 と は「サイト」レベルで観察されます。プロットはサイト内に階層的にネストされています。これは、共変量データが関連付けられたの 100 個の観測値です。x1x2yx1x2y

基本的にサイトごとに10回のサイトレベルの観測を複製する データフレームを使用して、ギザギザのy関数としてモデル化できますが、問題はx1ありません。x2plot_datax2

しかし、私が本当にやりたいことは、 がプロット レベルの観察を示し、サイトにインデックスを付けるよう に、モデルを階層的に適合させることです。これを行うには、以下の JAGS コードをどのように変更できますか?y[i] ~ x1[i] + x2[j][i][j]

パッケージ 'jagsUI' を使用して R でモデルを実行しています。モデルは少しだけ実行を開始しますが、メッセージが表示されます。

「mat [、「逸脱」] のエラー: 添え字が範囲外さらに: 警告メッセージ: 順序で.params(samples、parameters.to.save、DIC、verbose = verbose): JAGS は逸脱を監視しませんでした。

JAGS の DIC モジュールの問題である可能性があると考える前に、このエラーに遭遇したことは一度もなかったので、モデルを実行するときにモジュールをロードするようにしました。これを引き起こしている可能性のあるものについて何か提案はありますか?

これが私のコードです：