Thread Bioperl: remoteblast speichert Output file nicht (18 answers)
Opened by paidopoieo at 2005-12-02 02:27

paidopoieo
 2005-12-02 03:35
#60682 #60682
User since
2005-12-02
96 Artikel
BenutzerIn
[default_avatar]
Hallo esskar,

Title : saveoutput
Usage : my $saveoutput = $self->save_output($filename)
Function: Method to save the blast report
Returns : 1 (throws error otherwise)
Args : string [rid, filename]

das waere der beitrag zur methode save_output

und wenn ich mein programm starte dann laeuft es einwandfrei durch, speichert auch das file, nur leider ist es leer.....

das der Code zur methode save_output

Code: (dl )
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
sub save_output  {

my ($self, $filename) = @_;
if( ! defined $filename ) {
$self->throw("Can't save blast output. You must specify a filename to save to.");
}
#should be set when retrieving blast
my $blastfile = $self->file;
#open temp file and output file, have to filter out some HTML
open(TMP, $blastfile) or $self->throw("cannot open $blastfile");
open(SAVEOUT, ">$filename") or $self->throw("cannot open $filename");
my $seentop=0;
while(<TMP>) {
next if (/<pre>/);
if( /^(?:[T]?BLAST[NPX])\s*.+$/i ||
/^RPS-BLAST\s*.+$/i ) {
$seentop=1;
}
next if !$seentop;
if( $seentop ) {
print SAVEOUT;
}
}
close SAVEOUT;
close TMP;
return 1;

}


der output meines programms ist folgender:

entering blast....Blast entered successfully
submit Sequence...just do it....
RKKWWWRW
Sequence submitted successfully
entering while loop for saving Output....
retrieved Results successfully
1133487045-22411-147042628522.BLASTQ1
File saved successfully


submit Sequence...just do it....
RKKWWWWK
Sequence submitted successfully
entering while loop for saving Output....
retrieved Results successfully
1133487046-1295-107621241415.BLASTQ4
File saved successfully


vielen Dank fuer deine Bemuehungen
Dankeschoen

View full thread Bioperl: remoteblast speichert Output file nicht