Dear Tom, I called the macro but it did not work. I am sure I did some thing wrong but I dont know where it is. Could you please see the below codes and help me where is the problem? ? /* ********************************************************************************* */ /* ******************** W R D S R E S E A R C H M A C R O S ******************** */ /* ********************************************************************************* */ /* WRDS Macro: WINSORIZE */ /* Summary : Winsorizes or Trims Outliers */ /* Date : April 14, 2009 */ /* Author : Rabih Moussawi, WRDS */ /* Variables : - INSET and OUTSET are input and output datasets */ /* - SORTVAR: sort variable used in ranking */ /* - VARS: variables to trim and winsorize */ /* - PERC1: trimming and winsorization percent, each tail (default=1%) */ /* - TRIM: trimming=1/winsorization=0, default=0 */ /* ********************************************************************************* */ %MACRO WINSORIZE (INSET=,OUTSET=,SORTVAR=,VARS=,PERC1=1,TRIM=0); /* List of all variables */ %let vars = %sysfunc(compbl(&vars)); %let nvars = %nwords(&vars); /* Display Output */ %put ### START.; /* Trimming / Winsorization Options */ %if &trim=0 %then %put ### Winsorization; %else %put ### Trimming; %put ### Number of Variables: &nvars; %put ### List of Variables: &vars; options nonotes; /* Ranking within &sortvar levels */ %put ### Sorting... ; proc sort data=&inset; by &sortvar; run; /* 2-tail winsorization/trimming */ %let perc2 = %eval(100-&perc1); %let var2 = %sysfunc(tranwrd(&vars,%str( ),%str(__ )))__; %let var_p1 = %sysfunc(tranwrd(&vars,%str( ),%str(__&perc1 )))__&perc1 ; %let var_p2 = %sysfunc(tranwrd(&vars,%str( ),%str(__&perc2 )))__&perc2 ; /* Calculate upper and lower percentiles */ proc univariate data=&inset noprint; by &sortvar; var &vars; output out=_perc pctlpts=&perc1 &perc2 pctlpre=&var2; run; %if &trim=1 %then %let condition = %str(if myvars(i)>=perct2(i) or myvars(i)<=perct1(i) then myvars(i)=. ); %else %let condition = %str(myvars(i)=min(perct2(i),max(perct1(i),myvars(i))) ); %if &trim=0 %then %put ### Winsorizing at &perc1.%... ; %else %put ### Trimming at &perc1.%... ; /* Save output with trimmed/winsorized variables */ data &outset; merge &inset (in=a) _perc; by &sortvar; if a; array myvars {&nvars} &vars; array perct1 {&nvars} &var_p1; array perct2 {&nvars} &var_p2; do i = 1 to &nvars; if not missing(myvars(i)) then do; &condition; end; end; drop i &var_p1 &var_p2; run; /* House Cleaning */ proc sql; drop table _perc; quit; options notes; %put ### DONE . ; %put ; %MEND WINSORIZE; /* ********************************************************************************* */ /* ************* Material Copyright Wharton Research Data Services *************** */ /* ****************************** All Rights Reserved ****************************** */ /* ********************************************************************************* */ %MACRO WINSORIZE (INSET=clean.main11 ,OUTSET=clean.winsorized ,SORTVAR=none ,VARS=mdr, bdr, ebit_ta, fa_ta ,PERC1=1,TRIM=0);
... View more
