Output diagnostics

Source: R/output_diagnostics.R

output_diagnostics returns diagnostics for a fit MixSIAR model

output_diagnostics(
  jags.1,
  mix,
  source,
  output_options = list(summary_save = TRUE, summary_name = "summary_statistics",
    sup_post = FALSE, plot_post_save_pdf = TRUE, plot_post_name = "posterior_density",
    sup_pairs = FALSE, plot_pairs_save_pdf = TRUE, plot_pairs_name = "pairs_plot", sup_xy
    = TRUE, plot_xy_save_pdf = FALSE, plot_xy_name = "xy_plot", gelman = TRUE, heidel =
    FALSE, geweke = TRUE, diag_save = TRUE, diag_name = "diagnostics", indiv_effect =
    FALSE, plot_post_save_png = FALSE, plot_pairs_save_png = FALSE, plot_xy_save_png =
    FALSE, diag_save_ggmcmc = TRUE,      return_obj = FALSE)
)

Arguments

jags.1

rjags model object, output from run_model function
mix

output from load_mix_data
source

output from load_source_data
output_options

list containing options for plots and saving:

  • summary_save: Save the summary statistics as a txt file? Default = TRUE

  • summary_name: Summary statistics file name (.txt will be appended). Default = "summary_statistics"

  • sup_post: Suppress posterior density plot output in R? Default = FALSE

  • plot_post_save_pdf: Save posterior density plots as pdfs? Default = TRUE

  • plot_post_name: Posterior plot file name(s) (.pdf/.png will be appended) Default = "posterior_density"

  • sup_pairs: Suppress pairs plot output in R? Default = FALSE

  • plot_pairs_save_pdf: Save pairs plot as pdf? Default = TRUE

  • plot_pairs_name: Pairs plot file name (.pdf/.png will be appended) Default = "pairs_plot"

  • sup_xy: Suppress xy/trace plot output in R? Default = TRUE

  • plot_xy_save_pdf: Save xy/trace plot as pdf? Default = FALSE

  • plot_xy_name: XY/trace plot file name (.pdf/.png will be appended) Default = "xy_plot"

  • gelman: Calculate Gelman-Rubin diagnostic test? Default = TRUE

  • heidel: Calculate Heidelberg-Welch diagnostic test? Default = FALSE

  • geweke: Calculate Geweke diagnostic test? Default = TRUE

  • diag_save: Save the diagnostics as a .txt file? Default = TRUE

  • diag_name: Diagnostics file name (.txt will be appended) Default = "diagnostics"

  • indiv_effect: artifact, set to FALSE

  • plot_post_save_png: Save posterior density plots as pngs? Default = FALSE

  • plot_pairs_save_png: Save pairs plot as png? Default = FALSE

  • plot_xy_save_png: Save xy/trace plot as png? Default = FALSE

  • diag_save_ggmcmc: Save ggmcmc diagnostics as pdf? Default = TRUE

  • return_obj Return ggplot objects for later modification? Default = FALSE

Value

named list of three data frames (one each for Gelman, Heidelberg-Welch, and Geweke), but only if return_obj = TRUE