Are you blogging on PH? Get your free blog.

Matlab

Moderators: None (Apply to moderate this forum)
Number of threads: 630
Number of posts: 1072

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



 
Official Programmer's Heaven Blogs
Web Hosting | Browser and Social Games | Gadgets

Popular resources on Programmersheaven.com
Assembly | Basic | C | C# | C++ | Delphi | Java | JavaScript | Pascal | Perl | PHP | Python | Ruby | Visual Basic
© Copyright 2010 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.
bootstrapLabs Logo A BootstrapLabs project.