/** Approximate Score (Agresti-Coull) Confidence Interval for a Proportion **/ data ldl; input decrease_ldl @@; decrease=(decrease_ldl>0); cards; -7.6 13.9 27.1 14.8 6.3 42.0 41.7 18.5 31.2 24.0 ; run; proc sort data=ldl; by descending decrease; run; proc print;run; proc freq data=ldl order=data; tables decrease/binomial(ac) alpha=.05; ods select binomialCLs; run;