## ----setup, include = FALSE---------------------------------------------------
knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>"
)
library(gpciLindApproxProgII)

## ----example------------------------------------------------------------------
# Failure times and progressive removal scheme
x <- c(0.5, 1.2, 2.1, 3.4, 4.8)
r <- c(1, 0, 2, 0, 1)

# Fit model using Lindley approximation and chain generation
fit <- lindley_prog_gpci(
  x = x,
  r_removals = r,
  distribution = dist_weibull(),
  USL = 6, LSL = 0,
  chain_length = 200,
  burn_in = 50,
  thinning = 1,
  B = 50
)

# Print Summary Table
summary(fit)

## ----plot-fit, fig.width=7, fig.height=5--------------------------------------
plot(fit)