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) 
}