hdrscatterplot function

Scatterplot showing bivariate highest density regions