We use the Markov Chain Monte Carlo (MCMC) method to fit the model to the data, and adopt an adaptive Metropolis-Hastings (M-H) algorithm to carry out the MCMC procedure.