OpenBUGS: инициализация модели

Question

У меня проблема с инициализацией следующей модели в OpenBUGS

model
{
#likelihood
for (t in 1:n) { yisigma2[t] <- 1/exp(theta[t]);
y[t] ~ dnorm(0,yisigma2[t]);
}
#Priors
mu ~ dnorm(0,0.1);
phistar ~ dbeta(20,1.5);
itau2 ~ dgamma(2.5,0.025);
beta <- exp(mu/2);
phi <- 2*phistar-1;
tau <- sqrt(1/itau2);
theta0~dnorm(mu, itau2)
thmean[1] <- mu + phi*(theta0-mu);
theta[1] ~ dnorm(thmean[1],itau2);
for (t in 2:n) { thmean[t] <- mu + phi*(theta[t-1]-mu);
theta[t] ~ dnorm(thmean[t],itau2);
}
}

Это мои данные

list(y=c(-0.0383 , 0.0019 ,......-0.0094),n=945)

И это список моих инициалов

list(phistar= 0.98, mu=0, itau2=50)

Проверка модели, загрузка данных и этапы компиляции в порядке. При загрузке инициалов OpenBUGS говорит, что начальные значения загружены, но цепочка содержит неинициализированные переменные. Затем я попытался инициализировать theta0, но проблема не исчезла. Может ли кто-нибудь помочь мне с этим? Спасибо Халид

Answer 1

Я новичок в OpenBugs, но разве вы не должны указывать распределение для начальных единиц, а не однозначное значение? что-то вроде? inits <- function () {list (alpha = rnorm (1), beta = rnorm (1), sigma = rlnorm (1))}