Ejemplos de read_hmmout en Python

Ejemplo n.º 1

Archivo: pfam-domain-counts.py Proyecto: lierhan/bioinformatics

def did2count(hmmout, pids=False):
  pid2pfamdomains = pfam.read_hmmout(hmmout)
  if not pids: keys = pid2pfamdomains.keys()
  else: keys = pids
  did2count = {}
  for pid in keys:
    if not pid2pfamdomains.has_key(pid): continue
    for did in list(set([d.get_attr("hmm_name") for d in pid2pfamdomains[pid]])):
      if not did2count.has_key(did): did2count[did] = 0
      did2count[did] += 1
  return did2count

Ejemplo n.º 2

Archivo: pfam-domain-counts.py Proyecto: sjjose2/bioinformatics

def did2count(hmmout, pids=False):
    pid2pfamdomains = pfam.read_hmmout(hmmout)
    if not pids: keys = pid2pfamdomains.keys()
    else: keys = pids
    did2count = {}
    for pid in keys:
        if not pid2pfamdomains.has_key(pid): continue
        for did in list(
                set([d.get_attr("hmm_name") for d in pid2pfamdomains[pid]])):
            if not did2count.has_key(did): did2count[did] = 0
            did2count[did] += 1
    return did2count

Ejemplo n.º 3

def main(args):
    pid2domains = pfam.read_hmmout(args['hmmoutfile'])
    for pid, domains in pid2domains.iteritems():
        domains2clan(domains)

Ejemplo n.º 4

Archivo: pfam-pid2clan.py Proyecto: lierhan/bioinformatics

def main( args ):
  pid2domains = pfam.read_hmmout(args['hmmoutfile'])
  for pid, domains in pid2domains.iteritems():
    domains2clan(domains)

Ejemplo n.º 5

Archivo: pfam-pid2arrangement.py Proyecto: lierhan/bioinformatics

def main( args ):
  pid2domains = pfam.read_hmmout(args['hmmoutfile'])
  for pid, domains in pid2domains.iteritems():
    sys.stderr.write("\r" + pid + " "*20)
    arr = domains2arrangement(domains)
    print pid + "\t" + arr          

Ejemplo n.º 6

Archivo: pfam-pid2arrangement.py Proyecto: sjjose2/bioinformatics

def main(args):
    pid2domains = pfam.read_hmmout(args['hmmoutfile'])
    for pid, domains in pid2domains.iteritems():
        sys.stderr.write("\r" + pid + " " * 20)
        arr = domains2arrangement(domains)
        print pid + "\t" + arr