Added MLE calculation to get initial values in MCMC and plots