dnNbinomial
- Negative Binomial Distribution
size : | Natural (pass by const reference) |
Number of failures. | |
p : | Probability (pass by const reference) |
Probability of success. |
p ~ dnBeta(1.0,1.0)
x ~ dnNegativeBinomial(r=10,p)
x.clamp(8)
moves[1] = mvSlide(p, delta=0.1, weight=1.0)
monitors[1] = screenmonitor(printgen=1000, separator = " ", x)
mymodel = model(p)
mymcmc = mcmc(mymodel, monitors, moves)
mymcmc.burnin(generations=20000,tuningInterval=100)
mymcmc.run(generations=200000)