topic Re: Extract upper triangular elements using SAS IML in New SAS User

Extract upper triangular elements using SAS IML

Deva_123 — Thu, 18 Jul 2019 18:57:19 GMT

Hi All,

I am trying extract upper triangular elements using SAS/IML, I am getting the result incorrectly.

proc iml;
corr = {0 -26 44 34 -85
26 0 65 -56 43
-44 -65 0 21 7
-34 56 -21 0 21
85 -43 -7 -21 0};
r = row(corr);
c = col(corr);
upperTri = loc(r < c);
v = corr[upperTri];
print v;

Desired output :

-26

-85

-56

But I am getting output as :

-26
44
34
-85
26
0
65
-56
43
-44
-65
0
21
7
-34
56
-21
0
21
85
-43
-7
-21
0

Please help me!

Re: Extract upper triangular elements using SAS IML

PaigeMiller — Thu, 18 Jul 2019 19:10:36 GMT

I think this identifies the problem

proc iml;
corr = {0 -26 44 34 -85
26 0 65 -56 43
-44 -65 0 21 7
-34 56 -21 0 21
85 -43 -7 -21 0}; 
print corr;
show corr;

Re: Extract upper triangular elements using SAS IML

koyelghosh — Thu, 18 Jul 2019 19:11:44 GMT

@Deva_123I just reshaped it into a square matrix (5x5) and got the output as below.

Output

The code is as below.

proc iml;
corr = {0 -26 44 34 -85
26 0 65 -56 43
-44 -65 0 21 7
-34 56 -21 0 21
85 -43 -7 -21 0};
corr = shape(corr,5);
r = row(corr);
c = col(corr);
upperTri = loc(r < c);
v = corr[upperTri];
print v;
run;

Best wishes

Re: Extract upper triangular elements using SAS IML

Ksharp — Fri, 19 Jul 2019 12:45:32 GMT

Paige is right. You missed comma ',' .

And why not post it at IML forum ?


proc iml;
corr = {0 -26 44 34 -85 ,
26 0 65 -56 43 ,
-44 -65 0 21 7 ,
-34 56 -21 0 21 ,
85 -43 -7 -21 0};
r = row(corr);
c = col(corr);
upperTri = loc(r < c);
v = corr[upperTri];
print v;

quit;

OUTPUT:

Re: Extract upper triangular elements using SAS IML

Deva_123 — Mon, 22 Jul 2019 14:00:28 GMT

Thanks a lot for your quick reply! This is working perfectly!