[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;