model;
{
   a ~ dnorm( 0.0,1.0E-6)
   b ~ dnorm( 0.0,1.0E-6)
   for( i in 1 : N ) {
      Y[i] ~ dt(mu[i],tau,2)
   }
   for( i in 1 : N ) {
      mu[i] <- a + b * (X[i])
   }
   tau ~ dgamma(0.01,0.01)
   sd <- 1 / sqrt(tau)
}

data;
Y[]    X[]
3.4 5.5
3.8 5.9
9.1 6.5
2.2 3.3
2.6 3.6
2.9 4.6
2.0 2.9
2.7 3.6
1.9 3.1
3.4 4.9
END

list(N=10)

inits;
list(tau=1,a=0,b=0)