print_joints() now does sth promising...
Fixed in_bin to actually work if you have different number of bins to posterior samples....Also added some typically dodgy ggplot2 code to plot joint distro