[Bioperl-l] can't start 'water'
Vesko Baev
vesko_baev at abv.bg
Sat Sep 6 08:29:26 EDT 2003
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 - Направи си сайт! Лесно, бързо и професионално.
More information about the Bioperl-l
mailing list