[Bioperl-l] can't start 'water'
Jason Stajich
jason at cgt.duhs.duke.edu
Sat Sep 6 11:03:41 EDT 2003
do you have EMBOSS installed?
On Sat, 6 Sep 2003, Vesko Baev wrote:
> Hi,
> In my computer I have bioperl-pack & Bioperl-run-pack, and I HAVE water.pm,but:
> message in dos-prompt:
>
> -------------------- WARNING ---------------------
> MSG: Application [water] is not available!
> ---------------------------------------------------
> Can't call method "run" on an undefined value at water.pl line 22.
>
> Here is my script:
> #!usr/bin/perl
> use Bio::Factory::EMBOSS;
> use Bio::AlignIO;
> my $seqobj1 = Bio::Seq->new(-seq => "gacgggtatctttaggcggacttagg");
> my $seqobj2 = Bio::Seq->new(-seq => "gacgtgtatctttaggcggtcttacc");
> my @seqs_to_check; # this would be a list of seqs to compare
> push (@seqs_to_check,seqobj2);
> # get an EMBOSS application object from the EMBOSS factory
> $factory = new Bio::Factory::EMBOSS;
> # $application = $factory->program('embossversion');
> # run the application with an optional hash containing parameters
> # $result = $application->run(); # returns a string or creates a file
> # print $result . "\n";
>
> $water = $factory->program('water');
>
> # here is an example of running the application
> # water can compare 1 seq against 1->many sequences
> # in a database using Smith-Waterman
>
> my $wateroutfile = 'mirna.water';
> $water->run({ '-sequencea' => $seqobj1,
>
> '-seqall' => \@seqs_to_check,
> '-gapopen' => '10.0',
> '-gapextend' => '0.5',
> '-outfile' => $wateroutfile});
> # now you might want to get the alignment
> use Bio::AlignIO;
> my $alnin = new Bio::AlignIO(-format => 'emboss',
> -file => $wateroutfile);
>
> while( my $aln = $alnin->next_aln ) {
> # process the alignment -- these will be Bio::SimpleAlign objects
> }
>
> AND when I wrote 'use Bio::Tools::Run::PiseApplication::water;'
> the message was:
> Can't locate XML/Parser/PerlSAX.pm in @INC?!?!?!
>
> WHAT TO DO?
> Vesko
> Thenks!!
>
>
>
> -----------------------------------------------------------------
> http://nova.GBG.bg - Направи си сайт! Лесно, бързо и професионално.
> _______________________________________________
> Bioperl-l mailing list
> Bioperl-l at portal.open-bio.org
> http://portal.open-bio.org/mailman/listinfo/bioperl-l
>
--
Jason Stajich
Duke University
jason at cgt.mc.duke.edu
More information about the Bioperl-l
mailing list