Model
{
	theta0   ~ dnorm(0.0,0.1)
	LOWINC   ~ dnorm(-20,0.1)I(-10,0)
	TSAL     ~ dnorm(0.0,0.1)
	TEXP     ~ dnorm(0.0,0.1)
	PGIFT    ~ dnorm(0.0,0.1)
	STEAR    ~ dnorm(0.0,0.1)
	PCTSTATE ~ dnorm(0.0,0.1)
	INSTPUP  ~ dnorm(0.0,0.1)
	LAPASS   ~ dnorm(0.0,0.1) 
	LBURSD   ~ dnorm(0.0,0.1) 
	Y93      ~ dnorm(0.0,0.1)
	Y94      ~ dnorm(0.0,0.1)
	Y95      ~ dnorm(0.0,0.1) 
	Y96      ~ dnorm(0.0,0.1)
	Y97      ~ dnorm(0.0,0.1)
	Y98      ~ dnorm(0.0,0.1) 
	tau      ~ dgamma(16,6)

	for (i in 1 : N) {
                lambda[i] <- theta0 + LAPASS*X9[i] + LBURSD*X10[i] + LOWINC*X2[i] 
			   + TSAL*X3[i] + TEXP*X4[i] + PGIFT*X5[i] + STEAR*X6[i]
			   + PCTSTATE*X7[i] + INSTPUP*X8[i] + Y93*X11[i] + Y94*X12[i] 
			   + Y95*X13[i] + Y96*X14[i] + Y97*X15[i] + Y98*X16[i] + epsilon[i]
		Y[i]       ~ dnorm(lambda[i], tau)
		epsilon[i] ~ dnorm(0.0, tau)
	}
}