model

# y= Teneur en proteines
# alpha=paramtres 
# prec=1/variance 

{
     for (i in 1:N) {
          y[i] ~ dnorm(mu[i],prec)
          mu[i]<-alpha
                          }

   ##Prior##
 
          alpha~dnorm(0.0,1.0E-7)
	     var~dunif(0, 10000) 
		prec<-1/var						
}

