Next Article in Journal
A Definite Integral Involving the Logarithmic Function in Terms of the Lerch Function
Next Article in Special Issue
Moving Information Horizon Approach for Dynamic Game Models
Previous Article in Journal
On a Periodic Boundary Value Problem for a Fractional–Order Semilinear Functional Differential Inclusions in a Banach Space
Previous Article in Special Issue
Cournot Duopoly Games: Models and Investigations
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Robust Portfolio Optimization in an Illiquid Market in Discrete-Time

Financial Engineering and Risk Management laboratory, National Research University Higher School of Economics, 20 Myasnitskaya Ulitsa, 101000 Moscow, Russia
Mathematics 2019, 7(12), 1147; https://doi.org/10.3390/math7121147
Submission received: 27 October 2019 / Revised: 20 November 2019 / Accepted: 21 November 2019 / Published: 24 November 2019
(This article belongs to the Special Issue Mathematical Game Theory 2019)

Abstract

:
We present a robust dynamic programming approach to the general portfolio selection problem in the presence of transaction costs and trading limits. We formulate the problem as a dynamic infinite game against nature and obtain the corresponding Bellman-Isaacs equation. Under several additional assumptions, we get an alternative form of the equation, which is more feasible for a numerical solution. The framework covers a wide range of control problems, such as the estimation of the portfolio liquidation value, or portfolio selection in an adverse market. The results can be used in the presence of model errors, non-linear transaction costs and a price impact.
MSC:
91A25; 91G10
JEL Classification:
C61; C7; G11

1. Introduction

In this paper, we consider a robust approach to the constrained portfolio optimization problem in the presence of transaction costs and a price impact on the market. Optimal portfolio selection has become a popular topic in the mathematical finance literature since the pioneering works of Samuleson [1] and Merton [2]. Both papers consider a stochastic dynamic programming approach to the portfolio optimization problem on an efficient frictionless market with no trading limits. The continuous-time Samuelson model has been extended in various studies, see Reference [3] for an overview of the field. The works by Magill & Constantinides [4] and Davis & Norman [5], among others, paved the way for a new direction in research by considering a market with friction and introducing linear transaction costs, see Reference [6] for an overview. Portfolio selection for a discrete-time market, however, has not been given as much attention in the academic literature as its continuous-time counterpart, since in most cases it does not provide the Hamilton-Jacoby-Bellman equation (or quasi-variational inequality) to get the solution in a closed form, but requires a numerical solution instead. The discrete-time approach has been implemented in Bertsimas & Lo [7] and Almgren [8,9] for the optimal execution problem when an investor chooses the optimal strategy of the portfolio liquidation over a given number of consecutive periods.
The assumptions of the classic financial market model have undergone heavy criticism in view of the LTCM crisis in 1998 [10]. One of the main topics of discussion was the VAR risk metric, which assumed the Samuelson price model based on geometric Brownian motion and underestimated the risk of sudden price movements. Another question concerned the significance of market liquidity and transaction costs. Many works have been presented to deal with both problems. In References [8,9,11], the authors introduce a trade-off between profit and volatility in a portfolio optimization framework in the presence of costs. Reference [12] accounts for the sudden price movements by introducing jumps into the price model. Nevertheless, most frameworks are still based on the arithmetic or the geometric Brownian motion model. In this paper, we present a discrete-time portfolio optimization framework in the presence of ambiguity in an illiquid market. The approach can be used to account for the impact of model error, which might be significant in many cases and especially during a financial crisis, when the investor needs a conservative estimate of the portfolio liquidation value, or when the portfolio needs to be rebalanced in a short period of time.
We adapt the framework of Smirnov [13,14] and introduce a portfolio optimization problem as a game against nature where an investor tries to maximize the expected worst-case reward from a portfolio by making general and non-restrictive assumptions about the market. The robust optimization techniques have been introduced for a variety of economic problems recently. Deng et al. in Reference [15] combines the robust approach with the Markowitz framework to solve a one-period portfolio selection problem in a frictionless market, while Reference [16] solves a multi-stage asset allocation problem with linear costs by considering a series of linear programming problems with ambiguous constraints. In Reference [17], the authors present a game theoretic approach to the risk management control problem in the market with ambiguous coefficients of the Samuelson model. Robust optimization techniques have also been used to quantify the impact of the model error when assessing portfolio risks, see for example, Reference [18]. In this paper, we provide a framework for the multi-period portfolio optimization in the presence of model error by introducing the ambiguity in the market price distribution. We consider a general convex transaction costs function to account for the non-linear cost growth, and non-Markov dynamics of the system to account for a price impact.
The paper is structured as follows. In Section 1 we provide an overview of the problem. Section 2 provides a short overview of the research. In Section 3, we formulate the optimization framework, introduce the assumptions about the market, and present the corresponding Bellman-Isaacs equation. In Section 4, we study the sufficient conditions that allow the equation to be simplified, making it more feasible for a numerical solution. Section 5 discusses the results.

2. Overview of the Main Results

Consider an investor who tries to maximize the reward from the provided portfolio by the end of the investment horizon. Usually, the investor assumes some market price model, which allows them to estimate the distribution of the portfolio future wealth for any appropriate strategy, and make the correct investment decision. In this research, we assume that the investor does not possess the knowledge of the market price model. Instead, they have the estimates of the expected value and range of the future price movements. Both estimates are usually relevant for any investment process and can be interpreted in business terms. Below, we formulate the optimization problem in the presence of ambiguity and reduce it to the corresponding Bellman-Isaacs equation.
The rest of the research is devoted to simplifying the Bellman-Isaacs equation, since, in its original form, the equation requires estimating the worst expected value of the value function among all the appropriate candidate market models, each represented by the corresponding probability measure. Optimization over a set of measures does not provide an analytic solution in many practical cases, and the problem should be solved numerically. Nevertheless, some insight about the solution can be obtained in the most general case. In Theorem 2 below we show that the worst expected value is actually attained over a subset of atomic measures where the maximum number of atoms is defined by the dimension of the considered space. The result can be strengthened if (i) the value function remains concave with respect to the ambiguous parameters throughout the investment periods, and (ii) the problem can be formulated in terms of the position market values instead of volumes and market prices. In Theorem 4 we provide the corresponding sufficient conditions and prove that in this case the worst expected value is actually attained over a subset of atomic measures concentrated in the extreme points of the support. This allows to search for the optimal measure only on a subset of the support which might prove useful for a numerical implementation.

3. The Robust Portfolio Optimization Problem

Consider a market with a finite horizon and discrete time t = 0 , 1 , , N . The proposed framework describes a consecutive process of decision-making by an investor who manages a portfolio trying to attain maximum reward at the end of the strategy. Each decision is made at the beginning of the investment period, while the outcome is observed at the end of the period. Throughout the paper, t will denote the beginning of the corresponding period.
Assume that there are n risky assets and one risk-free asset on the market. Let X t i be the market price of the i-th risky asset at time t and let X t = ( X t 1 , , X t n ) . Let Y t be the market price of the risk-free asset. In the rest of the paper, the term “price” means the market price. The price of any risky asset at time t might be undetermined at t 1 while the exact value of Y t is assumed to be known at t 1 (such property will be referred to as predictability later on). An example of a risk-free asset is the price of money in the portfolio’s base currency. Another example is a sovereign bond account if the bonds are assumed risk-free and liquid, and if they are issued and traded in the portfolio’s base currency. For the rest of the paper, we assume that all prices are positive. Note that, while this is a classic assumption that is reasonable in many practical cases, it does not fully account for the risk of bankruptcy of the issuer. We also do not account for the price discretization and assume that X t 1 , , X t n , Y t R .
By portfolio at time t we mean a predictable vector of asset volumes H ^ t = ( H t 1 , , H t n , H t Y ) . Predictability means that the portfolio structure for the next investment period is decided at the beginning of that period. We assume that all the assets are infinitely divisible: H t 1 , , H t n , H t Y R . At the start of the strategy the investor is provided with the initial portfolio H ^ 0 . Let H t = ( H t 1 , , H t n ) R n denote the volumes of the risky assets at t.
To account for the friction in the market, we assume that each trade by the investor incurs transaction costs. Costs can be categorized as explicit, such as broker commissions and fees, and implicit, which are due to the insufficient liquidity of the market: a trade moves the corresponding market price against it, therefore the average price of the trade might be worse than the market price if the volume is sufficiently large. Thus, in the frictionless market, a portfolio can be rebalanced with no loss of market value. In the real market, however, any rebalancing might lead to the loss of portfolio market value, which can be thought of as carrying additional costs. We assume that there are no costs incurred by the trades in the risk-free asset.
We assume that all information about the market is contained in the prices of the assets. The system state can therefore be represented as a combination of the price values and the portfolio. Let W t Y = H t Y Y t be the market value of the risk-free position at time t. Since the risk-free asset is assumed to be absolutely liquid, it is easier to operate in terms of W t Y then work with H t Y and Y t separately. Thus, by the system state S t at time t, we would mean a combination of H t , X t and W t Y . We will also refer to the system history S 0 , , S t up to the moment t as S ¯ t . For the ease of notation, the functional dependence on the “·” argument will mean the dependence on the system’s past history throughout the paper. The dependence will sometimes be explicitly written out to clarify the dependence on the components of S ¯ t .
Let C t ( · , H t ) denote the costs carried when transitioning to the portfolio H t . The function C t ( · , H t ) is assumed predictable. We do not introduce any additional inflows and outflows of assets, hence the following budget equation should hold for all t:
Δ H t Y Y t 1 + Δ H t X t 1 = C t ( · , H t ) ,
where Δ is the backward difference operator, that is, Δ H t = H t H t 1 . Both in Equation (1) and in the rest of the paper, the product of vectors should be understood as the inner product. Equation (1) allows us to express H t Y as a function of H t and the past history as
H t Y ( · , H t ) = Y t 1 1 W t 1 Y ( H t H t 1 ) X t 1 C t ( · , H t ) ,
thus the optimization problem will be stated in terms of H t instead of H ^ t in the rest of the paper, and by the portfolio strategy we mean the sequence of H 1 , , H N .
We assume that at each time t the portfolio structure might be subject to the predictable constraints:
H t D t ( · ) .
The constraint sets D t ( · ) can be interpreted as the trading limits, which might serve the risk-management purposes or originate from the aims of the portfolio strategy. For example, consider a portfolio liquidation problem where an investor tries to close the risky positions over N consecutive periods. The problem can be formulated in several ways, one of which is gradually reducing the allowed volumes of risky assets via the constraints, that is,
D 1 ( · ) D 2 ( · ) D N ( · ) .
D t ( · ) = would have meant the inability to continue with the investment process and the early termination of the strategy. This case would have complicated the presentation of the results with technical details, so we assume the non-emptiness of D t ( · ) throughout the paper. As for the other properties, from the risk-management point of view, it is usually reasonable to assume the boundedness of the constraint sets, since the unbounded sets might encourage entering infinitely large risky positions that are subject to infinite risk.
At t 1 , the risk-free price Y t is assumed to be known by the investor. Consider the risk-free price dynamics in a multiplicative form as
Y t = R t Y t 1 , t = 1 , , N ,
where R t > 0 is assumed predictable. R t 1 can be interpreted as the risk-free rate of return. By combining (2) and (3), we can represent W t Y as a function of portfolio H t and the past history as
W t Y ( · , H t ) = H t Y Y t = R t H t Y Y t 1 = R t W t 1 Y ( H t H t 1 ) X t 1 C t ( · , H t ) .
As for the risky assets, the only assumption about the prices considers the range of returns. Let the risky price dynamics be represented in a multiplicative form as
X t i = ξ t i X t 1 i , ξ t i > 0 , i = 1 , , n , t = 1 , , N ,
where ξ t i 1 is the return of the i-th asset over the investment period starting at t 1 . Let ξ t = ( ξ t 1 , , ξ t n ) R n and let Λ ( ξ t ) be a diagonal matrix with ξ t on the main diagonal. Then we can write
X t = Λ ( ξ t ) X t 1 , t = 1 , , N .
The main assumption about market ambiguity is that we assume that ξ t belongs to a predefined set K t ( S 0 , , S t 1 ) R n for each t = 1 , , N , whereas the sets K 1 ( S ¯ 0 ) , , K N ( S ¯ N 1 ) are assumed convex and compact.
Below, we present the notion of optimality and the Bellman-Isaacs equation. Classic portfolio selection frameworks usually assume that the investor maximizes the expected reward function of the portfolio at the end of the strategy. for example, in Reference [7], the authors consider the negated value of the execution costs as a reward, while in Reference [8], the mixture of the execution costs and the market risk is optimized. Many frameworks are based on expected utility maximization, where the form of the utility reflects the risk-aversion of the investor. Our approach accounts for both the uncertainty and the ambiguity of the market—the former refers to the stochastic nature of the market given a specific market model, while the latter refers to the uncertainty about the market model itself. In the presence of ambiguity, the standard expected von Neumann-Morgenstern utility cannot be used to define the optimal strategy, see the Ellsberg paradox [19]. Instead, we say that the portfolio strategy is optimal if it maximizes the robust Savage representation as defined in Reference [20] (see Theorem 2.80). For this, we assume the reasonable additional axioms of uncertainty aversion and certainty independence presented in Reference [20], when describing the behavior of the investor. The form of the robust Savage representation resembles the worst-case expected reward, which makes it a natural extension of the classic framework.
Remark 1.
Optimization of the worst-case expected von Neumann-Morgenstern utility might still be appropriate when applied to a particular investment process. Assume, for example, that the utility function is defined by the risk-averse decision-maker based on the subjective estimates of the market. At the same time, the optimal strategy is obtained by another agent, the analyst, who tries to maximize the expected value of the provided utility in the worst case based on the analyst’s independent estimate. Thus, from the analyst’s point of view, the utility is an exogenous reward function which should be used in the framework.
We represent the ambiguity of the market via a set of probability measures describing price dynamics. For each t = 1 , , N , consider a set of probability measures P t ( · ) , where each measure describes the uncertainty of risky prices over the corresponding period (in terms of the distribution of ξ t , see (5)). Let P t , N ( · ) be the set of probability measures describing the risky prices across the investment periods starting at t 1 or later. Throughout the paper, we assume the following properties:
  • for each t = 1 , , N and every Q P t ( · ) , supp Q K t ( · ) and K t ( · ) has compact and convex values;
  • for each t = 1 , , N and every Q P t ( · ) , ξ Q ( d ξ ) = E t ( · ) , where E t ( · ) K t ( · ) ;
  • (The Rectangularity assumption) for each admissible portfolio strategy,
    P t , N ( S ¯ t 1 ) = P t ( S ¯ t 1 ) × × P N ( S ¯ N 1 ) .
Property 1 refers to the previously mentioned assumption about the price dynamics. The integral of the vector ξ in Property 2 should be interpreted element-wise as ξ i Q ( d ξ ) = E t i ( · ) for i = 1 , , n and E t ( · ) = E t 1 ( · ) , , E t n ( · ) . Property 2 means that at the beginning of each period the investor has certain expectations about the future price values, whereas E ( · ) is the expected value based on the investor’s estimate. The assumption is not restrictive since, in many practical cases, the investment process is based on the estimate of future price movements, which is provided by the analyst. Besides, dropping the assumption would lead to the degeneration of the strategy in some practical cases, which will be illustrated at the end of the Section. Property 3 is a generalization of the Rectangularity assumption of References [21,22], that can be interpreted as the independence of market ambiguity in the future from the present market model. As an example of ambiguity that does not comply with (6), consider a single-asset market with the multiplicative dynamics (5) and N = 2 . For each t, let K t ( · ) [ d , u ] , 0 < d < u , and let ξ t be either uniformly distributed on [ d , u + d 2 ] with probability p or uniformly distributed on ( u + d 2 , u ] with probability 1 p , where p [ 0 , 1 ] is unknown, which makes the market model ambiguous. Let Q t , p denote the corresponding probability measure at t. Then P t = Q t , p : p [ 0 , 1 ] for t = 1 , 2 , and
P 1 × P 2 = Q 1 , p × Q 2 , p : p , p [ 0 , 1 ] .
On the other hand,
P 1 , 2 = Q 1 , p × Q 2 , p : p [ 0 , 1 ] ,
hence P 1 , 2 P 1 × P 2 . The Rectangularity assumption means that the investor assumes that the market is “replenished” at each period in terms of ambiguity, which is reasonable in the context of the worst-case framework. See Reference [22] for a similar example for the Ellsberg paradox.
In this paper, we do not focus on the formal definition of the admissible strategies of the investor and the market, for a formal definition see for example, Reference [23]. Let J ( S ¯ N ) be the reward function which depends on the portfolio structure and the market state. We assume that J has values from R { } , where the value means the failure of the investment strategy. Let H be the set of the admissible portfolio strategies, and let H t , N ( · ) be the set of the admissible strategies over the periods starting at t 1 or later. Let P = P 1 , N ( S 0 ) be the set of probability measures describing the market dynamics across all periods. Then we can say that the portfolio strategy H * H is optimal if
sup H H inf Q P J ( S ¯ N ) d Q = inf Q P J ( S ¯ N * ) d Q ,
where S ¯ N * is the system history for the strategy H * . The supremum in (7) might not be achieved in general. However, it is usually the case in practice, if the problem is stated correctly, so throughout the paper we assume the existence of an optimal strategy. Note that inf Q P J ( · ) d Q is the robust Savage representation we mentioned earlier.
Remark 2.
Definition (7) requires integrability of the underlying reward function. In the presented practical framework, we do not focus on the question of measurability but assume the required form of measurability as needed. It will be shown in the next section, that we are mainly interested in the infimum over the subset of atomic measures from P t ( · ) that are concentrated in n + 1 or fewer points of K t ( · ) , t = 1 , , N , and every function is integrable with respect to any measure from that subset.
Let V t 1 * ( S ¯ t 1 ) be the supremum of the robust Savage representation as estimated at time t 1 , t = 1 , , N :
V t 1 * ( S ¯ t 1 ) = sup H t , N ( S ¯ t 1 ) inf Q t , N P t , N ( S ¯ t 1 ) J ( S ¯ N ) d Q t , N .
For t = N , let V N * ( S ¯ N ) J ( S ¯ N ) . We refer to V t * as the “value function” in the rest of the paper. The proposed portfolio optimization problem can be solved via the dynamic programming principle which can be applied to the proposed framework:
Theorem 1.
1. The value function V t * ( S ¯ t ) satisfies the following Bellman-Isaacs equation:
V t 1 * ( S ¯ t 1 ) = sup H D t ( S ¯ t 1 ) inf Q P t ( S ¯ t 1 ) V t * S ¯ t d Q , t = 1 , , N ,
V N * ( S ¯ N ) = J ( S ¯ N ) .
2. If for each t = 1 , , N the supremum in (9) is achieved at some H t * , then H * = { H 1 * , , H N * } is an optimal strategy.
Proof. 
The proof is technical and has been moved to the Appendix A. □
The dynamic programming principle makes the problem more feasible for the numerical solution by allowing us to reconstruct the value function for each investment period through an iterative process. By combining the formulae for the risk-free position value (4) and the price dynamics (5) with (9) and (10), we obtain the final form of the Bellman-Isaacs equation for the portfolio investment problem as
V t 1 * ( S ¯ t 1 ) = sup H D t ( S ¯ t 1 ) inf Q P t ( S ¯ t 1 ) V t * S ¯ t 1 , H , Λ ( ξ ) X t 1 , W t Y ( S ¯ t 1 , H ) Q ( d ξ ) ,
V N * ( S ¯ N ) = J ( S ¯ N ) .
The numerical solution of the provided equation might prove difficult due to the additional minimization problem in the recursive formula, which is not present in the classic stochastic framework. The problem is even more relevant due to the high dimension of the state space. In practical cases, the number of risky assets covered by the investment strategy can be quite high. For example, when the portfolio replicates a market index. In the next section, we provide sufficient conditions which allow the class of measures in the infimum problem to be narrowed down to atomic measures, which should make the Bellman-Isaacs equation more feasible for practical use.
Note that the value function V t * ( · ) might attain both and + values. While the case V t * ( · ) = arises naturally given the definition of J ( · ) , V t * ( · ) = + usually means an incorrect problem statement, since the possibility of infinite reward usually implies infinite risk of the optimal strategy. Therefore, in the rest of the paper we assume that V t * ( · ) < + for each t. The assumption is satisfied for example, for the bounded reward function:
Proposition 1.
If J ( S ¯ N ) a for some a R , then V t * ( S ¯ t ) a , for each t = 0 , , N .
Proof. 
We conduct the proof by backward induction. For t = N , the statement is true by definition of the value function. Assume that V t * ( · ) a for some t N . Then inf Q P t ( S ¯ t 1 ) V t * S ¯ t d Q a and from the Bellman-Isaacs equation we get the statement for V t 1 * ( · ) . □
If the reward function is unbounded, then proving that V t * ( · ) < + would require some additional continuity assumptions, which is out of the scope of the presented research.
At the end of the current Section, we return to the question of the robust strategy in the absence of any expectations about risky prices, that is when E t ( · ) , as defined above, is also ambiguous. For example, consider a one period problem ( N = 1 ) for a single risky asset ( n = 1 ) with no costs ( C t ( · , H ) 0 ) and phase constraints of the form D t ( · ) = [ H ̲ , H ¯ ] , where H ̲ < 0 and H ¯ > 0 , so that positive, negative and zero volumes are allowed. Let R t int K t ( · ) , which means that the risky price rate can either be higher or lower than the risk-free rate R t (the assumption is natural, since otherwise there would be obvious arbitrage in the market). Assume that X t = ξ t X t 1 , where the distribution of ξ t is ambiguous, and belongs to the family of distributions P t with support K t ( · ) . Note that for N = 1 , the Rectangularity property is automatic. Assume that P t contains all the Dirac measures, that is, the admissible strategies of the “market” include the deterministic strategies. Let the initial risky price be X 0 and let the initial portfolio have H 0 volume invested in the risky asset, while its risk-free position value is W 0 Y . Let the reward J ( · ) be a strictly monotone function of the portfolio market value H N X N + W N Y . For any admissible strategy H, the market value of the portfolio at t = 1 would be
W ( H ) = H X t + R t ( W t 1 Y ( H H t 1 ) X t 1 ) = R t ( W t 1 Y + H t 1 X t 1 ) + ( ξ R t ) H X t 1 ,
where ξ K t ( · ) . Note that if H > 0 and the “market” chooses a deterministic strategy δ ξ * , where ξ * < R t , then W ( H ) < R t ( W t 1 Y + H t 1 T X t 1 ) = W ( 0 ) ; similarly, if H < 0 , then we have W ( H ) < W ( 0 ) for the deterministic “market” strategy δ ξ * with ξ * > R t . Note that W ( 0 ) does not depend on the risky price (there are no risky investments in the portfolio). Therefore, the value of the robust Savage representation for any strategy H 0 is less than its value for H = 0 , which makes H = 0 the optimal strategy. The optimality of H = 0 means that the investor should liquidate the risky positions completely. One can easily extrapolate the example to the multiperiod problem and see that H = 0 is optimal for each period, which compromises the whole idea of the investment process. This example shows that to obtain a non-trivial strategy in this case, the investor needs to have some prior expectations about the future price values.

4. Simplified Forms of the Bellman-Isaacs Equation

This research provides a practical framework for the portfolio optimization problem. Therefore, we attempt to make economically reasonable assumptions while providing the interpreted solution via the Bellman-Isaacs equation. Still, a practical approach should allow for either an analytic or numerical solution. The Equations (11) and (12) does not allow for an analytic representation of the value function in most practical cases, while solving it numerically might be computationally difficult compared to the usual Bellman equation, since we need to solve an additional minimization problem at every step. Below, we consider particular cases of the framework that allow us to simplify the numerical solution of the internal minimization problem.
Let ext K t ( · ) be the set of extreme points of K t ( · ) . Let P t n ( · ) (respectively, P t n * ( · ) ) be the set of probability measures from P t ( · ) concentrated at n + 1 or fewer points (respectively, extreme points) of K t ( · ) . Below we show that in some cases, the infimum over P t in the Bellman-Isaacs equation is equivalent to the infimum over P t n or even P t n * . For this, we use the following result:
Lemma 1.
Let f ( x ) be a real-valued Borel measurable function on convex and compact K t ( · ) R n . Then (i) the infimum of f ( x ) Q ( d x ) over the set of probability measures P t ( · ) equals the infimum over P t n ( · ) ; (ii) if f ( x ) is concave, the infimum over P t ( · ) coincides with the infimum over P t n * ( · ) .
Proof. 
The proof is technical and has been moved to the Appendix A. A similar result for a continuous f ( x ) can be found in Smirnov [24]. The idea of their proof applies to our case as well. The proof of the Lemma 1 is based on their idea and is provided in this paper for consistency. □
Lemma 1 immediately yields the following result, which allows us to narrow the set of measures in the Bellman-Isaacs equation:
Theorem 2.
Assume that V t * ( · ) < + and the sets K t ( · ) are compact and convex for each t = 1 , , N . Then
V t 1 * ( S ¯ t 1 ) = sup H D t ( S ¯ t 1 ) inf Q P t n ( S ¯ t 1 ) V t * S ¯ t 1 , H , Λ ( ξ ) X t 1 , W t Y ( S ¯ t 1 , H ) Q ( d ξ ) ,
V N * ( S ¯ N ) = J ( S ¯ N ) .
Proof. 
Assume that for some t, S ¯ t 1 and H, H D t ( S ¯ t 1 ) , the subfunction ξ V t * S ¯ t 1 , H , Λ ( ξ ) X t 1 , W t Y ( S ¯ t 1 , H ) attains for some ξ K t ( S ¯ t 1 ) , ξ E t ( S ¯ t 1 ) . Then the infimum is attained at the atomic measure concentrated at the points ξ , ξ K t ( S ¯ t 1 ) , such that E t ( S ¯ t 1 ) belongs to the interval connecting ξ and ξ . If is attained exactly at E t ( S ¯ t 1 ) , then the corresponding atomic measure is concentrated at E t ( S ¯ t 1 ) . If is not attained on K t ( S ¯ t 1 ) , then we use Lemma 1 to complete the proof. □
Lemma 1 allows us to further narrow the class of measures when finding the infimum in the Bellman-Isaacs equation. For this, the subfunctions ξ V t * S ¯ t 1 , H , Λ ( ξ ) X t 1 , W t Y ( S ¯ t 1 , H ) in (11) should be concave. Below, we will prove the concavity via backward induction by imposing restrictions on the form of the cost function, the reward function, the constraint sets and the measure sets. Under the restrictions, the dependence on the risky position volumes H would only be allowed through the dependence on the market values H X , where “∘” means the element-wise (Hadamard) product.
For further convenience, we sometimes represent the dependence on the system history S ¯ t as the dependence on the arguments S ¯ t 1 , H t , X t and W t Y , thus decomposing the dependence on the latest system state and the history before that. We also uphold this notation for t = 0 , in which case the arguments should be interpreted as H 0 , X 0 and W 0 Y , meaning that there is no history at the beginning of the strategy. We also sometimes represent the system history S ¯ t as a combination of the history of the portfolio risky volumes H ¯ t , the history of the risk-free position value W ¯ t Y and the price history X ¯ t up to time t. For example, V t * ( S ¯ t ) might be equivalently written as V t * ( H ¯ t , X ¯ t , W ¯ t Y ) or V t * ( H ¯ t 1 , X ¯ t 1 , W ¯ t 1 Y , H t , X t , W t Y ) , whenever the notation is clear from the context. In the following proofs, we also make some monotonicity assumptions, in which case the comparison of vectors and arrays of arguments should be understood element-wise.
Assume that the constraint sets D t are monotonic in the following sense: let d t ( W ¯ t 1 Y ) be any of the subfunctions W ¯ t 1 Y D t ( S ¯ t 1 ) , then
d t ( W ¯ t 1 Y , 1 ) d t ( W ¯ t 1 Y , 2 ) if W ¯ t 1 Y , 1 W ¯ t 1 Y , 2 .
The assumption means that the increase in the portfolio risk-free position should provide additional trading options while keeping the previous ones available. The assumption might be inappropriate in some cases, for example, when the increase in risk-free funds leads to additional constraints on the less risky positions in favor of the more risky ones to stimulate the more risky and profitable strategies. We do not consider such cases, believing that the investor is reasonable enough to account for the increased risk-free reserves when choosing the strategy, thus the additional trading limits are not required.
Assume that the measure sets P t are monotonic in the following sense: let p t ( W ¯ t 1 Y ) be any of the subfunctions W ¯ t 1 Y P t ( S ¯ t 1 ) , then
p t ( W ¯ t 1 Y , 1 ) p t ( W ¯ t 1 Y , 2 ) if W ¯ t 1 Y , 1 W ¯ t 1 Y , 2 .
The assumption means that the increase in the risk-free funds does not lead to more ambiguity in the market, according to the investor’s subjective estimate. Otherwise, we would have assumed that the investor becomes more uncertain when provided with the additional risk-free liquid assets, which is unreasonable.
Below, we provide sufficient conditions for the concavity and monotonicity of the value function.
Theorem 3.
Let V t * ( · ) be the value function in the Bellman-Isaacs Equations (11) and (12). Assume that for each t = 1 , , N
  • the subfunctions H ¯ t Y , W ¯ t 1 Y C t ( S ¯ t 1 , H t ) are convex;
  • the subfunctions W ¯ t 1 Y C t ( S ¯ t 1 , H t ) are non-increasing;
  • the subfunctions H ¯ N , W ¯ N Y J ( S ¯ N ) are concave;
  • the subfunctions W ¯ N Y J ( S ¯ N ) are non-decreasing;
  • the constraint sets D t have convex values and monotone in the sense of (15);
  • the measure sets P t are monotone in the sense of (16);
  • V t * ( S ¯ t ) < + .
Then the subfunctions H ¯ t , W ¯ t Y V t * ( S ¯ t ) are concave and the subfunctions W ¯ t Y V t * ( S ¯ t ) are non-decreasing for t = 1 , , N .
Proof. 
The proof is technical and has been moved to the Appendix A. □
The convexity of the cost function is a common requirement in financial literature, which is justified in most practical cases due to the well-known price impact effect in the market. The assumption that the costs are non-increasing with respect to the risk-free position value is economically reasonable, since the additional risk-free and liquid funds should create more favorable conditions when choosing the execution strategy, while incurring no additional costs themselves. The concavity of the reward function is reasonable if the investor is risk-averse, while the monotonicity with respect to the risk-free position value is natural for many practical cases, when adding to the total portfolio value at no additional costs should increase the portfolio reward.
Let W t X = H t X t be the vector of risky position values at time t, and let W ¯ t X be the corresponding history up to time t, t = 0 , , N . Let S t W = ( W t X , W t Y ) be the state of the portfolio in terms of the position market values. S t W can be derived from S t and combines the market and portfolio states. Let S ¯ t W be the corresponding state history up to time t. Next, we show how the class of measures in the Bellman-Isaacs equation can be narrowed, if the value function can be represented as a function of S ¯ t W .
Theorem 4.
Assume that all the assumptions of Theorem 3 hold and, in addition,
  • the cost function has the following form:
    C t ( S ¯ t 1 , H t ) C ˜ t ( S ¯ t 1 W , H t X t 1 ) ;
  • the reward function has the following form:
    J ( S ¯ N ) J ˜ ( S ¯ N W ) ;
  • the constraint sets D t ( · ) are given in the form of constraints on the risky position market values:
    H D t ( S ¯ t 1 ) H X t 1 D ˜ t ( S ¯ t 1 W ) ;
  • the probability measure sets P t ( · ) have the following form:
    P t ( S ¯ t 1 ) P ˜ t ( S ¯ t 1 W ) .
Then
V t 1 * ( S ¯ t 1 ) = sup H D t ( S ¯ t 1 ) inf Q P t n * ( S ¯ t 1 ) V t * S ¯ t 1 , H , Λ ( ξ ) X t 1 , W t Y ( S ¯ t 1 , H ) Q ( d ξ ) ,
V N * ( S ¯ N ) = J ( S ¯ N ) .
Proof. 
The statement can be proven by using Lemma 1, so we only need to prove the concavity of the subfunctions ξ V t * S ¯ t 1 , H , Λ ( ξ ) X t , W t Y ( S ¯ t 1 , H ) for each t = 1 , , N . First, we will prove by backward induction that the value function V t * can be represented in the form
V t * ( S ¯ t ) V ˜ t * ( S ¯ t W ) .
The statement holds for t = N due to (18). Assume that it is true for some t N . From (17), we get that W t Y ( S ¯ t 1 , H t ) can be represented as a function of S ¯ t 1 W and H t X t 1 . Indeed, from (4) we have
W t Y ( S ¯ t 1 , H t ) = R t W t 1 Y H t X t 1 + H t 1 X t 1 C t ( S ¯ t 1 , H t )
= R t W t 1 Y 1 · ( H t X t 1 ) + 1 · ( H t 1 X t 1 ) C ˜ t ( S ¯ t 1 W , H t X t 1 ) W ˜ t Y ( S ¯ t 1 W , H X t 1 ) ,
where 1 means the vector ( 1 , 1 , , 1 ) R n . Since Λ ( ξ ) is a diagonal matrix, (23) implies that
V t * S ¯ t 1 , H , Λ ( ξ ) X t 1 , W t Y ( S ¯ t 1 , H ) V ˜ t * S ¯ t 1 W , H Λ ( ξ ) X t 1 , W ˜ t Y ( S ¯ t 1 W , H X t 1 )
V ˜ t * S ¯ t 1 W , Λ ( ξ ) ( H X t 1 ) , W ˜ t Y ( S ¯ t 1 W , H X t 1 ) .
Therefore, V t * S ¯ t 1 , H , Λ ( ξ ) X t 1 , W t Y ( S ¯ t 1 , H ) is, in fact, a function of S ¯ t 1 W , H X t 1 and ξ , which we will denote as ϕ t ( S ¯ t 1 W , H X t 1 , ξ ) . Then, by substituting (19) and (20) into the Bellman-Isaacs equation, we get that
V t 1 * ( S ¯ t 1 ) = sup H X t 1 D ˜ t ( S ¯ t 1 W ) inf Q P ˜ t ( S ¯ t 1 W ) ϕ t ( S ¯ t 1 W , H X t 1 , ξ ) Q ( ξ ) .
Therefore, V t 1 * ( S ¯ t 1 ) is a function of S ¯ t 1 W , which proves the induction step, thus proving that the value function admits the representation (23) for t = 0 , , N . Then, from (23) and the diagonality of Λ ( ξ ) , we get that, for any ξ K t ( S ¯ t 1 ) ,
V t * ( S ¯ t 1 , H , Λ ( ξ ) X t 1 , W t Y ) = V ˜ t * ( S ¯ t 1 W , H Λ ( ξ ) X t 1 , W t Y ) = V ˜ t * ( S ¯ t 1 W , ( Λ ( ξ ) H ) X t 1 , W t Y ) = V t * ( S ¯ t 1 , Λ ( ξ ) H , X t 1 , W t Y ) .
By Theorem 3, the subfunctions H V t * ( S ¯ t 1 , H , X t , W t Y ) are concave. Therefore, the subfunctions ξ V t * S ¯ t 1 , Λ ( ξ ) H , X t , W t Y are concave as compositions of the linear function ξ Λ ( ξ ) H and the concave function. Then, (25) implies that the subfunctions ξ V t * S ¯ t 1 , H , Λ ( ξ ) X t , W t Y are concave as well, which means the concavity of the subfunctions ξ V t * S ¯ t 1 , H , Λ ( ξ ) X t , W t Y ( S ¯ t 1 , H ) and proves the statement. □
The conditions of Theorem 4 might seem too restrictive. However, they are satisfied in a variety of practical cases. Below, we provide some examples covered by the Theorem.
The form of the cost function (17) includes the widely-used affine model of costs, see Reference [5]. At the same time, it allows us to introduce the non-linear cost growth, which might be relevant for large portfolios.
As the example of the concave and monotonic reward function of the form (18), consider a utility function which depends on the portfolio value, for example,
J ˜ ( S ¯ N W ) = U ( H N X N + W N Y )
or
J ˜ ( S ¯ N W ) = U ( H N X N + W N Y C N + 1 ( S ¯ N W ) ) ,
where U is a concave non-decreasing utility function and C N + 1 ( S ¯ N ) is the costs of liquidating the remaining risky assets in the portfolio after the end of the strategy (assume that H ¯ N , W ¯ N Y C N + 1 ( S ¯ N ) are convex and W ¯ N Y C N + 1 ( S ¯ N ) are non-increasing). In the first case, the investor wants to maximize the market value of the portfolio, and in the latter case, the liquidation value of the portfolio is being maximized.
As the example of the constraint sets covered by Theorem 4, consider a set of predictable functions g t i ( · , x ) for t = 1 , , N and i = 1 , , M , such that the corresponding subfunctions x g t i ( · , x ) are quasi-convex. Let β t i ( · ) 0 be predictable coefficients and let W t 1 = H t 1 X t 1 + W t 1 Y be the market value of the portfolio at t 1 . Consider the constraint sets of the form:
D t ( · ) = H : g t i ( · , H X t 1 ) β t i ( · ) W t 1 , i = 1 , , M .
Quasi-convexity of g t i implies the convexity of D t ( · ) values, while non-negativity of β t i ( · ) implies the monotonicity as required by Theorem 4. The introduced form of constraints includes several important trading limit types, for example:
  • g t i ( · , H X t 1 ) = H i X t 1 i , i = 1 , , n , means the constraints on the maximum value of each risky position in terms of the percentage of the portfolio market value;
  • g t i ( · , H X t 1 ) = H i X t 1 i , i = 1 , , n , means the constraints on the minimum value of each risky position, which limits short-selling;
  • g t i ( · , H X t 1 ) = j = 1 n a t i , j ( · ) H j X t 1 j means the constraints on the combinations of asset positions, which is useful for limiting investments in a group of products, which, for example, represent a sector of economy or constitute an index; a special case of this type of constraints is g t i ( · , H X t 1 ) = j = 1 n H j X t 1 j , which might be used to limit the risk-free asset short-selling.
Theorem 4 covers several price impact effects. For example, consider the measure sets P t with the expected value E t ( · ) = E t 1 ( · ) , , E t n ( · ) defined as
ln E t i ( · ) = k = 1 t 1 μ k i + k = 1 t 1 λ k i ( Δ H k X k ) + k = 1 t 1 γ k i ( Δ H k X k ) e ρ k i ( t k ) ,
where μ k i , λ k i ( w ) , γ k i ( w ) and ρ k i are assumed to be known, i = 1 , , n . For the i-th risky asset and time k, μ k i is the expected drift of the market price, λ k i ( w ) is the expected permanent price impact of the trade with market value w, γ k i ( w ) is the expected temporary price impact of the trade with market value w, and ρ k i characterizes the expected resilience rate of the market price. The logarithm in the formula guarantees that the expected value of the multiplicative coefficient of the price dynamics is positive. Let the support of the measure set be defined as
K t ( · ) = { E t ( · ) } + M t ,
where the “+” operator means the Minkowski sum, and M t characterizes both the uncertainty and the ambiguity of the market prices. Assume that 0 int M t to guarantee that E t ( · ) K t ( · ) , t = 1 , , N . The provided forms of E t ( · ) and K t ( · ) capture several aspects of the price impact which are relevant for portfolio optimization in a market with limited liquidity, see for example, References [8,11,25,26]. The structure of M t characterizes the investor’s estimate of the market ambiguity. For example, the investor might assume the rectangular structure to disregard any prior knowledge about the dependencies between the asset prices. On the other hand, by considering the elliptical form of M t , the investor assumes that the prices will not attain the respective extreme values all at the same time, thus introducing some prior knowledge about the market into the framework.

5. Conclusions

In this paper, we present a robust approach to a general portfolio optimization problem. We introduce the ambiguity in the market by formulating the problem as a game against nature, where the investor attempts to maximize the worst-case reward from the portfolio at the end of the strategy. The robust approach considers the model error when selecting an optimal portfolio strategy, which is relevant when the parameters of the model cannot be reliably estimated, or when the model-based forecasts are untrustworthy, for example, during a financial crisis. The framework might be used as a decision support system when comparing portfolio strategies based on the model-independent estimates about the future market price values. Our goal is to present a practical approach, therefore all the assumptions are economically reasonable and well-interpreted.
Calculation of the value function in the Bellman-Isaacs Equations (11) and (12) involves a complex problem of optimizing over a set of distributions with the common mean and support. A direct solution of the problem would have required some approximation and parametrization of the measure set which would have raised the question of approximation error. Theorem 2 shows that such approximation is unnecessary and the problem can be reduced to optimization over a finite set of points of the support, thus can be solved numerically similar to the classic Bellman equation. Indeed, given the n + 1 atoms of the distribution in R n , the probabilities of the atoms can be calculated directly as the corresponding barycentric coordinates of the distribution mean. Due to the recursive nature of the Equation (11), a derivative-free method, for example, a random search method, might be more appropriate for estimating the value of the robust Savage representation inf Q P t ( S ¯ t 1 ) V t * ( S ¯ t ) d Q , and for reconstructing the value function. Theorem 4 states that under some assumptions, the search space for the method can be reduced, since the atoms of the optimal measure are guaranteed to be the extreme points of the support.
The research covers a wide range of problem statements which might include non-linear transaction costs and price impact, and various forms of trading limits. The provided framework can be applied for such classic portfolio optimization problems as the optimal execution of trades and the portfolio selection, and might be used as an analytic tool for the strategy evaluation and risk-management purposes in an ambiguous market.

Funding

Support from the Basic Research Program of the National Research University Higher School of Economics is gratefully acknowledged.

Acknowledgments

The author would like to thank Sergey Smirnov, Faculty of Computational Mathematics and Cybernetics, Lomonosov Moscow State University, for helpful ideas and insight.

Conflicts of Interest

The author declares no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

Appendix A

Proof of Theorem 1
Proof. 
We begin with the proof of the first part of the Theorem. Equations (9) and (10) for t = N are true by definition of the value function, thus we only need to prove (9) for t = 1 , , N 1 in case N > 1 . For any t = 1 , , N 1 and any tail strategy h t , N H t , N ( S ¯ t 1 ) over the periods starting at t 1 or later, let
V t 1 ( S ¯ t 1 , h t , N ) = inf Q t , N P t , N ( S ¯ t 1 ) J ( S ¯ N ) d Q t , N ,
where S ¯ N is the history for strategy h t , N . Let H t be the part of strategy h t , N which covers the period starting at t 1 , and let h t + 1 , N be the rest of h t , N . Then for t < N , due to the Rectangularity property of P t , N , we have
V t 1 ( S ¯ t 1 , h t , N ) = inf Q t , N P t , N ( S ¯ t 1 ) J ( S ¯ N ) d Q t , N
= inf Q P t ( S ¯ t 1 ) d Q inf Q t + 1 , N P t + 1 , N ( S ¯ t ) J ( S ¯ N ) d Q t + 1 , N = inf Q P t ( S ¯ t 1 ) V t ( S ¯ t , h t + 1 , N ) d Q .
By definition (8), V t 1 * ( S ¯ t 1 ) = sup h t , N H t , N V t 1 ( S ¯ t 1 , h t , N ) . Therefore,
V t 1 * ( S ¯ t 1 ) = sup h t , N H t , N ( S ¯ t 1 ) inf Q P t ( S ¯ t 1 ) V t ( S ¯ t , h t + 1 , N ) d Q
sup h t , N H t , N ( S ¯ t 1 ) inf Q P t ( S ¯ t 1 ) V t * ( S ¯ t ) d Q = sup H t D t ( S ¯ t 1 ) inf Q P t ( S ¯ t 1 ) V t * ( S ¯ t ) d Q ,
because V t * ( S ¯ t ) does not depend on h t + 1 , N . On the other hand, since V t * ( S ¯ t ) is the supremum of V t ( S ¯ t , h t + 1 , N ) on H t + 1 , N ( S ¯ t ) , then for any ε > 0 there exists h t + 1 , N ε H t + 1 , N ( S ¯ t ) such that V t ( S ¯ t , h t + 1 , N ε ) V t * ( S ¯ t ) ε . Then we have
V t 1 * ( S ¯ t 1 ) = sup h t , N H t , N ( S ¯ t 1 ) inf Q P t ( S ¯ t 1 ) V t ( S ¯ t , h t + 1 , N ) d Q
= sup H t D t ( S ¯ t 1 ) sup h t + 1 , N H t + 1 , N ( S ¯ t ) inf Q P t ( S ¯ t 1 ) V t ( S ¯ t , h t + 1 , N ) d Q
sup H t D t ( S ¯ t 1 ) inf Q P t ( S ¯ t 1 ) V t ( S ¯ t , h t + 1 , N ε ) d Q
sup H t D t ( S ¯ t 1 ) inf Q P t ( S ¯ t 1 ) V t * ( S ¯ t ) d Q ε .
Since (A3) is true for any ε > 0 , from (A2) and (A3) we get Equation (9).
To prove the second part of the Theorem, note that H * is optimal iff V 0 ( S ¯ 0 , H * ) = V 0 * ( S ¯ 0 ) , which follows from the definition of the value function. Let h t , N * be the tail strategies of H * for t = 1 , , N . Next, we will prove that
V t 1 ( S ¯ t 1 , h t , N * ) = V t 1 * ( S ¯ t 1 ) , t = 1 , , N ,
by backward induction. Since h 1 , N * H * by definition, it will prove the optimality of H * . For t = N , (A4) follows from (9) and the definition of V t 1 ( S ¯ t 1 , h t , N ) . Assume that V t ( S ¯ t , h t + 1 , N * ) = V t * ( S ¯ t ) for some t N . Let S ¯ t * be the combination of the history S ¯ t 1 and the state S t * , induced by H t * . Then from (A1) we get that
V t 1 ( S ¯ t 1 , h t , N * ) = inf Q P t ( S ¯ t 1 ) V t ( S ¯ t * , h t + 1 , N * ) d Q = inf Q P t ( S ¯ t 1 ) V t * ( S ¯ t * ) d Q
= sup H t D t ( S ¯ t 1 ) inf Q P t ( S ¯ t 1 ) V t * ( S ¯ t ) d Q = V t 1 * ( S ¯ t 1 ) ,
which proves the induction step and, consequently, the second part of the Theorem. □
Proof of Lemma 1
Proof. 
The first statement follows from Theorem 3.2 of Reference [27] and implies that the infimum of f ( x ) Q ( d x ) over P t ( · ) equals the infimum of i = 1 n + 1 p i f ( x i ) over x i K t ( · ) and p i 0 , i = 1 , , n + 1 , such that i = 1 n + 1 p i x i = E t ( · ) and i = 1 n + 1 p i = 1 . The second part can be proven by using the first statement. From the Krein-Milman theorem [28] and the concavity of f ( x ) , we get that the infimum over measures concentrated at n + 1 or fewer points of K t ( · ) is greater or equal than the infimum over measures concentrated at ( n + 1 ) 2 or fewer extreme points of K t ( · ) , that is the infimum of i = 1 ( n + 1 ) 2 p i f ( x i ) over x i ext K t ( · ) and p i 0 , i = 1 , , ( n + 1 ) 2 , such that i = 1 ( n + 1 ) 2 p i x i = E t ( · ) and i = 1 ( n + 1 ) 2 p i = 1 . For any given x 1 , , x ( n + 1 ) 2 ext K t ( · ) , this problem is a linear programming problem in the standard form, which is known to have a basic feasible solution p 1 * , , p ( n + 1 ) 2 * with n + 1 or fewer nonzero elements. Therefore, the infimum actually coincides with the infimum over measures concentrated at n + 1 or fewer extreme points of K t ( · ) , which proves the result. □
Proof of Theorem 3
Proof. 
The statement is proven by backward induction. For t = N , it holds by the assumptions about J ( · ) . Let the subfunctions H ¯ t , W ¯ t Y V t * ( S ¯ t ) be concave and let the subfunctions W ¯ t Y V t * ( S ¯ t ) be non-decreasing for some t N .
Consider the function
ν t ( H ¯ t 1 , X ¯ t 1 , W ¯ t 1 Y , H t , ξ ) = V t * H ¯ t 1 , X ¯ t 1 , W ¯ t 1 Y , H t , Λ ( ξ ) X t 1 , W t Y ( H ¯ t 1 , X ¯ t 1 , W ¯ t 1 Y , H t ) .
Note that the subfunctions H ¯ t , W ¯ t 1 Y W t Y ( S ¯ t 1 , H t ) are concave and the subfunctions W ¯ t 1 Y W t Y ( S ¯ t 1 , H t ) are non-decreasing due to the assumptions about C t ( S ¯ t 1 , H t ) and the definition (4) of W t Y ( S ¯ t 1 , H ) . This, when combined with the concavity and monotonicity of V t * , implies that the subfunctions H ¯ t , W ¯ t 1 Y ν t ( H ¯ t 1 , X ¯ t 1 , W ¯ t 1 Y , H t , ξ ) are concave and that the subfunctions W ¯ t 1 Y ν t ( H ¯ t 1 , X ¯ t 1 , W ¯ t 1 Y , H t , ξ ) are non-decreasing. Therefore, the subfunctions H ¯ t , W ¯ t 1 Y ν t ( H ¯ t 1 , X ¯ t 1 , W ¯ t 1 Y , H t , ξ ) Q ( d ξ ) are concave and the subfunctions W ¯ t 1 Y ν t ( H ¯ t 1 , X ¯ t 1 , W ¯ t 1 Y , H t , ξ ) Q ( d ξ ) are non-decreasing by the properties of the integral. Consider
ϕ t ( S ¯ t 1 , H t ) = inf Q P t ( S ¯ t 1 ) ν t ( H ¯ t 1 , X ¯ t 1 , W ¯ t 1 Y , H t , ξ ) Q ( d ξ ) .
The subfunctions H ¯ t , W ¯ t 1 Y ϕ t ( S ¯ t 1 , H t ) are concave as the infima of the concave functions, and the subfunctions W ¯ t 1 Y ϕ t ( S ¯ t 1 , H t ) are non-decreasing due to the monotonicity of P t and the properties of infimum.
By definition, V t 1 * ( S ¯ t 1 ) = sup H t D t ( S ¯ t 1 ) ϕ t ( S ¯ t 1 , H t ) . For any admissible price history X ¯ t 1 and any ε > 0 , consider the admissible risky volume histories H ¯ t 1 1 , H ¯ t 1 2 and risk-free value histories W ¯ t 1 Y , 1 , W ¯ t 1 Y , 2 . By the definition of supremum, there exist H t 1 D t ( H ¯ t 1 1 , X ¯ t 1 , W ¯ t 1 Y , 1 ) and H t 2 D t ( H ¯ t 1 2 , X ¯ t 1 , W ¯ t 1 Y , 2 ) such that
ϕ t ( H ¯ t 1 1 , X ¯ t 1 , W ¯ t 1 Y , 1 , H t 1 ) V t 1 * ( H ¯ t 1 1 , X ¯ t 1 , W ¯ t 1 Y , 1 ) ε ,
ϕ t ( H ¯ t 1 2 , X ¯ t 1 , W ¯ t 1 Y , 2 , H t 2 ) V t 1 * ( H ¯ t 1 2 , X ¯ t 1 , W ¯ t 1 Y , 2 ) ε .
Consider some α 1 , α 2 0 , α 1 + α 2 = 1 . From the concavity of ϕ t , we get that
V t 1 * ( α 1 H ¯ t 1 1 + α 2 H ¯ t 1 2 , X ¯ t 1 , α 1 W ¯ t 1 Y , 1 + α 2 W ¯ t 1 Y , 2 )
ϕ ( α 1 H ¯ t 1 1 + α 2 H ¯ t 1 2 , X ¯ t 1 , α 1 W ¯ t 1 Y , 1 + α 2 W ¯ t 1 Y , 2 , α 1 H t 1 + α 2 H t 2 )
α 1 ϕ ( H ¯ t 1 1 , X ¯ t 1 , W ¯ t 1 Y , 1 , H t 1 ) + α 1 ϕ ( H ¯ t 1 2 , X ¯ t 1 , W ¯ t 1 Y , 2 , H t 2 )
α 1 V t 1 * ( H ¯ t 1 1 , X ¯ t 1 , W ¯ t 1 Y , 1 ) + α 2 V t 1 * ( H ¯ t 1 2 , X ¯ t 1 , W ¯ t 1 Y , 2 ) ε .
Since the above inequality holds for any ε > 0 , the subfunctions H ¯ t 1 , W ¯ t 1 Y V t 1 * ( S ¯ t 1 ) are concave.
To prove the monotonicity of W ¯ t Y V t * ( S ¯ t ) , for any admissible price history X ¯ t 1 and risky volume history H ¯ t 1 , consider the admissible risk-free value histories W ¯ t 1 Y , 1 and W ¯ t 1 Y , 2 , such that W ¯ t 1 Y , 1 W ¯ t 1 Y , 2 . Then, D t ( H ¯ t 1 , X ¯ t 1 , W ¯ t 1 Y , 1 ) D t ( H ¯ t 1 , X ¯ t 1 , W ¯ t 1 Y , 2 ) due to the monotonicity of D t ( · ) , and we use the monotonicity of W ¯ t 1 Y ϕ t ( S ¯ t 1 , H t ) to get that
V t 1 * ( H ¯ t 1 , X ¯ t 1 , W ¯ t 1 Y , 1 ) = sup H t D t ( H ¯ t 1 , X ¯ t 1 , W ¯ t 1 Y , 1 ) ϕ t ( H ¯ t 1 , X ¯ t 1 , W ¯ t 1 Y , 1 , H t )
sup H t D t ( H ¯ t 1 , X ¯ t 1 , W ¯ t 1 Y , 2 ) ϕ t ( H ¯ t 1 , X ¯ t 1 , W ¯ t 1 Y , 2 , H t ) = V t 1 * ( H ¯ t 1 , X ¯ t 1 , W ¯ t 1 Y , 2 ) ,
which proves that the subfunctions W ¯ t 1 Y V t 1 * ( S ¯ t 1 ) are non-decreasing. This concludes the proof of the induction step and proves the statement. □

References

  1. Samuelson, P. Lifetime portfolio selection by dynamic stochastic programming. Rev. Econ. Stat. 1969, 51, 239–246. [Google Scholar] [CrossRef]
  2. Merton, R. Lifetime portfolio selection under uncertainty: The continuous-time case. Rev. Econ. Stat. 1969, 51, 247–257. [Google Scholar] [CrossRef]
  3. Karatzas, I.; Shreve, S. Methods of Mathematical Finance; Springer: Berlin/Heidelberg, Germany, 1998; p. 415. [Google Scholar]
  4. Magill, M.; Constantinides, G. Portfolio selection and transactions costs. J. Econ. Theory 1976, 13, 245–263. [Google Scholar] [CrossRef]
  5. Davis, M.H.A.; Norman, A.R. Portfolio selection with transaction costs. Math. Oper. Res. 1990, 15, 676–713. [Google Scholar] [CrossRef]
  6. Kabanov, Y.; Safarian, M. Markets with Transaction Costs: Mathematical Theory; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2009. [Google Scholar]
  7. Bertsimas, D.; Lo, A.W. Optimal control of execution costs. J. Financ. Mark. 1998, 1, 1–50. [Google Scholar] [CrossRef]
  8. Almgren, R.; Chriss, N. Value under liquidation. Risk 1999, 12, 61–63. [Google Scholar]
  9. Almgren, R. Optimal execution with nonlinear impact functions and trading-enhanced risk. Appl. Math. Financ. 2003, 10, 18. [Google Scholar] [CrossRef]
  10. Jorion, P. Risk management lessons from Long-Term Capital Management. Eur. Financ. Manag. 2000, 6, 277–300. [Google Scholar] [CrossRef]
  11. Obizhaeva, A.; Wang, J. Optimal trading strategy and supply/demand dynamics. J. Financ. Mark. 2013, 16, 1–32. [Google Scholar] [CrossRef]
  12. Framstad, N.; ∅Ksendal, B.; Sulem, A. Optimal consumption and portfolio in a jump diffusion market with proportional transaction costs. J. Math. Econ. 2001, 35, 233–257. [Google Scholar] [CrossRef]
  13. Smirnov, S.N. Guaranteed deterministic approach to superhedging: Market model, trading constraints and Bellman-Isaacs equations. Math. Game Theory Appl. 2018, 10, 59–99. [Google Scholar]
  14. Smirnov, S.N. Guaranteed deterministic approach to superhedging: “No arbitrage” properties of the market. Math. Game Theory Appl. 2019, 11, 68–95. [Google Scholar]
  15. Deng, X.T.; Li, Z.F.; Wang, S.Y. A minimax portfolio selection strategy with equilibrium. Eur. J. Oper. Res. 2005, 166, 278–292. [Google Scholar] [CrossRef]
  16. Ben-Tal, A.; Margalit, T.; Nemirovski, A. Robust modeling of multi-stage portfolio problems. In High Performance Optimization; Springer: Berlin/Heidelberg, Germany, 2000; pp. 303–328. [Google Scholar]
  17. Talay, D.; Zheng, Z. Worst case model risk management. Financ. Stoch. 2002, 6, 517–537. [Google Scholar] [CrossRef]
  18. Glasserman, P.; Xu, X. Robust risk measurement and model risk. Quant. Financ. 2014, 14, 29–58. [Google Scholar] [CrossRef]
  19. Ellsberg, D. Risk, ambiguity, and the Savage axioms. Q. J. Econ. 1961, 75, 643–669. [Google Scholar] [CrossRef]
  20. Föllmer, H.; Schied, A. Stochastic Finance: An Introduction in Discrete Time; Walter de Gruyter: Berlin, Germany, 2004. [Google Scholar]
  21. Epstein, L.G.; Schneider, M. Recursive multiple-priors. J. Econ. Theory 2003, 113, 1–31. [Google Scholar] [CrossRef]
  22. Iyengar, G.N. Robust Dynamic Programming. Math. Oper. Res. 2005, 30, 257–280. [Google Scholar] [CrossRef]
  23. González-Trejo, J.; Hernández-Lerma, O.; Hoyos-Reyes, L.F. Minimax control of discrete-time stochastic systems. SIAM J. Control Optim. 2002, 41, 1626–1659. [Google Scholar] [CrossRef]
  24. Smirnov, S.N. A Guaranteed Deterministic Approach to Superhedging: The Case of Convex Payoff Functions on Options. Mathematics 2019. under review. [Google Scholar]
  25. Schied, A.; Schöneborn, T. Risk aversion and the dynamics of optimal liquidation strategies in illiquid markets. Financ. Stoch. 2009, 13, 181–204. [Google Scholar] [CrossRef]
  26. Vath, V.L.; Mnif, M.; Pham, H. A model of optimal portfolio selection under liquidity risk and price impact. Financ. Stoch. 2007, 11, 51–90. [Google Scholar] [CrossRef]
  27. Winkler, G. Extreme points of moment sets. Math. Oper. Res. 1988, 13, 581–587. [Google Scholar] [CrossRef]
  28. Krein, M.; Milman, D. On extreme points of regular convex sets. Stud. Math. 1940, 9, 133–138. [Google Scholar] [CrossRef] [Green Version]

Share and Cite

MDPI and ACS Style

Andreev, N. Robust Portfolio Optimization in an Illiquid Market in Discrete-Time. Mathematics 2019, 7, 1147. https://doi.org/10.3390/math7121147

AMA Style

Andreev N. Robust Portfolio Optimization in an Illiquid Market in Discrete-Time. Mathematics. 2019; 7(12):1147. https://doi.org/10.3390/math7121147

Chicago/Turabian Style

Andreev, Nikolay. 2019. "Robust Portfolio Optimization in an Illiquid Market in Discrete-Time" Mathematics 7, no. 12: 1147. https://doi.org/10.3390/math7121147

APA Style

Andreev, N. (2019). Robust Portfolio Optimization in an Illiquid Market in Discrete-Time. Mathematics, 7(12), 1147. https://doi.org/10.3390/math7121147

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop