Dear BUGS Friends,
I am new to winbugs and trying to run a case-control study using
bugs. I prepared my input data in SAS format with Y=1,0 disease indicator
and 5 covariates, x1-x5 (1,0). When I run the winbugs, I got lots 'trap
errors'. Would someone please help me to see if I input my data wrong or
specified my model wrong. Your help will be greatly appreciated.
Thanks
Quanhe yang
Centers for disease control and prevention
model
{
for (i in 1 :N) {
Y[i] ~ dbern(p[i])
x1[i]~dbern(q1)
x2[i]~dbern(q2)
x3[i]~dbern(q3)
x4[i]~dbern(q4)
x5[i]~dbern(q5)
b[i] ~ dnorm(mui[i],tau)
mui[i] <-mu+beta0 + beta1 * x1[i] + beta2 *
x2[i] + beta3*x3[i] +beta4*x4[i]+beta5*x5[i]
logit(p[i]) <- b[i]
}
q1~dunif(0.0, 1.0)
q2~dunif(0.0, 1.0)
q3~dunif(0.0, 1.0)
q4~dunif(0.0, 1.0)
q5~dunif(0.0, 1.0)
beta0 ~ dnorm(0.0,1.0E-6)
beta1 ~ dnorm(0.0,1.0E-6)
beta2 ~ dnorm(0.0,1.0E-6)
beta3 ~ dnorm(0.0,1.0E-6)
beta4 ~ dnorm(0.0,1.0E-6)
beta5 ~ dnorm(0.0,1.0E-6)
mu~dnorm(0.0, 1.0E-6)
tau ~ dgamma(0.001,0.1)
sigma <- 1 / sqrt(tau)
}
list( N=614)
Y[ ] x1[ ] x2[ ] x3[ ] x4[ ] x5[ ]
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 1 0 0
1 0 0 1 0 0
1 0 0 1 1 0
1 1 0 1 1 0
1 0 1 0 0 0
1 0 1 0 1 0
1 0 0 0 0 0
1 0 0 0 1 0
1 0 0 1 0 0
1 1 1 0 1 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 0 1 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 0 0 0
1 0 0 0 0 0
1 1 0 0 1 0
1 0 0 1 0 0
1 0 0 1 0 0
1 1 1 0 0 0
1 0 0 1 1 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 1 0 0
1 0 0 0 0 0
1 0 0 0 1 0
1 0 0 1 0 0
1 1 1 0 1 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 1 0 0
1 0 0 1 1 0
1 0 0 0 1 0
1 0 0 1 0 0
1 0 0 0 0 1
1 0 1 0 0 0
1 0 0 0 0 0
1 1 0 0 0 0
1 0 1 0 0 0
1 1 0 0 1 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 0 0 0
1 1 0 0 1 0
1 0 1 0 0 0
1 0 0 1 0 0
1 0 0 1 0 0
1 0 0 0 0 0
1 1 0 1 0 0
1 0 0 0 0 0
1 0 0 1 1 0
1 0 0 1 0 0
1 1 0 0 0 0
1 0 0 1 0 1
1 0 0 1 0 0
1 1 1 1 1 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 1 0 0
1 1 0 1 0 1
1 0 0 1 0 0
1 1 0 1 1 0
1 0 1 0 0 0
1 0 0 0 0 0
1 1 0 1 0 0
1 0 0 1 0 1
1 0 0 0 1 0
1 0 0 0 0 0
1 0 1 1 0 0
1 0 0 1 0 0
1 0 0 1 0 0
1 0 0 1 0 0
1 1 0 0 0 0
1 0 0 1 0 0
1 0 0 1 0 0
1 1 0 0 1 0
1 0 1 1 0 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 1 0 0
1 0 0 1 0 0
1 1 0 1 1 0
1 0 0 1 0 0
1 0 0 1 1 0
1 0 0 0 0 0
1 1 0 0 0 0
1 0 0 1 0 0
1 1 1 0 1 0
1 0 1 1 0 0
1 0 0 1 0 0
1 1 0 0 1 0
1 0 0 0 0 0
1 1 0 0 1 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 1 0 1 0
1 0 0 0 0 0
1 0 0 0 0 0
1 1 0 0 0 0
1 0 0 0 1 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 1 0 0
1 0 0 1 0 0
1 0 0 1 1 0
1 0 0 0 0 0
1 1 1 0 1 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 0 0 0
1 0 0 0 0 0
1 1 1 1 1 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 1 0 0
1 0 0 0 0 0
1 0 0 1 0 0
1 1 0 1 0 0
1 0 1 1 1 1
1 1 0 1 1 0
1 0 0 1 0 0
1 0 0 1 1 0
1 0 0 1 0 0
1 0 0 1 0 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 0 1 0
1 0 0 1 0 0
1 0 0 0 0 0
1 1 1 0 1 0
1 0 0 1 0 0
1 1 0 0 0 0
1 1 0 0 0 0
1 0 0 0 0 0
1 0 1 0 1 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 1 0 0
1 1 0 0 0 0
1 0 0 1 1 0
1 0 0 1 0 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 0 1 0
1 0 0 0 1 0
1 1 0 0 1 0
1 0 0 1 1 0
1 0 0 1 1 0
1 0 0 1 0 0
1 0 0 0 0 0
1 1 0 1 1 0
1 0 0 1 0 0
1 0 0 0 0 0
1 0 0 1 0 0
1 1 0 0 0 0
1 0 0 1 0 0
1 0 0 0 0 0
1 0 0 1 1 0
1 0 0 0 1 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 1 1 0
1 1 0 0 0 0
1 0 0 1 0 0
1 0 0 1 1 0
1 0 0 1 0 0
1 0 0 1 1 0
1 1 0 1 0 0
1 1 1 0 0 0
1 0 0 1 0 0
1 1 0 0 1 0
1 0 0 0 0 0
1 1 0 1 1 0
1 0 0 0 1 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 0 1 0
1 1 1 1 0 0
1 0 0 1 0 0
1 0 0 1 0 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 1 0 0
1 1 0 0 0 0
1 1 0 0 0 0
1 0 0 1 1 0
1 0 0 0 1 0
1 0 1 0 0 0
1 1 0 0 1 0
1 0 0 0 0 0
1 1 1 0 1 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 0 0 0
1 0 0 0 0 0
1 1 0 1 0 0
1 0 0 1 0 0
1 1 0 0 1 0
1 0 0 1 0 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 1 0 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 0 0 0
1 0 0 0 0 0
1 1 1 1 1 0
1 0 0 0 0 0
1 0 0 1 0 1
1 0 0 0 1 0
1 0 0 1 0 0
1 1 0 1 0 0
1 0 0 0 0 0
1 0 0 0 1 0
1 0 0 1 0 0
1 1 0 0 0 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 0 0 0 0
1 1 0 1 1 0
1 0 0 0 0 0
1 1 0 0 0 0
1 0 0 0 0 0
1 0 0 0 0 0
1 0 1 0 0 0
1 0 0 0 0 0
1 0 0 1 1 0
1 0 0 0 0 0
1 0 0 0 1 0
1 0 0 0 0 0
1 0 0 1 0 0
1 0 1 1 0 0
1 0 0 0 0 0
1 1 1 1 1 0
1 1 1 0 0 0
1 0 0 1 0 0
1 0 1 0 0 0
1 0 0 1 0 0
1 0 0 0 0 0
1 1 0 1 0 0
1 0 0 0 0 0
1 0 0 1 1 0
1 0 0 0 0 0
1 1 0 0 1 0
1 0 0 1 0 0
1 0 0 0 1 0
1 0 0 1 1 0
1 0 0 0 0 0
1 1 0 0 1 0
1 0 0 1 1 0
1 1 0 1 0 0
1 1 0 0 0 0
1 0 1 1 1 0
1 0 0 0 1 0
1 0 0 0 0 0
1 0 0 0 0 1
1 0 1 0 0 0
1 0 0 1 0 0
1 0 0 1 0 0
1 0 0 1 0 0
1 0 0 0 0 0
0 0 0 0 0 0
0 1 0 0 1 0
0 0 1 0 0 0
0 1 0 0 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 1 0 1 1 0
0 0 0 1 0 0
0 1 0 1 1 0
0 0 0 0 0 0
0 0 0 1 0 0
0 1 0 1 1 0
0 1 0 1 0 0
0 1 0 1 1 0
0 1 0 0 1 0
0 0 0 1 0 0
0 1 0 1 1 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 1 1 1 0 0
0 0 0 0 0 0
0 1 1 0 1 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 1 0 1 1 0
0 1 0 0 1 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 1 0 0 1 0
0 1 1 1 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 1 0
0 0 0 0 0 0
0 0 0 0 0 0
0 1 0 0 1 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 1 0 1 1 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 1 0
0 1 1 0 1 0
0 0 0 1 0 0
0 0 1 1 0 0
0 0 0 0 1 0
0 0 0 1 0 0
0 0 0 0 0 0
0 1 0 0 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 1 0 1 0 0
0 0 0 1 0 0
0 1 0 1 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 1 0 1 1 0
0 0 0 1 0 0
0 1 0 0 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 1 1 0
0 0 0 0 0 0
0 0 1 1 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 1 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 1 1 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 1 0 0 0 0
0 0 0 0 1 0
0 0 0 0 0 0
0 0 0 0 1 0
0 0 0 1 0 0
0 1 0 1 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 1 1 0
0 0 0 0 1 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 1 0 0 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 0 1 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 1 0 1 0 0
0 1 0 1 0 0
0 0 0 1 1 0
0 1 0 1 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 1 1 0
0 1 0 0 1 0
0 0 0 1 0 1
0 0 0 1 0 0
0 1 0 1 0 0
0 0 0 0 0 0
0 1 0 0 1 0
0 1 1 0 1 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 1 0 0 1 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 1 1 0 1 0
0 0 0 0 1 0
0 1 0 1 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 1 0 0 1 0
0 0 0 0 0 0
0 1 1 0 1 0
0 0 0 1 0 0
0 1 0 0 0 0
0 1 0 1 1 0
0 1 0 1 0 0
0 0 1 0 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 1 0 1 0 0
0 0 0 0 0 0
0 0 0 0 1 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 1 0 1 0 0
0 1 0 1 1 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 1 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 1 0 0 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 0 1 0
0 0 0 1 1 0
0 0 1 1 0 0
0 0 1 0 0 0
0 0 1 0 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 1 1 1 1 0
0 0 0 1 0 0
0 0 0 1 0 0
0 1 0 0 1 0
0 0 0 1 0 0
0 1 0 1 0 0
0 0 0 0 0 0
0 0 1 0 0 0
0 1 1 0 1 0
0 1 0 0 1 0
0 1 0 1 1 0
0 0 0 1 0 0
0 1 0 1 1 0
0 0 0 0 1 0
0 0 0 0 1 0
0 0 0 0 0 0
0 0 0 1 0 0
0 1 0 1 1 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 1 1 0
0 1 0 1 0 0
0 0 1 1 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 0 1 0
0 0 0 0 0 0
0 1 0 1 1 0
0 1 1 0 1 0
0 0 0 0 1 0
0 0 0 0 0 0
0 1 0 0 1 0
0 0 0 1 0 0
0 1 0 0 1 0
0 0 0 1 0 0
0 1 0 1 1 0
0 0 0 0 0 0
0 1 0 0 0 0
0 0 0 1 0 0
0 0 0 0 1 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 0 1 0
0 1 0 0 1 0
0 0 0 1 0 0
0 0 0 1 0 0
0 1 0 1 1 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 0 0 0
0 1 1 0 1 0
0 1 0 1 1 0
0 0 0 0 0 0
0 1 1 1 1 0
0 0 1 0 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 1 0 0 1 0
0 1 1 1 1 0
0 0 0 0 0 0
0 0 0 1 0 0
0 1 0 1 1 0
0 1 0 0 1 0
0 0 0 1 0 0
0 1 1 1 1 0
0 1 0 0 1 0
0 1 0 0 0 0
0 1 0 0 0 0
0 1 0 1 1 0
0 1 0 0 0 0
0 0 0 1 0 0
0 1 0 1 1 0
0 1 1 0 1 0
0 1 0 0 1 0
0 0 0 0 0 0
0 1 0 1 1 0
0 1 0 0 1 0
0 0 0 1 0 0
0 1 0 0 0 0
0 0 1 0 0 0
0 1 1 0 1 0
0 0 0 0 0 0
0 0 0 0 1 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 1 0 0
0 0 0 1 1 0
0 1 0 1 1 0
0 0 0 0 0 0
0 0 0 0 0 0
0 1 0 1 0 0
0 0 0 0 1 0
0 0 0 0 1 0
0 0 0 0 0 0
0 0 0 1 0 0
0 0 0 0 1 0
0 0 0 1 0 0
0 0 0 0 0 0
0 0 0 0 0 0)
list(beta0 = 0, beta1 = 0, beta2 = 0, beta3 = 0, beta4=0, beta5=2, tau = 1)
-------------------------------------------------------------------
This list is for discussion of modelling issues and the BUGS software.
For help with crashes and error messages, first mail [log in to unmask]
To mail the BUGS list, mail to [log in to unmask]
Before mailing, please check the archive at www.jiscmail.ac.uk/lists/bugs.html
Please do not mail attachments to the list.
To leave the BUGS list, send LEAVE BUGS to [log in to unmask]
If this fails, mail [log in to unmask], NOT the whole list
|