The optimal method for pricing Bermudan options by simulation

Least‐squares methods enable us to price Bermudan‐style options by Monte Carlo simulation. They are based on estimating the option continuation value by least‐squares. We show that the Bermudan price is maximized when this continuation value is estimated near the exercise boundary, which is equivalent to implicitly estimating the optimal exercise boundary by using the value‐matching condition. Localization is the key difference with respect to global regression methods, but is fundamental for optimal exercise decisions and requires estimation of the continuation value by iterating local least‐squares (because we estimate and localize the exercise boundary at the same time). In the numerical example, in agreement with this optimality, the new prices or lower bounds (i) improve upon the prices reported by other methods and (ii) are very close to the associated dual upper bounds. We also study the method's convergence.     

A FIRST‐ORDER BSPDE FOR SWING OPTION PRICING: CLASSICAL SOLUTIONS

In a companion paper, we studied a control problem related to swing option pricing in a general non‐Markovian setting. The main result there shows that the value process of this control problem can uniquely be characterized in terms of a first‐order backward stochastic partial differential equation (BSPDE) and a pathwise differential inclusion. In this paper, we additionally assume that the cash flow process of the swing option is left‐continuous in expectation. Under this assumption, we show that the value process is continuously differentiable in the space variable that represents the volume in which the holder of the option can still exercise until maturity. This gives rise to an existence and uniqueness result for the corresponding BSPDE in a classical sense. We also explicitly represent the space derivative of the value process in terms of a nonstandard optimal stopping problem over a subset of predictable stopping times. This representation can be applied to derive a dual minimization problem in terms of martingales.     

A Markov regime‐switching ARMA approach for hedging stock indices

This study considers the hedging effectiveness of applying the N‐state Markov regime‐switching autoregressive moving‐average (MRS‐ARMA) model to the S&P‐500 and FTSE‐100 markets. The distinguishingfeature of this study is to incorporate the observations of serially correlated stockreturns into the hedging analysis. To resolve the problem of N^T possible routes induced by the presence of MA parameters associated with the algorithm of Hamilton JD ( 1989 ) and a sample of size T, we propose an algorithm by combining the ideas of Hamilton JD ( 1989 ) and Gray SF ( 1996 ). We find that the hedging performances of the three proposed MRS‐MA(1) strategies herein are superior to their corresponding MRS counterparts considered in Alizadeh A and Nomikos N ( 2004 ) over the out‐of‐sample periods, even when we realistically track the transaction costs generated from rebalancing the hedged portfolios. © 2010 Wiley Periodicals, Inc. Jrl Fut Mark 31:165–191, 2011     

A copula‐based regime‐switching GARCH model for optimal futures hedging

The article develops a regime‐switching Gumbel–Clayton (RSGC) copula GARCH model for optimal futures hedging. There are three major contributions of RSGC. First, the dependence of spot and futures return series in RSGC is modeled using switching copula instead of assuming bivariate normality. Second, RSGC adopts an independent switching Generalized Autoregressive Conditional Heteroscedasticity (GARCH) process to avoid the path‐dependency problem. Third, based on the assumption of independent switching, a formula is derived for calculating the minimum variance hedge ratio. Empirical investigation in agricultural commodity markets reveals that RSGC provides good out‐of‐sample hedging effectiveness, illustrating importance of modeling regime shift and asymmetric dependence for futures hedging. © 2009 Wiley Periodicals, Inc. Jrl Fut Mark 29:946–972, 2009     

Double continuation regions for American and Swing options with negative discount rate in Lévy models

In this paper, we study perpetual American call and put options in an exponential Lévy model. We consider a negative effective discount rate that arises in a number of financial applications including stock loans and real options, where the strike price can potentially grow at a higher rate than the original discount factor. We show that in this case a double continuation region arises and we identify the two critical prices. We also generalize this result to multiple stopping problems of Swing type, that is, when successive exercise opportunities are separated by i.i.d. random refraction times. We conduct an extensive numerical analysis for the Black–Scholes model and the jump‐diffusion model with exponentially distributed jumps.     

A METHOD FOR PRICING AMERICAN OPTIONS USING SEMI‐INFINITE LINEAR PROGRAMMING

We introduce a new approach for the numerical pricing of American options. The main idea is to choose a finite number of suitable excessive functions (randomly) and to find the smallest majorant of the gain function in the span of these functions. The resulting problem is a linear semi‐infinite programming problem, that can be solved using standard algorithms. This leads to good upper bounds for the original problem. For our algorithms no discretization of space and time and no simulation is necessary. Furthermore it is applicable even for high‐dimensional problems. The algorithm provides an approximation of the value not only for one starting point, but for the complete value function on the continuation set, so that the optimal exercise region and, for example, the Greeks can be calculated. We apply the algorithm to (one‐ and) multidimensional diffusions and show it to be fast and accurate.     

CONVERGENCE OF A LEAST‐SQUARES MONTE CARLO ALGORITHM FOR AMERICAN OPTION PRICING WITH DEPENDENT SAMPLE DATA

We analyze the convergence of the Longstaff–Schwartz algorithm relying on only a single set of independent Monte Carlo sample paths that is repeatedly reused for all exercise time‐steps. We prove new estimates on the stochastic component of the error of this algorithm whenever the approximation architecture is any uniformly bounded set of L² functions of finite Vapnik–Chervonenkis dimension (VC‐dimension), but in particular need not necessarily be either convex or closed. We also establish new overall error estimates, incorporating bounds on the approximation error as well, for certain nonlinear, nonconvex sets of neural networks.     

A FIRST‐ORDER BSPDE FOR SWING OPTION PRICING

We study an optimal control problem related to swing option pricing in a general non‐Markovian setting in continuous time. As a main result we uniquely characterize the value process in terms of a first‐order nonlinear backward stochastic partial differential equation and a differential inclusion. Based on this result we also determine the set of optimal controls and derive a dual minimization problem.     

Optimal equilibria for time‐inconsistent stopping problems in continuous time

For an infinite‐horizon continuous‐time optimal stopping problem under nonexponential discounting, we look for an optimal equilibrium, which generates larger values than any other equilibrium does on the entire state space. When the discount function is log subadditive and the state process is one‐dimensional, an optimal equilibrium is constructed in a specific form, under appropriate regularity and integrability conditions. Although there may exist other optimal equilibria, we show that they can differ from the constructed one in very limited ways. This leads to a sufficient condition for the uniqueness of optimal equilibria, up to some closedness condition. To illustrate our theoretic results, a comprehensive analysis is carried out for three specific stopping problems, concerning asset liquidation and real options valuation. For each one of them, an optimal equilibrium is characterized through an explicit formula.     

Efficient trinomial trees for local-volatility models in pricing double-barrier options

A local-volatility (LV) model captures the volatility smile while retaining the preference freedom of the Black–Scholes model. Past attempts to construct a smile-consistent tree for the LV surface do not guarantee validity. This paper presents an efficient and valid smile-consistent tree for the LV model. The only assumption is that the LV surface be upper- and lower-bounded. With this tree, double-barrier options can be priced with fast convergence even in the presence of volatility smile. This is confirmed numerically. An implied tree is also presented. It recovers the LV surface reasonably well.     

PRICING OF HIGH‐DIMENSIONAL AMERICAN OPTIONS BY NEURAL NETWORKS

Pricing of American options in discrete time is considered, where the option is allowed to be based on several underlyings. It is assumed that the price processes of the underlyings are given Markov processes. We use the Monte Carlo approach to generate artificial sample paths of these price processes, and then we use the least squares neural networks regression estimates to estimate from this data the so‐called continuation values, which are defined as mean values of the American options for given values of the underlyings at time t subject to the constraint that the options are not exercised at time t. Results concerning consistency and rate of convergence of the estimates are presented, and the pricing of American options is illustrated by simulated data.     

The valuation of American options in a multidimensional exponential Lévy model

We consider the problem of valuation of American options written on dividend‐paying assets whose price dynamics follow a multidimensional exponential Lévy model. We carefully examine the relation between the option prices, related partial integro‐differential variational inequalities, and reflected backward stochastic differential equations. In particular, we prove regularity results for the value function and obtain the early exercise premium formula for a broad class of payoff functions.     

Equilibrium concepts for time‐inconsistent stopping problems in continuous time

A new notion of equilibrium, which we call strong equilibrium, is introduced for time‐inconsistent stopping problems in continuous time. Compared to the existing notions introduced in Huang, Y.‐J., & Nguyen‐Huu, A. (2018, Jan 01). Time‐consistent stopping under decreasing impatience. Finance and Stochastics, 22(1), 69–95 and Christensen, S., & Lindensjö, K. (2018). On finding equilibrium stopping times for time‐inconsistent markovian problems. SIAM Journal on Control and Optimization, 56(6), 4228–4255, which in this paper are called mild equilibrium and weak equilibrium, respectively, a strong equilibrium captures the idea of subgame perfect Nash equilibrium more accurately. When the state process is a continuous‐time Markov chain and the discount function is log subadditive, we show that an optimal mild equilibrium is always a strong equilibrium. Moreover, we provide a new iteration method that can directly construct an optimal mild equilibrium and thus also prove its existence.     

INTERNATIONAL RESERVE MANAGEMENT: A DRIFT‐SWITCHING REFLECTED JUMP‐DIFFUSION MODEL

We study the cost of shocks, that is, jump risk, with respect to reserve management when the reserve process is formulated as a drift‐switching jump diffusion with a reflecting barrier at 0. Inspired by the Brownian drift switching model, our model results in a more realistic dynamic behavior of international reserves than the buffer stock model. The new model can capture both the jump behavior in reserve dynamics and the leptokurtic feature of the increment distribution which has a higher peak and two asymmetric heavier tails than the normal distribution. Through the selection of an initial distribution that reflects certain steady state behaviors, the reserve process becomes a regenerative process. This selection enables us to derive a closed‐form expression for the total expected discounted cost of managing reserves, thus helping us to numerically find management strategies that minimize costs. The numerical results show that shocks at the reserve level have a significant effect on reserve management strategies and that model misspecification can result in nonnegligible additional costs.     

Brand switching in clothing: the role of variety‐seeking drive and product category‐level characteristics

The study explores the impact of internal factors such as variety‐seeking drive, loyalty, involvement, perceived risk, pleasure and similarity perceptions on brand‐switching propensity in clothing purchases using multiple regression. The model derived shows that five internal factors including involvement, perceived risk, loyalty and variety‐seeking drive have an impact on brand switching propensity. The findings highlight that highly involved and loyal consumers are likely to show brand‐switching propensities in their clothing purchases, thus buying from a repertoire of brands and stores. These findings imply that retailers should ensure that their clothing brands or stores are among those in the multi‐brand repertoire of the consumer.     

Azéma martingales for Bessel and CIR processes and the pricing of Parisian zero‐coupon bonds

In this paper, we study the excursions of Bessel and Cox–Ingersoll–Ross (CIR) processes with dimensions . We obtain densities for the last passage times and meanders of the processes. Using these results, we prove a variation of the Azéma martingale for the Bessel and CIR processes based on excursion theory. Furthermore, we study their Parisian excursions, and generalize previous results on the Parisian stopping time of Brownian motion to that of the Bessel and CIR processes. We obtain explicit formulas and asymptotic results for the densities of the Parisian stopping times, and develop exact simulation algorithms to sample the Parisian stopping times of Bessel and CIR processes. We introduce a new type of bond, the zero‐coupon Parisian bond. The buyer of such a bond is betting against zero interest rates, while the seller is effectively hedging against a period where interest rates fluctuate around 0. Using our results, we propose two methods for pricing these bonds and provide numerical examples.     

The robust pricing–hedging duality for American options in discrete time financial markets

We investigate the pricing–hedging duality for American options in discrete time financial models where some assets are traded dynamically and others, for example, a family of European options, only statically. In the first part of the paper, we consider an abstract setting, which includes the classical case with a fixed reference probability measure as well as the robust framework with a nondominated family of probability measures. Our first insight is that, by considering an enlargement of the space, we can see American options as European options and recover the pricing–hedging duality, which may fail in the original formulation. This can be seen as a weak formulation of the original problem. Our second insight is that a duality gap arises from the lack of dynamic consistency, and hence that a different enlargement, which reintroduces dynamic consistency is sufficient to recover the pricing–hedging duality: It is enough to consider fictitious extensions of the market in which all the assets are traded dynamically. In the second part of the paper, we study two important examples of the robust framework: the setup of Bouchard and Nutz and the martingale optimal transport setup of Beiglböck, Henry‐Labordère, and Penkner, and show that our general results apply in both cases and enable us to obtain the pricing–hedging duality for American options.