Three-Parameter Weibull Distribution

 /****************************************************************/
 /*          S A S   S A M P L E   L I B R A R Y                 */
 /*                                                              */
 /*    NAME: CAPW3A                                              */
 /*   TITLE: Three-Parameter Weibull Distribution                */
 /* PRODUCT: QC                                                  */
 /*  SYSTEM: ALL                                                 */
 /*    KEYS: Capability Analysis, Weibull Distribution,          */
 /*   PROCS: CAPABILITY IML                                      */
 /*    DATA:                                                     */
 /*                                                              */
 /*     REF: SAS/QC Software:  Usage and Reference, Version 6,   */
 /*          First Edition, Volume 1 and Volume 2                */
 /*    MISC:                                                     */
 /*                                                              */
 /*  Examples in the documentation were created using the        */
 /*    statement:  ods listing style=statistical;                */
 /*                                                              */
 /****************************************************************/

 options ps=60 ls=80 nodate;


 data Plastic;
    label Strength='Strength in psi';
    input Strength @@;
    cards;
 30.26  31.23  71.96  47.39  33.93  76.15  42.21
 81.37  78.48  72.65  61.63  34.90  24.83  68.93
 43.27  41.76  57.24  23.80  34.03  33.38  21.87
 31.29  32.48  51.54  44.06  42.66  47.98  33.73
 25.80  29.95  60.89  55.33  39.44  34.50  73.51
 43.41  54.67  99.43  50.76  48.81  31.86  33.88
 35.57  60.41  54.92  35.66  59.30  41.96  45.32
 ;

   goptions htext   = 3.0 pct
            htitle  = 3.5 pct;
   title 'Three-Parameter Weibull Fit';
   proc capability data=Plastic noprint;
      spec lsl=25 cleft=orange clsl=black;
      histogram Strength  / weibull(fill color = paoy
                                         theta = est)
                           cfill  = paoy
                           cframe = ligr
                           nolegend;
      inset lsl='LSL' lslpct / cfill = ywh pos=nw;
      inset weibull          / format=6.2 pos=ne cfill = ywh;
   run;



 goptions reset=all;