Looks to me like your first 4 steps are just these two steps. Plus I generated SIGN variables that can probably be used to reduce the PROC RANK steps to one instead of two.
data final ;
input Eink :$2. Bg He ;
cards;
01 55.487739 2.722371
02 127.485645 18.670678
03 182.95309 53.529988
04 354.69612 122.703093
05 450.777081 137.596942
06 461.617933 157.457164
07 429.38197 192.283645
08 370.644235 173.688658
09 281.256479 180.848677
10 170.426973 137.045263
11 118.22022 128.45571
12 93.656461 104.596442
13 68.561622 112.710123
14 42.443895 77.742469
15 50.933162 136.718895
16 35.285987 104.219131
17 23.123863 78.424041
18 16.724659 89.48717
19 9.271308 46.396113
20 4.870804 21.574738
21 4.83902 19.753305
22 3.487616 9.23105
23 0.860551 2.136266
24 1.202043 1.861731
;
data min_max ;
length Eink $2 Xmin Xmax 8;
do xmax = 250 to 3500 by 250,4000,4500,5000,6000,7000,8000,10000,15000,25000,. ;
i+1;
eink=put(i,z2.);
xmin=sum(0,lag(xmax));
output;
end;
drop i;
run;
proc summary data=FINAL ;
var Bg He ;
output out=sums(drop=_freq_ _type_) sum=sum_Bg sum_He ;
run;
data STEP4;
merge final min_max ;
by eink;
if _n_=1 then set sums;
cum_Bg + Bg;
Fx_Bg = sum_bg/2 - cum_Bg;
Sign_Bg = sign(fx_bg);
cum_He + He;
Fx_He = sum_he/2 - cum_He;
Sign_He = sign(fx_he);
run;
So if Bg and He represent the variables you need to process then you might use macro code to generate those blocks that use them.
...
sum=sum_Bg sum_He ;
...
cum_Bg + Bg;
Fx_Bg = sum_bg/2 - cum_Bg;
Sign_Bg = sign(fx_bg);
cum_He + He;
Fx_He = sum_he/2 - cum_He;
Sign_He = sign(fx_he);
...
But you could also transpose the data first. Then no need to generate code (other than the list of variables in the VAR statement of the TRANSPOSE step).
proc transpose data=FINAL out=TALL(rename=(col1=VALUE)) name=VARNAME ;
by Eink;
var Bg He ;
run;
proc sort data=TALL;
by varname eink;
run;
proc summary data=TALL nway;
by varname ;
var value;
output out=sums(drop=_freq_ _type_) sum=sum_value ;
run;
data step4_tall;
merge tall sums ;
by varname ;
if _n_=1 then do;
set min_max;
declare hash h(dataset:'min_max');
h.definekey('eink');
h.definedata('Xmin','Xmax');
h.definedone();
end;
if first.varnum then cum_value=.;
h.find();
cum_value + value ;
fx_value = sum_value/2 - cum_value ;
Sign_value = sign(fx_value);
run;
