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的分布而不是单点值吗？就像是？     inits＆lt; - function（）{list（alpha = rnorm（1），beta = rnorm（1），sigma = rlnorm（1））}