Reliability Growth Analysis.
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) orPiecewise 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%.
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"