hdr.cde function

Calculate highest density regions continuously over some conditioned variable.