Histogram with Superimposed Beta Curve

 /****************************************************************/
 /*          S A S   S A M P L E   L I B R A R Y                 */
 /*                                                              */
 /*    NAME: CAPBETA                                             */
 /*   TITLE: Histogram with Superimposed Beta Curve              */
 /* PRODUCT: QC                                                  */
 /*  SYSTEM: ALL                                                 */
 /*    KEYS: Capability Analysis, Beta Distribution,             */
 /*   PROCS: CAPABILITY                                          */
 /*    DATA:                                                     */
 /*                                                              */
 /*     REF: SAS/QC Software:  Examples                          */
 /*    MISC:                                                     */
 /*                                                              */
 /****************************************************************/

options ps=60 ls=80;


data plates;
   label gap='Plate Gap in cm';
   input gap @@;
   cards;
0.746  0.357  0.376  0.327  0.485
1.741  0.241  0.777  0.768  0.409
0.252  0.512  0.534  1.656  0.742
0.378  0.714  1.121  0.597  0.231
0.541  0.805  0.682  0.418  0.506
0.501  0.247  0.922  0.880  0.344
0.519  1.302  0.275  0.601  0.388
0.450  0.845  0.319  0.486  0.529
1.547  0.690  0.676  0.314  0.736
0.643  0.483  0.352  0.636  1.080
;

title1 'Fitted Beta Distribution of Plate Gaps';
title2 '50 Measurements';

proc capability data=plates noprint;
   var gap;
   specs lsl=0.3 clsl=white llsl=1  wlsl=2
         usl=0.8 cusl=white lusl=20 wusl=2;
   histogram /
      beta(theta=0 scale=1.9 color=yellow w=2)
      midpoints = 0.2 to 1.8 by 0.2
      legend    = legend1
      cframe    = gray
      cfill     = green;
   legend1 cframe=gray cborder=black;
run;

goptions reset=all;