topic Re: Obtain different values between two vectors in PROC IML in SAS/IML Software and Matrix Computations

Obtain different values between two vectors in PROC IML

Amy_W — Fri, 28 Feb 2014 04:02:05 GMT

Hi,

I have two vectors like

A = {6,4,3,2,5,1};

B = {2,3};

I want to get another vector which contains the different elements comparing A and B. So I tried using SETDIF function.

C = setdif(A,B);

which returns the vector C is {1 4 5 6}. But, what I want is the not-sorted vector {6 4 5 1}.

Does anyone know how I could get the vector I want?

Thanks very much!

Re: Obtain different values between two vectors in PROC IML

IanWakeling — Fri, 28 Feb 2014 14:47:17 GMT

I am thinking there should be a simpler solution than this:

C = remove( A, loc( ( (repeat(A, 1, nrow(B) )=t(B)) [ ,+] ) > 0) );

But I can't see any. This works for your example, but note that it will retain any repeated values in A, so A={6,4,3,2,5,1,4} would return {6 4 5 1 4}.

Re: Obtain different values between two vectors in PROC IML

Amy_W — Fri, 28 Feb 2014 17:50:50 GMT

Thank you, IanWakeling. I have searched for hours yesterday before I posted the question here and I didn't find anything simpler than your solution here. Actually I would like to retain repeated values in the vector C so thanks again.

Re: Obtain different values between two vectors in PROC IML

Rick_SAS — Fri, 28 Feb 2014 20:46:24 GMT

Use the ELEMENT function:

A = {6,4,3,2,5,1};

B = {2,3};

C = A[ loc(^element(A,B)) ];

Re: Obtain different values between two vectors in PROC IML

IanWakeling — Fri, 28 Feb 2014 22:19:07 GMT

That's perfection. Thank you Rick.

Re: Obtain different values between two vectors in PROC IML

Yao_W — Tue, 04 Mar 2014 17:27:24 GMT

Thank you, Rick. That is awesome.

Re: Obtain different values between two vectors in PROC IML

Rick_SAS — Wed, 19 Mar 2014 14:28:18 GMT

I blogged about this problem: Finding elements in one vector that are not in another vector - The DO Loop

Re: Obtain different values between two vectors in PROC IML

Yao_W — Mon, 14 Apr 2014 16:30:26 GMT

I saw it. I subscribed your blog. Very helpful.