This function computes the maximum likelyhood estimator of the given distribution parameters. Taken from the MASS package.
Pierre-Matthieu Pair, Régis Lebrun.