Bayesian Mode Inference R package BayesMultiMode version 1.0.0 A two-step Bayesian approach for mode inference (Baştürk et al., 2026) <doi:10.18637/jss.v116.i03>. First, a mixture distribution is fitted on the data using a sparse finite mixture (SFM) Markov chain Monte Carlo (MCMC) algorithm. The number of mixture components does not have to be known; the size of the mixture is estimated endogenously through the SFM approach. Second, the modes of the estimated mixture at each MCMC draw are retrieved using algorithms specifically tailored for mode detection. These estimates are then used to construct posterior probabilities for the number of modes, their locations and uncertainties, providing a powerful tool for mode inference. Software Depends: R (>= 3.5.0) Imports: assertthat, bayesplot, dplyr, ggplot2 (>= 3.4.0), ggpubr, gtools, magrittr, MCMCglmm, mvtnorm, posterior, sn, stringr, tidyr, Rdpack Suggests: testthat (>= 3.0.0) Paul Labonne <labonnepaul@gmail.com> Comprehensive R Archive Network (CRAN) Nalan Baştürk [aut], Jamie Cross [aut], Peter de Knijff [aut], Lennart Hoogerheide [aut], Paul Labonne [aut, cre], Herman van Dijk [aut] GPL (>= 3) 2026-06-01 application/tgz https://CRAN.R-project.org/package=BayesMultiMode doi:10.32614/CRAN.package.BayesMultiMode