diff --git a/feature_calculation/protein.py b/feature_calculation/protein.py index f31965e..b9d5363 100755 --- a/feature_calculation/protein.py +++ b/feature_calculation/protein.py @@ -83,12 +83,12 @@ def percent_thermolabile(AA_counts): def EKQH(AA_counts): #calculate the ratio of EK/QH in the proteome EK = 0.0 - QT=0.0 + QH=0.0 for x in ['E','K']: EK=EK+AA_counts[x] - for x in ['Q','T']: - QT=QT+AA_counts[x] - return EK/QT + for x in ['Q','H']: + QH=QH+AA_counts[x] + return EK/QH def EFMR(AA_counts): #calculate the fraction of EFMR in the proteome @@ -124,8 +124,6 @@ def percent_ERK(AA_counts): counts= sum(AA_counts.values()) for x in ['E','R','K']: ERK=ERK+AA_counts[x] - for x in AAs: - counts=counts+AA_counts[x] return ERK/counts def LKQ(AA_counts):