[SAS代码模板]抽样_surveyselect

SAS抽样代码模板

黄色部分为套用部分，红色部分为可选部分
——————————模板——————————
proc surveyselect data=总体数据 out=样本数据 method=抽样方法 n=抽取样本;
strata 分层变量;
run;
————————————————————
method指定抽样方法：
l srs：简单无重复随机抽样，可以用n=指定需要抽取的样本数，也可以用samprate=指定要抽取的样本占总体的比例。
l sys：系统抽样，需要指定样本（用语句sampsize=），从总体N中，每隔N/sampsize个抽取一个。
l 当指定srs方法，并使用strata语句时，则为分层抽样（需要先将总体按照分类变量排序sort），此时n（或者samprate）可以指定分别每一层的样本数（或比例），来实现不等比例抽样。比如，

proc surveyselect data=test3 out=results3 method=srs n=(30,20,50,40);
strata 分层变量;
run;

如果碰到这样的问题肿么办？？
总体可以根据变量A分成a1和a2两部分，总共打算抽取N个数据（比如50），想从a1中抽取n1个（比如30，不足n1=30个，则抽取全部，即有27个就抽27个），a2中抽取剩下的N-n1个，肿么抽？？？

给出如下抽样代码，建议有选择的套用（黄色部分是套用部分）：
假设总体数据名为x，有变量A（A=a1，a2），现在需要抽取50个样本：抽取A=a1的样本30个（不足30个抽取全部），剩下的抽取A=a2的数据，直到样本总数50个。

*将该中心的数据分成两部分：a1和a2，假设总体数据集为a;
data a1 a2;
set x;
if A=a1 then output a1;
if A=a2 then output a2;
run;

*对a1抽样：不足30抽取全部，大于30抽取30;
%let disd=%sysfunc(open(a1));
%let nobs=%sysfunc(attrn(&disd,nobs));
%let rc=%sysfunc(close(&dsid));
%if &nobs<=30 %then %do;
data samplea1;
set a1;
run;
%end;
%else %do;
proc surveyselect data=a1 out=samplea1 method=srs n=30;
run;
%end;
*对a2进行抽样:a1部分不足30的，a2补齐剩下的30-&nobs，a1抽取30的，part2抽取20个;
%let nobs2=%eval(50-&nobs);
%if &nobs<=30 %then %do;
proc surveyselect data=a2 out=samplea2 method=srs n=&nobs2;
run;
%end;
%else %do;
proc surveyselect data=a2 out=samplea2 method=srs n=20;
run;
%end;
*合并两表;
data sample;
set a1 a2;
run;
proc export data=sapmle outfile="D:sample.xls" dbms = excel2003;
run;

原文始发于微信公众号（PPV课数据科学社区）：[SAS代码模板]抽样_surveyselect