Borrar filtros
Borrar filtros

Error using the function gethmmprof

2 visualizaciones (últimos 30 días)
Arnold Hayer
Arnold Hayer el 3 de Oct. de 2015
Editada: per isakson el 5 de Oct. de 2015
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 comentarios
Walter Roberson
Walter Roberson el 3 de Oct. de 2015
Please show the error messages.
gethmmprof appears to be part of the Bioinformatics Toolbox
per isakson
per isakson el 3 de Oct. de 2015
Editada: per isakson el 5 de Oct. de 2015
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?

Iniciar sesión para comentar.

Respuesta aceptada

Paola Favaretto
Paola Favaretto el 5 de Oct. de 2015
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

Más respuestas (0)

Categorías

Más información sobre Genomics and Next Generation Sequencing en Help Center y File Exchange.

Etiquetas

Community Treasure Hunt

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

Start Hunting!

Translated by