aboutsummaryrefslogtreecommitdiff
path: root/src/score_mers.py
diff options
context:
space:
mode:
authorCalvin Morrison <mutantturkey@gmail.com>2014-03-03 14:59:42 -0500
committerCalvin Morrison <mutantturkey@gmail.com>2014-03-03 14:59:42 -0500
commitecb8a85eb5cf25e0aba4a3d2d7abf94f63f17297 (patch)
tree241afb5a58e6004fa6de7d9bc0c8c6072ab05b61 /src/score_mers.py
parent3206ee1ccdcfac65e6af7f0c9913854985b3ee58 (diff)
update defaults and align variable names with global env variables
Diffstat (limited to 'src/score_mers.py')
-rwxr-xr-xsrc/score_mers.py14
1 files changed, 7 insertions, 7 deletions
diff --git a/src/score_mers.py b/src/score_mers.py
index 1a54ba1..1c9d4ae 100755
--- a/src/score_mers.py
+++ b/src/score_mers.py
@@ -35,17 +35,17 @@ class Mer:
# import our variables
cpus = int(os.environ.get("cpus", cpu_count()));
min_mer_range = int(os.environ.get("min_mer_range", 6));
-max_mer_range = int(os.environ.get("max_mer_range", 10));
+max_mer_range = int(os.environ.get("max_mer_range", 12));
min_mer_count = int(os.environ.get("min_mer_count", 0));
max_select = int(os.environ.get("max_select", 15));
-max_check = int(os.environ.get("max_check", 35));
+max_check = int(os.environ.get("max_check", 35));
max_mer_distance = int(os.environ.get("max_mer_distance", 5000));
-nb_max_consecutive_binding = int(os.environ.get("max_consecutive_binding", 4));
+max_consecutive_binding = int(os.environ.get("max_consecutive_binding", 4));
binding = { 'A': 'T', 'T': 'A', 'C': 'G', 'G': 'C', '_': False }
-def max_consecutive_binding(mer1, mer2):
+def get_max_consecutive_binding(mer1, mer2):
if len(mer2) > len(mer1):
mer1, mer2 = mer2, mer1
@@ -198,9 +198,9 @@ def pop_bg(mer):
def load_heterodimer_dic(selected_mers):
for (mer1, mer2) in combinations(selected_mers, 2):
- res = max_consecutive_binding(mer1, mer2)
- heterodimer_dic[(mer1, mer2)] = res > nb_max_consecutive_binding
- heterodimer_dic[(mer2, mer1)] = res > nb_max_consecutive_binding
+ res = get_max_consecutive_binding(mer1, mer2)
+ heterodimer_dic[(mer1, mer2)] = res > max_consecutive_binding
+ heterodimer_dic[(mer2, mer1)] = res > max_consecutive_binding
# print res, heterodimer_dic[(mer1, mer2)]
def main():