I'm trying to find coefficients in theoretical function of uranium mining by minimizing distance between theretical function and actual values. And I'm getting the following:
NOTE: The problem has 201 variables (0 free, 0 fixed).
NOTE: The problem has 0 linear constraints (0 LE, 0 EQ, 0 GE, 0 range).
NOTE: The problem has 0 nonlinear constraints (0 LE, 0 EQ, 0 GE, 0 range).
WARNING: No objective has been specified. A constant zero objective will be used.
NOTE: The problem is a pure network instance. The ALGORITHM=NETWORK option is recommended for solving problems with this structure.
NOTE: The LP presolver value AUTOMATIC is applied.
NOTE: The LP presolver removed all variables and constraints.
NOTE: Optimal.
NOTE: Objective = 0.
Here is my code.
proc optmodel PRESOLVER=NONE;
* set initializing;
set <str> BLOCKS;
set PERIODS=1..100;
* data reading;
num reserves{BLOCKS};
num extraction{BLOCKS, PERIODS}, ls{BLOCKS, PERIODS};
read data libor.activeblocks_b into BLOCKS=[eqtypecode] reserves;
read data libor.active_b into [eqtypecode p_code] ls extraction;
var A{BLOCKS} >=0.1 <=3, V{BLOCKS} >=0.1 <=2, d{BLOCKS} >=0.1 <=1.5;
impvar _textraction{b in BLOCKS, p in PERIODS}=(reserves[b]/(1+(A[b]*(exp(-V[b]*(ls[b,p]-(d[b]/ls[b,p])))))));
min _dist{b in BLOCKS}=sum{p in PERIODS} ((_textraction[b,p]-extraction[b,p]*1000)^2);
solve;
for {b in BLOCKS, p in PERIODS} do;
textraction[b,p]=(reserves[b]/(1+(A[b]*(exp(-V[b]*(ls[b,p]-(d[b]/ls[b,p])))))));
end;
create data libor.coef from [block]=BLOCKS A V d;
create data libor.teor from [block p_code]={BLOCKS, PERIODS} textraction extraction ls;
quit;
Please help, I can't understand what's wrong with my code.
If you have some priorities among the objectives, you can minimize them one at a time sequentially, adding a new optimality constraint at each step. See the SECONDARY OBJECTIVE section of this SAS Global Forum 2014 paper.
If you want to perform multiobjective optimization and output a Pareto frontier, you can use PROC OPTLSO.
To avoid the failure and warning, you can either set the variables to an initial feasible solution before calling the solver or use the MULTISTART option in the SOLVE WITH NLP statement.
Your MIN statement declares several objectives (one for each b), but you have not specified which one you want to use. Either declare only one objective or use the OBJ clause in the SOLVE statement to pick one.
If you have some priorities among the objectives, you can minimize them one at a time sequentially, adding a new optimality constraint at each step. See the SECONDARY OBJECTIVE section of this SAS Global Forum 2014 paper.
If you want to perform multiobjective optimization and output a Pareto frontier, you can use PROC OPTLSO.
To avoid the failure and warning, you can either set the variables to an initial feasible solution before calling the solver or use the MULTISTART option in the SOLVE WITH NLP statement.
SAS Innovate 2025 is scheduled for May 6-9 in Orlando, FL. Sign up to be first to learn about the agenda and registration!
Learn how to run multiple linear regression models with and without interactions, presented by SAS user Alex Chaplin.
Find more tutorials on the SAS Users YouTube channel.