全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件
1828 1
2013-08-16
大家好,新人写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);
}

}

二维码

扫码加我 拉你入群

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

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

全部回复
2016-11-24 17:32:39
我也 遇到了这个问题,请问楼主解决了么?如能告知,不甚感激
二维码

扫码加我 拉你入群

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

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

相关推荐
栏目导航
热门文章
推荐文章

说点什么

分享

扫码加好友,拉您进群
各岗位、行业、专业交流群