estimate rho with ML

To increase numerical stability, on all thresholded models, first
estimate the rho term using ML, and then fix rho in the MCMC using ML.