FBMS (version 1.2)

diagn_plot: Plot convergence of best/median/mean/other summary log posteriors in time

Description

Plot convergence of best/median/mean/other summary log posteriors in time

Usage

diagn_plot(
  res,
  FUN = median,
  conf = 0.95,
  burnin = 0,
  window = 5,
  ylim = NULL,
  ...
)

Value

A list of summary statistics for checking convergence with given confidence intervals

Arguments

res

Object corresponding gmjmcmc output

FUN

The summary statistics to check convergence

conf

Which confidence intervals to plot

burnin

How many first populations to skip

window

Sliding window for computing the standard deviation

ylim

Limits for the plotting range; if unspecified, min and max of confidence intervals will be used

...

Additional graphical parameters passed to plot and lines functions, e.g. col, lwd, lty, main, xlab, ylab, ylim

Examples

result <- gmjmcmc(y = matrix(rnorm(100), 100),
                  x = matrix(rnorm(600), 100), 
                  P = 2, 
                  transforms =  c("p0", "exp_dbl"))
diagnstats <- diagn_plot(result)