kwc_wrapper: Robust depth-based changepoints for functional and...

In ggchangepoint: Combines Changepoint Analysis with 'ggplot2'

Robust depth-based changepoints for functional and multivariate data

Description

Wraps KWCChangepoint (Ramsay and Chenouri): the functional Kruskal-Wallis covariance test, which ranks observations by statistical depth and segments on the ranks. Because it never touches the values themselves it is insensitive to heavy tails and outlying curves, which is exactly where the moment-based functional tests degrade.

Usage

kwc_wrapper(
  x,
  algorithm = c("fkwc", "dwbs"),
  depth = NULL,
  change_in = c("covariance", "distribution"),
  seed = NULL,
  ...
)

Arguments

x	A numeric matrix or data frame with one row per observation (time point) and one column per grid location or coordinate.
algorithm	"fkwc" (default; pruned exact linear time over the depth ranks) or "dwbs" (depth-based wild binary segmentation).
depth	Depth function. For "fkwc" one of "RPD" (default), "FM", "LTR", "FMd", "RPDd"; for "dwbs" one of "spat", "hs", "mahal", "mahal75".
change_in	Reported change type: "covariance" (default) or "distribution". The test is sensitive to both; this only labels the result.
seed	Optional seed: the random-projection depths and the wild binary segmentation both randomise. 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.
...	Additional arguments passed to the engine.

Value

A ggcpt object. Multivariate input is reduced to one series per observation by taking the cross-sectional mean of the columns, and that is the series stored on the result: autoplot() draws it, tidy()'s cp_value reads it, and $segments$param_estimate and augment()'s .fitted/.resid are computed from it. It is not any one column of the input. The full input is kept in $data_wide for autoplot(type = "coordinates").

References

\insertRef

ramsay2024kwcggchangepoint

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(), mcp_wrapper(), mosum_wrapper(), network_wrapper(), not_wrapper(), npmojo_wrapper(), nsp_wrapper(), ocd_wrapper(), pilliat_wrapper(), segmented_wrapper(), smuce_wrapper(), sn_wrapper(), strucchange_wrapper(), taylor_wrapper(), tguh_wrapper(), trend_wrapper(), var_wrapper(), wbs2_wrapper(), wbs_wrapper(), wbsts_wrapper()

Examples

set.seed(2026)
X <- matrix(rnorm(100 * 20), nrow = 100)
X[51:100, ] <- X[51:100, ] * 3
kwc_wrapper(X, seed = 1)

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