0

私は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はサイトと時間のインデックスであり、GTはサイトの数と年数です。

エラーが出ますが

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 次元を持っている点が異なります。

4

1 に答える 1