私はJAGSが初めてで、R2jags
パッケージを介してRでモデルを実行しています。モデル コードは、Kéry & Schaub 2012 (「Bayesian Population Analysis using WinBUGS」)、pg 399 から取得したコードに基づいています。
カイ二乗不一致測度が計算されます
model {
....
for(g in 1:G) {
for (t in 1:T) {
...
E[g,t] <- pow((y[g,t] - eval[g,t]),2) / eval[g,t]
...
}#t
}#g
fit <- sum(E[,])
}#model
ここでg
、 とt
はサイトと時間のインデックスであり、G
とT
はサイトの数と年数です。
エラーが出ますが
Error in jags.model(model.file, data = data, inits = init.values, n.chains = n.chains, :
RUNTIME ERROR:
Compilation error on line 140.
Cannot evaluate subset expression for fit
WinBUGS と比較して JAGS で使用される構文が異なることが原因ですか? コードは本で使用されているものと同じですが、本の例のように 3 次元ではなく 2 次元を持っている点が異なります。