I want to manually estimate the highest posterior density (HPD) interval using R base functions. The main reason is that I conceptually understand the difference but want to
