gp.plot {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 )