how to model a bimodal distribution