The function hdrbw() for highest density region estimation bandwidth selection, within the R package:

Hyndman, R.J. (2018). hdrcde 3.2. Highest density regions and conditional density estimation. R package. cran.r-project.org