If a matched character is defined by a number of catagorical variables including binned variables, then SAS proc surveyselect with stratified sampling may help in your case.
楼主按照你的程序,我运行了一下,生成了4个文件a1 control case and hb,是吧case和control纵向合并吗,按照case的数目,把control中多余的删除是不是就是1:1匹配了,因为case和control加起来就是a1的数据,
我是用年龄和性别匹配,改了一下参数,不过老是出现一下错误:
proc sort data=data out=a1;
by sex age malb_yn;
data a1;set a1;
pp=sex; *班级匹配变量pp;
run;
data a_case a_control;set a1;
if malb_yn =1 then output a_case; *病例数据集a_case;
if malb_yn=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 #
%let k=%eval(&i-1);
proc sql noprint;
select caseid into: idx separated by ' ' from a_case; *病例ID号 到宏idx;
%let xx=%scan(&idx,&i,' ') ;
data case&i; set a_case;
if caseid=&xx then call symput("pp_v",pp); *相同匹配条件pp;
if caseid=&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 caseid having n(id)=1; *剔除已选择control样本;
%end;
data sample;set sample:;run; *生成匹配样本库;
data hb;set a_case sample;run; *生成总库;
proc datasets lib=work;save a1 data a_case a_control sample hb; run; *清除过程文件;
%mend;
%ss;