//Taux de natalité
tauxNatalite = data(:, 4);

//Taux de mortalité infantile
tauxMortaliteInfantile = data(:, 7);

//regression
[x, y] = reglin(tauxNatalite', tauxMortaliteInfantile');

//affichage
plot2d(tauxNatalite, tauxMortaliteInfantile, -1);
plot(tauxNatalite', x*tauxNatalite' + y);
xtitle("Taux de mortalité infantile en fonction du taux de natalité", "Taux de natalité (pour 1000 habitants)", "Taux de mortalité infantile (pour 1000 naissances)");
legend(["Nuage de points", "Régression linéaire"]);

//coefficient de correlation
r = correl(tauxNatalite, tauxMortaliteInfantile);
mprintf("Coefficient de correlation : %f\n", r);