Blame
|
1 | # 3. 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 | ## 3.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 | ## 3.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 | ## 3.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 |  |
||||||
|
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 | ## 3.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 |  |
||||||
|
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 | ## 3.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 | In supervised learning the model parameterizes the conditional $p(y \mid x; \theta)$, so the same principle applies to the conditional likelihood of the targets: |
||||||
| 67 | ||||||||
| 68 | $$\boxed{ \ell(\theta) = \sum_{i=1}^{m} \log p\!\left(y^{(i)} \mid x^{(i)}; \theta\right) }$$ |
|||||||
| 69 | ||||||||
| 70 | Maximizing $\ell$ is the same as minimizing the cost $J(\theta) = -\ell(\theta)$, so the likelihood view and the cost-minimization view of [General concepts](/en/Machine%20Learning/02%20General%20concepts) are two faces of one objective. |
|||||||
| 71 | ||||||||
|
72 | Maximum a posteriori instead maximizes the posterior, which multiplies the likelihood by a prior on $\theta$: |
||||||
| 73 | ||||||||
| 74 | $$\boxed{ \theta_{\mathrm{MAP}} = \arg\max_\theta \; p(D \mid \theta)\, p(\theta) }$$ |
|||||||
| 75 | ||||||||
|
76 | *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 the next module. With abundant data the prior washes out and the two estimators agree. |
||||||
|
77 | |||||||
|
78 | *The next module turns these principles into a first concrete model: linear regression, where maximum likelihood and maximum a posteriori both land on closed-form fits.* |
||||||
|
79 | |||||||
| 80 | --- |
|||||||
|
81 | Next: [Linear regression](/en/Machine%20Learning/04%20Linear%20regression) · [Course overview](/en/Machine%20Learning) |
||||||
