请教Winbugs做跳跃模型程序初始化问题！

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

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

立即领取

感谢您参与论坛问题回答

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

+2 论坛币

模型如下：
model
{ for( i in 1 : N ) {
za[i] ~ dnorm(muz,tauz)
}
for( i in 1 : N ) {
lama[i] ~ dbeta(2,1)
ja[i] ~ dbern(lama[i])
}
for( i in 1 : N ) {
jumpsize[i]<-za[i]*ja[i]
}
for( i in 1 : N ) {g[i]<-pow(r[i],2 * lata)
tau[i] <- isigmma2 / g[i]
}
for( i in 1 : N ) {
mu[i] <- r[i] + alpha0 + alpha1 * r[i] + alpha2 * pow(r[i],2) + alpha3 / r[i] + za[i] * ja[i]
}
for( i in 1 : N ) {
p[i] ~ dnorm(mu[i],tau[i])
}
lamada ~ dbeta(2,1)
alpha0 ~ dnorm( 0.0,1)
alpha1 ~ dnorm( 0.0,1)
alpha2 ~ dnorm( 0.0,1)
alpha3 ~ dnorm( 0.0,1)
lata ~ dnorm( 0.0,2)
isigmma2 ~ dgamma(4,2)
muz ~ dnorm( -0.01,20)
tauz ~ dgamma(4,2)
sigmma<-1/sqrt(isigmma2)
sigmmaz<-1/sqrt(tauz)
}
初始化如下：
#INITIAL VALUES
list(alpha0=0, alpha1 = 0,alpha2=0,alpha3=0,isigmma2 = 4,lata=0,muz=-0.01,tauz=4,lamada=2,sigmma=NA,
sigmmaz=NA)
然后初始化时，提示还有未初始化的变量，我感觉需要的都初始化了，请各位看看还有那个变量未初始化。

扫码加我 拉你入群

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

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

分享0 收藏0 回帖

关键词：winbugs WINBUG BUGS bug Win 模型 问题 程序