Numerical methods for stochastic programs with second order dominance constraints with applications to portfolio optimization

Inspired by the successful applications of the stochastic optimization with second order stochastic dominance (SSD) model in portfolio optimization, we study new numerical methods for a general SSD model where the underlying functions are not necessarily linear. Specifically, we penalize the SSD constraints to the objective under Slater’s constraint qualification and then apply the well known stochastic approximation (SA) method and the level function method to solve the penalized problem. Both methods are iterative: the former requires to calculate an approximate subgradient of the objective function of the penalized problem at each iterate while the latter requires to calculate a subgradient. Under some moderate conditions, we show that w.p.1 the sequence of approximated solutions generated by the SA method converges to an optimal solution of the true problem. As for the level function method, the convergence is deterministic and in some cases we are able to estimate the number of iterations for a given precision. Both methods are applied to portfolio optimization problem where the return functions are not necessarily linear and some numerical test results are reported.     

ALM models based on second order stochastic dominance

Computational Management Science - We propose asset and liability management models in which the risk of underfunding is modelled based on the concept of stochastic dominance. Investment decisions...     

Stability analysis of stochastic programs with second order dominance constraints

In this paper we present a stability analysis of a stochastic optimization problem with stochastic second order dominance constraints. We consider a perturbation of the underlying probability measure in the space of regular measures equipped with pseudometric discrepancy distance (Römisch in Stochastic Programming. Elsevier, Amsterdam, pp 483–554, 2003). By exploiting a result on error bounds in semi-infinite programming due to Gugat (Math Program Ser B 88:255–275, 2000), we show under the Slater constraint qualification that the optimal value function is Lipschitz continuous and the optimal solution set mapping is upper semicontinuous with respect to the perturbation of the probability measure. In particular, we consider the case when the probability measure is approximated by an empirical probability measure and show an exponential rate of convergence of the sequence of optimal solutions obtained from solving the approximation problem. The analysis is extended to the stationary points.     

Multi-stage portfolio selection problem with dynamic stochastic dominance constraints

We study the multi-stage portfolio selection problem where the utility function of an investor is ambiguous. The ambiguity is characterized by dynamic stochastic dominance constraints, which are able to capture the dynamics of the random return sequence during the investment process. We propose a multi-stage dynamic stochastic dominance constrained portfolio selection model, and use a mixed normal distribution with time-varying weights and the K-means clustering technique to generate a scenario tree for the transformation of the proposed model. Based on the scenario tree representation, we derive two linear programming approximation problems, using the sampling approach or the duality theory, which provide an upper bound approximation and a lower bound approximation for the original nonconvex problem. The upper bound is asymptotically tight with infinitely many samples. Numerical results illustrate the practicality and efficiency of the proposed new model and solution techniques.

     

Convergence analysis of stationary points in sample average approximation of stochastic programs with second order stochastic dominance constraints

Sample average approximation (SAA) method has recently been applied to solve stochastic programs with second order stochastic dominance (SSD) constraints. In particular, Hu et al. (Math Program 133:171–201, 2012) presented a detailed convergence analysis of $\epsilon $ -optimal values and $\epsilon $ -optimal solutions of sample average approximated stochastic programs with polyhedral SSD constraints. In this paper, we complement the existing research by presenting convergence analysis of stationary points when SAA is applied to a class of stochastic minimization problems with SSD constraints. Specifically, under some moderate conditions we prove that optimal solutions and stationary points obtained from solving sample average approximated problems converge with probability one to their true counterparts. Moreover, by exploiting some recent results on large deviation of random functions and sensitivity analysis of generalized equations, we derive exponential rate of convergence of stationary points.     

An incremental bundle method for portfolio selection problem under second-order stochastic dominance

Numerical Algorithms - We propose a preconditioner to accelerate the convergence of the GMRES iterative method for solving the system of linear equations obtained from discretize-then-optimize...     

Crossing points of distributions and a theorem that relates them to second order stochastic dominance

We state formal definitions for crossing points in pairs of distributions and give a detailed proof of a theorem that relates those points to the second order stochastic dominance (SSD). The theorem states that the fulfillment of the area balance conditions for SSD at the t values that correspond to crossing points, and at the limit t→∞, is a necessary and sufficient condition for its fulfillment at all t: {−∞<t<∞}, as required for the existence of SSD. We provide examples for the application of the theorem in the case of continuous distributions, including a continuous counter example to prove that the Mean-Variance criterion is not sufficient to state preferences under risk aversion.     

Second-order stochastic dominance constrained portfolio optimization: Theory and computational tests

Due to the definition of second-order stochastic dominance (SSD) in terms of utility theory, portfolio optimization with SSD constraints is of major practical interest. We contribute to the field in two ways: first, we present a self-contained theory with some new results and new proofs of known results; second, we perform a set of tests for computational efficiency. We provide new and simple arguments for the formulation of SSD constraints in a mathematical programming framework. For many individuals, an SSD constraint may seem too severe wherefore various relaxations (ASSD), have been proposed. We introduce yet another relaxation, directional SSD, where a candidate portfolio is admissible if a step from the benchmark in the direction of the candidate yields a dominating portfolio. Optimal step size depends on individual preferences reflected by the objective function. We compare computational efficiency of seven approaches for SD constrained portfolio problems, including SSD and ASSD constrained cases.     

Valid inequalities and restrictions for stochastic programming problems with first order stochastic dominance constraints

Stochastic dominance relations are well studied in statistics, decision theory and economics. Recently, there has been significant interest in introducing dominance relations into stochastic optimization problems as constraints. In the discrete case, stochastic optimization models involving second order stochastic dominance constraints can be solved by linear programming. However, problems involving first order stochastic dominance constraints are potentially hard due to the non-convexity of the associated feasible regions. In this paper we consider a mixed 0–1 linear programming formulation of a discrete first order constrained optimization model and present a relaxation based on second order constraints. We derive some valid inequalities and restrictions by employing the probabilistic structure of the problem. We also generate cuts that are valid inequalities for the disjunctive relaxations arising from the underlying combinatorial structure of the problem by applying the lift-and-project procedure. We describe three heuristic algorithms to construct feasible solutions, based on conditional second order constraints, variable fixing, and conditional value at risk. Finally, we present numerical results for several instances of a real world portfolio optimization problem. This research was supported by the NSF awards DMS-0603728 and DMI-0354678.     

Asymptotic behavior for second order stochastic evolution equations with memory

In this paper, a general second order integro-differential evolution equation with memory driven by multiplicative noise is considered. We prove the existence of global mild solution and asymptotic stability of the zero solution using Lyapunov function techniques. Moreover, we discuss three examples to show that the asymptotic stability results can be applied to various partial differential equations.     

Weak second order S-ROCK methods for Stratonovich stochastic differential equations

It is well known that the numerical solution of stiff stochastic ordinary differential equations leads to a step size reduction when explicit methods are used. This has led to a plethora of implicit or semi-implicit methods with a wide variety of stability properties. However, for stiff stochastic problems in which the eigenvalues of a drift term lie near the negative real axis, such as those arising from stochastic partial differential equations, explicit methods with extended stability regions can be very effective. In the present paper our aim is to derive explicit Runge–Kutta schemes for non-commutative Stratonovich stochastic differential equations, which are of weak order two and which have large stability regions. This will be achieved by the use of a technique in Chebyshev methods for ordinary differential equations.     

Relaxations of linear programming problems with first order stochastic dominance constraints

Linear stochastic programming problems with first order stochastic dominance (FSD) constraints are non-convex. For their mixed 0-1 linear programming formulation we present two convex relaxations based on second order stochastic dominance (SSD). We develop necessary and sufficient conditions for FSD, used to obtain a disjunctive programming formulation and to strengthen one of the SSD-based relaxations.     

Supplement: Efficient weak second order stochastic Runge-Kutta methods for non-commutative Stratonovich stochastic differential equations

This paper gives a modification of a class of stochastic Runge-Kutta methods proposed in a paper by Komori (2007). The slight modification can reduce the computational costs of the methods significantly.     

Portfolio construction based on stochastic dominance and target return distributions

Mean-risk models have been widely used in portfolio optimization. However, such models may produce portfolios that are dominated with respect to second order stochastic dominance and therefore not optimal for rational and risk-averse investors. This paper considers the problem of constructing a portfolio which is non-dominated with respect to second order stochastic dominance and whose return distribution has specified desirable properties. The problem is multi-objective and is transformed into a single objective problem by using the reference point method, in which target levels, known as aspiration points, are specified for the objective functions. A model is proposed in which the aspiration points relate to ordered outcomes for the portfolio return. This concept is extended by additionally specifying reservation points, which act pre-emptively in the optimization model. The theoretical properties of the models are studied. The performance of the models on real data drawn from the Hang Seng index is also investigated.     

Processing second-order stochastic dominance models using cutting-plane representations

Second-order stochastic dominance (SSD) is widely recognised as an important decision criterion in portfolio selection. Unfortunately, stochastic dominance models are known to be very demanding from a computational point of view. In this paper we consider two classes of models which use SSD as a choice criterion. The first, proposed by Dentcheva and Ruszczyński (J Bank Finance 30:433–451, 2006), uses a SSD constraint, which can be expressed as integrated chance constraints (ICCs). The second, proposed by Roman et al. (Math Program, Ser B 108:541–569, 2006) uses SSD through a multi-objective formulation with CVaR objectives. Cutting plane representations and algorithms were proposed by Klein Haneveld and Van der Vlerk (Comput Manage Sci 3:245–269, 2006) for ICCs, and by Künzi-Bay and Mayer (Comput Manage Sci 3:3–27, 2006) for CVaR minimization. These concepts are taken into consideration to propose representations and solution methods for the above class of SSD based models. We describe a cutting plane based solution algorithm and outline implementation details. A computational study is presented, which demonstrates the effectiveness and the scale-up properties of the solution algorithm, as applied to the SSD model of Roman et al. (Math Program, Ser B 108:541–569, 2006).     

A new representation for second order stochastic integral-differential operators and its applications

In this paper, we’ll prove new representation theorems for a kind of second order stochastic integraldifferential operator by stochastic differential equations(SDEs) and backward stochastic differential equations(BSDEs) with jumps, and give some applications.     

Multi-objective stochastic programming for portfolio selection

Generally, in the portfolio selection problem the Decision Maker (DM) considers simultaneously conflicting objectives such as rate of return, liquidity and risk. Multi-objective programming techniques such as goal programming (GP) and compromise programming (CP) are used to choose the portfolio best satisfying the DM’s aspirations and preferences. In this article, we assume that the parameters associated with the objectives are random and normally distributed. We propose a chance constrained compromise programming model (CCCP) as a deterministic transformation to multi-objective stochastic programming portfolio model. CCCP is based on CP and chance constrained programming (CCP) models. The proposed program is illustrated by means of a portfolio selection problem from the Tunisian stock exchange market.     

On the relative efficiency of nth order and DARA stochastic dominance rules

It is known that third order stochastic dominance implies DARA dominance while no implications exist between higher orders and DARA dominance. A recent contribution points out that, with regard to the problem of determining lower and upper bounds for the price of a financial option, the DARA rule turns out to improve the stochastic dominance criteria of any order. In this paper the relative efficiency of the ordinary stochastic dominance and DARA criteria for alternatives with discrete distributions are compared, in order to see if the better performance of DARA criterion is also suitable for other practical applications. Moreover, the operational use of the stochastic dominance techniques for financial choices is deepened.     

On the numerical expansion of a second order stochastic process

The utility of the classical Karhunen-Loève expansion of a second order process is limited to its practical derivation, because it depends upon the solution of a Fredholm integral equation associated with it whose kernel is the covariance function of the process. So, in this paper we study two numerical procedures for solving such equations, the Rayleigh-Ritz and the collocation methods, and also essay two different bases of L²-orthogonal functions in order to perform the algorithms, Legendre polynomials and trigonometric functions, on two well-known processes, the Wiener-Lévy process and the Brownian-bridge. The accuracy of the numerical results in relation to the real ones, as well as comparative studies among both procedures are also included.     

Representation of stochastic processes of second order and linear operations

Series representations are obtained for the entire class of measurable, second order stochastic processes defined on any interval of the real line. They include as particular cases all earlier representations; they suggest a notion of “smoothness” that generalizes well-known continuity notions; and they decompose the stochastic process into two orthogonal parts, the smooth part and a strongly discontinuous part. Also linear operations on measurable, second order processes are studied; it is shown that all “smooth” linear operations on a process, and, in particular, all linear operations on a “smooth” process, can be approximated arbitrarily closely by linear operations on the sample paths of the process.