model
{
   for(i in 1:nObs) 
   {
      mu[i] <- gamma[idnum[i]] + beta1*sweeks[i]
      weight[i] ~ dnorm(mu[i],tauEps)	
   }	

   for (iSubj in 1:nSubj)
   {
       gamma[iSubj] ~ dnorm(beta0,tauU)
   }  
   beta0  ~ dnorm(0,1.0E-8)
   beta1  ~ dnorm(0,1.0E-8)
   tauU   ~ dgamma(0.01,0.01)
   tauEps ~ dgamma(0.01,0.01)
   sigU <- 1/sqrt(tauU)
   sigEps <- 1/sqrt(tauEps)
}