Not a member of Pastebin yet?
Sign Up,
it unlocks many cool features!
- def list_locus(dico): #Trie les locus. Une liste pour tous les locus présents, une pour les locus à 1 transcrit
- #Et une pour les locus à plusieurs transcrits
- li=[]
- for i in dico:
- li.append(i["locus_name"])
- list_all_locus=sorted(list(set(li))) #set(liste) renvoie la même liste débarassée des redondances
- locus_one_transcript=[]
- locus_multi_transcript=[]
- for i in list_all_locus:
- a=[];b=[]
- for dico in data:
- if dico["locus_name"]==i:
- a.append(dico["transcript_name"]) #Stocke le nom de tous les transcrits pour un locus
- b=list(set(a))
- if len(b)==1: #Liste des noms non-redondante. Si il y a un élément, c'est un locus à un seul transcrit
- locus_one_transcript.append(i)
- else:
- locus_multi_transcript.append(i)
- if len(locus_one_transcript)+len(locus_multi_transcript)!=len(list_all_locus):
- print "Erreur dans la gestion des locus"
- print len(locus_one_transcript)
- print len(locus_multi_transcript)
Add Comment
Please, Sign In to add comment