Sure; but the question is how to plot contours nicely without binning and/or sampling ![]()
I realized that to determine HPD cutoffs, I need to sample anyway. Then I can make a really dense grid, and just use contour.
There are some hints here about what I want, but ATM just using contour on a very dense grid is probably the easiest solution.