TAMO.MD.AlignAce
index
/home/David_Gordon/docs/TAMO/MD/AlignAce.py

AlignAce.py -- Interface / Wrapper for the AlignACE program.
 
Download AlignACE separately from author's website (see TAMO/paths.py for instructions).
We have succesfully built AlignACE on linux and Win32 platforms (via cygwin)
 
This program is executable from the command line, and it runs the AlignACE several times
with different number seeds.  Type "$TAMO/MD/AlignAce.py" for options.
 
Copyright (2005) Whitehead Institute for Biomedical Research (except as noted below)
All Rights Reserved
 
Author: David Benjamin Gordon

 
Modules
       
TAMO
math
os
random
re
sys
tempfile

 
Classes
       
AlignAce
MetaAce

 
class AlignAce
    ################################################################################
AlignAce                                                                     #
# Execute, Load, and analyze AlignACE output                                   #
################################################################################
 
Usage:
theAce = AlignAce.AlignAce('filename.fsa',width=10)
 
-or-
 
theAce = AlignAce.AlginAce('filename.ace')
 
The first form executes AlignACE with the specificied input sequences, and the second
is used to parse alignace output from a previous run.
 
Typically, the user is only interested in the ".motifs" member variable, which contains
a list of all the discovered Motifs (as Motif objects).
 
  Methods defined here:
__init__(self, file='', width=10, args='', seed='')
__repr__(self)

 
class MetaAce
    Usage: theMetaAce = AlignAce.MetaAce(fastafile, width=10, iterations = 5, gcback=0.38)
 
Class to invoke AlignACE many times and combine the results.  A random number
generator is used to ensure that AlignACE is seeded differently each time.
 
Typically, the user is only interested in the ".results" member variable, which contains
a list of all the discovered Motifs (as Motif objects).
 
  Methods defined here:
__init__(self, fastafile='', width=10, iterations=5, gcback=0.38)
flatten(self, thresh=0.5)
flatten(thresh=0.5)
This function does some messy reductions in metaACE lists, to cut the number of
identical motifs reported.  A better solution is to use one of the clustering
methods elsewhere in the TAMO package.

 
Functions
       
flatten(motifs, thresh=0.5, prop='church', secondary_prop='pvalue')
flatten(motifs,thresh=0.5,prop='church',secondary_prop='pvalue'):
This function does some messy reductions in metaACE lists, to cut the number of
identical motifs reported.  A better solution is to use one of the clustering
methods elsewhere in the TAMO package.
 
"prop" and "secondary_prop" are used to settle ties between otherwise similar motifs.
"church" refers to the AlignACE "group specificity score," and "pvalue" refers to
the Enrichment score.
main()