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 ?
[1] n’écrire que ce qui est après le >
2004-2024 © Mathazay - Tous droits réservés
Ce site est géré sous SPIP 3.0.17 [21515] et utilise le squelette EVA-Web 4.2
Dernière mise à jour : lundi 7 septembre 2015