model { for (i in 1:n) { t[i] ~ dweib(2, mu[i]) I(t.cens[i], ) mu[i] <- exp(-beta0 - beta1*x[i]) } beta0 ~ dnorm( 7.53, 25) beta1 ~ dnorm(0,.0001) # reference (non-informative) prior } # end of BUGS code