weibulltools (version 0.5.4)

confint_betabinom: Beta Binomial Confidence Bounds for Quantiles and/or Probabilities

Description

This non-parametric approach calculates confidence bounds for quantiles and/or failure probabilities using a procedure that is similar to that used in calculating median ranks. The location-scale parameters estimated by rank regression are needed.

Usage

confint_betabinom(x, event, loc_sc_params, distribution = c("weibull",
  "lognormal", "loglogistic"), bounds = c("two_sided", "lower", "upper"),
  conf_level = 0.95, direction = c("y", "x"))

Arguments

x

a numeric vector which consists of lifetime data. x is used to specify the range of confidence region(s).

event

a vector of binary data (0 or 1) indicating whether unit i is a right censored observation (= 0) or a failure (= 1).

loc_sc_params

a (named) numeric vector of estimated location and scale parameters for a specified distribution. The order of elements is important. First entry needs to be the location parameter \(\mu\) and the second element needs to be the scale parameter \(\sigma\).

distribution

supposed distribution of the random variable. The value can be "weibull", "lognormal" or "loglogistic". Other distributions have not been implemented yet.

bounds

a character string specifying the interval(s) which has/have to be computed. Must be one of "two_sided" (default), "lower" or "upper".

conf_level

confidence level of the interval. The default value is conf_level = 0.95.

direction

a character string specifying the direction of the computed interval(s). Must be either "y" (failure probabilities) or "x" (quantiles).

Value

A data frame containing the lifetime characteristic, interpolated ranks as a function of probabilities, the probabilities which are used to compute the ranks and computed values for the specified confidence bound(s).

Examples

# NOT RUN {
obs   <- seq(10000, 100000, 10000)
state <- c(0, 1, 1, 0, 0, 0, 1, 0, 1, 0)

df_john <- johnson_method(x = obs, event = state)
mrr <- rank_regression(x = df_john$characteristic,
                       y = df_john$prob,
                       event = df_john$status,
                       distribution = "weibull",
                       conf_level = .95)
conf_betabin <- confint_betabinom(x = df_john$characteristic,
                                  event = df_john$status,
                                  loc_sc_params = mrr$loc_sc_coefficients,
                                  distribution = "weibull",
                                  bounds = "two_sided",
                                  conf_level = 0.95,
                                  direction = "y")
# }