Сообщение об ошибке 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)
Любая помощь будет оценена! Я схожу с ума от этого!