gp.plot: Gumbel Probability Plot

In weibullness: Goodness-of-Fit Test for Weibull Distribution (Weibullness)

Gumbel Probability Plot

Description

gp.plot produces a Gumbel probability plot.

Usage

gp.plot(x, plot.it=TRUE, a, col.line="black", lty.line=1, 
        xlim=NULL, ylim=NULL, main=NULL, sub=NULL, xlab=NULL, ylab="Probability", ...) 

Arguments

x	a numeric vector of data values. Missing values are allowed.
plot.it	logical. Should the result be plotted?
a	the offset fraction to be used; typically in (0,1). See ppoints.
col.line	the color of the straight line.
lty.line	the line type of the straight line.
xlim	the x limits of the plot.
ylim	the y limits of the plot.
main	a main title for the plot, see also title.
sub	a sub title for the plot.
xlab	a label for the x axis, defaults to a description of x.
ylab	a label for the y axis, defaults to "Probability".
...	graphical parameters.

Details

The Gumbel probability plot is based on taking the logarithm of the Gumbel cumulative distribution function twice.

Value

A list with the following components:

x	The sorted data
y	-log(-log(ppoints(n,a=a)))

Author(s)

Chanseok Park

See Also

plot, qqnorm, qqplot, wp.plot, iwp.plot, ep.plot.
bs.plot for the Birnbaum-Saunders probability plot in package bsgof.

Examples

 
x = c(-3.16, -3.07, -2.24, -1.8, -1.48, -0.92, -0.87, -0.41, -0.06, 1.15)
# With cosmetic lines 
gp.plot(x, main="Gumbel Probability Plot", col.line="red", 
        xlab="Lifetimes of bearings", lty.line=1, pch=3)
hline = -log(-log(c( (1:5)/100, (1:9)/10)))
abline( h=hline, col=gray(0.1), lty=3, lwd=0.5 )
abline( v=seq(-4, 2,by=0.5), col=gray(0.1), lty=3, lwd=0.5 )

weibullness documentation built on May 29, 2024, 1:27 a.m.