DATA Step, Macro, Functions and more

Randomization algorithm

Accepted Solution Solved
Reply
Frequent Contributor
Posts: 81
Accepted Solution

Randomization algorithm

Hi,

 

I have a rolled up data set that contains a metric that counts the number of members and Dollars spent:

AgeGroup Region Members Spend

20-24 West 1 50

20-24 East 4 250

20-24 South 2 80

...

For any cell where Members LE 3 I randomize to 0 or 3 and adjust the spend accordingly:

AgeGroup Region Members Spend

20-24 West 0 0

20-24 East 4 250

20-24 South 3 120

 

At the end I evaluate the overall % difference to calclate the possible bias from new to old:

MemberBias = 7/7-1 = 0%

SpendBias = 370/380-1=-2.6%

 

Here's where i'm stuck.

IF MemberBias OR SpendBias is GE +/- 1% re-run the randomization until both is under +-1% then keep that final dataset.

 

I'm not sure how to structure the logic in my code to do this, any help greatly appreciated.

Thanks!


Accepted Solutions
Solution
‎09-28-2015 09:16 AM
SAS Employee
Posts: 24

Re: Randomization algorithm

I suggest using a macro approach. If you create a macro with a %do %until loop you can repeat your data confidentialized & proc sql steps until your bias criteria is met. Something like the following:

 

%macro bias;

%let ok_bias=0;

%do %until (&ok_bias=1);

 

data confidentialized; set rollup(Rename=(Members=tmpMembers Spend=tmpSpend)) ;
If tmpMembers LE 3 then do ;

...

run;

 

*** BIAS Check ;
Proc SQL ;
Create table

...

quit;

 

data _null_;

set BiasCheck;

if memberbias lt .01 and spendbias lt .01 then call symput('ok_bias','1');

run;

 

%end;

%mend;

 

%bias;

View solution in original post


All Replies
Super User
Posts: 10,474

Re: Randomization algorithm

How are you 'randomizing'?

 

It almost sounds as if you might want to look into Proc Surveyselect and the amount spent might be considered a size measure.

 

 

Frequent Contributor
Posts: 81

Re: Randomization algorithm

Sorry 'Randomizing' may not be the right title for this thread.

 

I am randomizing using Rand("uniform")

 

Data Confidentialized ; Set Rollup(Rename=(Members=tmpMembers Spend=tmpSpend)) ;
If tmpMembers LE 3 then do ;
    Random=rand("Uniform");
    Members=0 ;
    IF Random GT 0.5 then Members=3 ;

 

    Dollars=(Members/tmpMembers)*tmpSpend;
end ;

IF tmpMembers GT 3 then do ; Members =tmpMembers ; Spend=tmpSpend ; end ;
 run ;

*** BIAS Check ;
Proc SQL ;
Create table BiasCheck as SElect
    sum(Members)/sum(tmpMembers)-1 as MemberBias format Percent8.2,
    Sum(Spend)/Sum(tmpSpend)-1 as SpendBias format Percent8.2
From Confidentialized ;
quit;

proc print data=biascheck ; run ;

 

 

Essentially, based on the results of the Bias check, i'd like to re-iterate the confidentialization code so I can 'promote' the version where both bias is less than 1% difference.



Solution
‎09-28-2015 09:16 AM
SAS Employee
Posts: 24

Re: Randomization algorithm

I suggest using a macro approach. If you create a macro with a %do %until loop you can repeat your data confidentialized & proc sql steps until your bias criteria is met. Something like the following:

 

%macro bias;

%let ok_bias=0;

%do %until (&ok_bias=1);

 

data confidentialized; set rollup(Rename=(Members=tmpMembers Spend=tmpSpend)) ;
If tmpMembers LE 3 then do ;

...

run;

 

*** BIAS Check ;
Proc SQL ;
Create table

...

quit;

 

data _null_;

set BiasCheck;

if memberbias lt .01 and spendbias lt .01 then call symput('ok_bias','1');

run;

 

%end;

%mend;

 

%bias;

Frequent Contributor
Posts: 81

Re: Randomization algorithm

Thank you for your suggestion, it works great.

How would I alter the code to have a counter for the number of iterations it took to finally be promoted?

 

 

Super User
Posts: 10,474

Re: Randomization algorithm

Before the %do intialize your counter and then increment within the loop

%let counter=0;

%do .....;

%let counter = %eval(&counter + 1);

 

 

IF you want to see the value of counter outside the loop you will need to make it GLOBAL variable though.

Frequent Contributor
Posts: 81

Re: Randomization algorithm

Thank you ballardw.

I'm trying to alter the code so the loop will break after 100 iterations OR bias is below the threshold and somehow print out which one was it that kicked it out.

Thanks again for your help.
SAS Employee
Posts: 24

Re: Randomization algorithm

Give this a try.

 

%macro bias;

%let counter=1;

%let ok_bias=0;

%do %until (&ok_bias=1 or counter gt 100);

 

data confidentialized; set rollup(Rename=(Members=tmpMembers Spend=tmpSpend)) ;
If tmpMembers LE 3 then do ;

...

run;

 

*** BIAS Check ;
Proc SQL ;
Create table

...

quit;

 

data _null_;

set BiasCheck;

if memberbias lt .01 and spendbias lt .01 then call symput('ok_bias','1');

run;

 

%let counter=%eval(&counter + 1);

 

%end;

%put &ok_bias was the ok_bias value;

%put &counter iterations occurred.;

%mend;

 

%bias;

☑ This topic is SOLVED.

Need further help from the community? Please ask a new question.

Discussion stats
  • 7 replies
  • 256 views
  • 0 likes
  • 3 in conversation