dnBivariatePoisson
- Bivariate Poisson Distribution
theta1 : | RealPos (pass by const reference) |
theta2 : | RealPos (pass by const reference) |
theta0 : | RealPos (pass by const reference) |
th0 ~ dnUniform(0.0, 10.0)
th1 ~ dnUniform(0.0, 10.0)
th2 ~ dnUniform(0.0, 10.0)
x ~ dnBivariatePoisson(th0, th1, th2)
x.clamp([3, 3, 3])
moves[1] = mvSlide(th0, delta=0.01, weight=1.0)
moves[2] = mvSlide(th1, delta=0.01, weight=1.0)
moves[3] = mvSlide(th2, delta=0.01, weight=1.0)
monitors[1] = mnScreen(printgen=20000, th0)
mymodel = model(th1)
mymcmc = mcmc(mymodel, monitors, moves)
mymcmc.burnin(generations=20000, tuningInterval=100)
mymcmc.run(generations=200000)