/* Estimation d'une densité multivariée */

data densite_m;
do i = 1 to 400;
x = rannor(1);
y = abs(x)**0.5 + rannor(2);
output;
end;run;

data densite_m2;
do i = 1 to 80;
x = rannor(1);
y = abs(x)**0.5 + rannor(2);
output;
end;run;

   ods html;
   ods graphics on;
   proc kde data=densite_m;
      bivar x y/ bivstats plots=surface contour scatterplot;
   run;
  proc kde data=densite_m2;
      bivar x y/ bivstats plots=surface contour scatterplot;
   run;

   ods graphics off;
   ods html close; 








/* Refaire tourner le bloc avec moins des données pour mettre en évidence le besoin d'un
minimum
de données avec les méthodes non paramétriques */