load('model.mat');
changeCobraSolver('ibm_cplex', 'all');
%% Select reactions to be explorer
[transRxns,nonTransRxns] = findTransRxns(model,true);
[tmp,ATPMnumber] = ismember(rxn_atpm, nonTransRxns);
[tmp,BioMassnumber] = ismember(rxn_biomass, nonTransRxns);
nonTransRxnsLength = length(nonTransRxns);
%% Candidate reactions to be knocked out
selectedRxns = {nonTransRxns{[1:ATPMnumber-1,ATPMnumber+1:BioMassnumber-1,BioMassnumber+1:nonTransRxnsLength]}}:
%% Optnock options
options.targetRxn = rxn_target;
options.vMax = 1000;
options.numDel = 3;
options.numDelSense = 'L';
constrOpt.rxnList = {rxn_biomass};
constrOpt.values = [0.2];
constrOpt.sense = 'GE';
%% Launch Optknock
optKnockSol = OptKnock(model, selectedRxns, options, constrOpt);
>> optKnockSol = OptKnock(model, selectedRxns, options, constrOpt);
0
Error using cplexlink1280
CPLEX Error 1003: Bad argument to Callable Library routine.
Error in Cplex/subsasgn
Error in solveCobraMILP (line 326)
cplexlp.Param.output.writelevel.Cur = cobraParams.printLevel;
Error in OptKnock (line 207)
optKnockSol = solveCobraMILP(bilevelMILPproblem,'printLevel',verbFlag);
>> cplexlp.Param.output.writelevel.Cur = cobraParams.printLevel
ans=
Cplex handle
Properties:
Model: [1x1 struct]
Param: [1x1 struct]
DisplayFunc: @disp
--
---
You received this message because you are subscribed to the Google Groups "COBRA Toolbox" group.
To unsubscribe from this group and stop receiving emails from it, send an email to cobra-toolbo...@googlegroups.com.
For more options, visit https://groups.google.com/d/optout.