Implementation of sampler described in Kemp et al. (2007)