Introduction to the Global Optimization Toolbox
Overview of Features
Returns Data and Minimum Acceptable Return
Maple has partnered with Noesis Solutions to develop a new version of the Maple Global Optimization Toolbox that is powered by Optimus® technology, first released with the Global Optimization Toolbox 17.
The new global solver offers two solver methods, a differential evolution algorithm and a hybrid algorithm that uses interpolating response surface models, as well as new options to guide the search. In addition, the updated toolbox is fully compatible with existing code.
Highlights of some of the new options available for the two new solver methods:
Differential Evolution Algorithm
Average step width — the variation in design variables used as stopping criterion
Population size — the number of designs evaluated during each iteration
Inverse crossover probability — the probability that variables might be the same as predecessors without adjustment
Target weighting factor — applies weighted differences of the variables of the randomly selected designs of the previous generation, allowing you to set the balance between computation speed and probability of success
Adaptive Stochastic Search Methods
Theta method — methods for maximizing the predictive quality of the model, maximizing the likelihood of the set of points with respect to the model, and minimizing the semi-norm of the correlation matrix
Number of sigma — how much attention is given to the estimated standard deviation when selecting the next points
Nugget — the smoothing factor for the kriging response surface model
Optimum search method — different methods for locating new promising points in the design space
The following sample application solves an investment optimization problem in two ways.
Traditional investment performance benchmarks, like the Sharpe Ratio, approximate the returns distribution with mean and standard deviation. This, however, assumes the distribution is normal. Many modern investments vehicles, like hedge funds, display fat tails, and skew and kurtosis in the returns distribution. Hence, they cannot be adequately benchmarked with traditional approaches.
One solution, proposed by Shadwick and Keating in 2002 is the Omega Ratio. This divides the returns distribution into two halves – the area below a target return, and the above a target return. The Omega Ratio is simply the former divided by the latter. A higher value is better.
For a set of discrete returns, the Omega Ratio is given by
where L is a target return and R is a vector of returns.
This application finds the asset weights that maximize the Omega Ratio of a portfolio of ten investments, given their simulated monthly returns and a target return.
This is a non-convex problem, and requires global optimizers for a rigorous solution. However, a transformation of the variables (only valid for Omega Ratios of over 1) converts the optimization into a linear program.
This application implements both approaches, the former using Maple's Global Optimization Toolbox, and the latter using Maple's linear programming features. For the data set provided in this application, both approaches give comparable results.
Monthly hedge fund returns
Number of funds
Number of returns for each fund
OmegaRatio≔procL,returns,weights local weightedReturns,above,below,N,S:N≔LinearAlgebraColumnDimensionreturns:S≔LinearAlgebraRowDimensionreturns: weightedReturns≔seq⁡add⁡returnsi,j⋅weightsj,j=1..N,i=1..S:below≔selectx→evalb⁡x<L,weightedReturns:above≔selectx→evalb⁡x>L,weightedReturns: return add⁡a − L,ainaboveadd⁡L − a,ainbelowend proc:
"Strawman" portfolio of equal weights at the target return
The new Global Optimization Toolbox finds the optimum asset weights easily.
Optimized Omega Ratio
Optimized investment weights
Finding the optimum asset weights using the linear programming method:
Optimized Omega Ratio
Download Help Document