Model {
		for (i in 1:N) {
			MMSE[i] ~ dbeta(alpha[i],beta[i])
			alpha[i] <- mu[i]*phi
			beta[i] <- phi*(1-mu[i])
			mu[i] <- ilogit(b0+b1*Age[i]+b2*MdV[i]+b3*AD[i]+b4*Anx[i])
		}
		
		# priors
		b0 ~ dnorm(0,0.001)
		b1 ~ dnorm(0,0.001)
		b2 ~ dnorm(0,0.001)
		b3 ~ dnorm(0,0.001)
		b4 ~ dnorm(0,0.001)
		phi ~ dunif(0,10)
		#phi ~ dgamma(0.001,0.0001)
}

# initials
list(b0=1,b1=-0.1,b2=0.5,b3=0.5,b4=0.5,phi=3)
list(b0=2,b1=-0.5,b2=-0.5,b3=-0.5,b4=-0.5,phi=5)