model { a ~ dnorm(alpha, tau.alpha); b ~ dnorm(beta, tau.beta); g ~ dnorm(gamma, tau.gamma); alpha ~ dnorm(1.0, 1.0E-4); beta ~ dnorm(1.0, 1.0E-4); gamma <- 180 - alpha - beta; }