Summarize JAGS Output for Hierarchical MPT Models
Provide clean and readable summary statistics tailored to MPT models based on the JAGS output.
summarizeMPT(mcmc, mptInfo, probs = c(0.025, 0.5, 0.975), summ = NULL)
mcmc |
the actual mcmc.list output of the sampler of a fitted MPT model (accesible via |
mptInfo |
the internally stored information about the fitted MPT model (accesible via |
probs |
quantile probabilities used to compute credibility intervals |
summ |
optional argument for internal use |
The MPT-specific summary is computed directly after fitting a model. However, this function might be used manually after removing MCMC samples (e.g., extending the burnin period).
# Remove additional burnin samples and recompute MPT summary ## Not run: # start later or thin (see ?window) mcmc.subsamp <- window(fittedModel$runjags$mcmc, start=3001, thin=2) new.mpt.summary <- summarizeMPT(mcmc.subsamp, fittedModel$mptInfo) new.mpt.summary ## End(Not run)
Please choose more modern alternatives, such as Google Chrome or Mozilla Firefox.