Two-Parameter Weibull Q-Q Plot

 /****************************************************************/
 /*          S A S   S A M P L E   L I B R A R Y                 */
 /*                                                              */
 /*    NAME: CAPWEIB1                                            */
 /*   TITLE: Two-Parameter Weibull Q-Q Plot                      */
 /* PRODUCT: QC                                                  */
 /*  SYSTEM: ALL                                                 */
 /*    KEYS: Capability Analysis, Quantile-Quantile Plots,       */
 /*   PROCS: CAPABILITY                                          */
 /*    DATA:                                                     */
 /*                                                              */
 /*    MISC:                                                     */
 /*                                                              */
 /****************************************************************/
  options ps=60 ls=80;

  data failures;
     input time @@;
     cards;
  27.96  45.78  50.90  23.63  18.22
  44.91   9.44  45.20  43.80  46.81
  33.66  18.07  16.08  29.40  32.40
  52.92  29.35  18.05  16.75  21.37
  23.83  37.33  27.36  48.92  39.56
  28.02  19.27   6.12  65.80  21.07
  39.77   1.01  50.95  71.46  27.67
  20.24  30.56  14.98  45.79   5.95
  46.85  48.30  36.74  16.06  25.49
  40.16  40.59  51.31  20.43  57.33
   2.75  53.67  33.60  36.64  24.55
  33.47  28.05  32.46  33.39  31.35
  49.39  51.84  28.93  38.73   7.59
  47.98  60.80  45.20  56.51  34.16
  40.67  29.86  40.31  38.46  38.98
  41.66  43.29  40.75  52.92  33.84
  37.92  34.51  50.72   6.54  45.33
  34.99  23.18  70.52  46.94  16.39
  20.42  38.53  43.28  25.85   4.65
  56.43  51.78  51.15  12.05  28.54
  ;

  symbol1 c=white v=star h=1.5;

  title 'Two-Parameter Weibull Q-Q Plot of Failure Times';

  proc capability data=failures noprint;
     var time;
     qqplot / weibull2
              square
              href   = -5 to 1
              lhref  = 33
              vref   = 1 to 4
              lvref  = 33
              cframe = gray
              haxis  = axis1
              vaxis  = axis2;

     axis1 label=('log-Weibull Quantile');
     axis2 label=('log of Failure Time ');

  run;

  goptions reset=all;