model{

      ## Likelihood ##
      for( i in 1 : N ) {
        M[i] ~ dnorm(mu[i], prec.obs[i])
        mu[i] ~ dnorm(theta, prec.date[i])
        prec.obs[i] <- 1/(s[i]*s[i])
        
        # Introduction of an individual random effect : use of the Shrinkage distribution
        u[i] ~ dunif(0,1) 
        prec.date[i] <- 1/s02 * u[i]/(1-u[i])     
        sd[i] <- 1/sqrt(prec.date[i])
      }
      
      ## Prior ##
      theta ~ dunif(ta,tb)
      
      }
