model {

  itodist[1,1:4] ~ ddirch(alphaito1[])
  itodist[2,1:4] ~ ddirch(alphaito2[])

  uetzdist[1,1:2] ~ ddirch(alphabin1[])
  uetzdist[2,1:2] ~ ddirch(alphabin2[])

  litdist[1,1:2] ~ ddirch(alphabin1[])
  litdist[2,1:2] ~ ddirch(alphabin2[])

  ppiprior[1:2] ~ ddirch(alphappi[])

  for( i in 1 : M ) {
       ppi[i] ~ dcat(ppiprior[])

       ito[i] ~ dcat(itodist[ppi[i],])
       uetz[i] ~ dcat(uetzdist[ppi[i],])
       lit[i] ~ dcat(litdist[ppi[i],])
  }

}