model
{
	for(i in 1:N)
	{
		prece[i]<-1/desvar[i]
		Y[i] ~ dnorm(mu[i], prece[i])
  mu[i] <- alpha + bhhpers*HHPERS[i]+bhhemp1*HHEMP1[i]+bedterhh*EDTERHH[i]+bagehh*AGEHH[i] + bsexhh*SEXHH[i]+ u[i]

		u[i] ~ dnorm(0, precu)

		rankmu[i]<-rank(mu[1:N], i)
                rankmu10[i]<- (1-step(rankmu[i]-29))
		rankmu20[i]<- (1-step(rankmu[i]-58))
		rankpovline[i]<-step(mu[i]-693.695)


	}

	precu ~ dgamma (a0,b0)

	alpha ~ dflat()

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

	sigmau<-1/precu
}