powered by
Plot convergence of best/median/mean/other summary log posteriors in time
diagn_plot(res, FUN = median, conf = 0.95, burnin = 0, window = 5, ylim = NULL)
A list of summary statistics for checking convergence with given confidence intervals
Object corresponding gmjmcmc output
The summary statistics to check convergence
which confidence intervals to plot
how many first populations to skip
sliding window for computing the standard deviation
limits for the plotting range, if unspecified, min and max of confidence intervals will be used
result <- gmjmcmc(matrix(rnorm(600), 100), P = 2, gaussian.loglik, NULL, c("p0", "exp_dbl")) diagnstats <- diagn_plot(result)
Run the code above in your browser using DataLab