model
{
   for(i in 1:nObs) 
   {
      log(mu[i]) <- beta0 + U[idnum[i]] + beta1*progabide[i] 
                  + beta2*sszcntbas[i] + beta3*sage[i]
      szcnt[i] ~ dpois(mu[i])	
   }	
  
   for (iSubj in 1:nSubj)
   {
      U[iSubj] ~ dnorm(0,tauU)
   }

   beta0  ~ dnorm(0,1.0E-8)
   beta1  ~ dnorm(0,1.0E-8)
   beta2  ~ dnorm(0,1.0E-8)
   beta3  ~ dnorm(0,1.0E-8)
   tauU   ~ dgamma(0.01,0.01)
   sigU <- 1/sqrt(tauU)
}