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;