Probability
Foundations
Maximum Likelihood (MLE)

Maximum Likelihood Estimation (MLE)

🎯

Probability is about predicting data given a model. Statistics is about predicting the model (parameters) given the data. MLE is the primary tool for this reversal.

The Core Concept: The Likelihood Function

Imagine you have a coin. You don't know if it's fair. You flip it 10 times and get 8 Heads. What is the most "likely" probability of heads (pp) for this coin?

In probability, we define the Probability Function P(x∣θ)P(x \mid \theta), which tells us the chance of seeing data xx given parameter θ\theta. In statistics, we flip this into the Likelihood Function L(θ∣x)L(\theta \mid x):

L(θ∣x)=P(x∣θ)L(\theta \mid x) = P(x \mid \theta)

The math is the same, but the perspective changes: the data xx is now fixed (the 8 heads you saw), and the parameter θ\theta is the variable we are trying to find.

The Recipe for MLE

Finding the maximum of a function usually involves calculus. However, because most probability functions involve products (due to Independence), the derivatives get messy. We use a trick: the Log-Likelihood.

  1. Define the Likelihood: Write the joint probability of your data points.
  2. Take the Natural Log: ℓ(θ)=ln⁡(L(θ))\ell(\theta) = \ln(L(\theta)). (Since ln⁡\ln is a monotonically increasing function, the maximum of the log is the same as the maximum of the original).
  3. Differentiate: Find ddθℓ(θ)\frac{d}{d\theta} \ell(\theta).
  4. Solve for Zero: Set the derivative to zero and solve for θ^\hat{\theta}.

Why "Maximum" Likelihood?

If we assume our coin has p=0.5p=0.5, the chance of getting 8/10 heads is quite low (~4%). If we assume p=0.8p=0.8, the chance is much higher (~30%). MLE says: "The best estimate for the world is the one that makes our observed reality most probable."

Real-World Connection: The Normal Distribution

If you take a set of measurements x1,x2,…,xnx_1, x_2, \dots, x_n and assume they come from a Normal Distribution, the MLE for the mean (μ\mu) is simply the Sample Average:

μ^MLE=1n∑xi\hat{\mu}_{MLE} = \frac{1}{n} \sum x_i

This is why we use the average so often—it is mathematically the most likely center of a Gaussian world.

Test Your Knowledge

Example: MLE for a Bernoulli Trial

You flip a coin nn times and observe kk successes. Find the MLE for the probability of success pp.

View Step-by-Step Solution
  1. Likelihood: L(p)=pk(1−p)n−kL(p) = p^k (1-p)^{n-k}
  2. Log-Likelihood: ℓ(p)=kln⁡(p)+(n−k)ln⁡(1−p)\ell(p) = k \ln(p) + (n-k) \ln(1-p)
  3. Derivative: ddpℓ(p)=kp−n−k1−p\frac{d}{dp} \ell(p) = \frac{k}{p} - \frac{n-k}{1-p}
  4. Solve: Set to 0: kp=n−k1−p  ⟹  k(1−p)=p(n−k)\frac{k}{p} = \frac{n-k}{1-p} \implies k(1-p) = p(n-k) k−kp=np−kp  ⟹  k=np  ⟹  p^=knk - kp = np - kp \implies k = np \implies \hat{p} = \frac{k}{n}

The MLE for pp is simply the proportion of successes observed.