Last chance! 50% off unlimited learning
Sale ends in
plot method for class gsPLMIX. It builds a suite of plots, visual convergence diagnostics and credible intervals for the MCMC samples of a Bayesian mixture of Plackett-Luce models. Graphics can be plotted directly into the current working device or stored into an external file placed into the current working directory.
# S3 method for gsPLMIX
plot(x, file = "ggmcmc-output.pdf", family = NA,
plot = NULL, param_page = 5, width = 7, height = 10,
dev_type_html = "png", post_est = "mean", max_scale_radar = NULL,
...)Object of class gsPLMIX returned by the gibbsPLMIX function.
Character vector with the name of the file to be created in the current working directory. Defaults is "ggmcmc-output.pdf". When NULL, plots are directly returned into the current working device (not recommended). This option allows also the user to work with an opened pdf (or other) device. When the file has an html file extension, the output is an Rmarkdown report with the figures embedded in the html file.
Character string indicating the name of the family of parameters to be plotted. A family of parameters is considered to be any group of parameters with the same name but different numerical values (for example w[1], w[2], etc). Default is NA meaning that all the parameters in the chain are plotted. Alternatively, one can choose "w", "p", "log_lik", "deviance" or "objective".
Character vector containing the names of the desired plots. Default is NULL meaning that all the plots and convergence diagnostics are built (see 'Details').
Number of parameters to be plotted in each page. Defaults is 5.
Numeric scalar indicating the width of the pdf display in inches. Defaults is 7.
Numeric scalar indicating the height of the pdf display in inches. Defaults is 10.
Character vector indicating the type of graphical device for the html output. Default is "png". Alternatively, one can choose "svg".
Character string indicating the point estimates of the Plackett-Luce mixture parameters to be computed from the gsPLMIX class object and then plotted in the current working device. Default is "mean". Alternatively, one can choose "median".
Numeric scalar indicating the maximum value on each axis of the radar plot for the support parameter point estimates. Default is NULL meaning that the maximum of the estimated support parameters is used.
Further arguments passed to or from other methods (not used).
Plots of the MCMC samples include histograms, densities, traceplots, running means plots, overlapped densities comparing the complete and partial samples, autocorrelation functions, crosscorrelation plots and caterpillar plots of the 90 and 95% equal-tails credible intervals. Note that the latter are created for the support parameters (when either family=NA or family="p"), for the mixture weights in the case family=NA or family="w"), for the log-likelihood values (when family="log_lik"), for the deviance values (when family="deviance"). Convergence tools include the potential scale reduction factor and the Geweke z-score. These functionalities are implemented with a call to the ggs and ggmcmc functions of the ggmcmc package (see 'Examples' for the specification of the plot argument) and for the objective function values (when family="objective").
By recalling the chartJSRadar function from the radarchart package and the routines of the ggplot2 package, plot.gsPLMIX additionally produces a radar plot of the support parameters and, when
Ashton, D. and Porter, S. (2016). radarchart: Radar Chart from 'Chart.js'. R package version 0.3.1. https://CRAN.R-project.org/package=radarchart
Wickham, H. (2009). ggplot2: Elegant Graphics for Data Analysis. Springer-Verlag New York.
Fernandez-i-Marin, X. (2006). ggmcmc: Analysis of MCMC Samples and Bayesian Inference, Journal of Statistical Software, 70(9), pages 1--20, DOI: 10.18637/jss.v070.i09.
ggs, ggmcmc, chartJSRadar and ggplot
# NOT RUN {
# Not run:
data(d_carconf)
GIBBS <- gibbsPLMIX(pi_inv=d_carconf, K=ncol(d_carconf), G=5, n_iter=30, n_burn=10)
# Not run:
# Plot posterior samples supplied as an gsPLMIX class object
# plot(GIBBS)
# Selected plots of the posterior samples of the support parameters
# plot(GIBBS, family="p", plot=c("compare_partial","Rhat","caterpillar"), param_page=6)
# Selected plots of the posterior samples of the mixture weights
# plot(GIBBS, family="w", plot=c("histogram","running","crosscorrelation","caterpillar"))
# Selected plots of the posterior log-likelihood values
# plot(GIBBS, family="log_lik", plot=c("autocorrelation","geweke"), param_page=1)
# Selected plots of the posterior deviance values
# plot(GIBBS, family="deviance", plot=c("traceplot","density"), param_page=1)
# }
Run the code above in your browser using DataLab