CALCUL DES COEFFICIENTS (PENTE ET ORDONNEE A L'ORIGINE)
DE LA DROITE DES MOINDRES CARRES PONDERES
L'équation de la Droite des Moindres Carrés Pondérés
(DMCP) est :
y* = a + b.x
où "b" est la pente, et "a" l'ordonnée à l'origine.
Le résidu de l'observation i est :
yi* - yi = (a + b.xi ) - yi
La DMCP minimise la Somme des Carrés des Résidus pondérés
:
SSRw =
i
wi.(yi* - yi)²
=
i
wi.(a + b.xi - yi)²
Les coefficients "a" et "b" sont
obtenus en annulant les dérivées partielles de SSRw par
rapport à "a" et à "b" :
*
SSRw
/
a
= 2.(a.
i
wi + b.
i
wi.xi -
i
wi.yi ) = 0
*
SSRw
/
b
= 2.(a.
i
wi.xi + b.
i
wi.xi² -
i
wi.xi.yi ) = 0
Ce système de deux équations linéaires se résoud facilement :
|
b = [( a = [( |
Ces équations ne sont utiles que si la valeur de chaque poids wi est connue individuellement, c'est à dire si la variance var(yi) est connue en chaque point.
Remarques :
1) Si tous les poids ont la même valeur (homoscedasticité), ces équations redonnent les équations normales de la Régression Linéaire Simple :
b =
i(xi
-
).(yi
-
)
/
i(xi
-
)²
a =
- b.![]()
2) Si les poids sont inversement proportionnels à x, ces équations se réduisent à :
|
b = (I.Sy - J.n) / (I.Sx - n²) Pente a = (J.Sx - n.Sy ) / (I.Sx - n²) Ord. à l'origine |
avec les notations suivantes :
* I =
i
(1/xi) ,
* J =
i
(yi /xi).
* Sx =
i
xi
* Sy =
i
yi
I et J ne sont pas définis si un xi,
disons x0, est nul. Nous laissons en exercice le calcul des
expressions de "a"
de"b" quand w0 ![]()
.