sas如何进行分层抽样

已有一个总体数据集rawdata，现想从该总体里分层抽取一个样本。按地区、年龄段、性别分层抽取，每层抽取的数量为指定数字，如下图所示。用sas如何实现？

是有点复杂，没想到什么具体的过程步可以直接实现。理论上可以通过数据步给每个记录随机数，然后根据年龄、性别、分组按需要的数量取随机数最大/小的相应个记录。

孤单的我们 发表于 2020-10-10 09:08
是有点复杂，没想到什么具体的过程步可以直接实现。理论上可以通过数据步给每个记录随机数，然后根据年龄、 ...

最先想到的也是这杨操作，但是这个分组表只展示了一部分用作举例，完整的分组有七百多个，一个个设置基本上没有操作性。

hi11 发表于 2020-10-12 22:48
最先想到的也是这杨操作，但是这个分组表只展示了一部分用作举例，完整的分组有七百多个，一个个设置基本 ...

不论多少个，只要有完整的列表，数据集读取后，代码自动就做完了。700个和7个代码都一样

1. /*generate the test dataset*/
   
2. data raw;
   
3.         do age=1 to 5;
   
4.                 do sex=1 to 2;
   
5.                         do grp=1 to 9;
   
6.                                 do i=1 to 100;
   
7.                                         x=rannor(201013);
   
8.                                         output;
   
9.                                 end;
   
10.                         end;
    
11.                 end;
    
12.         end;
    
13. run;
    
14. 
    
15. /*generate the sample rule set*/
    
16. data rule;
    
17.         do age=1 to 5;
    
18.                 do sex=1 to 2;
    
19.                         do grp=1 to 9;
    
20. /*                                n is the sample size for each strata*/
    
21.                                         n=ceil(uniform(201013)*10);
    
22.                                         output;
    
23.                         end;
    
24.                 end;
    
25.         end;
    
26. run;
    
27. 
    
28. /*add the random sequence for the raw data*/
    
29. data raw1;
    
30.         set raw;
    
31.         randseq=rannor(201013);
    
32. /*        add the record ID*/
    
33.         id=_n_;
    
34. run;
    
35. 
    
36. proc sql noprint;
    
37.         create table raw2 as
    
38.         select a.*,b.n from raw1 as a
    
39.         left join rule as b
    
40.         on a.age eq b.age and a.sex eq b.sex and a.grp eq b.grp
    
41.         order by age,sex,grp,randseq;
    
42. quit;
    
43. 
    
44. data wanted;
    
45.         set raw2;
    
46.         by age sex grp;
    
47.         retain k 0;
    
48.         if first.grp then k=0;
    
49.         k+1;
    
50.         if k<=n ;
    
51. run;
    
52.        

复制代码