IdentifiantMot de passe
Loading...
Mot de passe oublié ?Je m'inscris ! (gratuit)
Navigation

Inscrivez-vous gratuitement
pour pouvoir participer, suivre les réponses en temps réel, voter pour les messages, poser vos propres questions et recevoir la newsletter

Bioinformatique Perl Discussion :

problème pour lancer un blast avec bioperl


Sujet :

Bioinformatique Perl

Vue hybride

Message précédent Message précédent   Message suivant Message suivant
  1. #1
    Membre éclairé
    Profil pro
    Inscrit en
    Mars 2011
    Messages
    299
    Détails du profil
    Informations personnelles :
    Localisation : France

    Informations forums :
    Inscription : Mars 2011
    Messages : 299
    Par défaut problème pour lancer un blast avec bioperl
    Bonjour,

    j'aimerai lancer un blast avec bioperl avec le code suivant :
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
     
    #!/usr/bin/perl -w
     
     
    use Bio::Seq;
    use Bio::Tools::Run::StandAloneBlast;
    #use strict;
     
    my($blast_obj) = Bio::Tools::Run::StandAloneBlast->new(-program=>'blastn', -database=>'At.fa');
    my($seq_obj)= Bio::Seq->new(-id => "test_query", -seq=>"AATATATGCGCGAT");
     
    my($blast_report) = $blast_obj->blastall($seq_obj);
    mais j'obtiens cette erreur :
    Use of uninitialized value $_[0] in join or string at /usr/lib/perl/5.12/File/Spec/Unix.pm line 41.

    ------------- EXCEPTION: Bio::Root::Exception -------------
    MSG: blastall call crashed: 256 | Aucun fichier ou dossier de ce type | /usr/bin/blastall -d "/At.fa" -i /tmp/4P9BDVACuV -p blastn -o \/tmp\/cV_T77dD_7

    STACK: Error::throw
    STACK: Bio::Root::Root::throw /usr/share/perl5/Bio/Root/Root.pm:472
    STACK: Bio::Tools::Run::StandAloneNCBIBlast::_runblast /usr/share/perl5/Bio/Tools/Run/StandAloneNCBIBlast.pm:437
    STACK: Bio::Tools::Run::StandAloneNCBIBlast::_generic_local_blast /usr/share/perl5/Bio/Tools/Run/StandAloneNCBIBlast.pm:411
    STACK: Bio::Tools::Run::StandAloneNCBIBlast::blastall /usr/share/perl5/Bio/Tools/Run/StandAloneNCBIBlast.pm:263
    STACK: ./B_4.pl:11
    ----------------------

    je cherche sur internet mais je ne trouve pas de bonne réponse.

    J'espère que vous pourrez m'aider.

    Alaninho

  2. #2
    Responsable Perl et Outils

    Avatar de djibril
    Homme Profil pro
    Inscrit en
    Avril 2004
    Messages
    19 822
    Détails du profil
    Informations personnelles :
    Sexe : Homme
    Localisation : France

    Informations forums :
    Inscription : Avril 2004
    Messages : 19 822
    Par défaut
    Déjà sans tester quoi que ce soit, commencez par mieux écrire votre code :
    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
    6
    7
    8
    9
    10
    11
    12
    13
    14
    15
    16
    17
    18
    19
    20
    21
    22
    23
    #!/usr/bin/perl
    #===============================================================================
    # Author : 
    # Date   :
    # Main   :
    #===============================================================================
    use Carp;
    use strict;
    use warnings;
     
    use Bio::Seq;
    use Bio::Tools::Run::StandAloneBlast;
     
    my $blast_obj = Bio::Tools::Run::StandAloneBlast->new(
      -program  => 'blastn', 
      -database => 'At.fa',
    );
    my $seq_obj = Bio::Seq->new(
      -id  => 'test_query',
      -seq => 'AATATATGCGCGAT',
    );
     
    my $blast_report = $blast_obj->blastall($seq_obj);

  3. #3
    Membre éclairé
    Profil pro
    Inscrit en
    Mars 2011
    Messages
    299
    Détails du profil
    Informations personnelles :
    Localisation : France

    Informations forums :
    Inscription : Mars 2011
    Messages : 299
    Par défaut
    Merci, mais ça marche toujours pas.

    Alaninho

  4. #4
    Membre averti
    Profil pro
    Inscrit en
    Juin 2010
    Messages
    22
    Détails du profil
    Informations personnelles :
    Localisation : France

    Informations forums :
    Inscription : Juin 2010
    Messages : 22
    Par défaut
    Vu le message d'erreur, il semblerait que blastall, l'exécutable de blast, ne soit pas dans le répertoire /usr/bin/

    Avant de faire l'appel à par l'intermédiaire d'un script, as-tu essayé d'exécuter blastall dans un terminal?

+ Répondre à la discussion
Cette discussion est résolue.

Discussions similaires

  1. [JDOM] Problème pour lire les noeuds avec JDOM
    Par Eclix dans le forum Format d'échange (XML, JSON...)
    Réponses: 7
    Dernier message: 26/05/2014, 18h11
  2. Problème pour lancer une application avec Java
    Par bintij dans le forum Débuter avec Java
    Réponses: 4
    Dernier message: 23/10/2011, 20h49
  3. [Système]Problème pour lancer l'appli au démarrage
    Par CVERRIER3 dans le forum API standards et tierces
    Réponses: 2
    Dernier message: 23/11/2005, 16h07
  4. problème pour faire un insert avec une image
    Par vbcasimir dans le forum Modules
    Réponses: 3
    Dernier message: 02/11/2005, 09h21
  5. Réponses: 1
    Dernier message: 17/05/2004, 11h29

Partager

Partager
  • Envoyer la discussion sur Viadeo
  • Envoyer la discussion sur Twitter
  • Envoyer la discussion sur Google
  • Envoyer la discussion sur Facebook
  • Envoyer la discussion sur Digg
  • Envoyer la discussion sur Delicious
  • Envoyer la discussion sur MySpace
  • Envoyer la discussion sur Yahoo