@@ -42,8 +42,8 @@ def synthetic_peptides_by_subsequence(
42
42
if not isinstance (lengths , dict ):
43
43
lengths = dict ((length , 1.0 ) for length in lengths )
44
44
45
- lengths = pandas .Series (lengths )
46
- lengths /= len (lengths )
45
+ lengths_series = pandas .Series (lengths )
46
+ lengths_series /= len (lengths )
47
47
48
48
num_binders = int (round (num_peptides * fraction_binders ))
49
49
num_non_binders = num_peptides - num_binders
@@ -52,15 +52,20 @@ def synthetic_peptides_by_subsequence(
52
52
peptides = []
53
53
54
54
# Generate non-binders
55
- for (length , weight ) in lengths .items ():
55
+ for (length , weight ) in lengths_series .items ():
56
56
peptides .extend (
57
57
random_peptides (round (weight * num_non_binders ), round (length )))
58
58
59
59
for binding_core in binding_subsequences :
60
60
# Generate binders
61
- lengths_binders = lengths .ix [lengths .index >= len (binding_core )]
61
+ lengths_binders = lengths_series .ix [
62
+ lengths_series .index >= len (binding_core )
63
+ ]
62
64
normalized_lengths_binders = (
63
- lengths_binders / lengths_binders .sum () / len (binding_subsequences ))
65
+ lengths_binders /
66
+ lengths_binders .sum () /
67
+ len (binding_subsequences ))
68
+
64
69
for (length , weight ) in normalized_lengths_binders .items ():
65
70
if length >= len (binding_core ):
66
71
num_peptides_to_make = int (round (weight * num_binders ))
0 commit comments