plot.hdrconf function

Plot HDRs with confidence intervals