proc means data=Clean.Patients noprint ; var HR; where HR ne (900); output out=Mean_Std(drop=_type_ _freq_) mean= std= / autoname; run; title "Outliers for HR Based on 2 Standard Deviations"; data _null_; file print; set Clean.Patients(keep=Patno HR ); ***bring in the means and standard deviations; if _n_ = 1 then set Mean_Std; **_n_ is the index**; if HR lt HR_Mean - 2*HR_StdDev and not missing(HR) or HR gt HR_Mean + 2*HR_StdDev then put Patno= HR=; run