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!

Author

David LeBauer, Shawn Serbin

Examples

if (FALSE) { # \dontrun{
summary <- pecan.ma.summary(
 trait.mcmc,
 settings$pfts$pft,
 settings$outdir,
 settings$meta.analysis$threshold)
} # }