data a;
drop bias;
do i = 1 to 100000;
bias = 0.5;
if ranuni(0) > bias then a1 = 1; else a1 = 0;
if ranuni(0) > bias then a2 = 1; else a2 = 0;
if ranuni(0) > bias then a3 = 1; else a3 = 0;
if ranuni(0) > bias then a4 = 1; else a4 = 0;
if ranuni(0) > bias then a5 = 1; else a5 = 0;
if ranuni(0) > bias then a6 = 1; else a6 = 0;
if ranuni(0) > bias then a7 = 1; else a7 = 0;

value = a1 + 2 * a2 + 4 * a3 + 8 * a4 + 16 * a5 + 32 * a6 + 64 * a7 ;


* a5 = (a1 and not a2) or (not a1 and a2);
* a6 = (a3 and not a4) or (not a3 and a4);
* a7 = (a5 and not a6) or (not a5 and a6);

output;
end;
run;
proc univariate data=a plot;
var value;
run;
data b; set a;
if a1 = 1 then a1 = 0; else a1 = 1;
if a3 = 1 then a3 = 0; else a3 = 1;
if a5 = 1 then a5 = 0; else a5 = 1;
if a7 = 1 then a7 = 0; else a7 = 1;

value = a1 + 2 * a2 + 4 * a3 + 8 * a4 + 16 * a5 + 32 * a6 + 64 * a7 ;
value1 = value * MOD(i,14);
run;
proc univariate data=b plot;
var value;
run;

proc arima data=b;
     identify var = value1;
  run;