发帖
雷达卡
各位，大家好！遇到一个棘手而又紧迫的问题。
我使用WinBUGS编写了一个二水平随机效应模型。但是update后，只能迭代几百次。很是郁闷！调整先验分布、初始值，均不能改善。所以，请各位帮助一下。看调整哪里比较合适。怎么才能收敛。因为我不是学数理统计的，我是学医学的。所以，请大家务必说的详细一些！如果留下联系方式，电话或者EMAIL请教的话更好！我的QQ：55996194
看其他文章和帖子都是迭代几万次。我的数据量比较大。一共4212个记录，5个变量。其中包括2个随机变量。
model{
for (i in 1:m)
{
for (t in 1:T)
{
#Binomial likelihood for observed counts
bin[i,t] ~ dbern(prop[i,t])
logit(prop[i,t])<-a+b*year1[t]+c[1]*lst_rank[i,t]+c[2]*m3_3[i,t]
}
a<-a1+a2+c[3]*km_rank
b<-b1+b2
a2 ~ dnorm(0,tau.a2)
b2 ~ dnorm(0,tau.b2)
}
#Other priors
for (i in 1:3){
c~ dnorm(0,1.0E-6)
}
a1 ~ dnorm(0,1.0E-6)
b1 ~ dnorm(0,1.0E-6)
tau.a2 ~ dgamma(0.01,0.01)
tau.b2 ~ dgamma(0.01,0.01)
}
京公网安备 11010802022788号
