## ----include = FALSE----------------------------------------------------------
knitr::opts_chunk$set(collapse = TRUE, comment = "#>",
                      fig.width = 7, fig.height = 4.2)

## ----setup, message = FALSE---------------------------------------------------
library(shewhartr)
library(ggplot2)
library(dplyr)

## -----------------------------------------------------------------------------
data(cvd_recife)
head(cvd_recife)

## -----------------------------------------------------------------------------
fit_imr <- shewhart_i_mr(cvd_recife,
                         value = new_deaths,
                         index = date)
broom::glance(fit_imr)

## -----------------------------------------------------------------------------
fit <- shewhart_regression(
  cvd_recife,
  value      = new_deaths,
  index      = .t,
  model      = "loglog",
  phase_rule = "we_seven_same",   # legacy WE rule used in the original analysis
  rules      = c("nelson_1_beyond_3s", "we_seven_same")
)
broom::glance(fit)
length(fit$fits)        # number of phases

## -----------------------------------------------------------------------------
autoplot(fit)

## -----------------------------------------------------------------------------
shewhart_box_cox(cvd_recife$new_deaths + 1)$lambda_hat