model{
 for(i in 1:n){
   Death[i] ~ dbern(p[i])
   logit(p[i]) <- beta[1] + beta[2]*ISS[i]
                    + beta[3]*RTS[i]
                    + beta[4]*AGE[i]
                    + beta[5]*TI[i]
                    + beta[6]*AGE[i]*TI[i]
 }
 for(i in 1:6){
   ptilde[i] ~ dbeta(a[i],b[i])
   v[i] <- log(ptilde[i]/(1-ptilde[i]))} 
for(i in 1:6){
   beta[i] <- inprod(Xtinv[i,1:6], v[1:6])
 }
}