model
{
   for(i in 1:n) 
   {
      mu[i] <- beta0 + beta1*sage[i] + inprod(u[],Z[i,])
      sstrontium.ratio[i] ~ dnorm(mu[i],tauEps)	
   }	
   for (k in 1:numKnots)
   {
      u[k] ~  dnorm(0,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)
}