topic Re: How to find frequency of i in each bin of data? in SAS/IML Software and Matrix Computations

How to find frequency of i in each bin of data?

ervinodsingh — Thu, 11 May 2017 11:21:20 GMT

Hi i have data with two variables one is score and other is binary variable "label".

score label

245 1

261 0

250 1

300 1

270 0

I am interested to find counts of 1 in label variable for the given manual bins. I have written code but it is giving only frequency of squares in the bins as.

proc iml;
use work.prtf_data;
read all var {Score} into x;
close work.prtf_data;
 

cutpts = {.M 261 273 283 291 298 305 312 330 341 342 .I};
r = bin(x, cutpts);
 
call tabulate(BinNumber, Freq, r);
lbls = {"< 261" "262-273" "274-283" "284-291" "292-298" "299-305" "306-312" "313-330" "331-341" "> 342"};
print Freq[colname=lbls];

i want my results as

"score" "" 1 freq"

<261 2

262-273 0

.....

....

Please suggest me any solution. Thanks

Re: How to find frequency of i in each bin of data?

Rick_SAS — Thu, 11 May 2017 13:26:53 GMT

So you only want to count the data for which label=1? If so, use a WHERE clause to subset the data when you read it in:

proc iml;
use work.prtf_data where(label=1);
read all var {Score} into x;
close work.prtf_data;

I can't tell from your example if the LABEL variable is numeric or character. If it is character, then use

WHERE(label='1')

For the general question "how do I count the number of observations in uneven bins," see "Bin observations by using custom cut points and unevenly spaced bins."

Re: How to find frequency of i in each bin of data?

Ksharp — Thu, 11 May 2017 13:22:25 GMT

data class;
input score    label;
cards;
245       1
261       0
250       1
300       1
270       0
200       0
;
run;
proc iml;
use class;
read all var {score label} ;
close ;
 

cutpts = {.M 262 .I};
r = bin(score, cutpts);

do i=1 to max(r);
 idx=loc(r=i);

 group=group//i;
 freq=freq//sum(label[idx]);
end;

lbls = t( {"< 261" "262-273"} );
print lbls group freq;
quit;

Re: How to find frequency of i in each bin of data?

ervinodsingh — Mon, 15 May 2017 05:58:11 GMT

Hi Rick your solution also worked for me. Thanks