示例#1
0
def dbEntropy(inFileName,wordsize):
	"""
	deprecated
	"""
	vec_len=[]
	stat=Stat()
	file_db=open(inFileName)
	seq=Bioseq()
	numseq=0
	while 1:
		seq.read(file_db)
		if seq.sequence==None:
			break
		i=seq.entropy(wordsize)
		stat.add(i)

		numseq=numseq+1
		print 'sequence #',numseq,'=',seq.getLength(),'[',seq.header[0:40],'...] entropy',i
		vec_len.append((-i,numseq,seq.header))

	file_db.close()
	vec_len.sort()
	for s in vec_len:
		print 'I=',-s[0],'=> #',s[1], s[2]
	print stat.string()
	return vec_len