The question about function relaxModelWithDGoSlacks

[pwzation]

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

[agmpo]

Hello! I also have this same problem...though it is via testing the scripts associated with the i-At-LSPHERE repository on the model provided there (i.e. Leaf1). When it gets to the TFAmodel fitting step, it states:

It is possible to obtain TFA-feasible solutions with the original DGo values calculated by GCM!
 
Gapfilling...Unrecognized function or variable 'GFModelMaster'.

Error in gapFillModelTFA (line 273)
[ActRxnsAll,foundSolution] = performGapfillingTFA(nutrients,growthNutrients,GFModelMaster,excRxns,minMed,vitamins,modelOrig,maxNumAlt);

Which if you try to figure out what went wrong by running the line:
[GFModelMaster, conflict] = PrepareForGapFilling(sourceModel,{modelDB},'',0,flagTFA,{},[],thermoData);

I get:

Error using relaxModelWithDGoSlacks (line 175)
retry

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

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

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

So it seems to me like the relaxModelWithDGoSlacks script is not spitting out something that it should in this case. I am new to both matlab and metabolic modeling and am trying to learn more about them using these scripts, and any help would be deeply appreciated. Thanks!