“…Numerous R packages offer functionalities to estimate conventional parametric and semiparametric survival models via maximum likelihood estimation (MLE) or Bayesian MCMC methods (Diez, 2013;Therneau, 2019;Wang, Chen, Wang, & Yan, 2019;Zhou, Hanson, & Zhang, 2020). Other R packages focus on estimation of parametric or semi-parametric cure survival models using MLE (Amdahl, 2019;Beger, Hill, Metternich, Minhas, & Ward, 2017;Cai, Zou, Peng, & Zhang, 2012;Han, Zhang, & Shao, 2017). To our knowledge, there is no R package that fits parametric mixture cure models, including the MF survival model, via Bayesian MCMC (e.g., slice sampling) methods that offer a powerful yet flexible tool for estimating such models.…”