epfl-lcsb / mattfa Goto Github PK

Hello!
I'm now trying the matTFA on matlab. When I running the code: [GFModelMaster, conflict] = PrepareForGapFilling(sourceModel,{modelDB},'',0,flagTFA,{},[],thermoData); it gives me this error:
It is possible to obtain TFA-feasible solutions with the original DGo values calculated by GCM!
Wrong use of relaxModelWithDGoSlacks
retry

Error in convToTFA (line 646)
[model, relaxedDGoVarsValues, modelwDGoSlackVars] = relaxModelWithDGoSlacks(model, minObjSolVal, []);

Error in convToTFAStru4gf (line 29)
model = convToTFA(model, DBThermo, rxnsNoThermo,'DGo',[], 0.95*sol.f, [], [], 0);

Error PrepareForGapFilling (line 94)
model4 = convToTFAStru4gf(model3,DBThermo,flagTFA,rxnsNoThermo,metabData);

To be honest, I'm fresh on this function. So how can I check and maintain my DGo value?

Best,
Peng

I cannot open any of the models in the tutorials or the thermodynamic data file. The error is the same for all
% Not a binary MAT-file. Try load -ASCII to read as text.
and when I use ascii I get this error:
%Unknown text on line number 1 of ASCII file small_ecoli.mat "version".
It seems like I am missing something. could asist me on this?

Hi

I have run the matTFA for 'GSmodel_Ecoli.mat'. It fails because of infeasibility. I was wondering if the codes works fine on your systems? if not, please let me know what might be the cause of this error?

Thank you

Relaxation of DGo_MECDPS : 28.391 +- 2.256 (sigma) ---> 28.391 +- 72.256
Relaxation of DGo_REPHACCOAI : 16.774 +- 4.121 (sigma) ---> 16.774 +- 24.121
It was not possible to obtain TFA-feasible solutions with the original DGo values calculated by GCM!
We had to relax the following DGo ranges as follows:

ans =

'varName'           'LB-before'    'UB-before'    'LB-relaxed'    'UB-relaxed'
'DGo_MECDPS'        [  26.1357]    [  30.6469]    [   12.6021]    [   44.1805]
'DGo_REPHACCOAI'    [  12.6529]    [  20.8957]    [    8.5316]    [   25.0171]

Error using relaxModelWithDGoSlacks (line 171)
Although we relaxed the DGo variables it was not possible to find a feasible solution!

Error in convToTFA (line 645)
[model, relaxedDGoVarsValues, modelwDGoSlackVars] = relaxModelWithDGoSlacks(model, minObjSolVal,
rxnNameListNoDGoRelax);

Error in tutorial_GEM (line 53)
tmp = convToTFA(prepped_m,ReactionDB,[],'DGo');