plot.abrem: Abrem Object Plotting

In abremout: Output functions for the Abernethy Reliability Methods package (abrem), including plotting.

Description

This function adds the .abrem method to plot from the graphics package. Currently, plot.abrem() behaves as follows:

• Displays the (life)time observations and fits (if any)of an abrem object. The calculated Weibull or Lognormal cumulative distribution is plotted on Weibull probability paper. Two-parameter Weibull plots are displayed as a straight line.

• (right-)Censored data (suspensions) are supported, where some observations are not failures.

• The mark() function draws a point on the plot with dotted lines extending towards the axes, for marking points of interest.

Usage

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

Arguments

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

Details

The code... argument can be any graphical parameter that can be supplied to plot.default, and any option that is set by the function options.abremplot.

Value

Currently, plot.abrem returns no value. In the future it will return the abrem object, with any additional plotting options appended.

To Do/Known issues/Bugs

• The legend positions are slightly to the right, overlapping the right Y axis (especially for small horizontal ranges).

• Plotting confidence bounds with blives=c() and conf.n=0 produces errors and "numeric(0)" in the legend.

• "numeric(0)" in the legend when supplying custom blives values.

• return arguments, with additional plotting information attached.

• Legends are still shown even when no fits were calculated: more specific when survival is not installed, i still attempts to show legends for method.fit="surv".

Author(s)

Jurgen Symynck [e-mail]

See Also

Abrem, abrem.fit, options.abremplot

Examples

data(p2.2,package="abrem")
da <- abrem.fit(Abrem(p2.2),dist="weibull",
    method.fit=c("mrr","xony","qbeta"),pch=3)
da <- abrem.conf(da,method.conf.blives="mcpivotal")
da <- abrem.fit(da,dist="weibull",
    method.fit=c("mrr","xony","bernard"),
    col="red",pch=4,lty=2)






plot(da,xlim=c(1e-2,1e4))

abremout documentation built on May 2, 2019, 6:09 p.m.