model {
y <- 10
y ~ dbin(0.5, N)
N ~ dcat(p[])
for (i in 1:100) {p[i] <- 1/100}
}

node mean sd MC error 2.5% median 97.5% start sample
N 21.01 4.702 0.01445 13.0 20.0 32.0 1 100000