Vous êtes ici : Accueil > Salle du prof > Il est libre Math > Ajustement linéaire avec R : résumé
Publié : 31 mars 2011

Ajustement linéaire avec R : résumé

Avec un exercice de bac pour STG-CGRH [1] :

Deux séries :

> rang<-c(1,6,10,11,12,13)
> conso<-c(8.6,10.4,12.2,11.9,12.1,12.2)
> gx<-mean(rang)
> gy<-mean(conso)
> gx;gy
[1] 8.833333
[1] 11.23333

Puis l’ajustement :

> ajust<-lm(conso~rang) # la fonction d'ajustement y en x
> ajust$coefficients # b et a
(Intercept)        rang
 8.4559157   0.3144246

Et la fonction affine :

> a<-ajust$coefficients[2]
> b<-ajust$coefficients[1]
> y<-function(x){b+a*x} # la droite d'ajustement
> y(16) # une estimation
Intercept
13.48671
> x<-function(y){(y-b)/a}#le calcul d'antécédent
> x(10)
Intercept
4.910825

Enfin la représentation graphique :

> plot(conso~rang) # le nuage y/x
> abline(ajust) # la droite
> points(gx,gy,col="blue",cex=1.5) # le point G, cex=1,5 pour agrandir
> text(gx,gy,"G",pos=3,col="blue") # le nom du point G

Quoi de plus simple ?

Notes

[1n’écrire que ce qui est après le >