model{ for(i in 1:n){ e[i]~dt(mu,tau,nu) } mu ~dnorm(0, tau2) tau ~dgamma(5,4) tau2 <- 10*tau nu ~dgamma(2, 0.5) }