I think best option is what you need.
[pre]
data reg;
array v y x1-x14;
do i = 1 to 100;
do j = 1 to dim(v);
v =ranuni(1234);
end;
output;
end;
run;
proc reg outest=est;
model Y = X1-X4 /* 4 confounders */
X5-x14
/ include=4 best=1024 selection=adjrsq aic sbc;
run;
quit;
proc print n;
run;
[/pre]
... View more