1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
|
Author: Laszlo Kajan <lkajan@rostlab.org>
Description: move code outside of loop
Forwarded: http://lists.alioth.debian.org/pipermail/debian-med-packaging/2012-October/017448.html
--- a/score_conservation.py
+++ b/score_conservation.py
@@ -136,8 +136,9 @@
aa_num += 1
+ freqsum = (sum(seq_weights) + len(amino_acids) * pc_amount)
for j in range(len(freq_counts)):
- freq_counts[j] = freq_counts[j] / (sum(seq_weights) + len(amino_acids) * pc_amount)
+ freq_counts[j] = freq_counts[j] / freqsum
return freq_counts
|
