Generate summary statistics and diagnostics for PEcAn meta.analysis

Generate summary statistics and diagnostics for PEcAn meta.analysis

pecan.ma.summary(mcmc.object, pft, outdir, threshold = 1.2, gg = FALSE)

Arguments

mcmc.object

JAGS mcmc output from pecan.ma
pft

plant functional type
outdir

output directory
threshold

Gelman-Rubin convergence diagnostic (MGPRF) default = 1.2 following Bolker 2008 Ecological Models and Data in R
gg

produce extra diagnostic plots using the "ggmcmc" package? Caution: very slow!

Examples

# NOT RUN {
summary <- pecan.ma.summary(
 trait.mcmc,
 settings$pfts$pft,
 settings$outdir,
 settings$meta.analysis$threshold)
# }