If you have a PH account, you can customize your PH profile.

Matlab

Moderators: None (Apply to moderate this forum)
Number of threads: 597
Number of posts: 1018

This Forum Only
Post New Thread

Report
bioinformatics toolbox- command from turoial won't work Posted by markmh on 13 Nov 2009 at 2:26 PM
I am trying to do an HMM alignment, and I am talking the commands straight out of the tutorial, and I am getting errors. It is funny, because I can do other tutorials, its not like the whole thing is broken.
This is what happens:

hmmmodel = gethmmprof('PF00002')

is the command.

error:
??? Error using ==> pfamhmmread at 70
Input is not a valid PFAM file.

Error in ==> getsangerdata at 180
out = pfamhmmread(s);

Error in ==> gethmmprof at 82
model = getsangerdata(accessnum,'database','hmm',varargin{:});

Does anyone know what is going on?
Thanks.
Mark
Report
Re: bioinformatics toolbox- command from turoial won't work Posted by miac on 7 Dec 2009 at 4:45 PM
Hi Mark:

Try this instead:

<for example Pfam Id= PF01596>

seed_seqs = gethmmalignment(1596,'type','seed');
seqdisp(seed_seqs(),'row',80)

-Mia
Report
Re: bioinformatics toolbox- command from turoial won't work Posted by miac on 7 Dec 2009 at 4:46 PM
Hi Mark:

<for example Pfam Id= PF01596>
Try this instead:

seed_seqs = gethmmalignment(1596,'type','seed');
seqdisp(seed_seqs(),'row',80)

-Mia



 
Popular resources and forums for programmers on Programmersheaven.com
Assembly, Basic, C, C#, C++, Delphi, Java, JavaScript, Pascal, Perl, PHP, Python, Ruby, Visual Basic
© Copyright 2009 Programmersheaven.com - All rights reserved.
Reproduction in whole or in part, in any form or medium without express written permission is prohibited.
Violators of this policy may be subject to legal action. Please read our Terms Of Use and Privacy Statement for more information.
Publisher: Lars Hagelin. Read the latest words from the publisher here.
Be the first to sign up for Lars Hagelin’s In-depth Outsourcing Newsletter here.
bootstrapLabs Logo A BootstrapLabs project.