Statistiques à deux variables avec Xcas

Voici comment on peut résoudre un exercice classique sur les statistiques à deux variables.

Enoncé


Soit la série statistique suivante :

$X$ 0 5 10 15 20 25 30
$Y$ 18 21 25 30 36 42 50

Donner les principaux paramètres de cette série (moyenne, écart-type, …)

Résolution à l’aide de Xcas


X:=[0,5,10,15,20,25,30]

$$[0,5,10,15,20,25,30]$$

Y:=[18,21,25,30,36,42,50]

$$[18,21,25,30,36,42,50]$$

Donner les principaux paramètres de cette série (moyenne, écart-type, …)

var2stat_e(X,Y):= local x,y,S; S:=zip((x,y)->x*y,X,Y); print("- N : "+size(X)+";"); print("------------- première série x ------------------"); print("- Moyenne des x: "+normal(moyenne(X))+";"); print("- Variance des x: "+normal(variance(X))+";"); print("- Ecart-type des x: "+normal(ecart_type(X))+";"); print("- Somme des x : "+normal(sum(X))+";"); print("- Somme des x^2 : "+normal(sum(X^2))+";"); print("- Médiane des x : "+normal(median(X))+";"); print("- Q1 des x: "+normal(quartile1(X))+";"); print("- Q3 des x: "+normal(quartile3(X))+";"); print("- Min des x: "+normal(min(X))+";"); print("- Max des x: "+normal(max(X))+";"); print("------------- Deuxième série y ------------------"); print("- Moyenne des y: "+normal(moyenne(Y))+";"); print("- Variance des y: "+normal(variance(Y))+";"); print("- Ecart-type des y: "+normal(ecart_type(Y))+";"); print("- Somme des y : "+normal(sum(Y))+";"); print("- Somme des y^2 : "+normal(sum(Y^2))+";"); print("- Médiane des y: "+normal(median(Y))+";"); print("- Q1 des y : "+normal(quartile1(Y))+";"); print("- Q3 des y: "+normal(quartile3(Y))+";"); print("- Min des y: "+normal(min(Y))+";"); print("- Max des y: "+normal(max(Y))+";"); print("------------- Paramètres communs ------------------"); print("- Somme des xy : "+normal(sum(S))+";"); print("- Covariance des xy : "+normal(covariance(X,Y))+";"); :;

// Parsing var2stat_e

// Success compiling var2stat_e

var2stat_e(X,Y)

“- N : 7;”

“————- première série x ——————”

“- Moyenne des x: 15;”

“- Variance des x: 100;”

“- Ecart-type des x: 10;”

“- Somme des x : 105;”

“- Somme des x^2 : 2275;”

“- Médiane des x : 15.0;”

“- Q1 des x: 5.0;”

“- Q3 des x: 20.0;”

“- Min des x: 0;”

“- Max des x: 30;”

“————- Deuxième série y ——————”

“- Moyenne des y: 222/7;”

“- Variance des y: 5666/49;”

“- Ecart-type des y: sqrt(5666)/7;”

“- Somme des y : 222;”

“- Somme des y^2 : 7850;”

“- Médiane des y: 30.0;”

“- Q1 des y : 21.0;”

“- Q3 des y: 36.0;”

“- Min des y: 18;”

“- Max des y: 50;”

“————- Paramètres communs ——————”

“- Somme des xy : 4075;”

“- Covariance des xy : 745/7;”

Si on veut les valeurs approchées :

var2stat_a(X,Y):= local x,y,S; S:=zip((x,y)->x*y,X,Y); print("- N : "+size(X)+";"); print("------------- première série x ------------------"); print("- Moyenne des x: "+evalf(moyenne(X))+";"); print("- Variance des x: "+evalf(variance(X))+";"); print("- Ecart-type des x: "+evalf(ecart_type(X))+";"); print("- Somme des x : "+evalf(sum(X))+";"); print("- Somme des x^2 : "+evalf(sum(X^2))+";"); print("- Médiane des x : "+eval(median(X))+";"); print("- Q1 des x: "+evalf(quartile1(X))+";"); print("- Q3 des x: "+evalf(quartile3(X))+";"); print("- Min des x: "+evalf(min(X))+";"); print("- Max des x: "+evalf(max(X))+";"); print("------------- Deuxième série y ------------------"); print("- Moyenne des y: "+evalf(moyenne(Y))+";"); print("- Variance des y: "+evalf(variance(Y))+";"); print("- Ecart-type des y: "+evalf(ecart_type(Y))+";"); print("- Somme des y : "+evalf(sum(Y))+";"); print("- Somme des y^2 : "+evalf(sum(Y^2))+";"); print("- Médiane des y: "+eval(median(Y))+";"); print("- Q1 des y : "+evalf(quartile1(Y))+";"); print("- Q3 des y: "+evalf(quartile3(Y))+";"); print("- Min des y: "+evalf(min(Y))+";"); print("- Max des y: "+evalf(max(Y))+";"); print("------------- Paramètres communs ------------------"); print("- Somme des xy : "+evalf(sum(S))+";"); print("- Covariance des xy : "+evalf(covariance(X,Y))+";"); :;

// Parsing var2stat_a

// Success compiling var2stat_a

var2stat_a(X,Y)

“————- première série x ——————”

“- Moyenne des x: 15.0;”

“- Variance des x: 100.0;”

“- Ecart-type des x: 10.0;”

“- Somme des x : 105.0;”

“- Somme des x^2 : 2275.0;”

“- Médiane des x : 15.0;”

“- Q1 des x: 5.0;”

“- Q3 des x: 20.0;”

“- Min des x: 0.0;”

“- Max des x: 30.0;”

“————- Deuxième série y ——————”

“- Moyenne des y: 31.7142857143;”

“- Variance des y: 115.632653061;”

“- Ecart-type des y: 10.753262438;”

“- Somme des y : 222.0;”

“- Somme des y^2 : 7850.0;”

“- Médiane des y: 30.0;”

“- Q1 des y : 21.0;”

“- Q3 des y: 36.0;”

“- Min des y: 18.0;”

“- Max des y: 50.0;”

“————- Paramètres communs ——————”

“- Somme des xy : 4075.0;”

“- Covariance des xy : 106.428571429;”