Я пытаюсь понять, как JAGS рассчитывает отклонение и информационный критерий отклонения (DIC).
Функция dic.samples
в rjags
в R выдает ошибку, если у вас есть только одна цепочка.
stop("2 or more parallel chains required")
Формула расчета DIC включает расчет ожидаемого отклонения и отклонения при ожидаемых значениях параметров в модели. Я предполагаю, что такие оценки могут быть получены из одной цепочки (хотя иногда могут возникать проблемы сходимости, которые можно идентифицировать только с более чем одной цепочкой).
В руководстве указано, что:
Монитор pD оценивает вклад наблюдаемого стохастического узла в эффективное число параметров (pD) [3] путем сравнения отклонения отклонения по нескольким цепочкам [1]. Он создается с использованием типа опциона (pD). Если модель имеет только одну цепочку, то монитор pD не может быть определен.
Итак, мой вопрос:
- Почему JAGS требует не менее двух цепочек для расчета DIC?
- Или, точнее, почему для монитора pD требуются две цепи?