* some code on Gaussian dichotomisation