model
{
for (i in 1:I) {
cases[i] ~ dpois(mu[i]);
log(mu[i]) <- log(pyr[i])+alphac[age[i]]+gammac[period[i]]+betac[cohort[i]];
}
alphamean[1] <- 0.0;
alphaprec[1] <- taua*1.0E-6;
alphamean[2] <- 0.0;
alphaprec[2] <- taua*1.0E-6;}
大家交流一下


雷达卡



京公网安备 11010802022788号







