AbstractComparison of appropriate models to describe observational data is a fundamental task of science. The Bayesian model evidence, or marginal likelihood, is a computationally challenging, yet crucial, quantity to estimate to perform Bayesian model comparison. We introduce a methodology to compute the Bayesian model evidence in simulation-based inference (SBI) scenarios (often called likelihood-free inference). In particular, we leverage the recently proposed learned harmonic mean estimator and exploit the fact that it is decoupled from the method used to generate posterior samples, i.e. it requires posterior samples only, which may be generated by any approach. This flexibility, which is lacking in many alternative methods for computing the model evidence, allows us to develop SBI model comparison techniques for the three main neural density estimation approaches, including neural posterior estimation, neural likelihood estimation, and neural ratio estimation. We demonstrate and validate our SBI evidence calculation techniques on a range of inference problems, including a gravitational wave example. Moreover, we further validate the accuracy of the learned harmonic mean estimator, implemented in the harmonic software, in likelihood-based settings. These results highlight the potential of harmonic as a sampler-agnostic method to estimate the model evidence in both likelihood-based and simulation-based scenarios.