sen‘s slope检验

 其中  和  是时间  和   的数据值。这 N 个  值的中位数就是Sen的斜率估计值。

import pandas as pd
import numpy as np
#修改至当前文件夹工作路径
data = pd.read_csv('1.csv',header = 0)
print(data)
data_length = data.shape[0]
for i in data.columns[1:]:
	empty_list = []
	c = 0
	for k in range(data_length):
		c += 1
		for j in range(data_length-c):
			data_k = data[i][j]
			data_j = data[i][j+c]
			result = (data_j-data_k)/c
			empty_list.append(result)
	sen_slope = np.median(empty_list)
	print(sen_slope)