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