model{
 for(i in 1:n){
   Death[i] ~ dbern(p[i])
   logit(p[i]) <- gamma[1] + gamma[2]*(ISS[i]-14.3)/11
                    + gamma[3]*(RTS[i]-7.29)/1.25
                    + gamma[4]*(AGE[i]-31.4)/17
                    + gamma[5]*TI[i]
                    + gamma[6]*(AGE[i]-31.4)*TI[i]/17
 }
 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){   gamma[i] <- inprod(Ztinv[i,1:6], v[1:6])
 }
}