fix: calculateFfactor can take protData input

SysBioChalmers · May 26, 2024 · 5f3c14d · 5f3c14d
1 parent 5e9ca54
commit 5f3c14d
Showing 1 changed file with 8 additions and 9 deletions.
diff --git a/src/geckomat/limit_proteins/calculateFfactor.m b/src/geckomat/limit_proteins/calculateFfactor.m
@@ -61,20 +61,19 @@
     uniprot     = uniprotDB.ID(b(a));
     level(~a)   = [];
     clear protData
-    protData.uniprot = uniprot;
+    protData.uniprotIDs = uniprot;
     protData.level   = level;
+    % Get MW and abundance (unit does not matter, f is fraction)
+    [~,idx] = ismember(protData.uniprotIDs,uniprotDB.ID);
+    protData.MW = uniprotDB.MW(idx);
+    protData.abundances = protData.level .* protData.MW;
 end
 
-% Get MW and abundance (unit does not matter, f is fraction)
-[~,idx] = ismember(protData.uniprot,uniprotDB.ID);
-protData.MW = uniprotDB.MW(idx);
-protData.abundance = protData.level .* protData.MW;
-
-totalProt = sum(protData.abundance);
+totalProt = sum(protData.abundances);
 
 % Get enzymes in model
-enzymesInModel = ismember(protData.uniprot,enzymes);
-totalEnz = sum(protData.abundance(enzymesInModel));
+enzymesInModel = ismember(protData.uniprotIDs,enzymes);
+totalEnz = sum(protData.abundances(enzymesInModel));
 
 f = totalEnz/totalProt;
 end