openbugs不能generate initial value

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

大家好，新人写ARSV模型，不过遇到问题不能运行，显示 unable to generate initial value for node <theta[1]> of type UpdaterAMblock.UpdaterHetro, 实在是不知道什么原因了，model syntex is correct， load data 也没问题，大家能帮我找找原因么？万分感谢，模型代码如下

model{

for (i in 1:N) {
    Yisigma2[i] <- 1/exp(theta[i]);
    Y[i]~ dnorm(0,Yisigma2[i]);
}

###########################

phistar ~ dbeta(20,1.5);
phi <- 2*phistar - 1;
psistar ~ dbeta(20,1.5);
psi <- 2*psistar - 1;
mu ~ dnorm(0,0.1);
beta <- exp(mu/2);
itau2 ~ dgamma(2.5,0.025);
tau <- sqrt(1/itau2);

###########################

isigma2<-(1-phi*phi-psi*psi)*itau2;
theta0 ~ dnorm(mu,isigma2);
thetamean[1] <- mu + phi*(theta0-mu);
theta[1] ~ dnorm(thetamean[1],itau2);
thetamean[2] <- mu + phi*(theta0-mu) + psi*(theta[1]-mu);
theta[2] ~ dnorm(thetamean[2],itau2);
for (i in 3:N) {
    thetamean[i] <- mu + phi*(theta[i-1]-mu) + psi*(theta[i-2]-mu);
    theta[i] ~ dnorm(thetamean[i],itau2);
}

}

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏0 回帖

关键词：openbugs generate Initial value ERAT generate initial

我也 遇到了这个问题，请问楼主解决了么？如能告知，不甚感激