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

paidopoieo
 2005-12-02 02:27
#60680 #60680
User since
2005-12-02
96 Artikel
BenutzerIn
[default_avatar]
Hi,
Ich bin ziemlich verzweifelt, da ich seit zwei Tagen versuche meinen Output in einem file zu speichern, das file wird zwar angelegt, ist jedoch leer...Blast ist ein Bioinformatiktool, das ich uebers Netz anspreche.....da geht auch alles nur die files sind alle leer.....
hier der code......


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
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
#!/usr/bin/perl -w

use strict;
use warnings;
use Bio::SeqIO;
use Bio::Tools::Run::RemoteBlast;
use Bio::Seq;
use IO::String;
use Bio::SearchIO;



my $prog = 'blastp';
my $db = 'swissprot';
my $e_val= '20000';
my $matrix = 'PAM30';
#my $outfile = 'Output';

my @data;
my $line_dataArray;
my $rid;
my $count = 1;


my @params = ( '-prog' => $prog,
'-data' => $db,
'-expect' => $e_val,
'-matrix' => $matrix);


my $seqio_obj = Bio::SeqIO->new(-file => "Perm.txt",
-format => "raw",
);


print "entering blast....";
my $factory = Bio::Tools::Run::RemoteBlast->new (@params);
print "Blast entered successfully \n";




while( my $query = $seqio_obj->next_seq ) {
print "submit Sequence...just do it....\n";
my $r = $factory->submit_blast($query);
print $query->seq;
print "\n";

# Wait for the reply and save the output file
print "entering while loop for saving Output.... \n";
while ( my @rids = $factory->each_rid ) {
foreach my $rid ( @rids ) {
my $rc = $factory->retrieve_blast($rid);
print "retrieved Results successfully \n";
if( !ref($rc) ) {
if( $rc < 0 ) {
$factory->remove_rid($rid);
}
sleep 5;
} else {
#my $result = $rc->next_result;
print $rid;
print "\n";
my $filename = " $count .out";
$factory->save_output($filename);
print "File saved successfully \n";
$count++;
$factory->remove_rid($rid);
# }
}
}
print "\n";
print "\n";


}
}


Ich hoffe es kann mir jemand helfen, dankeschoen im voraus...

mfg
Hubert

View full thread Bioperl: remoteblast speichert Output file nicht