Dear all,
I have a problem initializing a model in WinBugs version 1.1.1.
This code is no problem for check model, data loaded and compile. However,
when I inits the list, it appears "trap type guard". If I use gen
inits, it appears "trap undefined real result".
Can somebody tell me what I am doing wrong? I look forward to your reply.
Thanks in Advance.
Kuo-Szu Chiang
-----------------------------------------------------------------------------
model
{
for (i in 1:N){
theta[i]~dnorm(mu,tau)
lambda[i]<-n[i]*p[i]
y[i]~dpois(lambda[i])
logit(p[i])<-theta[i]
}
mu~dnorm(0.0,1.0E-6)
sigma<-1/sqrt(tau)
u~dunif(0,1)
tau<-1/(((u/1-u)*s0)*((u/1-u)*s0))
s0<-0.00896
}
list(n=c(1019,1512,1424,54155,447,1907,1755,5756,509,350,473,329,7137,430,304,163,163,159,281,154,889,260,371,232,228,343,454,323,311,784,426,184,181,177,177,291,170,158,274,150,265,257,254,28937,445,447,329,206,313,314,314,202,198,183,292,178,287,282,164
,164,1923,3672,261,581,550,431,399,286,592,246,547,438,202,790,648,354,730,144,1093,384,278,596,1889,22514),
y=c(2,8,8,402,1,12,11,42,2,1,2,1,55,2,1,0,0,0,1,0,6,1,2,1,1,2,3,2,2,6,3,1,1,1,1,2,1,1,2,1,2,2,2,251,4,4,3,2,3,3,3,2,2,2,3,2,3,3,2,2,18,34,3,6,6,5,5,4,7,4,7,6,4,10,9,6,10,4,14,7,6,10,28,334),
N=84)
list(theta=c(0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.
1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1,0.1),tau=1,mu=0)
------------------------------------------------------------------------------
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|