关于MK突变监测的突变点读取

curtainsky

各路大神，想请教在利用时间序列数据做MK突变时，可以得到一张图，两条曲线相交即为突变点。
那想问下是否有办法读取具体的相交点横坐标？毕竟有大量数据，不可能一张一张图去看。
或者想请教是否有其他方法读取一段时间序列数据中的拐点？谢谢大家。

田小溪

这个问题我也很想知道

lqouc

假设x=1,n-1
(a(x)-b(x))*(a(x+1)-b(x+1))<0时
x或x+1可以认为是交叉点

curtainsky

lqouc 发表于 2017-12-27 15:31
假设x=1,n-1
(a(x)-b(x))*(a(x+1)-b(x+1))

您好！感谢回答。
想请问，这样的话不是要交叉点在X=0点附件交叉才可以满足吗？
我可能水平太低，不太理解，还希望你能具体一些说明。
我是想知道如何由一段时间序列数据中得出拐点的具体坐标。

curtainsky

lqouc 发表于 2017-12-27 15:31
假设x=1,n-1
(a(x)-b(x))*(a(x+1)-b(x+1))

您好！感谢回答。
想请问，这样的话不是要交叉点在X=0点附件交叉才可以满足吗？
我可能水平太低，不太理解，还希望你能具体一些说明。
我是想知道如何由一段时间序列数据中得出拐点的具体坐标。

lqouc

curtainsky 发表于 2017-12-27 15:57
您好！感谢回答。
想请问，这样的话不是要交叉点在X=0点附件交叉才可以满足吗？
我可能水平太低，不太 ...

a和b式mk得到的正反序列，序列长度为n，上面的表达式基本思想就是a-b在相邻两个点符号相反。
也就是在x点a<b，x+1点a>b，这样就算是交叉了。
上面写的那个表达式加上for循环和if块直接就能用了。

当然这里需要提醒你一点，mk检验有可能存在多个交点，所以还是需要你自己看图判断的，如果有特别特别多的图，那么只能说你要做的东西不适合直接用mk检验，需要对数据预处理。