Most of my posteriors are turning out to be bi-modal and using the mean (MAP) estimate is not possible. Instead, I want to use the mode estimate of the posterior samples. Is there an easy way to achieve this? I can plot the kde of posteriors using sns.kdeplot() but can’t seem to find an easy way to get the actual mode estimates. Thanks!