implemented MCMC method that follows Ma-Liu paper, i.e. does not rely on data augmentation technique; this method is however (in its current implementation) significantly slower