Thank you! I modified my code and it worked. I appreciate all the help from the forum. ... View more

Thank you! I just had one more question. I understand that when we want to print results within iml we can use colname and rowname to give row and column headings. But that is not possible when we are defining a user defined function using return statement? ... View more

Thank you! I am new to iml and so I am having a bit trouble understanding the syntax. I looked at the documentation and modified the code accordingly but it gave me this error.  proc iml; start mlr(dsname,xNames,yName); use (dsName); read all var xNames into x; read all var yName into y; close; n = nrow(x); /* number of observations */ m= ncol(x); /* number of variables */ x=j(n,1,1)||x; /* adding a column of 1 corresponding to intercept*/ /* Compute Xinv, the inverse of X’X */ return ( result ); finish mlr; r = mlr("sashelp.class", {"Height" "Age"}, "Weight"); print r; quit; ... View more

I tried print r to print all the components of r but it gave this error.   ... View more

Thank you! It gave me a result but it only gave me the output for first matrix Analysis_of_variance. Is it possible to get results for all the matrices using return? They are in my result matrix but they didnot get output. This is the output I got.   ... View more

Thanks! I modified the code but it still gives me this error. I am trying to create a function where I can input any dataset and if I call the function it would still give the output. I am confused as to what is wrong with my code. proc iml; INPUT: dsName = {dataset}; xNames = {'x1' 'x2' 'x3'}; yName ={'y'}; start mlr(dsname,xNames,yName); use (dsName); read all var xNames into x; read all var yName into y; close; n = nrow(x); /* number of observations */ m= ncol(x); /* number of variables */ x=j(n,1,1)||x; /* adding a column of 1 corresponding to intercept*/ /* Compute Xinv, the inverse of X’X and the vector of coefficient estimates Beta. */ xinv=inv(x`*x); beta= xinv*x`*y; finish mlr; run mlr("sashelp.class", {"Height" "Age"}, "Weight"); quit;   ... View more

I am trying to learn iml 🙂 ... View more

proc iml; start mlr(x,y); use dataset; read all var {'x1' 'x2' 'x3'} into x; read all var {'y'} into y; n = nrow(x); /* number of observations */ m= ncol(x); /* number of variables */ x=j(n,1,1)||x; /* adding a column of 1 corresponding to intercept*/ /* Compute Xinv, the inverse of X’X and the vector of coefficient estimates Beta. */ Computations finish mlr;  This is my code. I cannot understand what is wrong with the code. Any help would be appreciated. ... View more