model
{ #likelihood for(i in 1:N)
{ for(j in 1: M) { x[i,j] ~ dbinom(f[z[i],j], 2)

#prior z[i] ~ dcat(pi) f[z[i],j] ~ dunif(0,1)
}
}
}