ATLAS Ouvrir dans l'atlas →

ExerciceMaths Sup – MPSI, chapitre 21 — Dénombrement et probabilités · Variance et covariance

Exercice 84

Calcul

Soit XX et YY deux variables aléatoires réelles, avec V(X)>0\mathrm{V}(X) > 0. Déterminer les réels aa et bb qui minimisent la quantité E((Y−(aX+b))2)\mathrm{E}\Big(\big(Y - (aX + b)\big)^{2}\Big), et comparer avec la droite des moindres carrés évoquée dans le cours après le théorème 12.

Corrigé

Pour une variable ZZ et un réel bb, la linéarité de l'espérance donne E((Z−b)2)=E(Z2)−2bE(Z)+b2=V(Z)+(E(Z)−b)2,\mathrm{E}\big((Z - b)^{2}\big) = \mathrm{E}(Z^{2}) - 2b\mathrm{E}(Z) + b^{2} = \mathrm{V}(Z) + \big(\mathrm{E}(Z) - b\big)^{2}, par la formule de König-Huygens. On l'applique à Z=Y−aXZ = Y - aX, dont la variance se développe par bilinéarité (proposition 23) : V(Y−aX)=V(Y)−2aCov⁡(X,Y)+a2V(X).\mathrm{V}(Y - aX) = \mathrm{V}(Y) - 2a\operatorname{Cov}(X, Y) + a^{2}\mathrm{V}(X). On obtient, en posant a∗=Cov⁡(X,Y)V(X)a^{\ast} = \frac{\operatorname{Cov}(X, Y)}{\mathrm{V}(X)} et en complétant le carré, E((Y−(aX+b))2)=V(X) (a−a∗)2+(b−E(Y)+aE(X))2+V(Y)−Cov⁡(X,Y)2V(X).\mathrm{E}\Big(\big(Y - (aX + b)\big)^{2}\Big) = \mathrm{V}(X)\,(a - a^{\ast})^{2} + \big(b - \mathrm{E}(Y) + a\mathrm{E}(X)\big)^{2} + \mathrm{V}(Y) - \frac{\operatorname{Cov}(X, Y)^{2}}{\mathrm{V}(X)}. Les deux premiers termes sont positifs, puisque V(X)>0\mathrm{V}(X) > 0, et ils sont nuls simultanément si et seulement si a=a∗a = a^{\ast} et b=E(Y)−a∗E(X)b = \mathrm{E}(Y) - a^{\ast}\mathrm{E}(X). Ainsi le minimum est atteint en un unique couple : a=Cov⁡(X,Y)V(X),b=E(Y)−a E(X),min⁡=V(Y)−Cov⁡(X,Y)2V(X)\boxed{a = \frac{\operatorname{Cov}(X, Y)}{\mathrm{V}(X)}, \qquad b = \mathrm{E}(Y) - a\,\mathrm{E}(X), \qquad \min = \mathrm{V}(Y) - \frac{\operatorname{Cov}(X, Y)^{2}}{\mathrm{V}(X)}} Ce minimum est positif : on retrouve l'inégalité de Cauchy-Schwarz Cov⁡(X,Y)2⩽V(X)V(Y)\operatorname{Cov}(X, Y)^{2} \leqslant \mathrm{V}(X)\mathrm{V}(Y) du théorème 12, et il est nul si et seulement si P(Y=aX+b)=1\mathrm{P}(Y = aX + b) = 1 (théorème 11 (iii)). Si V(Y)>0\mathrm{V}(Y) > 0, il s'écrit V(Y)(1−ρ2)\mathrm{V}(Y)(1 - \rho^{2}), où ρ\rho est le coefficient de corrélation : plus ∣ρ∣|\rho| est proche de 11, mieux YY est approchée par une fonction affine de XX.

Si l'on munit un nuage de points (x1,y1)(x_{1}, y_{1}), …, (xm,ym)(x_{m}, y_{m}) de la probabilité uniforme, et que XX et YY sont l'abscisse et l'ordonnée d'un point tiré au hasard, la quantité minimisée est 1m∑i(yi−(axi+b))2\frac{1}{m}\sum_{i} \big(y_{i} - (ax_{i} + b)\big)^{2}, et l'on retrouve exactement la droite des moindres carrés de la proposition 14 du chapitre 19 : pente a=Cov⁡(X,Y)V(X)a = \frac{\operatorname{Cov}(X, Y)}{\mathrm{V}(X)}, comme le dit la remarque qui suit le théorème 12, et droite passant par le point moyen (xˉ,yˉ)=(E(X),E(Y))(\bar{x}, \bar{y}) = \big(\mathrm{E}(X), \mathrm{E}(Y)\big).

À retenir. La minimisation se fait sans dériver, en complétant deux carrés : d'abord en bb, à aa fixé, puis en aa. C'est la méthode du chapitre 19, où la droite des moindres carrés était une projection orthogonale ; ici, aX+baX + b est la meilleure approximation de YY par une fonction affine de XX, au sens de l'écart quadratique moyen.

Ouvrir dans l'atlas interactif →