import statistics as st def CalcCorr(s1,s2): n = min(len(s1),len(s2)) mean1 = st.mean(s1) mean2 = st.mean(s2) cov = sum((s1i - mean1) * (s2i - mean2) for s1i, s2i in zip(s1, s2)) / n stdS1 = st.stdev(s1) stdS2 = st.stdev(s2) if stdS1 == 0 or stdS2 == 0: return 0 return cov/(stdS1*stdS2)