Proc IML;
a=1;
print a;
A={1};
print A;
B={1 2 3 4 5};
print B;
C={1, 2, 3, 4, 5};
print C;
D=B*C;
print D;
E=a*C;
print E;

A= {1 0, 1 0};
B = {0 1, 0 1};
C = A + B;
Print C;

A= {1 1, 1 1};
B = {2 3, 2 3};
C = A * B;
Print C;
D = A * B * A;
Print D;
E = {2 3, 5 6};
Print E;
F = t(E);
Print F;
E= {2 3, 5 6};
Print E;
F = t(E);
Print F;
G = det(E);
Print G;
H = inv(E);
Print H;
quit;


data example;
infile datalines;
input X;
datalines;
1
2
3
0
1
2
;
run;
proc means data=example; run;



proc iml;
X = {1, 2, 3, 0, 1, 2};
print X;
/* counting the number of rows in the matrix = number of n */
n=nrow(X);
print n;
xbar = sum(X)/(n);
print xbar;
Y = X - xbar;
print Y;
VAR = t(Y)*Y;
print VAR;
SD = sqrt(VAR/(n-1));
print SD;
quit;



data multivariate;
infile datalines;
input X1 X2;
datalines;
42 4
52 5
48 4
58 3
;
run;
PROC CORR DATA=multivariate COV OUTP=cov;
  VAR X1 X2;
  RUN;


proc iml;
X = {42 4, 52 5, 48 4, 58 3};
print X;
n=nrow(X);
print n;
/* Generating an unit matrix */
U = J(1, 4, 1);
print U; 
S = U*X;
print S;
Xbar = S/n;
print Xbar;
/* Generating an unit matrix */
Q = J(4, 1, 1);
print Q;
D = X-(Q*Xbar);
print D;
Dt = t(D);
print Dt;
Sumsquare = t(D)*D;
print Sumsquare;
S = Sumsquare/(n-1);
print S;
quit;