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 =   [(iwi ).(iwi.xi.yi) - (iwi.xi).(i wi.yi)] / [(i wi).(iwi.xi²) - (iwi.x)²]           Pente

a = [(i wi.xi²).(iwi.yi) - (iwi.xi).(iwi.xi.yi)] / [(i wi).(iwi.xi²) - (iwi.x)²]     Ord. à l'Origine

 

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 .