5 meme分析后显示error at: maxw

进行meme分析,输入命令:/biosoft/meme/meme-v4.12.0/bin/meme ../pep_id-seq-F.fa -protein -oc ./ -nostatus -time 18000 -maxsize 6000000 -mod anr -nmotifs 10 -minw 6 maxw 100

显示:error at: maxw及以下信息,反复查看找不到原因,老师能否解决以下。

Usage:meme<dataset> [optional arguments]

<dataset> file containing sequences in FASTA format
[-h]print this message
[-o <output dir>]name of directory for output files
will not replace existing directory
[-oc <output dir>]name of directory for output files
will replace existing directory
[-text]output in text format (default is HTML)
[-dna]sequences use DNA alphabet
[-rna]sequences use RNA alphabet
[-protein]sequences use protein alphabet
[-alph <alph file>]sequences use custom alphabet
[-mod oops|zoops|anr]distribution of motifs
[-nmotifs <nmotifs>]maximum number of motifs to find
[-evt <ev>]stop if motif E-value greater than <evt>
[-nsites <sites>]number of sites for each motif
[-minsites <minsites>]minimum number of sites for each motif
[-maxsites <maxsites>]maximum number of sites for each motif
[-wnsites <wnsites>]weight on expected number of sites
[-w <w>]motif width
[-minw <minw>]minimum motif width
[-maxw <maxw>]maximum motif width
[-nomatrim]do not adjust motif width using multiple
alignment
[-wg <wg>]gap opening cost for multiple alignments
[-ws <ws>]gap extension cost for multiple alignments
[-noendgaps]do not count end gaps in multiple alignments
[-bfile <bfile>]name of background Markov model file
[-revcomp]allow sites on + or - DNA strands
[-pal]force palindromes (requires -dna)
[-maxiter <maxiter>]maximum EM iterations to run
[-distance <distance>]EM convergence criterion
[-psp <pspfile>]name of positional priors file
[-prior dirichlet|dmix|mega|megap|addone]
type of prior to use
[-b <b>]strength of the prior
[-plib <plib>]name of Dirichlet prior file
[-maxwords <maxwords>]maximum number of words to test as EM starts
[-spfuzz <spfuzz>]fuzziness of sequence to theta mapping
[-spmap uni|pam]starting point seq to theta mapping type
[-cons <cons>]consensus sequence to start EM from
[-heapsize <hs>]size of heaps for widths where substring 
search occurs
[-x_branch]perform x-branching
[-w_branch]perform width branching
[-allw]include all motif widths from min to max
[-bfactor <bf>]branching factor for branching search
[-maxsize <maxsize>]maximum dataset size in characters
[-nostatus]do not print progress reports to terminal
[-p <np>]use parallel version with <np> processors
[-time <t>]quit before <t> CPU seconds consumed
[-sf <sf>]print <sf> as name of sequence file
[-V]verbose mode
[-version]display the version number and exit
请先 登录 后评论

最佳答案 2019-04-04 17:47

缺失一个 -  :

attachments-2019-04-CSffKquK5ca55c095dc45.jpg

细心看看错误提示,一般也能找到错误;


linux基础不好建议学习:linux系统使用

请先 登录 后评论

其它 0 个回答