病例对照、匹配（配对）抽样 SAS 程序（原创）-经管之家官网！

匹配抽样+源代码+参考

1、匹配抽样

---

最近需要用到sas来对我的一个数据进行1：1抽样，我在网上找了些资料看，没怎么看懂，于是自己想了一个思路来编写这个抽样程序，现分享给大家。
首先，数据库js如下，含有唯一变量ID，需要匹配的变量 shcool、grade、class和age，case是定义变量，1是病例，0是对照。抽样要求school、grade、class要相同，age相差不超过2岁。
先说我的思路，首先是生成病例库和对照库， 然后依次从病例库中抽选1条，然后按照条件从对照库选择符合的样本组成待抽选库，然后从这个待抽选库中随机抽选1条，最后把抽出的这条从对照库中剔除，然后重复上述过程直到每条病例都匹配完成，最后合并抽出来的样本。
2、源代码

---

代码如下：

1. proc sort data=js out=a1;
   
2. by school grade class case age;
   
3. data a1;set a1;
   
4. pp=compress(school||grade||class); *班级匹配变量pp;
   
5. run;
   
6. data a_case a_control;set a1;
   
7. if case=1 then output a_case; *病例数据集a_case;
   
8. if case=0 then output a_control; *对照数据集a_control;
   
9. run;
   
10. 
    
11. %macro ss;
    
12. proc datasets lib=work;delete sample;run;*清除前次抽选样本，方便重复运行宏抽样程序;
    
13. proc sql noprint;
    
14. select count(*) into: num from a_case; *统计病例个数;
    
15. %do i=1 %to #
    
16. %let k=%eval(&i-1);
    
17. proc sql noprint;
    
18. select id into: idx separated by ' ' from a_case;*病例ID号 到宏idx;
    
19. %let xx=%scan(&idx,&i,' ') ;
    
20. data case&i;set a_case;
    
21. if id=&xx then call symput("pp_v",pp); *相同匹配条件pp;
    
22. if id=&xx then call symput("age_v",age); run; *范围匹配条件age;
    
23. data control&i;set a_control;
    
24. if pp="&pp_v" and abs(age-&age_v)<3; run; *生成符合匹配条件control样本;
    
25. proc surveyselect noprint data=control&i method=srs n=1 out=sample&i seed=1000; *n=1 ,1：1匹配;
    
26. data a_control;set a_control sample&i; run;
    
27. proc sql noprint;
    
28. create table a_control asselect * from a_controlgroup by id having n(id)=1; *剔除已选择control样本;
    
29. %end;
    
30. data sample;set sample:;run; *生成匹配样本库;
    
31. data hb;set a_case sample;run; *生成总库;
    
32. proc datasets lib=work;save a1 js a_case a_control sample hb; run; *清除过程文件;
    
33. %mend;
    
34. %ss;

复制代码
最后，说说程序缺陷，首先是运行效率低，因为是循环语句反复生成数据库，病例较多时，运行肯定慢； 其次未考虑到匹配未成功的情况；第三，程序还是有点啰嗦。 仅供参考。
3、参考

---

[程序分享] 病例对照、匹配（配对）抽样 SAS 程序（原创）
本文来自: 人大经济论坛 SAS专版 版，详细出处参考： https://bbs.pinggu.org/forum.php?mod=viewthread&tid=2845169&page=1&fromuid=3495972