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 :

Coverage, Blast, et BioPerl


Sujet :

Bioinformatique Perl

Vue hybride

Message précédent Message précédent   Message suivant Message suivant
  1. #1
    Membre averti
    Femme Profil pro
    Etudiante
    Inscrit en
    Avril 2011
    Messages
    10
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Localisation : France

    Informations professionnelles :
    Activité : Etudiante

    Informations forums :
    Inscription : Avril 2011
    Messages : 10
    Par défaut Coverage, Blast, et BioPerl
    Bonjour les développeurs !


    Alors voilà, je suis actuellement stagiaire dans un laboratoire suite à mes deux années en IUT Génie Bio (opt Bioinfo). Je dois faire un programme capable d'utiliser un certain nombre d'informations depuis un tblastn en local. Je travaille donc beaucoup avec les modules BioPerl (Bio::SeqIO, Bio::Seq, Bio::SearchIO), mais n'ayant jamais étudié ces modules ou même le Perl d'un point de vue POO, il y a quelques points où j'ai du mal, notamment celui qui suit, sur lequel je m'arrache les cheveux depuis quelques jours :

    Je veux récupérer le taux de couverture (coverage) entre ma séquence prot et ma séquence nucléique trouvée. J'ai bien trouvé quelques fonctions (frac_aligned_hit(), frac_aligned_query() ) mais il ne semble pas que ce soit ça, où alors je les utilise mal, car je n'obtiens jamais de résultats interessants (<10% en général). De plus, si c'est bien l'une de ces deux fonction, je n'ai absolument aucune idée de la bonne...


    Pensez vous pouvoir m'aider? Je vous remercie d'ors et déjà pour l'attention que vous porterez à ma question.

  2. #2
    Membre averti
    Femme Profil pro
    Etudiante
    Inscrit en
    Avril 2011
    Messages
    10
    Détails du profil
    Informations personnelles :
    Sexe : Femme
    Localisation : France

    Informations professionnelles :
    Activité : Etudiante

    Informations forums :
    Inscription : Avril 2011
    Messages : 10
    Par défaut
    Bon, alternative trouvée :


    Code : Sélectionner tout - Visualiser dans une fenêtre à part
    1
    2
    3
    4
    5
     
    $cov = ($hsp->end('sbjct')-$hsp->start('sbjct')+1)/3;
     
    	if ($cov < ($cover * $hit->logical_length('query'))/100) {return 0;}
    	else {return 1;}
    Woilà.

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

Discussions similaires

  1. Blast via Bioperl => erreur
    Par Kawaccino dans le forum Bioinformatique
    Réponses: 3
    Dernier message: 13/03/2013, 10h51
  2. problème pour lancer un blast avec bioperl
    Par alaninho dans le forum Bioinformatique
    Réponses: 3
    Dernier message: 24/11/2011, 14h59
  3. erreur code, Blast Bioperl
    Par Bioinfoe dans le forum Bioinformatique
    Réponses: 4
    Dernier message: 14/10/2009, 16h59
  4. Blast avec Bioperl
    Par Petitchaoui dans le forum Bioinformatique
    Réponses: 21
    Dernier message: 12/01/2009, 08h03
  5. [bioperl]parser un blast xml
    Par sohnic dans le forum Bioinformatique
    Réponses: 1
    Dernier message: 24/04/2007, 13h31

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