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