Xinyuan
September 25, 2016, 3:32pm
1
Dear all,
I was trying to minimize a nonnegative object. But the result shows unbounded. Could u help me find out what’s wrong here ?
The code is to perform the algorithm shown in the following formula:
load('data1.mat');
[N,d] = size (Xtrain);
beta = ones(d,1);
solution = zeros(10,d);
for k = 1:10
cvx_begin
variable x(d);
minimize ((Ytrain - Xtrain*x)'*(Ytrain - Xtrain*x)/N + 0.5* sum(abs(x)./(sqrt(abs(beta))+ eps)));
cvx_end
solution(k, 1:d) = x';
beta = x;
clear x;
end
Can you show the output from running it? Which iteration (value of k) did your get the unbounded message? Did you look at the beta which was used in that iteration? Did it have any very large magnitude components?
We haven’t seen data1.mat, so we don’;t know whether your input data might be very bad numerically (spanning many orders of magnitude). Perhaps the solver is having numerical difficulties with challenging input data, to include value of beta somwhere for k=2:10… Have you tried more than one solver?
Xinyuan
September 26, 2016, 12:01am
3
Thank you for such a quick reply.
The dataset has been uploaded to google drive at the following link.https://drive.google.com/drive/folders/0BwJJX4_LUeQNSFM0SDdoSl9iczg?usp=sharing
I set beta all ones at the beginning.
The “unbounded” first appear when k = 3, and the outcome of the first 3 steps are shown as below:
Calling SDPT3 4.0: 622 variables, 221 equality constraints
num. of constraints = 221
SDPT3: Infeasible path-following algorithms
version predcorr gam expon scale_data
0|0.000|0.000|3.6e+00|7.6e+05|2.5e+08| 0.000000e+00 0.000000e+00| 0:0:01| chol 1 1
number of iterations = 14
Status: Solved
Calling SDPT3 4.0: 622 variables, 221 equality constraints
num. of constraints = 221
SDPT3: Infeasible path-following algorithms
version predcorr gam expon scale_data
0|0.000|0.000|2.1e+00|7.6e+05|1.1e+19| 0.000000e+00 0.000000e+00| 0:0:00| chol 2 2
number of iterations = 17
Status: Solved
Calling SDPT3 4.0: 622 variables, 221 equality constraints
num. of constraints = 221
SDPT3: Infeasible path-following algorithms
version predcorr gam expon scale_data
0|0.000|0.000|2.1e+00|7.6e+05|1.9e+23| 0.000000e+00 0.000000e+00| 0:0:00| chol 2 2
number of iterations = 17
Status: Unbounded
I replicated your results using SDPT3. On k = 2, cvx_optval = -7.53569e+07. On k = 3, the solver states “sqlp stop: primal problem is suspected of being infeasible” Notice that it doesn’t state that the problem is definitely infeasible (of course, that’s the same message I get when I try it on an actual infeasible problem, such as 1 <= x <= 0). I think it is having numerical difficulty because some of the terms in the sum are abs(x) times 1/sqrt(eps) = abs(x) * 6.7e7… I guess this is giving SDPT3 numerical troubles.
SeDuMi did make it through k = 10, and converged to cvx_optval = 742.242, and was already getting close by k = 3, but was pretty wild on k = 2 (cvx_optbal = 295519). I have no idea whether this converged to the correct answer, but it converged. I’ll let you worry about how to fix things up.
Xinyuan
September 26, 2016, 1:29pm
5
Thank you so much. It works after I switch to SeDumi.
Then I was trying to replace the linear expression with logit withminimize (-sum(labelTrain.*log(exp(trainData*x)./(1+exp(trainData*x)))+(1-labelTrain).*log(1-exp(trainData*x)./(1+exp(trainData*x))))+ lambda* sum(abs(x)./(sqrt(abs(beta))+ eps)));
But it turns out an error
Error using + (line 83)
Error in - (line 21)
Error in test (line 23)log(exp(trainData x)./(1+exp(trainDatax)))+(1-labelTrain).log(1-exp(trainData x)./(1+exp(trainData x))))+
I do not understand this error information and how can I make it work?
You are in violation of CVX’s rules regarding log-concave functions.
You can find the rules in the CVX 3.0Beta Users Guide at http://web.cvxr.com/cvx/beta/doc/advanced.html#log-convexity . The rules stated there also apply to CVX 2,.1 even though they are not mentioned in the current 2.1 Users’ Guide. They are also discussed in the first answer by mcg in Log of sigmoid function .
“Regular” logistic regression can be done in CVX. See http://cvxr.com/cvx/examples/cvxbook/Ch07_statistical_estim/html/logistics.html for an example.
