全部版块 我的主页
论坛 计量经济学与统计论坛 五区 计量经济学与统计软件 winbugs及其他软件专版
1209 0
2018-03-20
求助大神们~!初学者,代码如下,操作到update步骤时,点击update左下方出现sorry something went wrong in procedure Updater.Mode in module UpdaterRejection.想请教下大神们是由于什么原因造成的,应该怎么修改,不吝赐教,灰常感谢~


Model
{
     #Likelihood
     for(i in 1:N)
{
     O[i] ~ dpois(mu[i]);
     log(mu[i])<- log(E[i]) + log.RR[i];
     log.RR[i] <- alpha + beta1*D[i]+ beta2*P[i] + beta3*C[i] + beta4*M[i]+U[i] + S[i];
     U[i] ~ dnorm(0,tau.u);
     RR[i]<-exp(log.RR[i]);
}
    #Priors
    S[1:N]~car.normal(adj[],weights[],num[],tau.s);
   for (k in 1:SumNumNeigh)
  {
    weights[k]<- 1;
    }

   tau.u~dgamma(0.01,0.001);
    tau.s~dgamma(0.5,0.0005);
    alpha~dflat();
    beta1~dnorm(0,0.0000001);
    beta2~dnorm(0,0.0000001);
    beta3~dnorm(0,0.0000001);   
        beta4~dnorm(0,0.0000001);
    #Functions
    sigma.u<-sqrt(1/tau.u);
    sigma.s<-sqrt(1/tau.s);
}

list(N=7,O=c(225,898,1086,1712,933,347,558),E=c(119,1227,682,1184,2797,457,2216),D=c(544,2136,4419,5941,2507,1621,2211),P=c(25496,46344,50133,45816,20031,28865,23530),C=c(12,200
,119,173,105,56,200),M=c(0.7144,5.2924,3.6764,3.3041,8.9979,1.797,7.4244),num=c(0,1,3,3,3,3,1),adj=c(3,2,4,6,3,5,6,4,6,7,3,4,5,5), SumNumNeigh=14)

#INITIAL VALUES
list(tau.u=1,tau.s=1,alpha=0,beta1=0,beta2=0,beta3=0,beta4=0,S=c(0,0,0,0,0,0,0))


二维码

扫码加我 拉你入群

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

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

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

说点什么

分享

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