model {

  mu ~ dnorm( 0.0, 1.0E-6 )
  tau ~ dgamma( 0.001, 0.001 )
  nu ~ dunif( 2.0, 12.0 )
 
  for ( i in 1:n ) {

    y[ i ] ~ dt( mu, tau, nu )
 
  }

  sigma <- 1.0 / sqrt( tau )
 
 y.new ~ dt( mu, tau, nu )
 
}