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) #non-informative }