JAGS - すべてのパラメータノードに関数を適用する

Question

私はJAGSを初めて使い、R2jagsパッケージ経由でRでモデルを実行しています。 は、モデルコードがKéry&シャウブ2012（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 

JAGSがWinBUGSとは異なる構文で使用されていますか？このコードは、本の例のように3つではなく2つの次元がある点を除いて、本で使用されているものと同じです。

Answer 1

あなたの質問の最後の部分に答えるために、JAGSの構文が違うのでエラーは発生しません（ただし、エラーメッセージはバグで異なって見えるかもしれません）。

は、実際には私はあなたが投稿したコードスニペットに何かを見ることができないし、次の再現性の例では、Yとのevalをデータに与えられたとき、それは、少なくとも動作することを示しています

m <- '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[,]) 

#data# G, T, y, eval 
#monitor# fit 

}#model 
' 

library('runjags') 

G=T <- 10 
y <- matrix(rnorm(100), nrow=G, ncol=T) 
eval <- matrix(rnorm(100), nrow=G, ncol=T) 

results <- run.jags(m) 

てきました140行目が何を参照しているかを検証しましたか？どちらの線140も、あなたが示していないものですか、あるいは次元数の異なるモデルのどこかにfitまたはEのいずれかを指定していますか？

これが当てはまらない場合でもエラーが発生した場合は、質問に問題を示す最小限の再現可能な例を追加してください（既に書いたものの下の--- EDIT ---の下にあることが好ましい）私たちはそれを手助けしようとすることができます。

マット