Now calculates MLE using R