Skip to contents

Reliability Growth Analysis.

Usage

rga(
  times,
  failures,
  model_type = "Crow-AMSAA",
  breakpoints = NULL,
  conf_level = 0.95
)

Arguments

times

A vector of cumulative times at which failures occurred.

failures

A vector of the number of failures at each corresponding time in times.

model_type

The model type. Either Crow-AMSAA (default) or Piecewise Weibull NHPP with change point detection.

breakpoints

An optional vector of breakpoints for the Piecewise Weibull NHPP model.

conf_level

The desired confidence level, which defaults to 95%.

Value

The function returns an object of class rga that contains the results for the model.

Examples

times <- c(100, 200, 300, 400, 500)
failures <- c(1, 2, 1, 3, 2)
result <- rga(times, failures)
print(result)
#> $model
#> 
#> Call:
#> stats::lm(formula = log_cum_failures ~ log_times)
#> 
#> Coefficients:
#> (Intercept)    log_times  
#>     -3.6145       0.7987  
#> 
#> 
#> $AIC
#> [1] -3.551005
#> 
#> $BIC
#> [1] -4.722691
#> 
#> $breakpoints
#> NULL
#> 
#> $fitted_values
#>        1        2        3        4        5 
#> 1.065602 2.562478 4.457440 6.703022 9.266388 
#> 
#> $lower_bounds
#>         1         2         3         4         5 
#> 0.7656742 2.1157687 3.7410494 5.4093372 7.1120671 
#> 
#> $upper_bounds
#>         1         2         3         4         5 
#>  1.483017  3.103503  5.311015  8.306102 12.073275 
#> 
#> $betas
#>     Estimate   Std. Error      t value     Pr(>|t|) 
#> 7.986756e-01 5.624037e-02 1.420111e+01 7.564937e-04 
#> 
#> $lambdas
#>     Estimate   Std. Error      t value     Pr(>|t|) 
#> 7.986756e-01 5.624037e-02 1.420111e+01 7.564937e-04 
#> 
#> attr(,"class")
#> [1] "rga"