model 

{for (i in 1:N)
   {for (j in 1:5)
      {mean[i,j]<-alpha[i]+beta[i]*(t[j]-22)
       y[i,j] ~ dnorm(mean[i,j],tau.y)
       }
    alpha[i] ~ dnorm(mu.alpha,tau.alpha)
    betamean[i]<-mu.beta+gamma*(alpha[i]-mu.alpha)
    beta[i] ~ dnorm(betamean[i],tau.beta)
    }

 mu.alpha ~ dnorm(0,0.0001)
 mu.beta ~ dnorm(0,0.0001)
 gamma ~ dnorm(0,0.00025)
 tau.y ~ dgamma(0.001,0.001)
 tau.alpha ~ dgamma(0.001,0.001)
 tau.beta ~ dgamma(0.001,0.001)
 }