data a;
infile 'C:\Users\gabriela\M2PRO\an2013\andon\card.txt';
input ident age med1$ syst50 diast50 haut50 poids50 chol50 status$ clinique$ med2$ syst62
diast62 chol62 poids62 diagn$ annee$ survie$;
run;


/*   ACP */
proc means data=a;
title 'Statistiques elementaires';
  var age syst50 diast50 poids50 haut50 chol50 syst62 diast62 chol62 poids62; run;
  ODS PDF file="C:\Users\gabriela\M2PRO\an2013\andon\card_sorties_photocopie.pdf" startpage=no;   /* c'est pour mettre dans un fichier pdf les sorties */
  ods graphics on;
  proc princomp data=a out=b plots(ncomp=4 flip)=(pattern(circles= 1.0) score);
    var age syst50 diast50 haut50 poids50 chol50 syst62 diast62 chol62 poids62;
        title 'ACP sur toutes les donnees quantitatives'; run;
         ods graphics off;
		 proc corr data=b;
		   var age syst50 diast50 haut50 poids50 chol50 syst62 diast62 chol62 poids62 prin1-prin4;
		   run;
        proc plot data=b;
        title "Premiere comp fonction de la 2eme, et STATUS SOCIO-ECON";
        plot prin2*prin1=status; run;
                proc plot data=b;
        title "Premiere comp fonction de la 2eme, et STATUS CLINIQUE";
        plot prin2*prin1=clinique; run;
		 ods pdf close;

        proc plot data=b;
        title "Premiere comp fonction de la 2eme, et DIAGNOSTIQUE 1962";
        plot prin2*prin1=med2; run;
proc plot data=b;
        title "Premiere comp fonction de la 2eme, et DIAGNOSTIQUE ISCHEMIC";
        plot prin2*prin1=diagn;
        proc plot data=b;
        title "Premiere comp fonction de la 2eme, et SURVIE";
        plot prin2*prin1=survie; run;
        proc plot data=b;
        title "Premiere comp fonction de la 3eme, et SURVIE";
        plot prin3*prin4=survie; run;
        proc plot data=b;
        title "Premiere comp fonction de la 3eme, et STATUS SOCIO-ECON";
        plot prin3*prin4=status;
                proc plot data=b;
        title "Premiere comp fonction de la 3eme, et STATUS CLINIQUE";
        plot prin3*prin4=clinique;
        proc plot data=b;
        title "Premiere comp fonction de la 3eme, et DIAGNOSTIQUE 1962";
        plot prin3*prin4=diag62;
proc plot data=b;
        title "Premiere comp fonction de la 3eme, et DIAGNOSTIQUE ISCHEMIC";
        plot prin3*prin4=diagn;
run;

/* AFD */

  data a;
infile 'C:\Users\gabriela\M2PRO\an2013\andon\card.txt';
input ident age med1$ syst50 diast50 haut50 poids50 chol50 status$ clinique$ med2$ syst62
diast62 chol62 poids62 diagn$ annee$ survie$;
run;
  proc discrim data=a anova;
  class survie;
  var age syst50 diast50 haut50 poids50 chol50 syst62 diast62 chol62 poids62 ; run;
proc discrim data=a crossvalidate;
  class survie;
  var age  chol62 poids62 ; run;
proc discrim data=a method=npar r=1.75 crossvalidate;
title " Methode nonparamétrique des plus proches voisins " ;
  class survie;
  var age  chol62 poids62 ; run;
  proc discrim data=a method=npar k=7 crossvalidate;
  class survie;
  var age  chol62 poids62 ; run;


/* ACC*/
proc cancorr data=a out=sortie red;
  var syst50 diast50  poids50 chol50;
  with syst62 diast62 chol62 poids62;
run;
proc plot data=sortie;
  plot v1*w1 ="*"$survie;
  run;