model extend;
{
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)

}

#initial values

list(mu=404.59, tau=0.04, nu=5.0)

#data

list(y= c(375, 392, 393, 397, 398, 398, 399, 399, 399, 399, 399, 399, 399, 400, 400, 400, 400, 401,
          401, 401, 401, 401, 401, 401, 401, 401, 401, 401, 401, 402, 402, 402, 402, 402, 402, 402,
          402, 403, 403, 403, 403, 403, 403, 404, 404, 404, 404, 404, 404, 404, 404, 404, 405, 405,
          405, 405, 405, 406, 406, 406, 406, 406, 406, 406, 406, 406, 406, 406, 406, 407, 407, 407,
          407, 407, 407, 407, 407, 408, 408, 408, 408, 408, 409, 409, 409, 409, 409, 410, 410, 410,
          410, 411, 412, 412, 412, 413, 415, 418, 423, 437), n=100)