model tonsils

{
 y[1:6]~dmulti(p[],1398)

 for (k in 1:6)
     {p[k]<-phi[k]/sum(phi[])
      phi[k]<-exp(eta[k])
           }

  eta[1]<-  beta[1]-2*beta[2]        -2*beta[4]
  eta[2]<-  beta[1]+  beta[2]-beta[3]+  beta[4]-beta[5]
  eta[3]<-  beta[1]+  beta[2]+beta[3]+  beta[4]+beta[5]
  eta[4]<- -beta[1]-2*beta[2]        +2*beta[4]
  eta[5]<- -beta[1]+  beta[2]-beta[3]-  beta[4]+beta[5]
  eta[6]<- -beta[1]+  beta[2]+beta[3]-  beta[4]-beta[5]
      
  beta[1]~dnorm(-1.5,0.4)
  beta[2]~dnorm(0,0.625)
  beta[3]~dnorm(0,0.625)
  beta[4]~dnorm(0,1.0)
  beta[5]~dnorm(0,1.0)

  rnormal<-p[1]/(p[1]+p[4])
  rlarge<-p[2]/(p[2]+p[5])
  rvlarge<-p[3]/(p[3]+p[6])
 
  lrrlarge<-log(rlarge/rnormal)
  lrrvlarge<-log(rvlarge/rnormal)

 }