发帖
雷达卡
我一定犯了什么很低级的错误是不是....呜呜....
/********************************************************************************/
/* Program: PSMatch_Multi.sas
/*
/* Platform: SAS 9.1.3
/*
/* Drug/Protocol: Generalized SAS Macro
/*
/* Description: Does N:1 optimized propensity score matching within specified
/* absolute differences of propensity score
/********************************************************************************/
%macro psmatch_multi (pat_dsn =prop_score_treated, /* Name of data set with patient data */
pat_idvar=did, /* Name of Patient ID variable in data set &PAT_DSN */
pat_psvar=prob_treat, /* Name of Propensity Score variable in data set &PAT_DSN */
cntl_dsn= prop_score_untreated, /* Name of data set with control data */
cntl_idvar=did, /* Name of Control ID variable in data set &CNTL_DSN */
cntl_psvar=prob_treat, /* Name of Propensity Score variable in data set &CNTL_DSN */
match_dsn=matched_pairs, /* Name of output data set with N:1 matches */
match_ratio=1, /* Number of control matches per patient */
score_diff=0.001, /* Maximum allowable absolute differences between propensity scores*/
) /* Optional input seed for random number generator */
;
/**********************************************************************/
/* Delete final matched pairs dataset, if it exists from a prior run
/**********************************************************************/
proc datasets nolist;
delete __final_matched_pairs;
run;
quit;
/********************************************/
/* Make all internal macro variables local
/********************************************/
%local __dsid __varnum __cntl_type __rc __num;
/***************************************************************************/
/* Determine characterisitcs of Control ID variable (numeric or character)
/***************************************************************************/
%let __dsid = %sysfunc(open(&cntl_dsn,i));
%let __varnum = %sysfunc(varnum(&__dsid, &cntl_idvar));
%let __cntl_type = %sysfunc(vartype(&__dsid, &__varnum));
%let __rc = %sysfunc(close(&__dsid));
%put &__cntl_type;
/**************************/
/* New director Matching Data
/**************************/
data __patmatch (keep = &pat_idvar &pat_psvar);
set &pat_dsn;
run;
/**************************/
/* Control Matching Data
/**************************/
data __contmatch (keep = &cntl_idvar &cntl_psvar);
set &cntl_dsn;
run;
/************************************************************/
/* Find all possible matches between directors and controls
/* Propensity scores must match within +/- &match
/************************************************************/
proc sql;
create table __matches0 as
select
p.&pat_idvar as pat_idvar,
c.&cntl_idvar as cntl_idvar,
p.&pat_psvar as pat_score,
c.&cntl_psvar as cntl_score
from __patmatch p left join __contmatch c
on abs(p.&pat_psvar - c.&cntl_psvar) <= &score_diff
order by pat_idvar;
quit;
/*************************************/
/* Data set of all possible matches
/*************************************/
data __possible_matches;
set __matches0;
/*-----------------------------------------*/
/* Create a random number for each match
/*-----------------------------------------*/
rand_num = ranuni(&seed);
/*-----------------------------------------------*/
/* Remove patients who had no possible matches
/*-----------------------------------------------*/
%if &__cntl_type = C %then %do;
if cntl_idvar ^= '';
%end;
%else %if &__cntl_type = N %then %do;
if cntl_idvar ^= .;
%end;
/*---------------------------*/
/* Create a dummy variable
/*---------------------------*/
n = 1;
run;
/******************************************************************/
/* Find the number of potential control matches for each patient
/******************************************************************/
proc freq data=__possible_matches noprint;
tables pat_idvar / out=__matchfreq (keep = pat_idvar count);
run;
/****************************************************************************/
/* Optimize control matching for directors based on number of possible matches
/* Pick matches for directors with the fewest number of possible matches first
/****************************************************************************/
data __matches_freq0;
merge __possible_matches
__matchfreq;
by pat_idvar;
run;
/********************************************************************/
/* Find the number of potiential directors matches for each control
/********************************************************************/
proc freq data=__possible_matches noprint;
tables cntl_idvar / out=__cntlfreq (keep = cntl_idvar count rename = (count = cntl_count));
run;
proc sort data=__matches_freq0;
by cntl_idvar;
run;
data __matches_freq;
merge __matches_freq0
__cntlfreq;
by cntl_idvar;
/*------------------------------------------------------*/
/* Take out directors with less than number of matches
/*------------------------------------------------------*/
if count >= &match_ratio;
run;
proc datasets nolist;
delete __matches0;
run;
quit;
/*****************************************************************/
/* Count the number of entries in the file of possible matches
/*****************************************************************/
%let __dsid = %sysfunc(open(__matches_freq,i));
%let __num = %sysfunc(attrn(&__dsid,nobs));
%let __rc = %sysfunc(close(&__dsid));
%do %while (&__num >= 1);
proc sort data=__matches_freq;
by count cntl_count rand_num pat_idvar;
run;
/*********************************************************/
/* Get first randomly selected patient with the minimum
/* number of matches
/*********************************************************/
data __first_pat_idvar (keep = pat_idvar);
set __matches_freq;
by n;
if first.n;
run;
/******************************************************************/
/* Get all matches for that patient
/* Select the first randomly selected for the number of matches
/******************************************************************/
proc sort data=__matches_freq;
by pat_idvar count cntl_count rand_num;
run;
data __all_first_id;
merge __matches_freq
__first_pat_idvar (in=i);
by pat_idvar;
if i;
num + 1;
run;
data __new_matched_pairs (keep = pat_idvar cntl_idvar pat_score cntl_score);
set __all_first_id;
label pat_idvar = "Director ID, original variable name &dir_idvar"
cntl_idvar = "Matched Control ID, original variable name &cntl_idvar"
pat_score = "Director Propensity Score, original var name &pat_psvar"
cntl_score = "Matched Control Propensity Score, orig var &cntl_psvar"
;
if num <= &match_ratio;
run;
/******************************************/
/* Remove directors with matched controls
/******************************************/
proc sort data=__new_matched_pairs (keep = pat_idvar)
out=__new_matched_pats nodupkey;
by pat_idvar;
run;
data __match_remove_pat;
merge __possible_matches
__new_matched_pats (in=id);
by pat_idvar;
if ^id;
run;
/************************************************************/
/* Remove all matched pairs that include selected controls
/************************************************************/
proc sort data=__new_matched_pairs (keep = cntl_idvar) out=__remove_cont;
by cntl_idvar;
run;
proc sort data=__match_remove_pat;
by cntl_idvar;
run;
data __match_remove_cont;
merge __match_remove_pat
__remove_cont (in=id);
by cntl_idvar;
if ^id;
run;
proc sort data=__match_remove_cont out=__possible_matches;
by pat_idvar;
run;
/********************************************************/
/* Add new matched pairs to set of final matched pairs
/********************************************************/
proc append base=__final_matched_pairs data=__new_matched_pairs;
run;
/******************************************************************/
/* Find the number of potential control matches for each patient
/******************************************************************/
proc freq data=__possible_matches noprint;
tables pat_idvar / out=__matchfreq (keep = pat_idvar count);
run;
/***************************************************************************/
/* Optimize control matching for patients based on number of possible matches
/* Pick matches for patients with the fewest number of possible matches first
/***************************************************************************/
data __matches_freq0;
merge __possible_matches
__matchfreq;
by pat_idvar;
run;
/********************************************************************/
/* Find the number of potential patient matches for each control
/********************************************************************/
proc freq data=__possible_matches noprint;
tables cntl_idvar / out=__cntlfreq (keep = cntl_idvar count rename = (count = cntl_count));
run;
proc sort data=__matches_freq0;
by cntl_idvar;
run;
data __matches_freq;
merge __matches_freq0
__cntlfreq;
by cntl_idvar;
/*------------------------------------------------------*/
/* Take out directors with less than number of matches
/*------------------------------------------------------*/
if count >= &match_ratio;
run;
/********************************************************/
/* Determine number of remaining possible matched pairs
/********************************************************/
%let __dsid = %sysfunc(open(__matches_freq,i));
%let __num = %sysfunc(attrn(&__dsid,nobs));
%let __rc = %sysfunc(close(&__dsid));
%end; /* of " %do %while (&__num >= 1); */
/********************************************************************************/
/* Create final output data set with one observation for each original patient
/* ID Variable names in output data set are PAT_IDVAR, PAT_SCORE, CNTL_IDVAR,
/* CNTL_SCORE
/* If no match for patient ID (PAT_IDVAR), then corresponding CNTL variables
/* (CNTL_IDVAR, CNTL_SCORE) are missing.
/********************************************************************************/
proc sort data=__final_matched_pairs;
by pat_idvar pat_score;
run;
data __patmatch_orig;
set __patmatch (rename= (&pat_idvar = pat_idvar &pat_psvar = pat_score));
run;
proc sort data=__patmatch_orig;
by pat_idvar;
run;
data &match_dsn (label = "Final Matched Pairs for Propensity Score Matching");
merge __final_matched_pairs
__patmatch_orig;
by pat_idvar pat_score;
run;
/***************************************************/
/* Delete all temporary datasets created by macro
/***************************************************/
proc datasets nolist;
delete __contmatch __final_matched_pairs __matches_freq0 __matches_freq
__match_pair0 __matchfreq __match_remove_cont __match_remove_pat
__new_matched_pairs __patmatch __patmatch_orig __possible_matches
__remove_cont __cntlfreq __first_pat_idvar __all_first_id
__new_matched_pats;
run;
quit;
%mend psmatch_multi;
提升卡
置顶卡
沉默卡
变色卡
抢沙发
千斤顶
显身卡
京公网安备 11010802022788号
