“…A number of methods and software packages have been proposed to aid researchers in capturing and characterizing unknown changes in time series data. However, they often fall short in practice because they are restricted to univariate data (Erdman & Emerson, 2007;Killick & Eckley, 2014;Ross, 2015) restricted to a specific statistic or statistical model (Barnett & Onnela, 2016;Dürre et al, 2015;Galeano & Wied, 2017;Zeileis et al, 2002) or impose other stringent assumptions on the data (Chen & Gupta, 2012;Davis et al, 2006;Lavielle, 2005). There are multivariate, non-parametric methods available (Arlot, Celisse, & Harchaoui, 2012;Bulteel et al, 2014;Lung-Yut-Fong, Lévy-Leduc, & Cappé, 2012;Matteson & James, 2014) that can potentially detect all types of changes, but they often lack power to detect changes in statistics that do not pertain to the mean level (e.g., correlation ;Cabrieto, Tuerlinckx, Kuppens, Wilhelm, et al, 2018) and do not uncover which statistics are exactly involved in the change.…”