Archive

Ma version du ¼

crouet64pau

Bonjour,
après vérifié manuellement qu’il n’existait pas de caractère de terminaison de la chaine (pas de #), voici ma version :
def count_gc_at(dna):

renvoie une fois pour toutes la longueur de la chaine

totalcar=len(dna)

#initialisation des variables
nb_GC=nb_AT=0
#lancement de la boucle
for i in range (0,len(dna)):

on incrémente que les génes pairés

        if dna[i]=="A"or dna[i]=="T": nb_AT+=1
        elif dna[i]=="G"or dna[i]=="C": nb_GC+=1    
#return "votre code"
# votre code

#calcul du ratio
GC=nb_GC/totalcar
AT=nb_AT/totalcar
#on renvoie les valeurs à la suite du script
return(GC,AT)


Je voudrais juste me permettre de préciser que demander les GC et AT n’est pas identique dans la formulation (selon moi) aux ratios de G+C et A+T dans la chaîne, je veux dire qu’on peut aussi comprendre qu’on cherche le nombre de fois où les bi-gram “AT” et "CG " sont trouvés dans la séquence.

boulbi777

Au début de l’exercice, il me semble qu’on a bien spécifié que l’exercice fait référence aux éléments vus dans la séquence 6 de la semaine 1. Partant de là, je pense qu’il n’y a pas trop d’ambiguité, puisqu’il n’y a pas de notion de bigramme dans cette dite séquence…

Sinon, j’apprécie ton algorithme, il faudrait juste éviter de déclarer beaucoup de variables si ce n’est pas nécessaire :wink: