Python计算两个时间序列的互相关

import pandas as pd
#刚开始做一次相关系数和P值分析
#把第二个列表复制一下
#原来的那个数组一直做右移,每移一次,记录一下右移的步长,相关系数和P值
#另外一个复制过的数组一直做左移,每移一次,记录一下左移的步长,相关系数和P值
from scipy.stats import pearsonr
A=[1,2,3,4,5,6,7,8,9,1,2,3]
B=[1.3,3.2,5.6,7,1,3,3,1,5,6,7,9]
lag=[]#滞后的步长
pearsonrs=[]
pvalues=[]

result=pearsonr(A,B)
lag.append(0)
pearsonrs.append(result[0])
pvalues.append(result[1])

###############不断进行右移,并记录所有参数
C=B.copy()
for i in range(20):
    C.insert(0,C.pop())
    result=pearsonr(A,C)
    lag.append(i+1)
    pearsonrs.append(result[0])
    pvalues.append(result[1])

###############不断进行左移,并记录所有参数
D=B.copy()
for i in range(20):
    D.insert(len(D),D[0])
    D.remove(D[0])
    result=pearsonr(A,D)
    lag.append(-(i+1))
    pearsonrs.append(result[0])
    pvalues.append(result[1])
  • 3
    点赞
  • 7
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值