Implement functions to calculate maximum likelihood estimates for parameters