From d715f6202073ec7481f75a9218511720f2c4d82e Mon Sep 17 00:00:00 2001
From: davidhcsic <davidhenriques@Davids-MacBook-Pro.local>
Date: Sat, 24 Aug 2024 12:21:26 +0200
Subject: [PATCH] Updated code for flux predictions in anaerobic conditions.
 Loads txt isntead of xlsx.

---
 code/modelTests/anaerobic_flux_predictions.m | 15 ++++++++-------
 1 file changed, 8 insertions(+), 7 deletions(-)

diff --git a/code/modelTests/anaerobic_flux_predictions.m b/code/modelTests/anaerobic_flux_predictions.m
index 53f2fb12..f68d93b1 100644
--- a/code/modelTests/anaerobic_flux_predictions.m
+++ b/code/modelTests/anaerobic_flux_predictions.m
@@ -1,9 +1,9 @@
 function R2=anaerobic_flux_predictions(model)
 
-[vals_flux text_flux]=xlsread('../../data/physiology/anaerobic_flux.xlsx','flux_data');
-
-% %% This is temporary; remove when rxn id is attributed
-% text_flux(strcmp(text_flux(:,7),'rxxS'),7)={'r_1239'};
+fluxTable=  readtable('../../data/physiology/flux_data_anaerobic.txt');
+fluxTable = table2cell(fluxTable);
+vals_flux=fluxTable;
+text_flux=fluxTable;
 
 sim_vals=[];
 
@@ -18,9 +18,10 @@
 
 for i=1:length(data_sets)
 
-    index_data_set=find(strcmp(text_flux(:,6),data_sets(i)));
+    index_data_set=find(strcmp(text_flux(:,6),data_sets(i,1)));
     index_glx=strcmp(text_flux(index_data_set,7),'r_1714');
-    model = setParam(model,'eq','r_1714',-vals_flux(index_data_set(index_glx),2) );
+    
+    model = setParam(model,'eq','r_1714',-cell2mat(vals_flux(index_data_set(find(index_glx)),4) ));
 
     %% Solve the LP problem
     res=solveLP(model,1);
@@ -32,7 +33,7 @@
 
     scaled_sim=abs(-100.*res.x(index_model)./res.x(findRxnIDs(model,'r_1714')));
 
-    data_vals=abs(vals_flux(index_data_set,3));
+    data_vals=abs(cell2mat(vals_flux(index_data_set,5)));
     merged_data=[merged_data data_vals'];
     merged_sim=[merged_sim scaled_sim'];
     merged_names=[merged_names text_flux(index_data_set,7)'];