Learn everything about Analytics

Next step in the world of SAS IML

SHARE
, / 0

In the last article on IML (here) , we introduced you to the world of Matrix language on SAS. We also talked about some basic commands on IML. This article will take you a step further in the same world. After reading this article, you can take up problem solving using SAS IML. In some of the coming articles on IML, we will also take some examples of cases where IML will be very handy.

A Quick Recap 

Following code will help you recall all the commands we discussed in the last article. Try guessing the output of the commands before looking at the output below.

proc iml;

reset log print;

A = {1 2,3 4};

B=t(A);

C= inv(A);

D = A*C;

run;

iml1 iml2Following are the description of each matrix defined in this code :

A : Original Matrix

B : Transpose of Matrix A

C : Inverse of Matrix A

D : Matrix multiplication of A and B which should ideally be an identity matrix (Using property of Inverse).

Now let’s add a few more types of command on IML in our dictionary.

Replacing elements of a Matrix

In many cases we need to manipulate individual datapoint to fix outliers or missing values. To do this operation in IML, you can use the following method (B has one value replaced while C has an entire column replaced) :

proc iml;

reset log print;

A = {1 2,3 4};

A[1,2] = 0;

B=A;

A[2, ] = 3;

C = A;

run;

iml3 iml4Summarize operators

Operation like normalization of probabilities, needs you to collapse one dimension of the matrix.  Following code can be used to collapse row or column :

proc iml;

reset log print;

A = {1 2,3 4};

c = a[+,];

c=a[##,];

d=a[,:];

run;

iml5

 

iml6

 

Control statements

For making any logic, we always need control statements/conditional statements. IML supports all IF-THEN, DO-WHILE, DO-UNTIL statements. Following are a commands which will help you write conditional statements :

proc iml;

reset log print;

A = {1 2,3 4};

B= {2 3,4 5};

if A < B then

msg = “all (A<B)”;

else msg = “some elements of A is greater than or equal to the corresponding element of B”;

run;

iml7 iml8

Note that in IML, the condition is applied on every element. If element wise operation at all row and column gives a TRUE output, the final output is TRUE. Following is an example of a do-end statement :

proc iml;

x=1:5;

sum1=0;

do i = 1 to ncol(x);

sum1=sum1+x[i];

end;

sum2 = sum(x);

print sum1 sum2;

run;

iml9 iml10

Concatenation Operators

Very often we try merging or appending the datasets. IML provides very easy way to do these operations. Following commands can be used for merging (horizontal concatenate) and appending (vertical concatenate) :

proc iml;

a = {1,2,3,4,5};

b={3,5,4,1,3};

c={0,1,0,0,1};

x=a||b||c;

y=a//b;

print x;

print y;

run;

iml11 iml12

Logical Operators

Multiple conditions can be used together to form complex conditions on IML. Following is a simple code which demonstrates the same :

proc iml;

do x=-3 to 3;

if (x>=-2) & (x<=2) then

y = sqrt(4-x##2);

else

y=”Unable to compute the function. “;

print y;

end;

run;

iml13 iml14

End Notes

This command list was not exhaustive but will definitely enable you to fulfill most of the purposes. In our following articles we will pick up some iterative problems like formulation of Newton-Raphson, Solving Linear Regression, Linear Programming etc. to demonstrate the usage of IML programming.

Did you find the article useful? Have you used IML before?  Do you plan to use IML in any of your business problems? If yes, share with us how IML programming will simplify or enable your analysis.

If you like what you just read & want to continue your analytics learningsubscribe to our emailsfollow us on twitter or like our facebook page.

 

Leave A Reply

Your email address will not be published.

Join world’s fastest growing Analytics Community
Receive awesome tips, guides, infographics and become expert at: