1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36
| #MATRICE DE COMPTAGE
#on lit toutes les lignes du fichier de references et on les met en clé dans un dictionnaire
#avec comme valeurs [0,0,0,0]
import numpy as np
def matricep (repertoire, reference):
ref=open(reference, "r")
m={'Sequence': ['a','b','c','d']}
def traitement(line):
return(line.strip().lstrip(">"))
for line in ref:
line = traitement(line)
m[line] = [0 for i in repertoire]
i=0
for fichier in repertoire:
with open(fichier, "r") as f:
j=0
for line in f:
line = traitement(line)
if line in m.keys():
j+=1
m[line][i]= j
i+=1
a=(np.array(list(m.items())))
print(a)
np.savetxt(r'C:\Users\33652\Desktop\out.txt', a, fmt='%s')
matricep([r'C:\Users\33652\Desktop\a.txt',r'C:\Users\33652\Desktop\b.txt',r'C:\Users\33652\Desktop\c.txt',r'C:\Users\33652\Desktop\d.txt'],r'C:\Users\33652\Desktop\ref.txt') |
Partager