.
. Starting Step 3.1 - Generate PPCC plot for Weibull
.
title automatic
title case upper
character x
line blank
y1label PPCC
multiplot 2 2
multiplot corner coordinates 2 2 98 98
set minmax 1
x1label gamma
ylimits 0.90 1.0
weibull ppcc plot y
justification left
move 22 31
let maxppcc = round(maxppcc,4)
text max ppcc = ^maxppcc
move 22 24
text shape = ^shape
let gamma1 = 1.5
let gamma2 = 2.5
ylimits 0.95 1.0
weibull ppcc plot y
delete gamma1 gamma2
ylimits
let maxppcc = round(maxppcc,4)
let shape2 = round(shape,3)
move 22 31
text max ppcc = ^maxppcc
move 22 24
text shape = ^shape2
let gamma = shape
y1label DATA
x1label THEORETICAL
weibull probability plot y
justification left
let ppa02 = round(ppa0,3)
let ppa12 = round(ppa1,3)
move 22 85
text location = ^ppa02
move 22 78
text slope = ^ppa12
set relative histogram area
title histogram with overlaid pdf
y1label AREA
x1label POLISHED WINDOW STRENGTH
xlimits 10 50
xtic offset 0 5
relative histogram y
multiplot 2 2 4
limits freeze
line solid
char BLANK
plot weipdf(x,gamma,ppa0,ppa1) for x = 20 0.1 45
limits
end of multiplot
title case asis
title automatic
x1label
y1label
line solid
character blank
xtic offset 5 5
LET X = DATA .01 .025 .05 .95 .975 .99
LET Y2 = WEIPPF(X,GAMMA,PPA0,PPA1)
SET WRITE REWIND OFF
WRITE WEIBPPF.OUT "Estimated percent points using Weibull Distribution"
WRITE WEIBPPF.OUT " "
WRITE WEIBPPF.OUT "PERCENT POINT        FAILURE TIME "
SET WRITE FORMAT F5.2,15X,F7.2
WRITE WEIBPPF.OUT X Y2
SET WRITE FORMAT
SET WRITE REWIND ON
. Finished Step 3.1: Page Back for Output