Error using the function gethmmprof

1 次查看(过去 30 天)
I am trying to use the "gethmmprof" function. Using the example given in the documentation, hmm_7tm = gethmmprof(2), I only get a series of error messages. Is the function up to date with the current database format? Thanks!
  2 个评论
Walter Roberson
Walter Roberson 2015-10-3
Please show the error messages.
gethmmprof appears to be part of the Bioinformatics Toolbox
per isakson
per isakson 2015-10-3
编辑:per isakson 2015-10-5
With R2013a
Error using pfamhmmread (line 85)
Record 1 in the inputfile is not a valid PFAM file. PFAMHMMREAD supports HMMMER2.0 and
HMMER3/b formatted files.
Error in getsangerdata (line 175)
out = pfamhmmread(s);
Error in gethmmprof (line 79)
model = getsangerdata(accessnum,'database','hmm',varargin{:});
and the two first lines of the downloaded string are
HMMER3/f [3.1b1 | May 2013]
NAME 7tm_2
Thus, with R2013a it seems to be a problem with the file format.
Which Matlab release?

请先登录,再进行评论。

采纳的回答

Paola Favaretto
Paola Favaretto 2015-10-5
The Bioinformatics Toolbox version you are running (R2013a) supports only HMMER2.0 and HMMER3/b. The format of the current database is HMMER3/f. The difference in the content is not significant, but the format is not recognized and therefore the function throws an error. If you contact Mathworks customer support they will send you a patch.
-Paola

更多回答(0 个)

类别

Help CenterFile Exchange 中查找有关 Genomics and Next Generation Sequencing 的更多信息

标签

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

Translated by