Ледяной покров реки кажется цельным, но под ним есть трещины, которые меняют течение. Так и временной ряд: за ровной поверхностью чисел скрываются моменты, когда правила игры меняются. Эти моменты называют точками смены, и их поиск похож на простукивание льда в поисках слабых мест.
Классические методы ищут трещины, предполагая, что лёд однороден, то есть данные независимы. Но в реальности соседние значения ряда связаны, как соседние льдины: сегодняшняя температура зависит от вчерашней. Авторы предложили иначе: они разбивают ряд на окна, оценивают в каждом окне коэффициенты авторегрессии (зависимости от прошлого), а затем получают последовательность почти независимых многомерных точек. На этих точках уже можно применять быстрые алгоритмы поиска переломов.
Чтобы понять, сколько трещин действительно есть, а не примерещилось, используют штрафной критерий. Он похож на судью, который штрафует за каждую лишнюю трещину: сумма квадратичных потерь внутри сегментов плюс штраф за каждую точку смены. Авторы доказали, что если штраф растёт как логарифм от логарифма размера выборки, то метод почти наверняка выбирает истинное число точек смены. А вот критерий типа AIC с постоянным штрафом всегда даёт ненулевую вероятность переобучения, то есть находит лишние трещины.
На реальных данных метод сработал точно. В температуре востока США за 1895–2015 годы он нашёл точки смены в 1901, 1929, 1944, 2009 и 2012 годах, совпадающие с фазами Атлантической мультидекадной осцилляции. В индексе Эль-Ниньо за 1854–2015 годы обнаружен сдвиг около 1979 года, связанный с Тихоокеанской декадной осцилляцией. Лёд треснул именно там, где ожидали климатологи.
Сегодня метод можно использовать в климатологии для поиска смен климатических режимов, в финансах для обнаружения изменений волатильности, в сетевой безопасности для выявления аномалий трафика. Быстрое и надёжное простукивание льда экономит время и позволяет заметить трещину до того, как она станет полыньёй.
Температура востока США за последние 120 лет меняла свой режим всего пять раз, и эти смены совпадают с фазами Атлантической мультидекадной осцилляции, которая длится по 20–40 лет.