我的程序如下:
model
{
y2 ~ dbin(p1,x1)
x1 <- n-y1
p1 <- 1-npv
y4 ~ dbin(ppv,x0)
x0 <- y1*p0
ppv <- (se*p)/(se*p+(1-p)*(1-sp))
npv <- (sp*(1-p))/(sp*(1-p)+(1-se)*p)
se ~ dbeta(13.32,6.28)
sp ~ dbeta(21.20,2.06)
y2 <- x1*(1-npv)
y4 <- x0*ppv
p <- (y2+y3+y4)/n
p0 ~ dunif(0,0.2)
p ~ dbeta(1,1)
}
list(n=90000,y1=1061,y3=88)
该程序是关于贝叶斯法估计患病率的,运行到compile时就死了。然后我把程序分成几个小程序运行时就提示multiple definition of node y2等问题。本人初学winbugs,目前做课题急需要,请这里的高手帮忙,谢谢!