Don't like ads? PRO users don't see any ads ;-)
Guest

Untitled

By: a guest on Jun 22nd, 2012  |  syntax: None  |  size: 0.41 KB  |  hits: 7  |  expires: Never
download  |  raw  |  embed  |  report abuse  |  print
Text below is selected. Please press Ctrl+C to copy to your clipboard. (⌘+C on Mac)
  1. def COFM(comat,normal=True):
  2.         CoFreq = []
  3.         for sp1 in xrange(len(comat)-1):
  4.                 for sp2 in xrange(sp1,len(comat)):
  5.                         NCo = 0
  6.                         for nsite in xrange(len(comat[0])):
  7.                                 if comat[sp1][nsite]*comat[sp2][nsite] > 0:
  8.                                         NCo += 1
  9.                         if normal:
  10.                             NCo /= len(comat[0])
  11.                         CoFreq.append(NCo)
  12.         return CoFreq
  13.        
  14. CoMat = [[0,0,1,0],[0,2,1,0],[1,1,0,1],[1,0,0,0]]
  15.  
  16. print COFM(CoMat)