- def COFM(comat,normal=True):
- CoFreq = []
- for sp1 in xrange(len(comat)-1):
- for sp2 in xrange(sp1,len(comat)):
- NCo = 0
- for nsite in xrange(len(comat[0])):
- if comat[sp1][nsite]*comat[sp2][nsite] > 0:
- NCo += 1
- if normal:
- NCo /= len(comat[0])
- CoFreq.append(NCo)
- return CoFreq
- CoMat = [[0,0,1,0],[0,2,1,0],[1,1,0,1],[1,0,0,0]]
- print COFM(CoMat)