一个关于选取性别、年龄匹配的SAS程序问题

是 否 +2 论坛币

k人 参与回答

经管之家送您一份

应届毕业生专属福利!

求职就业群

赵安豆老师微信：zhaoandou666

经管之家联合CDA

送您一个全额奖学金名额~ !

立即领取

感谢您参与论坛问题回答

经管之家送您两个论坛币！

+2 论坛币

我想从数据库中选择性别相同，年龄相差不超过2岁的病例与对照，我参考了论坛帖子https://bbs.pinggu.org/thread-2845169-1-1.html，自己的程序如下，但总是失败，却不知道问题出在哪里，恳请各位大侠指点迷津，不胜感激！
data match;
input id sex age case;
cards;
1 1 34 1
2 1 35 0
3 2 65 1
4 2 66 0
5 1 11 1
run;
data a1;
set match;
pp=sex; *性别匹配变量pp;
run;
data a_case a_control;set a1;
if case=1 then output a_case; *病例数据集a_case;
if case=0 then output a_control; *对照数据集a_control;
run;

%macro ss;
proc datasets lib=work;delete sample;run;*清除前次抽选样本，方便重复运行宏抽样程序;
proc sql noprint;
select count(*) into: num from a_case; *统计病例个数;
%do i=1 %to &num
%let k=%eval(&i-1);
proc sql noprint;
select id into: idx separated by ' ' from a_case;  *病例ID号 到宏idx;
%let xx=%scan(&idx,&i,' ') ;
data case&i;set a_case;
if id=&xx then call symput("pp_v",pp); *相同匹配条件pp;
if id=&xx then call symput("age_v",age); run; *范围匹配条件age;
data control&i;set a_control;
if pp="&pp_v" and abs(age-&age_v)<3; run; *生成符合匹配条件control样本;
proc surveyselect noprint data=control&i method=srs n=1 out=sample&i seed=1000;        *n=1 ,1：1匹配;
data a_control;set a_control sample&i; run;
proc sql noprint;
create table a_control as  select * from a_control  group by id having n(id)=1;        *剔除已选择control样本;
%end;
data sample;set sample:;run; *生成匹配样本库;
data hb;set a_case sample;run; *生成总库;
proc datasets lib=work;save a1 js a_case a_control sample hb; run; *清除过程文件;
%mend;
%ss;

扫码加我 拉你入群

请注明：姓名-公司-职位

以便审核进群资格，未注明则拒绝

分享0 收藏2 回帖

关键词：sas程序 proc sql separate datasets control 程序

少了个分号吧

data match;
input id sex age case;
cards;
1 1 34 1
2 1 35 0
3 2 65 1
4 2 66 0
5 1 11 1
run;
data a1;
set match;
pp=sex; *性别匹配变量pp;
run;
data a_case a_control;set a1;
if case=1 then output a_case; *病例数据集a_case;
if case=0 then output a_control; *对照数据集a_control;
run;

%macro ss;
proc datasets lib=work;delete sample;run;*清除前次抽选样本，方便重复运行宏抽样程序;
proc sql noprint;
select count(*) into: num from a_case; *统计病例个数;
%do i=1 %to &num
%let k=%eval(&i-1);
proc sql noprint;
select id into: idx separated by ' ' from a_case;  *病例ID号 到宏idx;
%let xx=%scan(&idx,&i,' ') ;
data case&i;set a_case;
if id=&xx then call symput("pp_v",pp); *相同匹配条件pp;
if id=&xx then call symput("age_v",age); run; *范围匹配条件age;
data control&i;set a_control;
if pp="&pp_v" and abs(age-&age_v)<3; run; *生成符合匹配条件control样本;
proc surveyselect noprint data=control&i method=srs n=1 out=sample&i seed=1000;        *n=1 ,1：1匹配;
data a_control;set a_control sample&i; run;
proc sql noprint;
create table a_control as  select * from a_control  group by id having n(id)=1;        *剔除已选择control样本;
%end;
data sample;set sample:;run; *生成匹配样本库;
data hb;set a_case sample;run; *生成总库;
proc datasets lib=work;save a1 js a_case a_control sample hb; run; *清除过程文件;
%mend;
%ss;

谢谢您的解答，可是我加了分号后，程序虽然可以运行，但生成的总库hb并没有达到匹配的效果，里面还是5条记录，这是怎么回事啊，是不是程序有错啊，恳请大侠帮着改一改啊，从而能够达到匹配的效果，不胜感激!

我也参考了这个帖子的程序，结果匹配出来的总库也不对，不知道楼主有没有解决这个问题？请教

病例对照的sas貌似很棘手耶

1. /* try this */
   
2. 
   
3. data a;
   
4. input id sex age case;
   
5. cards;
   
6. 1 1 34 1
   
7. 2 1 35 0
   
8. 3 2 65 1
   
9. 4 2 66 0
   
10. 5 1 11 1
    
11. ;
    
12. run;
    
13. 
    
14. proc sort data= a; by sex; run;
    
15. 
    
16. data matched;
    
17.     merge a(where = (case=1))  
    
18.           a(rename = (id=id1 age=age1  case=case1) where=(case1=0));
    
19.         by sex;
    
20.         if abs(age-age1) < 3;
    
21. run;

复制代码

请问这个问题有没有得到解决呢？
请指教

这个问题似乎没有最终最标准的抽样方法