Turn on suggestions

Auto-suggest helps you quickly narrow down your search results by suggesting possible matches as you type.

Showing results for

Options

- RSS Feed
- Mark Topic as New
- Mark Topic as Read
- Float this Topic for Current User
- Bookmark
- Subscribe
- Mute
- Printer Friendly Page

🔒 This topic is **solved** and **locked**.
Need further help from the community? Please
sign in and ask a **new** question.

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Posted 02-23-2017 10:04 AM
(1311 views)

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.

1 ACCEPTED SOLUTION

Accepted Solutions

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

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.

4 REPLIES 4

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Hi, Rob Pratt, thanks for quick reply!

What should I do if I wanna minimize all that objectives? Is there some procs?

And when I picked one objective using OBJ, it says solution status failed 😞 and this WARNING: Objective function cannot be evaluated at starting point. in the log.

What should I do if I wanna minimize all that objectives? Is there some procs?

And when I picked one objective using OBJ, it says solution status failed 😞 and this WARNING: Objective function cannot be evaluated at starting point. in the log.

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

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.

- Mark as New
- Bookmark
- Subscribe
- Mute
- RSS Feed
- Permalink
- Report Inappropriate Content

Thank you very much Rob!

http://support.sas.com/documentation/cdl/en/ormpex/67518/HTML/default/viewer.htm#ormpex_ex11_toc.htm this example helped me to avoid failure and I've found optiamal solution for one objective.

thanks for suggesting PROC OPTLSO, I'll try to use it!

http://support.sas.com/documentation/cdl/en/ormpex/67518/HTML/default/viewer.htm#ormpex_ex11_toc.htm this example helped me to avoid failure and I've found optiamal solution for one objective.

thanks for suggesting PROC OPTLSO, I'll try to use it!

Registration is open! SAS is returning to Vegas for an AI and analytics experience like no other! Whether you're an executive, manager, end user or SAS partner, SAS Innovate is designed for everyone on your team. Register for just $495 by 12/31/2023.

**If you are interested in speaking, there is still time to submit a session idea. More details are posted on the website. **

Multiple Linear Regression in SAS

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.