Sugi-89-265 Bretheim.pdf - sasCommunity

More documents

Recommendations

Info

EldUbit 2SAS* System statementsoptions nodate:libname in 'd;\sas':* compute logarithm for each value *:data one:set in. raw:logchg",logt charge.l ;proc sort data~one:by drg;'* compute mean of log values for each group *:proc means data=one mean noprint:var logchg~by drg:output out~stats m€an~mchg;* merge each observation with its group mean *:data twa:merge one stats:by drg;if first.drg then do:sum_sd""O :end:* compute geometric mean *:geo_mean ~ exp(mchg):* accumulate squared deviations ~:squardev ~ (charge - geo_mean)**2:sum_sd + squardev;* ,compute the standard deviation and threshold *.if last.drg an~ _freq_ > 1 then do:std_dev '" sqrt(sum_sd/Lfreq_ - 1)):thres = geo_mean ~ !strl_dev * 3}:output:end:procvarfQrmatlabeltitleltitle2run~print data~twolabel:drg _freq_ geo_mean std_dev thres:geo_mean std_dev thres dollarl0.2:_freq_'~'geo mean ~ 'GEO. MEAN'std~dev 'STD. DEV.'thres ' THRESHOLD' ;'OUTLIER iHRESHOLDS BASED ON'~'THREE STANDARD DEVIATIONS FROM THE GEOMETRIC MEA~;:1431
Exhibit 3Sample OutpntOUTLIER THRESHOLDS BASED ONTHREE STANDARD DEVIATIONS FROM THE GEOMETRIC MEANQ!lSllil.G.NGFA).MEANSTD.DEV.THRESHOLD123105201379433$18,623.466,162.691,356.76$2,729.59812.30157.90$26,812.238,599.591,830.48•1432
Page 1: COMPUTING THE STANDARD DEVIATION FR

Sugi-89-265 Bretheim.pdf - sasCommunity

You also want an ePaper? Increase the reach of your titles

Delete template?

Save as template?