taylor_wrapper: Taylor's change point analyzer

In ggchangepoint: Combines Changepoint Analysis with 'ggplot2'

Taylor's change point analyzer

Description

Wraps ChangePointTaylor::change_point_analyzer(): the bootstrap-and-recursion procedure of Wayne Taylor that the quality-control and Six Sigma community uses as its default. Each candidate is scored by the bootstrap probability that a change occurred there, which gives a confidence level per changepoint and a confidence interval for its location, both carried onto the result.

Usage

taylor_wrapper(
  x,
  n_bootstraps = 1000,
  min_candidate_conf = 0.5,
  min_conf = 0.9,
  conf_level = 0.95,
  seed = NULL
)

Arguments

x	A numeric vector.
n_bootstraps	Bootstrap samples per candidate. Defaults to 1000; the engine accepts 100 to 1,000,000.
min_candidate_conf	Minimum confidence for a candidate to be considered, between 0.3 and 1. Defaults to 0.5.
min_conf	Minimum confidence for a changepoint to be reported, between 0.5 and 1. Defaults to 0.9.
conf_level	Confidence level of the reported location intervals, between 0.9 and 0.999 (the engine's range). Defaults to 0.95.
seed	Optional seed (the procedure is bootstrap-based). The seed is scoped to this call: .Random.seed is saved and restored, so a seeded call inside a simulation loop does not pin the loop's own stream.

Value

A ggcpt object with ci_lower/ci_upper (so autoplot(show_ci = TRUE) works) and a confidence column.

Series length, and why you cannot interrupt it

This engine is written for the series lengths quality control sees (hundreds to low thousands), and it does not scale. At n = 10{,}000 with the default n_bootstraps = 1000 it runs for minutes, and more importantly it runs where R cannot look: a setTimeLimit() of 45 seconds was still not honoured after 170, and one of 125 seconds after 200, so the call had to be killed from outside the session. R checks elapsed-time limits and keyboard interrupts at the same points, which means Ctrl-C will not stop it either.

So size the call before starting it rather than after. n_bootstraps is the knob (the cost is roughly linear in it), and the “Benchmarks” article lists the methods that do scale to long series. That page is web-only, because the sweep behind it takes over twenty minutes: it is published at https://pursuitofdatascience.github.io/ggchangepoint/articles/benchmarks.html rather than built into the package, so vignette() will not find it.

References

\insertRef

taylor2000changepointggchangepoint

See Also

Other changepoint engines: bcp_wrapper(), beast_wrapper(), bfast_wrapper(), binsegrcpp_wrapper(), bocpd_wrapper(), capa_wrapper(), cpm_wrapper(), cpop_wrapper(), cpt_wrapper(), decafs_wrapper(), ecp_wrapper(), envcpt_wrapper(), esac_wrapper(), fabisearch_wrapper(), fastcpd_wrapper(), fcov_wrapper(), fmean_wrapper(), fpop_wrapper(), geomcp_wrapper(), hdcov_wrapper(), hdreg_wrapper(), idetect_wrapper(), inspect_wrapper(), kcp_wrapper(), kwc_wrapper(), mcp_wrapper(), mosum_wrapper(), network_wrapper(), not_wrapper(), npmojo_wrapper(), nsp_wrapper(), ocd_wrapper(), pilliat_wrapper(), segmented_wrapper(), smuce_wrapper(), sn_wrapper(), strucchange_wrapper(), tguh_wrapper(), trend_wrapper(), var_wrapper(), wbs2_wrapper(), wbs_wrapper(), wbsts_wrapper()

Examples

set.seed(2026)
taylor_wrapper(c(rnorm(60), rnorm(60, 3)), n_bootstraps = 200, seed = 1)

ggchangepoint documentation built on Nov. 5, 2026, 9:06 a.m.