model{
	for (i in 1:N){
		pheno[i]~dbern.aux(q[i])   #latent class model
		q[i]~dnorm(mu[i],1)
		mu[i]<-int+inprod(geno[i,1:102]*b[1:102])
	}

	for(i in 1:102){
		b[i]~dnorm(0, prec[i])
		prec[i]<-1/var[i]
		var[i] <-(1-gamma[i])*tau2 + gamma[i]*tau2*c2
		gamma[i]~dbern(p)
	}

	int~dnorm(0, 0.01) 
	p~dbeta(1,1)
}