weibullness (version 1.24.1)

gp.plot: 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", ...)

Value

A list with the following components:

x

The sorted data

y

-log(-log(ppoints(n,a=a)))

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.

Author

Chanseok Park

Details

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

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 )