crossplot function

crossplot: Plot the median and HPD interval bars for pairs of distribution