How do I use a minimization function in scipy with constraints How do I use a minimization function in scipy with constraints numpy numpy

How do I use a minimization function in scipy with constraints


You can do a constrained optimization with COBYLA or SLSQP as it says in the docs.

from scipy.optimize import minimizestart_pos = np.ones(6)*(1/6.) #or whatever#Says one minus the sum of all variables must be zerocons = ({'type': 'eq', 'fun': lambda x:  1 - sum(x)})#Required to have non negative valuesbnds = tuple((0,1) for x in start_pos)

Combine these into the minimization function.

res = minimize(getSharpe, start_pos, method='SLSQP', bounds=bnds ,constraints=cons)


Check .minimize docstring:

scipy.optimize.minimize(fun, x0, args=(), method='BFGS', jac=None, hess=None, hessp=None, \              bounds=None, constraints=(), tol=None, callback=None, options=None)

What matters the most in your case will be the bounds. When you want to constrain your parameter in [0,1] (or (0,1)?) You need to define it for each variable, such as:

bounds=((0,1), (0,1).....)

Now, the other part, sum(x)==1. There may be more elegant ways to do it, but consider this: instead of minimizing f(x), you minimize h=lambda x: f(x)+g(x), a new function essential f(x)+g(x) where g(x) is a function reaches it minimum when sum(x)=1. Such as g=lambda x: (sum(x)-1)**2.

The minimum of h(x) is reached when both f(x) and g(x) are at their minimum. Sort of a case of Lagrange multiplier method http://en.wikipedia.org/wiki/Lagrange_multiplier