Blame

e0287e lugonthier 2026-07-02 15:52:15
Add new content and images for machine learning and mathematics modules - Added images for regularization and high-dimensional inference. - Introduced Support Vector Machines (SVM) module with detailed explanations and images. - Created Decision Trees and Ensemble Methods module with comprehensive content and illustrations. - Added a Mathematics overview module and a refresher on mathematical concepts essential for machine learning. - Included SVG diagrams for Bayes' rule and multivariate Gaussian distribution.
1
# 4. Probabilistic formulation
2
3
Probability is the language machine learning uses to handle uncertainty. This module sets out the rules for discrete and continuous variables, takes a first look at information theory, shows the Bayesian way of turning probabilities into decisions, and defines the two estimation principles the course returns to again and again: maximum likelihood and maximum a posteriori.
4
5
**Objectives**
6
- State the rules of probability for discrete and continuous variables.
7
- Relate joint, conditional, and marginal probabilities by the sum and product rules and by Bayes' rule.
8
- Measure uncertainty with entropy, cross-entropy, and the Kullback-Leibler divergence.
9
- Make the decision that minimizes expected loss, and recover the maximum-a-posteriori classifier.
10
- Define the maximum-likelihood and maximum-a-posteriori estimators.
11
12
## 4.1 Probability, discrete and continuous
13
14
A random variable takes values with probabilities that are non-negative and sum or integrate to one. A discrete variable has a probability mass function, a continuous one a probability density function:
15
16
$$\boxed{ \sum_x p(x) = 1 \qquad \int p(x)\, dx = 1, \quad p(x) \ge 0 }$$
17
18
For a continuous variable, probability attaches to intervals through an integral, $P(a \le X \le b) = \int_a^b p(x)\, dx$, not to single points.
19
20
## 4.2 Joint, conditional, and Bayes
21
22
Two variables have a joint distribution $p(x, y)$. Summing (or integrating) out one variable gives the marginal, the sum rule, and the joint factors into a conditional times a marginal, the product rule:
23
24
$$\boxed{ p(x) = \sum_y p(x, y) \qquad p(x, y) = p(y \mid x)\, p(x) }$$
25
26
Rearranging the product rule both ways gives Bayes' rule, which flips a conditional:
27
28
$$\boxed{ p(y \mid x) = \frac{p(x \mid y)\, p(y)}{p(x)} }$$
29
30
Two variables are independent when the joint is the product of the marginals, $p(x, y) = p(x)\, p(y)$.
31
32
## 4.3 A little information theory
33
34
The entropy of a distribution measures its uncertainty, the average number of bits needed to describe an outcome:
35
36
$$\boxed{ H(X) = -\sum_x p(x)\log p(x) }$$
37
38
![Binary entropy](/en/Machine%20Learning/04%20Probabilistic%20formulation/a/entropy.png)
39
40
*For a two-outcome variable the entropy is largest at $p = 0.5$, where the outcome is hardest to predict, and zero when one outcome is certain.*
41
42
The cross-entropy measures the cost of using a model $q$ when the truth is $p$, and the Kullback-Leibler divergence measures how far $q$ sits from $p$:
43
44
$$\boxed{ H(p, q) = -\sum_x p(x)\log q(x) \qquad D_{\mathrm{KL}}(p \,\|\, q) = \sum_x p(x)\log\frac{p(x)}{q(x)} \ge 0 }$$
45
46
*Remark:* minimizing the cross-entropy between the true labels and a model's predictions is the same as maximizing the likelihood of those labels. This is why classification networks minimize cross-entropy, a thread picked up in later modules.
47
48
## 4.4 Bayesian decision theory
49
50
To classify an input $x$, the Bayesian rule uses the posterior over classes. Under the 0-1 loss, the decision that minimizes the expected loss is simply the most probable class, and because the posterior is proportional to the class-conditional density times the prior, it can be computed either way:
51
52
$$\boxed{ \hat{y} = \arg\max_y \; p(y \mid x) = \arg\max_y \; p(x \mid y)\, p(y) }$$
53
54
![Bayesian decision between two classes](/en/Machine%20Learning/04%20Probabilistic%20formulation/a/bayes-decision.png)
55
56
*Each class contributes its density scaled by its prior, and the decision boundary falls where the two are equal. On each side the class with the larger posterior wins.*
57
58
*Remark:* this is the optimal classifier, called the Bayes classifier. Every method later in the course is, in effect, an attempt to approximate these posteriors from data.
59
60
## 4.5 Maximum likelihood and maximum a posteriori
61
62
We rarely know the true distribution, so we estimate its parameters $\theta$ from data. Maximum likelihood picks the $\theta$ that makes the observed data most probable, usually maximized as a sum of log-likelihoods over the $m$ examples:
63
64
$$\boxed{ \theta_{\mathrm{MLE}} = \arg\max_\theta \sum_{i=1}^{m} \log p(x^{(i)} \mid \theta) }$$
65
66
Maximum a posteriori instead maximizes the posterior, which multiplies the likelihood by a prior on $\theta$:
67
68
$$\boxed{ \theta_{\mathrm{MAP}} = \arg\max_\theta \; p(D \mid \theta)\, p(\theta) }$$
69
70
*Remark:* maximum a posteriori is maximum likelihood plus a prior. A Gaussian prior on $\theta$ becomes an L2 penalty and a Laplace prior an L1 penalty, which is exactly the regularization of a later module. With abundant data the prior washes out and the two estimators agree.
71
72
*The next module turns these principles into concrete loss functions and the gradient descent that minimizes them.*
73
74
---
75
Next: [Linear regression](/en/Machine%20Learning/05%20Linear%20regression) · [Course overview](/en/Machine%20Learning)