multimodal density plots

[AndrewLJackson]

siberdensityplot() only works for unimodal distributions. At the very least it should check and warn if multimodal, but it should be easy enough to implement the multimodal version. My memory is that hdrcde::hdr.boxplot() does not pass on the argument new = FALSE to plot, but i could be wrong.