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] ~ ddirich(prior[1 : 2])
error[3,1 : 3] ~ ddirich(prior[1 : 3])
error[4,1 : 4] ~ ddirich(prior[1 : 4])
p[1 : 4] ~ ddirich(prior[]); # prior for p
}