Statistics
| Revision:

root / branches / pthreads-737 / examples / admb-re / logistic / bugs.txt @ 765

History | View | Annotate | Download (361 Bytes)

1
model
2
		{
3
			for( j in 1 : q ) {
4
				u[j] ~ dnorm(0.0,tau)
5
			}
6
			for( j in 1 : p ) {
7
#				b[j] ~ dnorm(0.0,1.0E-6)
8
				b[j] ~ dunif(-10,10)
9
			}			
10
			for( i in 1 : N ) {			 
11
			  logit(P[i]) <- inprod(x[i,],b[])+inprod(z[i,],u[])
12
			  y[i] ~ dbin(P[i],1)
13
			}
14
			tau ~ dgamma(0.001,0.001)I(.003,20000)
15
			sigma <- 1 / sqrt(tau)
16
			logsigma <- log(sigma)
17
		}