opencobra / cobratoolbox

The COnstraint-Based Reconstruction and Analysis Toolbox. Documentation:
https://opencobra.github.io/cobratoolbox
Other
240 stars 305 forks source link

`optimizeCbModel` returning `NaN` for solution #2255

Open nleroy917 opened 2 months ago

nleroy917 commented 2 months ago

Hello!

I'm having issues running singleGeneDeletion. Speficially, the call to the WT optimizeCbModel returns a solution with NaN, consequently leading to all my growth rate ratios also being NaN. I've believe narrowed it down to line 119 inside singleGeneDeletion.

I've confirmed that changing the line from:

solWT = optimizeCbModel(model,'max','one'); % by default uses the min manhattan distance norm FBA solution.

To:

solWT = optimizeCbModel(model,'max'); % by default uses the min manhattan distance norm FBA solution.

Seems to give something other than NaN, but I don't have the expertise to figure out a.) why this works, and b.) if the solution is still valid.

Any help or insight is appreciated!

I hereby confirm that I have:

System specs: