Сообщение об ошибке Trap для экспоненциального распределения для WinBUGS

Я работаю с WinBUGS и пытаюсь найти экспоненциальное распределение набора данных, но продолжаю сталкиваться с:

Ловушка - неопределенный реальный результат

Я пытался изменить мой dunif(), но я продолжаю сталкиваться с той же проблемой. Мои данные раньше работали нормально, но сейчас не работают. Кто-нибудь может мне с этим помочь?? Моя модель выглядит следующим образом:

модель

{

    for (i in 1:N)
    {
        zeros [i] <- 0
        phi [i] <- -log(L[i])
        zeros [i] ~ dpois(phi[i])
        L[i] <- exp(rho1-exp(rho1)*y[i])
    }

    rho1 ~ dunif (-.01, .01)
    lambda <- exp(rho1)

}

Данные

list(y = c(78, 2281, 963, 322, 1349, 777, 7482, 6328, 1728, 4489, 256, 613, 569, 3562, 1902, 2789, 2305, 361, 456, 1199), N=20)

это окно содержит начальные значения

list(rho1=-10, lambda=10)

Любая помощь будет оценена! Я схожу с ума от этого!

0 ответов

Другие вопросы по тегам