model
{
	for(i in 1:nin)
	{
		for(j in 1:n[i])
		{
		
		y[cumn[i]+j] ~ dnorm(mu[cumn[i]+j], prece)

		mu[cumn[i]+j] <- alpha + bhhpers*hhpers[cumn[i]+j] + bhhemp1*hhemp1[cumn[i]+j] +bedterhh*edterhh[cumn[i]+j] +bagehh*agehh[cumn[i]+j] +bsexhh*sexhh[cumn[i]+j] +u[s[i]]+v[s[i]]

		}
		u[s[i]] ~ dnorm(0, precu)
	}


	for(i in 1:N)
	{	
		mug[i]<-alpha+u[i]+v[i]+bhhpers*HHPERS[i]+bhhemp1*HHEMP1[i]+bedterhh*EDTERHH[i]+bagehh*AGEHH[i]+bsexhh*SEXHH[i]
	}

	v[1:N] ~ car.normal(adj[], weights[], num[], precv)

	precu ~ dgamma (a0,b0)
	precv ~ dgamma (a1,b1)
	prece ~ dgamma (a2,b2)

	alpha ~ dflat()
	bhhpers ~ dflat()
	bhhemp1 ~ dflat()
	bedterhh ~ dflat()
	bagehh ~ dflat()
	bsexhh ~ dflat()

	sigmau<-1/precu
	sigmav<-1/precv
	sigmae<-1/prece
}