Updated Holmes and Held, most of the important changes have been incorporated into the hhProbit function in bayesRegress