//Revenu
revenu = data(find(data(:,13) >=0), 13);

//Taux de mortalité infantile
tauxMortaliteInfantile = data(find(data(:,13) >=0), 7);

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

[a, b] = reglin(revenu', tauxMortaliteInfantile');

//affichage
plot2d(revenu, tauxMortaliteInfantile, -1);
plot(revenu', a*revenu'+b);
xtitle("Taux de mortalité infantile en fonction du revenu", "Revenu national brut ppa par habitant (en $US)", "Taux de mortalité infantile (pour 1000 naissances)");
legend(["Nuage de points", "Régression linéaire"]);