Mercurial > hg > cc > cirrus_work

comparison bin/spearman.py @ 26:5c5440e7854a

Find changesets by keywords (author, files, the commit message), revision number or hash, or revset expression.
a bit more
author Henry S. Thompson <ht@inf.ed.ac.uk>
date Tue, 15 Nov 2022 19:37:28 +0000
parents 50337cd1d16f
children 21da4d6521db
comparison
equal deleted inserted replaced
25:50337cd1d16f 26:5c5440e7854a
1 #!/usr/bin/env python3 1 #!/usr/bin/env python3
2 import numpy as np 2 import numpy as np
3 from numpy import loadtxt 3 from numpy import loadtxt
4 from scipy import stats 4 from scipy import stats
5 import statsmodels.api as sm 5 import statsmodels.api as sm
6 import matplotlib.pyplot as plt
6 import pylab 7 import pylab
7 8
8 import sys 9 import sys
9 10
10 cc19=loadtxt(sys.argv[1],delimiter=',') 11 cc19=loadtxt(sys.argv[1],delimiter=',')
11 cc19s_o=stats.spearmanr(cc19,nan_policy='omit') 12 cc19s_o=stats.spearmanr(cc19,nan_policy='omit')
12 cc19s_x=np.array([np.concatenate((cc19s_o.correlation[i][1:i],cc19s_o.correlation[i][i+1:])) for i in range(1,101)]) 13 cc19s_x=np.array([np.concatenate((cc19s_o.correlation[i][1:i],cc19s_o.correlation[i][i+1:])) for i in range(1,101)])
13 14 cc19s_xd=[stats.describe(cc19s_x[i]) for i in range(100)]