Ex 7.2 IVF: institutional comparison and ranking. Independent and exchangeable analysis done in single run model { for( i in 1 : 52 ) { r.fix[i] <- n[i] * est[i]/100 # effective number of successes r.fix[i] ~ dbin(p.fix[i], n[i]) p.fix[i] ~ dunif(0, 1) ranks.fix[i]<-rank(p.fix[],i) # rank of i'th clinic r.rand[i] <- n[i] * est[i]/100 r.rand[i] ~ dbin(p.rand[i], n[i]) logit(p.rand[i]) <- theta[i] theta[i] ~ dnorm(mu,tau2.inv) ranks.rand[i]<-rank(p.rand[],i) } mu ~dunif(-10,10) tau2.inv <- 1/(tau*tau) tau ~dunif(0,10) }