plot.abrem: 'Abrem' Object Plotting

In abrem: Abernethy Reliability Methods

Description

This function adds the .abrem method to plot from the graphics package.

Currently, the function plots the (life-)time observations, fits (if any) and confidence bounds for B-lives (if any) of an abrem object or a list of abrem objects on Weibull or Lognormal probability paper.

For each fit in the (list of) abrem object(s), legends are added to the plot, displaying the fit parameters and (if available) goodness-of-fit indicators and confidence information.

Usage

## S3 method for class 'abrem'
plot(x, ...)

Arguments

x	Object of class "abrem" or a list of abrem objects.
...	Options for plotting the abrem object; see section "Details".

Details

The ... argument can be any graphical parameter that can be supplied to plot.default, and any option that can be set by the function options.abrem. The options set in this way are applied to all graphical elements of the plot, overriding any previously supplied options.

One can pass a list of abrem objects to plot.abrem; in that case it is mandatory to use the full method name: plot.abrem(...) and not plot(...).

The calculated Weibull or Lognormal distribution fits are plotted on Weibull probability paper by default, but by passing the argument log = "xy" to the function, Lognormal paper is used (see options.abrem).

When a list of abrem objects is passed, the plot window is generated with the options of the first abrem object in the list.

Value

Currently, the function returns no value.

Author(s)

Jurgen Symynck [e-mail]

See Also

options.abrem

Examples

da1 <- Abrem(rweibull(5,3,1000))
da1 <- abrem.fit(da1,dist="weibull",  method.fit=c("rr","xony"),pch=3)
da1 <- abrem.fit(da1,dist="weibull",  method.fit="mle",col="red1")
da1 <- abrem.fit(da1,dist="weibull",  method.fit="mle-rba",col="orange",lty=2)
da1 <- abrem.fit(da1,dist="lognormal",method.fit=c("rr","xony"),
    col="steelblue3",pch=8)
da1 <- abrem.conf(da1)
plot(da1,col="black") # plot all in black
plot(da1,main="Weibull distributed observations") # plot with specified colors
    
### usage of lists ###
options.abrem(blives=0.1) # make the legend boxes a bit shorter...
da2 <- abrem.conf(abrem.fit(Abrem(runif(5,10,100),col="red")))
da3 <- abrem.conf(abrem.fit(Abrem(rweibull(5,2,1000),col="green4",pch=3)))
da4 <- abrem.conf(abrem.fit(Abrem(rlnorm(5,log(500),log(2)),col="blue3",pch=8),
    dist="lognormal"))

plot.abrem(list(da2,da3,da4),xlim=c(1,1e6),
    main="Uniformly distributed observations")

abrem documentation built on May 2, 2019, 4:49 p.m.