Legofit
infers population history from nucleotide site patterns.
Loading...
Searching...
No Matches
diverg.py File Reference

Functions

 diverg.usage ()
 diverg.openInput (fname)

Variables

bool diverg.dokl = False
 diverg.fname2 = None
 diverg.fname1 = sys.argv[1]
 diverg.f1 = openInput(fname1)
 diverg.f2 = openInput(fname2)
list diverg.pat1 = []
list diverg.prob1 = []
list diverg.pat2 = []
list diverg.prob2 = []
 diverg.line = line.strip()
 diverg.wid1 = max(9, len(fname1))
 diverg.wid2 = max(9, len(fname2))
int diverg.widpat = 7
 diverg.s1 = sum(prob1)
 diverg.s2 = sum(prob2)
str diverg.fmt = "%%%ds %%%ds" % (widpat, wid1)
float diverg.KLsum = 0.0
list diverg.kl = prob2[i]*log(prob2[i]/prob1[i])
 diverg.end