model {
n1 ~ dbin(p,N)
m2 ~ dbin(p,n1)
c2 ~ dbin(p,N-n1)
p ~ dunif(0,1)
for (i in 1:90)
{
pi[i] <- 0
}
for (i in 91:400)
{
pi[i] <- 1
}
N ~ dcat(pi)
}
