Bayesian Analysis for a Fractional Population Growth Model

Box 6

set.seed(123)

fit.movil<-jags(data=movil.data, inits = movil.inits, parameters.to.save=

     movil.params, n.chains =2, n.iter = 500000, n.burnin=250000,

           model.file="movil.model.jags")

  print(fit.movil)