Précédent   Forum des professionnels en informatique > Autres langages > Perl > Bioinformatique
Bioinformatique Toutes vos questions sur les scripts Perl associés à la bioinformatique, modules bioperl, projets bioinformatiques, etc ... Avant de poster, veuillez consulter les cours Perl et les critiques de livres.
Partagez cette discussion sur d'autres réseaux sociaux : Viadeo Twitter Google Facebook Digg Delicious MySpace Yahoo
Réponse Proposer ce sujet en actualité
 
Outils de la discussion
Publicité
'
Vieux 23/11/2011, 11h13   #1
Nouveau Membre du Club
 
Inscription : mars 2011
Messages : 101
Détails du profil
Informations forums :
Inscription : mars 2011
Messages : 101
Points : 33
Points : 33
Par défaut Undefined subroutine &Bio::SearchIO called

Bonjour,

lorsque je lance mon programme, j'ai cette erreur et je ne sais pas comment la résoudre.

j'aimerai lire le fichier d'un blast.

Merci

Alaninho
alaninho est déconnecté   Envoyer un message privé Réponse avec citation 00
Vieux 23/11/2011, 11h19   #2
Membre éprouvé
 
Avatar de Beniou
 
Homme
Inscription : novembre 2009
Messages : 348
Détails du profil
Informations personnelles :
Sexe : Homme
Âge : 32
Localisation : France, Nord (Nord Pas de Calais)

Informations forums :
Inscription : novembre 2009
Messages : 348
Points : 492
Points : 492
Bonjour,

Sans code à l'appui qui donne cette erreur, cela va être difficile de répondre à ton problème...
Beniou est déconnecté   Envoyer un message privé Réponse avec citation 00
Vieux 23/11/2011, 11h24   #3
Nouveau Membre du Club
 
Inscription : mars 2011
Messages : 101
Détails du profil
Informations forums :
Inscription : mars 2011
Messages : 101
Points : 33
Points : 33
désolé, je pensais l'avoir mis, le voici :
Code :
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
 
#!/usr/bin/perl -w
 
use Bio::SearchIO;
 
 
sub filtre_blast {
	$in=Bio::SearchIO(-format=>'blast', -file=>"out_blast.txt");
	while($result=$in->next_result){
		while($hit=$result->next_hit){
			while($hsp=$hit->next_hsp){
				next if($hsp->length('total')>100);
				next if($hsp->percent_identity>80);
				next if($hsp->evalue<0.05);
				print $result->query_name."\t".$hsp->start('query')."\t".$hsp->end('query')."\n";
			}
		}
	}
}
 
filtre_blast();
alaninho est déconnecté   Envoyer un message privé Réponse avec citation 00
Vieux 23/11/2011, 15h35   #4
Membre éprouvé
 
Avatar de Beniou
 
Homme
Inscription : novembre 2009
Messages : 348
Détails du profil
Informations personnelles :
Sexe : Homme
Âge : 32
Localisation : France, Nord (Nord Pas de Calais)

Informations forums :
Inscription : novembre 2009
Messages : 348
Points : 492
Points : 492
Bio::SearchIO est le module.

Pour lire ton fasta tu dois utiliser la méthode new de ce module:
Code :
1
2
 
$in=Bio::SearchIO->new(-format=>'blast', -file=>"out_blast.txt");
Beniou est déconnecté   Envoyer un message privé Réponse avec citation 00
Réponse Proposer ce sujet en actualité Cette discussion est résolue.
Outils de la discussion



Fuseau horaire GMT +2. Il est actuellement 14h36.


 
 
 
 
Partenaires

Hébergement Web