diagn_plot {FBMS}

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 = 10000)

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

Value

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

Examples

result <- gmjmcmc(matrix(rnorm(600), 100), P = 2, gaussian.loglik, NULL, c("p0", "exp_dbl"))
diagnstats <- diagn_plot(result)