时间序列|change point detection

2024-09-08 09:32

本文主要是介绍时间序列|change point detection,希望对大家解决编程问题提供一定的参考价值,需要的开发者们随着小编来一起学习吧!

change point detection 被称为变点检测,其基本定义是在一个序列或过程中,当某个统计特性(分布类型、分布参数)在某时间点受系统性因素而非偶然因素影响发生变化,我们就称该时间点为变点。变点识别即利用统计量或统计方法或机器学习方法将该变点位置估计出来。

Change Point Detection的类型

  • online
    指连续观察某一随机过程,监测到变点时停止检验,不运用到未来数据,主要用于事件预警。
  • offline
    从已经获得的时序数列中检测过去的变点位置,主要用作历史检验。

1. 控制图方法CUSUM

#####原画与介绍
根据累计数据微小的偏差,来探测数据分布是否发生了变化。这也是最古老最原始的,被业界广泛应用于工业质量检测、自动监测、金融方面。比如被阿里用于全景业务平台的监控解决方案。

基本模型流程

param: threshold
input:time series
for 0---->n:计算累计偏差量
如果累计偏差量大于阈值threshold,则视为这段time series与前段time series有明显不同,分割出两个区间windows,从开始出现偏差的点,可以推出change point
output:windows, change points

优缺点

  • 缺点:需要调参threshold,太小模型太敏感,太大模型太粗糙。
    解决方案:动态阈值预测:样本选择---->异常样本筛除—>样本截取—>预测基准值(一阶导数的状态值)。

  • 优点:比较简单,比较容易实施。业界研究与应用比较多,可以找到很相关材料与优化方法。

2. Probability Density Estimation

#####想法与原理
对于一组时间序列,出现change point前与出现change point后的概率密度分布会不同。

基本模型流程

  1. 用前n个点,probability density models—》estimate probability density function
  2. 然后用score来衡量,加入这个点后,概率密度分布的差异(加入这个点后,概率密度分布的变化大小 )
  3. score越高,这个点是change point的概率越高。

优缺点

  • 缺点:概率密度估计模型的确立,带参数。并且,probability density function很难做的准确,需要很大的数据量。

3. Direct Compute

想法

因为概率密度分布很难做的准确,因此衍生出了不计算评估概率密度分布,而是直接一个point前后的概率目的分布的差异,对于一个点前的数据和之后的数据,用一些模型/算法可以衡量他前后分布的差异。

可以比较多种差异:

  • 均值
  • 方差
  • 均值与方差

相关的模型:量化两种概率分布P和Q之间差异

  • Kernel Mean
  • non-parametric Gaussian kernel model
  • Kulback-Leibler Importance Estimation Procedure
  • Kernel FDA

如何分割time series?

  • BinSegment
  • bottom-up
  • Window-based
  • Segment Neighbourhood

优缺点
该方法很容易受data noise的影响,高维度可能效果很差。

4. Probability Method

对比change point前后分布的不同,这部分侧重于直接预测某个点是否是change point.

(1)Gaussian Process

想法:用t之前的N点构建一个time series prediction 的模型(这里就是GP模型),如果t点的值和预测的值偏差很大,就记一个potential Alarm,当连续几个值偏差很大,potential Alarm超过某个阈值时,就Raise Alarm,说明出现新的分布。

优缺点

这个方法对one-day-event不敏感。

###(2)Bayesian
一开始和GP差不多。

想法:对于一个点,给定他在上一个change point之后的信息,估计这个点是change point的概率。

优缺点
先验函数很难定义,定的不好,结果不准 。

一般来说offline的方法比online的更准确,然后GP会比Bayesian更准确。

5. Clustering Method

把change point detection看成,将time series分成很多的clusters。

  • (1)层次聚类:将许多time series进行聚类,在同一类中,如果有一个time series的行为与同聚类内其它成员差异较大,则视作发生了change。
  • (2)图形聚类:
  • (3)基于局部图形聚类与分割

6. 其它

基于图论、基于控制论、系统辩识方法等。

参考文献

[1]. Aminikhanghahi, Samaneh and Diane J. Cook. “A survey of methods for time series change point detection.” Knowledge and Information Systems 51 (2016): 339-367.

[2]. Liu, Song, Makoto Yamada, Nigel Collier and Masashi Sugiyama. “Change-Point Detection in Time-Series Data by Relative Density-Ratio Estimation.” Neural networks : the official journal of the International Neural Network Society 43 (2012): 72-83.

[3]. Itoh, Naoki and Juergen Kurths. “Change-Point Detection of Climate Time Series by Nonparametric Method.” (2010).

[4]. Cho, Haeran and Piotr Fryzlewicz. “Multiple change-point detection for high-dimensional time series via Sparsified Binary Segmentation.” (2013).

[5]. Saatci, Yunus, Ryan D. Turner and Carl E. Rasmussen. “Gaussian Process Change Point Models.” ICML (2010).

[6]. Lacasa, Lucas, Bartolo Luque, Fernando J Ballesteros, Jordi Luque and Juan Carlos Nuño. “From time series to complex networks: the visibility graph.” Proceedings of the National Academy of Sciences of the United States of America 105 13 (2008): 4972-5.

[7]. Harchaoui, Zaïd, Francis R. Bach and Eric Moulines. “Kernel Change-point Analysis.” NIPS (2008).

[8]. Jeske, Daniel R., Veronica Montes De Oca, Wolfgang Bischoff and Mazda Marvasti. “Cusum techniques for timeslot sequences with applications to network surveillance.” Computational Statistics & Data Analysis 53 (2009): 4332-4344.

[9]. Chib, Siddhartha and John M. Olin. “Estimation and comparison of multiple change-point models.” (1997).

这篇关于时间序列|change point detection的文章就介绍到这儿,希望我们推荐的文章对编程师们有所帮助!



http://www.chinasem.cn/article/1147738

相关文章

服务器集群同步时间手记

1.时间服务器配置(必须root用户) (1)检查ntp是否安装 [root@node1 桌面]# rpm -qa|grep ntpntp-4.2.6p5-10.el6.centos.x86_64fontpackages-filesystem-1.41-1.1.el6.noarchntpdate-4.2.6p5-10.el6.centos.x86_64 (2)修改ntp配置文件 [r

uva 10131 最长子序列

题意: 给大象的体重和智商,求体重按从大到小,智商从高到低的最长子序列,并输出路径。 代码: #include <iostream>#include <cstdio>#include <cstdlib>#include <algorithm>#include <cstring>#include <cmath>#include <stack>#include <vect

POJ1631最长单调递增子序列

最长单调递增子序列 import java.io.BufferedReader;import java.io.InputStream;import java.io.InputStreamReader;import java.io.PrintWriter;import java.math.BigInteger;import java.util.StringTokenizer;publ

fzu 2277 Change 线段树

Problem 2277 Change Time Limit: 2000 mSec    Memory Limit : 262144 KB  Problem Description There is a rooted tree with n nodes, number from 1-n. Root’s number is 1.Each node has a value ai.

MiniGPT-3D, 首个高效的3D点云大语言模型,仅需一张RTX3090显卡,训练一天时间,已开源

项目主页:https://tangyuan96.github.io/minigpt_3d_project_page/ 代码:https://github.com/TangYuan96/MiniGPT-3D 论文:https://arxiv.org/pdf/2405.01413 MiniGPT-3D在多个任务上取得了SoTA,被ACM MM2024接收,只拥有47.8M的可训练参数,在一张RTX

leetcode105 从前序与中序遍历序列构造二叉树

根据一棵树的前序遍历与中序遍历构造二叉树。 注意: 你可以假设树中没有重复的元素。 例如,给出 前序遍历 preorder = [3,9,20,15,7]中序遍历 inorder = [9,3,15,20,7] 返回如下的二叉树: 3/ \9 20/ \15 7   class Solution {public TreeNode buildTree(int[] pr

批处理以当前时间为文件名创建文件

批处理以当前时间为文件名创建文件 批处理创建空文件 有时候,需要创建以当前时间命名的文件,手动输入当然可以,但是有更省心的方法吗? 假设我是 windows 操作系统,打开命令行。 输入以下命令试试: echo %date:~0,4%_%date:~5,2%_%date:~8,2%_%time:~0,2%_%time:~3,2%_%time:~6,2% 输出类似: 2019_06

【MRI基础】TR 和 TE 时间概念

重复时间 (TR) 磁共振成像 (MRI) 中的 TR(重复时间,repetition time)是施加于同一切片的连续脉冲序列之间的时间间隔。具体而言,TR 是施加一个 RF(射频)脉冲与施加下一个 RF 脉冲之间的持续时间。TR 以毫秒 (ms) 为单位,主要控制后续脉冲之前的纵向弛豫程度(T1 弛豫),使其成为显著影响 MRI 中的图像对比度和信号特性的重要参数。 回声时间 (TE)

LeetCode:64. 最大正方形 动态规划 时间复杂度O(nm)

64. 最大正方形 题目链接 题目描述 给定一个由 0 和 1 组成的二维矩阵,找出只包含 1 的最大正方形,并返回其面积。 示例1: 输入: 1 0 1 0 01 0 1 1 11 1 1 1 11 0 0 1 0输出: 4 示例2: 输入: 0 1 1 0 01 1 1 1 11 1 1 1 11 1 1 1 1输出: 9 解题思路 这道题的思路是使用动态规划

O(n)时间内对[0..n^-1]之间的n个数排序

题目 如何在O(n)时间内,对0到n^2-1之间的n个整数进行排序 思路 把整数转换为n进制再排序,每个数有两位,每位的取值范围是[0..n-1],再进行基数排序 代码 #include <iostream>#include <cmath>using namespace std;int n, radix, length_A, digit = 2;void Print(int *A,