各位大侠帮帮忙,我是openbugs初学者,现有程序如下:
model
{
   alpha ~ dnorm(5,50)
   blpha ~ dnorm(5,50)
   for( i in 1 : N ) {
      p[i] ~ dbeta(alpha,blpha)
   }
   for( i in 1 : N ) {
      r[i] ~ dbin(p[i],n[i])
   }
   feta <- 1 / (1 + alpha + blpha)
}
list(N=47,n=c(1,9,10,6,17,12,9,5,5,2,1,10,1,15,9,12,4,7,8,8,9,2,12,8,4,5,8,9,13,3,2,9,5,3,6,2,3,4,7,5,6,8,6,5,4,4,1),r=c(0,0,0,1,2,1,1,0,1,1,1,2,1,2,1,2,1,2,5,1,4,1,0,1,2,4,2,4,6,1,2,1,1,3,4,2,1,3,3,2,3,5,3,3,3,3,1),p=c(0,0,0,0.167,0.118,0.083,0.111,0,0.2,0.5,1,0.2,1,0.133,0.111,0.167,0.25,0.286,0.625,0.125,0.444,0.5,0,0.125,0.5,0.8,0.25,0.444,0.462,0.333,1,0.111,0.2,1,0.667,1,0.333,0.75,0.429,0.4,0.5,0.625,0.5,0.6,0.75,0.75,1))
list(alpha=1,blpha=1)
前面所有的步骤都OK 但是最后updating时出现了错误:“sorry something went wrong in procedure UpdaterStd.Sample”,也就是通常说的最后无法updating。求哪位大侠帮帮小女子,怎么可以修改正确,小女子不胜感激