wblr.conf {WeibullR}R Documentation

Add Confidence Interval Bounds to wblr Objects

Description

This function adds confidence calculations to the last fit entity in a wblr object.

Usage

wblr.conf(x,...)

Arguments

x

Object of class "wblr".

...

Options for calculating confidence, and for plotting the results.

Details

This function adds confidence calculations to the last fit entity in a wblr object and adds the result to the object alongside any pre-existing confidence calculations.

Certain confidence calculations pertain to specific fit methods. The pivotal.rr confidence method can only apply to a rank regression fit, while a likelihood ratio calculation can only be performed on an mle fit.

Option arguments suitable for passing into wblr.conf:

method.conf

A character string describing the techniques used for calculating confidence interval bounds. Implemented methods are "pivotal-rr" (default), "bbb", "fm" "fmbounds", and "lrb". Methods must conform to the method.fit in the wblr.fit call immediately preceding the wblr.conf call. Method "pivotal-rr" requires a rank regression fit method. Methods "fm", "fmbounds, and "lrb" require a mle based fit.

dq

A named series of quantiles at which confidence interval bounds will be calculated.

num_dq

The number of points used for the "abrem" dq determination.

user_dq

A vector of quantiles set by user. Default c(seq(.01,.09,by=.01),seq(.10,.90,by=.10),seq(.91,.99, by=.01)).

ci

The double-sided confidence interval, also chi sq confidence level for likelihood ratio. Must be in a range <1 && >0, default is 0.9.

blife.pts

The probability points at which to report Blife on legend.

Specific controls for pivotal analysis only
Specific controls for likelihood ratio contour and bounds only

Additionally, graphical control options for the confidence interval bounds can be passed in with the following options:

Subsequent calls to wblr.fit and wblr.conf will inherit these options.

Value

The function returns its argument x, extended with the confidence calculations and any optional graphical and calculation arguments as passed to the function.

References

William Q. Meeker and Luis A. Escobar, (1998) "Statistical Methods for Reliability Data", Wiley-Interscience, New York

Robert B. Abernethy, (2008) "The New Weibull Handbook, Fifth Edition"

John I. McCool, (2012) "Using the Weibull Distribution: Reliability, Modeling and Inference"

Examples

## full dataset ##
da1 <- wblr(runif(10,100,1e4),label="Complete data")
da1 <- wblr.fit(da1)
da1 <- wblr.conf(da1,method.conf="pivotal-rr",col="red")
da1 <- wblr.conf(da1,method.conf="bbb",col="orange")
da1 <- wblr.fit(da1, method.fit="mle")
da1 <- wblr.conf(da1,method.conf="lrb",col="yellow3")

## plot datasets ##
## Not run: 
plot(da1,main="Comparison between different bound types.")

## End(Not run)

## with right-censored data
da2 <- runif(8,100,1e4)
da2 <- wblr(fail=da2,susp=rep(max(da2),2),label="Type II right-censored data")
da2 <- wblr.fit(da2)
da2 <- wblr.conf(da2,method.conf="pivotal-rr",col="blue1")
da2 <- wblr.conf(da2,method.conf="bbb",col="steelblue")
da2 <- wblr.fit(da2, method.fit="mle")
da2 <- wblr.conf(da2,method.conf="lrb",col="cyan3")
## Not run: 
plot(da2,main="Comparison between different bound types with right-censored data.")

## End(Not run)

[Package WeibullR version 1.0.10 Index]