TP IML

Options formdlim='*';
Title1 'tab';
data TAB;
INPUT x y;
   CARDS;
4 6
2 1
;
RUN;

PROC IML;
   USE tab;
   READ ALL VAR{x y} INTO X;
   PRINT x;
quit;

PROC IML;
   USE tab;
   READ ALL VAR{x y} INTO X;
   PRINT x;

   CREATE tab2 FROM X;
   APPEND FROM X;
   CLOSE tab2;
QUIT;
Title1 'tab2';
PROC PRINT data=tab2;
RUN;

Options formdlim='*';
title1 'matrices X Y, et produit matriciel';
PROC IML;
   X={1 2,3 4};
   Y=T(X);
   Z=X*Y;
   PRINT X, Y, Z;
QUIT;

title1 'matrices X Y, et produit element par element';
PROC IML;
   X={1 2,3 4};
   Y=T(X);
   Z=X#Y;
   PRINT X, Y, Z;
QUIT;

title1 'matrices X Y, et juxtaposition';
PROC IML;
   X={1 2,3 4};
   Y=T(X);
   Z=X||Y;
   PRINT X, Y, Z;
QUIT;

title1 'matrice X et acces à la valeur de la 1ere ligne, 2e colonne';
PROC IML;
   X={1 2,3 4};
   Z=X[1,2];
   PRINT X, Z;
QUIT;

title1 'matrice X et instructions';
PROC IML;
   X={1 2,3 4};
   Sommelign=X[,+];
   Sommecol=X[+,];
   multlign=X[,#];
   multcol=X[#,];
   sum2lig=X[,##];
   sum2col=X[##,];

   PRINT X, Sommelign, Sommecol, multlign, multcol, sum2lig, sum2col;
QUIT;

title1 'matrices X, Y, Id, et instructions';
PROC IML;
   X={1 2,3 4};
   Y=J(3,3,5);
   I=I(4);
   diago=DIAG(X);
   vectdiag=VECDIAG(X);
   eigenvect=EIGVEC(X);
   eigenval=EIGVAL(X);
   latrace_y=TRACE(Y);
   determ_x=DET(X);
   inv_x=INV(x);
   PRINT X, Y, I, diago, vectdiag, eigenvect, eigenval, latrace_y,determ_x,inv_x;
QUIT;

title1 'matrice X : operations et instructions';
PROC IML;
   X={1 2,3 4};
   Xcarre=X**2;
   Xcube=X**3;
   nbcol_X=ncol(X);
   nblig_X=nrow(X);
   PRINT X, Xcarre, Xcube, nbcol_X, nblig_X;
QUIT;

Options formdlim='*' ps=105;
title1 'boucles';
PROC IML;
   X={1 2,3 4};
   nbcol_X=ncol(X);
  nblig_X=nrow(X);
   do i=1 to nbcol_X;
     do j=1 to nblig_X;
       termemat=X[i,j];
       print termemat;
     end;
   end;
QUIT;

Options formdlim='*' ps=105;
title1 'IF...THEN...ELSE';
PROC IML;
  X={1 2,3 4};
  Y={9 8, 7 6, 5 4};
  nbcol_X=ncol(X);
  nblig_X=nrow(X);
  nbcol_Y=ncol(Y);
  nblig_Y=nrow(Y);
  print x;
  print y;
  if ((nblig_X=nblig_Y)&(nblig_X=nblig_Y))
  then Z=X*Y;
  else Z=0;
  print Z;
QUIT;

Options formdlim='*' ps=105;
title1 'IF...THEN...ELSE';
PROC IML;
  X={1 2,3 4};
  Y={9 8, 7 6};
  nbcol_X=ncol(X);
  nblig_X=nrow(X);
  nbcol_Y=ncol(Y);
  nblig_Y=nrow(Y);
  print x;
  print y;
  if ((nblig_X=nblig_Y)&(nblig_X=nblig_Y))
  then Z=X*Y;
  else Z=0;
  print Z;
QUIT;

Options formdlim='*' ps=105;
title1 'generation de var. aleat gaussienne';
PROC IML;
  u=repeat(0,100,1);
  x=normal(u);
  y=normal(u);
  z=x||y;
    PRINT z;
    CREATE tabout FROM z;
   APPEND FROM z;
   CLOSE tabout;
QUIT;

title1 'z en SAS';
PROC PRINT data=tabout;
RUN;

PROC GPLOT data=tabout;
PLOT col1*col2;
RUN;

Options formdlim='*' ps=105;
title1 'generation de var. iid gaussienne';
PROC IML;
   matr0={1 0, 0 1};
   u=repeat(0,100,1);
    x=normal(u);
    y=normal(u);
    z=x||y;
   z2=T(z);
   PRINT z2;
   Zr0=matr0*z2;
   Zr0=T(Zr0);
   PRINT Zr0;
   CREATE tabro0 FROM Z;
    APPEND FROM Z ;
    CLOSE tabro0;
QUIT;

title1 'z en SAS';
PROC PRINT data=tabro0;
RUN;

PROC GPLOT data=tabro0;
PLOT col1*col2;
RUN;