plot_convdiags: Checking convergence visually

Description Usage Arguments Author(s) See Also Examples

View source: R/mcmc_plots.R

Description

Based on a 'damcmc_res' object, this function will produce many graphs to help assess convergence visually, including running mean plots and autocorrelation plots for all the parameters. This function calls plot_runmean and plot_autocorr for all parameters so we do not have to it individually.

For examples see

http://faculty.missouri.edu/~micheasa/sppmix/sppmix_all_examples.html#plot_convdiags

Usage

1
2
plot_convdiags(fit, burnin = floor(fit$L/10), open_new_window = FALSE,
  maxlag = 100)

Arguments

fit

Object of class damcmc_res or bdmcmc_res.

burnin

Number of initial realizations to discard. By default, it is 1/10 of the total number of iterations.

open_new_window

Open a new window for the plot.

maxlag

The maximum lag value to consider. Default is 100.

Author(s)

Sakis Micheas

See Also

est_mix_damcmc, rmixsurf, plot_runmean, plot_autocorr, rsppmix

Examples

1
2
3
4
5
truemix_surf <- rmixsurf(m = 3, lambda=100, xlim = c(-3,3), ylim = c(-3,3))
plot(truemix_surf)
genPPP=rsppmix(intsurf = truemix_surf, truncate = FALSE)
fit = est_mix_damcmc(pp = genPPP, m = 3)
plot_convdiags(fit)

sppmix documentation built on Jan. 13, 2021, 10:04 p.m.