## ----include = FALSE----------------------------------------------------------
has_data <- requireNamespace("shewhartr", quietly = TRUE)
knitr::opts_chunk$set(collapse = TRUE, comment = "#>", fig.width = 7, fig.height = 4.5,
                      eval = has_data)

## ----setup--------------------------------------------------------------------
library(edetect)

## ----claims-------------------------------------------------------------------
claims <- shewhartr::claims_p
head(claims, 3)

p_chart <- edetect_p(x = claims$defects[16:30], n = claims$n[16:30],
                     phase1 = data.frame(x = claims$defects[1:15], n = claims$n[1:15]),
                     arl = 370, direction = "up")
edetect_alarm(p_chart)
autoplot(p_chart)

## ----claims-report------------------------------------------------------------
edetect_report(p_chart)

## ----pcb----------------------------------------------------------------------
defects <- shewhartr::pcb_solder$defects
c_chart <- edetect_c(defects[21:50], phase1 = defects[1:20], arl = 370, direction = "both")
edetect_alarm(c_chart)
max(tidy(c_chart)$evidence)
autoplot(c_chart)

## ----u------------------------------------------------------------------------
set.seed(1)
sizes <- sample(2:4, 20, replace = TRUE)
sizes <- sizes[cumsum(sizes) <= 50]
sample_id <- rep(seq_along(sizes), sizes)
samples <- data.frame(n = sizes,
                      x = as.integer(tapply(defects[seq_along(sample_id)], sample_id, sum)))
head(samples, 4)

u_chart <- edetect_u(samples$x, n = samples$n, center = 6, arl = 370, direction = "both")
edetect_alarm(u_chart)
autoplot(u_chart)