JAGS が逸脱度と逸脱情報基準(DIC)を計算する方法を理解しようとしています。
チェーンが 1 つしかない場合、Rのdic.samples関数 inはエラーをスローします。rjags
stop("2 or more parallel chains required")
DIC を計算するための式には、モデル内のパラメーターの期待値での期待逸脱度と逸脱度の計算が含まれます。そのような推定値は単一のチェーンから取得できると思います (複数のチェーンでのみ識別できる収束の問題が時々あるかもしれませんが)。
マニュアルには次のように記載されています。
pD モニターは、複数のチェーン [1] にわたる逸脱度を比較することにより、観測された確率的ノードから有効なパラメーター数 (pD) [3] への寄与を推定します。オプションタイプ(pD)を使用して作成されます。モデルにチェーンが 1 つしかない場合、pD モニターを定義できません。
だから私の質問は：
- JAGS が DIC を計算するために少なくとも 2 つのチェーンを必要とするのはなぜですか?
- より具体的には、なぜ pD モニターは 2 つのチェーンを必要とするのですか?