I want to remove outliers using median +/- 1.5 IQR (Qrange in SAS).
The proc univariate can generate median and Qrange, but how do I use these values in another proc or data step?
Another way is to use proc sql, but it seems proc sql summary function does not have qrange
or proc boxplot can also generate median and qrange, but I still need to use these values in another step.
You could save the median and iqr as macro variables, apply them in a Data step:
proc univariate data = sashelp.iris;
var petallength;
output out=boxStats median=median qrange = iqr;
run;
data _null_;
set boxStats;
call symput ('median',median);
call symput ('iqr', iqr);
run;
%put &median;
%put &iqr;
data trimmed;
set sashelp.iris;
if (petallength le &median + 1.5 * &iqr) and (petallength ge &median - 1.5 * &iqr);
run;
proc print data = trimmed;
run;
Use proc means, capture the data in a dataset and merge that into the main data set.
If you have SAS 9.4 proc SQL supports the median function otherwise it doesn't.
You could save the median and iqr as macro variables, apply them in a Data step:
proc univariate data = sashelp.iris;
var petallength;
output out=boxStats median=median qrange = iqr;
run;
data _null_;
set boxStats;
call symput ('median',median);
call symput ('iqr', iqr);
run;
%put &median;
%put &iqr;
data trimmed;
set sashelp.iris;
if (petallength le &median + 1.5 * &iqr) and (petallength ge &median - 1.5 * &iqr);
run;
proc print data = trimmed;
run;
Hi,
Please try this:
proc univariate data=sashelp.class;
var age;
output out=stats median=median Qrange=Qrange;
run;
data want(drop=median Qrange);
if _n_=1 then set stats;
set sashelp.class;
if age<median - 1.5*Qrange or age>median + 1.5*Qrange then delete;
run;
Available on demand!
Missed SAS Innovate Las Vegas? Watch all the action for free! View the keynotes, general sessions and 22 breakouts on demand.
Learn how use the CAT functions in SAS to join values from multiple variables into a single value.
Find more tutorials on the SAS Users YouTube channel.