{

  mu ~ dnorm( 0.0, 1.0E-6 )
  sigma ~ dunif( 0.0, 7.0 )
  nu ~ dunif( 2.0, 12.0 )

  for ( i in 1:n ) {

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

  }

  tau <- 1.0 / ( sigma * sigma )

  y.new ~ dt( mu, tau, nu )

}