Simulated maximum likelihood estimation applied to region-censored multi-dimensional thresholding data