-
Notifications
You must be signed in to change notification settings - Fork 16
Expand file tree
/
Copy pathdataset_gen.py
More file actions
55 lines (43 loc) · 1.59 KB
/
dataset_gen.py
File metadata and controls
55 lines (43 loc) · 1.59 KB
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
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
import numpy as np
import os
import random
import argparse
import pandas as pd
import statistics
import csv
parser = argparse.ArgumentParser()
parser.add_argument("--iteration_num", type=int, required=True)
parser.add_argument("--label", type=str, required=True)
parser.add_argument("--model_dir", type=str, required=True)
args = parser.parse_args()
def append_to_csv(name, sequence, iteration_num, output_file):
file_exists = os.path.exists(output_file) and os.stat(output_file).st_size > 0
with open(output_file, "a", newline="") as csvfile:
fieldnames = [
"name",
"sequence",
"lenght",
"iteration_num",
]
writer = csv.DictWriter(csvfile, fieldnames=fieldnames)
if not file_exists:
writer.writeheader()
writer.writerow({
"name" : name,
"sequence" : sequence,
"lenght" : len(sequence),
"iteration_num": iteration_num,
})
def generate_dataset(iteration_num, ec_label):
output_file = "logs.csv"
with open(f"seq_gen_{ec_label}_iteration{iteration_num}.fasta", "r") as f:
rep_seq = f.readlines()
sequences_rep = dict()
for line in rep_seq:
if ">" in line:
name = line.split("\t")[0].replace(">", "").replace("\n", "")
emb_identifier = line.replace(">", "").replace("\n", "")
else:
sequence = line.strip()
append_to_csv(name, sequence, iteration_num, output_file)
generate_dataset(args.iteration_num, args.label)