model{ for(i in 1:n){ y[i] ~ dnorm(theta[i], tau) theta[i] <- beta[1] - exp(beta[2] + beta[3]*xstd[i]) } beta[1:3] ~ dmnorm(beta0[1:3],prec0[1:3,1:3]) tau ~ dgamma(a,b) for(i in 1:K){ thetastar[i] <- beta[1] - exp(beta[2] + beta[3]*xstar[i]) } }