model{
   for(i in 1:n){
      muy[i]<-b0+b1*x1[i]+b2*x2[i]
      y[i]~dt(muy[i],1/sigma^2,tau)
      predy[i]~dt(muy[i],1/sigma^2,tau)
   }
   # Priors
   b0~dnorm(0,0.0001)
   b1~dnorm(0,0.0001)
   b2~dnorm(0,0.0001)
   sigma~dunif(0.001,1000)
   tau<-ntau+1
   ntau~dexp(1/29)
}