meme.adh.zoops

EM算法的改进

********************************************************************************
MEME - Motif discovery tool
********************************************************************************
MEME version 3.5.3 (Release date: 2006-03-08 09:30:25 +1000 (Wed, 08 Mar 2006))

For further information on how to interpret these results or to get
a copy of the MEME software please access http://meme.nbcr.net.

This file may be used as input to the MAST algorithm for searching
sequence databases for matches to groups of motifs.  MAST is available
for interactive use and downloading at http://meme.nbcr.net.
********************************************************************************


********************************************************************************
REFERENCE
********************************************************************************
If you use this program in your research, please cite:

Timothy L. Bailey and Charles Elkan,
"Fitting a mixture model by expectation maximization to discover
motifs in biopolymers", Proceedings of the Second International
Conference on Intelligent Systems for Molecular Biology, pp. 28-36,
AAAI Press, Menlo Park, California, 1994.
********************************************************************************


********************************************************************************
TRAINING SET
********************************************************************************
DATAFILE= /home/tbailey/tests/adh.s
ALPHABET= ACDEFGHIKLMNPQRSTVWY
Sequence name            Weight Length  Sequence name            Weight Length  
-------------            ------ ------  -------------            ------ ------  
2BHD_STREX               1.0000    255  3BHD_COMTE               1.0000    253  
ADH_DROME                1.0000    255  AP27_MOUSE               1.0000    244  
BA72_EUBSP               1.0000    249  BDH_HUMAN                1.0000    343  
BPHB_PSEPS               1.0000    275  BUDC_KLETE               1.0000    241  
DHES_HUMAN               1.0000    327  DHGB_BACME               1.0000    262  
DHII_HUMAN               1.0000    292  DHMA_FLAS1               1.0000    270  
ENTA_ECOLI               1.0000    248  FIXR_BRAJA               1.0000    278  
GUTD_ECOLI               1.0000    259  HDE_CANTR                1.0000    906  
HDHA_ECOLI               1.0000    255  LIGD_PSEPA               1.0000    305  
NODG_RHIME               1.0000    245  RIDH_KLEAE               1.0000    249  
YINL_LISMO               1.0000    248  YRTP_BACSU               1.0000    238  
CSGA_MYXXA               1.0000    166  DHB2_HUMAN               1.0000    387  
DHB3_HUMAN               1.0000    310  DHCA_HUMAN               1.0000    276  
FABI_ECOLI               1.0000    262  FVT1_HUMAN               1.0000    332  
HMTR_LEIMA               1.0000    287  MAS1_AGRRA               1.0000    476  
PCR_PEA                  1.0000    399  RFBB_NEIGO               1.0000    346  
YURA_MYXXA               1.0000    258  
********************************************************************************

********************************************************************************
COMMAND LINE SUMMARY
********************************************************************************
This information can also be useful in the event you wish to report a
problem with the MEME software.

command: meme /home/tbailey/tests/adh.s -mod zoops -protein -nostatus -nmotifs 2 

model:  mod=         zoops    nmotifs=         2    evt=           inf
object function=  E-value of product of p-values
width:  minw=            8    maxw=           50    minic=        0.00
width:  wg=             11    ws=              1    endgaps=       yes
nsites: minsites=        2    maxsites=       33    wnsites=       0.8
theta:  prob=            1    spmap=         pam    spfuzz=        120
em:     prior=       megap    b=           49980    maxiter=        50
        distance=    1e-05
data:   n=            9996    N=              33

sample: seed=            0    seqfrac=         1
Dirichlet mixture priors file: prior30.plib
Letter frequencies in dataset:
A 0.111 C 0.012 D 0.050 E 0.055 F 0.036 G 0.090 H 0.018 I 0.057 K 0.052 
L 0.092 M 0.027 N 0.041 P 0.041 Q 0.029 R 0.049 S 0.064 T 0.057 V 0.083 
W 0.010 Y 0.027 
Background letter frequencies (from dataset with add-one prior applied):
A 0.111 C 0.012 D 0.050 E 0.055 F 0.036 G 0.090 H 0.018 I 0.057 K 0.052 
L 0.092 M 0.027 N 0.041 P 0.041 Q 0.029 R 0.049 S 0.064 T 0.057 V 0.083 
W 0.010 Y 0.027 
********************************************************************************


********************************************************************************
MOTIF  1	width =   27   sites =  33   llr = 1084   E-value = 2.8e-167
********************************************************************************
--------------------------------------------------------------------------------
	Motif 1 Description
--------------------------------------------------------------------------------
Simplified        A  :162:56112:1:215:::4::::11:
pos.-specific     C  :1:::::::::::::::::::::::::
probability       D  ::::::::1:::1::::2:::::::::
matrix            E  ::::::::1:::21:::5::11:::::
                  F  :::::21:::3:::::1::::::::::
                  G  :2:::12::4:::1:1::111:5:1:1
                  H  :::::::::1::1:::1::::1:::::
                  I  :::::::11:::::::::1::::5:2:
                  K  ::::a:::1:::2:::1::121::1::
                  L  :::::::31:6:::4:215::1::1::
                  M  ::1::::1:1:1::2:1:1::::::::
                  N  :::::::::1:::1::::::::2:1:4
                  P  ::::::::::::::::::::4::::::
                  Q  ::::::::::::1::11:::11:::::
                  R  ::::::::2:::3::22::2111:4:1
                  S  :426::1:11:213:::1:1:1::::2
                  T  :1:2::::1::4:21:::::11::1:1
                  V  :::::::41::2:1::1:1::::3161
                  W  :::::1:::::::::::::::::::::
                  Y  9:::::::::::::::::2::2:::::

         bits    6.7                            
                 6.0                            
                 5.4                            
                 4.7 *                          
Information      4.0 *   *                      
content          3.4 *   *                      
(47.4 bits)      2.7 *  **     *                
                 2.0 *  ****   ***  * *    ** * 
                 1.3 ******** **** ** ** *******
                 0.7 ***************************
                 0.0 ---------------------------

Multilevel           YSASKAAVxGLTRSLAxELAPxGIRVN
consensus                 FGL  FSE MR D     V   
sequence                        V               
                                                
                                                
--------------------------------------------------------------------------------

--------------------------------------------------------------------------------
	Motif 1 sites sorted by position p-value
--------------------------------------------------------------------------------
Sequence name             Start   P-value                       Site          
-------------             ----- ---------            ---------------------------
YRTP_BACSU                  155  2.53e-22 GQRGAAVTSA YSASKFAVLGLTESLMQEVRKHNIRVS ALTPSTVASD
BUDC_KLETE                  152  4.81e-22 GHVGNPELAV YSSSKFAVRGLTQTAARDLAPLGITVN GFCPGIVKTP
AP27_MOUSE                  149  5.51e-21 AHVTFPNLIT YSSTKGAMTMLTKAMAMELGPHKIRVN SVNPTVVLTD
NODG_RHIME                  152  4.01e-20 GAIGNPGQTN YCASKAGMIGFSKSLAQEIATRNITVN CVAPGFIESA
FIXR_BRAJA                  189  1.18e-19 SRVHPFAGSA YATSKAALASLTRELAHDYAPHGIRVN AIAPGEIRTD
HMTR_LEIMA                  193  5.57e-19 TNQPLLGYTI YTMAKGALEGLTRSAALELAPLQIRVN GVGPGLSVLV
DHGB_BACME                  160  9.19e-19 WKIPWPLFVH YAASKGGMKLMTETLALEYAPKGIRVN NIGPGAINTP
GUTD_ECOLI                  154  2.44e-18 GKVGSKHNSG YSAAKFGGVGLTQSLALDLAEYGITVH SLMLGNLLKS
HDHA_ECOLI                  159  7.97e-18 AENKNINMTS YASSKAAASHLVRNMAFDLGEKNIRVN GIAPGAILTD
2BHD_STREX                  152  8.95e-18 GLMGLALTSS YGASKWGVRGLSKLAAVELGTDRIRVN SVHPGMTYTP
HDE_CANTR                   467  1.59e-17 GIYGNFGQAN YSSSKAGILGLSKTMAIEGAKNNIKVN IVAPHAETAM
FVT1_HUMAN                  186  1.78e-17 GQLGLFGFTA YSASKFAIRGLAEALQMEVKPYNVYIT VAYPPDTDTP
DHB2_HUMAN                  232  3.12e-17 GGAPMERLAS YGSSKAAVTMFSSVMRLELSKWGIKVA SIQPGGFLTN
YURA_MYXXA                  160  4.34e-17 AGFRGLPATR YSASKAFLSTFMESLRVDLRGTGVRVT CIYPGFVKSE
DHMA_FLAS1                  165  4.34e-17 SFMAEPEAAA YVAAKGGVAMLTRAMAVDLARHGILVN MIAPGPVDVT
RIDH_KLEAE                  160  6.03e-17 GVVPVIWEPV YTASKFAVQAFVHTTRRQVAQYGVRVG AVLPGPVVTA
DHES_HUMAN                  155  7.44e-16 GLMGLPFNDV YCASKFALEGLCESLAVLLLPFGVHLS LIECGPVHTA
BA72_EUBSP                  157  1.11e-15 GIFGSLSGVG YPASKASVIGLTHGLGREIIRKNIRVV GVAPGVVNTD
DHII_HUMAN                  183  1.23e-15 GKVAYPMVAA YSASKFALDGFFSSIRKEYSVSRVNVS ITLCVLGLID
BDH_HUMAN                   208  1.35e-15 GRMANPARSP YCITKFGVEAFSDCLRYEMYPLGVKVS VVEPGNFIAA