Python如何實(shí)現(xiàn)MK檢驗(yàn)

本篇內(nèi)容介紹了“Python如何實(shí)現(xiàn)MK檢驗(yàn)”的有關(guān)知識(shí)，在實(shí)際案例的操作過程中，不少人都會(huì)遇到這樣的困境，接下來就讓小編帶領(lǐng)大家學(xué)習(xí)一下如何處理這些情況吧！希望大家仔細(xì)閱讀，能夠?qū)W有所成！

MK檢驗(yàn)：時(shí)間序列進(jìn)行檢測(cè)，并找出突變點(diǎn)，本文參考網(wǎng)上的matlab程序改寫為python代碼如下：

import numpy as np
import pandas as pd
import matplotlib.pyplot as plt

#讀取時(shí)間序列數(shù)據(jù)
data = pd.read_csv('')
#定義時(shí)間和徑流數(shù)據(jù)列
x = list(range(len(data)))
y = data.to_list
#獲取樣本數(shù)據(jù)
n = len(y)

#正序計(jì)算
#定義累計(jì)量序列Sk，長度n，初始值為0
Sk = np.zeros(n)
UFk = np.zeros(n)
#定義Sk序列元素s
s = 0

#i從2開始，根據(jù)統(tǒng)計(jì)量UFk公式，i=1時(shí)，Sk（1）、E（1）、Var（1）均為0，此時(shí)UFk無意義，因此公式中，令UFk（1） = 0
for i in range(2,n):
    for j in range(1,i):
        if y[i]>y[j]:
            s += 1
    Sk[i] = s
    E = i * (i - 1)/4
    Var = i * (i - 1) * (2 * i + 5)/72
    UFk[i] = (Sk[i] - E)/np.sqrt(Var)

#逆序計(jì)算，構(gòu)造逆序列y2，長度為n，初值為0
#定義逆累計(jì)量序列Sk2，長度n，初始值為0
#定義逆統(tǒng)計(jì)量序列Sk2，長度n，初始值為0
y2 = np.zeros(n)
Sk2 = np.zeros(n)
UBk = np.zeros(n)

#s歸零
s = 0
#按時(shí)間序列逆轉(zhuǎn)樣本
y2 = y[::-1]

# i從2開始，根據(jù)統(tǒng)計(jì)量UBk公式，i=1時(shí)，Sk（1）、E（1）、Var（1）均為0，此時(shí)UFk無意義，因此公式中，令UBk（1） = 0
for i in range(2, n):
    for j in range(1, i):
        if y2[i] > y2[j]:
            s += 1
    Sk2[i] = s
    E = i * (i - 1) / 4
    Var = i * (i - 1) * (2 * i + 5) / 72
    UBk[i] = -(Sk2[i] - E) / np.sqrt(Var)

#逆轉(zhuǎn)逆序列
UBk2 = UBk[::-1]

#畫圖

plt.figure(figsize=(10,5))
plt.plot(range(1 ,n+1),UFk,label = 'UFk',color = 'orange')
plt.plot(range(1 ,n+1),UBk2,label = 'UBk',color = 'cornflowerblue')
plt.ylabel('UFk-UBk')
x_lim = plt.xlim()
plt.plot(x_lim,[-1.96,-1.96],'m--',color = 'r')
plt.plot(x_lim, [0,0],'m--')
plt.plot(x_lim,[1.96,1.96],'m--',color = 'r')
plt.show()

“Python如何實(shí)現(xiàn)MK檢驗(yàn)”的內(nèi)容就介紹到這里了，感謝大家的閱讀。如果想了解更多行業(yè)相關(guān)的知識(shí)可以關(guān)注億速云網(wǎng)站，小編將為大家輸出更多高質(zhì)量的實(shí)用文章！