Does a metropolis hastings for the Erlang distribution
mcmc.erlang(
dat,
prior.par1,
prior.par2,
init.pars,
verbose,
burnin,
n.samples,
sds = c(1, 1)
)a matrix of n.samples X 2 parameters, on the estimation scale
the data to fit
mean of priors. A negative binomial (for shape) and a normal for log(scale)
dispersion parameters for priors, dispersion for negative binomial, log scale sd for normal
the starting parameters on the reporting scale
how often to print an update
how many burnin iterations to do
the number of samples to keep and report back
the standard deviations for the proposal distribution