model {
yT ~ dbin(pT, nT)
pT ~ dunif(0, 1)
for (i in 1:8) {
sP[i] ~ dexp(theta)
}
theta ~ dgamma(0.001, 0.001)
surv.t <- pT/theta # expected survival with transplant
Is <- surv.t - 2
}


Data:
list(yT=8, nT=10, sP=c(2,3,4,4,6,7,10,12))

Inits:
list(theta=1)

node mean sd MC error 2.5% median 97.5% start sample
Is 3.136 2.271 0.003046 0.2014 2.669 8.836 1001 500000
pT 0.7495 0.1202 1.702E-4 0.482 0.7639 0.9396 1001 500000
surv.t 5.136 2.271 0.003046 2.201 4.669 10.84 1001 500000