data a;
do i=1 to 29;
x=ranuni(29);
output;
end;
format x f8.4;
run;
proc surveyselect data=a out=b method=sys sampsize=4 seed=10;
run;
proc means data=c mean std ;
var x;
output out=d (drop=_type_) mean=avg std=std;
run;

data a;
do i=1 to 29;
x=ranuni(29);
output;
end;
format x f8.4;
run;
proc surveyselect data=a out=b method=sys sampsize=4 seed=10;
run;
proc means data=c mean std ;
var x;
output out=d (drop=_type_) mean=avg std=std;
run;

data crackman;
do i=1 to 29;
x=ranuni(29);
output;
end;
format x f8.4;
run;
%crackman;
%macro crackman;
%do m=1 %to 30;
%a(k=&&m.);
%end;
%mend crackman;
%macro crackman1(k=);
%do n=1 %to 29;
proc surveyselect data=a out=b&n. method=sys sampsize=&n. seed=%eval(&k.*&n.);
run;
proc means data=b&n. mean std ;
var x;
output out=d&n. (drop=_type_) mean=avg std=std;
run;
%end;
%do n=1 %to 29;
proc append base=d1 data=d&n. force;
run;
%end;
data m&k.;
set d1;
run;
%mend crackman1;
结果里面的M1-M30为结果
_freq_为每一次抽样的样本个数 AVG为均值 STD为标准差
楼主看哈
我这个程序电脑上跑了几分钟，效率还是有点低 希望高手来补充，优化一下

data crackman;
do i=1 to 29;
x=ranuni(29);
output;
end;
format x f8.4;
run;
%crackman;
%macro crackman;
%do m=1 %to 30;
%crackman1(k=&&m.);
%end;
%mend crackman;
%macro crackman1(k=);
%do n=1 %to 29;
proc surveyselect data=a out=b&n. method=sys sampsize=&n. seed=%eval(&k.*&n.);
run;
proc means data=b&n. mean std ;
var x;
output out=d&n. (drop=_type_) mean=avg std=std;
run;
%end;
%do n=1 %to 29;
proc append base=d1 data=d&n. force;
run;
%end;
data m&k.;
set d1;
run;
%mend crackman1;
结果里面的M1-M30为结果
_freq_为每一次抽样的样本个数 AVG为均值 STD为标准差
楼主看哈
我这个程序电脑上跑了几分钟，效率还是有点低 希望高手来补充，优化一下

5# crackman

So good。