model
    {
       for (i in 1 : ns){
          nbiops[i] <- sum(biopsies[i, ]) 
          true[i] ~ dcat(p[])
          biopsies[i, 1 : 4] ~ dmulti(error[true[i], ], nbiops[i])
       }
       error[2,1 : 2] ~ ddirch(prior[1 : 2])
       error[3,1 : 3] ~ ddirch(prior[1 : 3])
       error[4,1 : 4] ~ ddirch(prior[1 : 4])
       p[1 : 4] ~ ddirch(prior[]); # prior for p
    }