The pricing of Asian options under stochastic interest rates

The purpose of this paper is to analyse the effect of stochastic interest rates on the pricing of Asian options. It is shown that a stochastic, in contrast to a deterministic, development of the term structure of interest rates has a significant influence. The price of the underlying asset, e.g. a stock or oil, and the prices of bonds are assumed to follow correlated two-dimensional Itô processes. The averages considered in the Asian options are calculated on a discrete time grid, e.g. all closing prices on Wednesdays during the lifetime of the contract. The value of an Asian option will be obtained through the application of Monte Carlo simulation, and for this purpose the stochastic processes for the basic assets need not be severely restricted. However, to make comparison with published results originating from models with deterministic interest rates, we will stay within the setting of a Gaussian framework.     

Real R&;D options with time-to-learn and learning-by-doing

We model R&D efforts to enhance the value of a product or technology before final development. Such efforts may be directed towards improving quality, adding new features, or adopting technological innovations. They are implemented as optional, costly and interacting control actions expected to enhance value but with uncertain outcome. We examine the interesting issues of the optimal timing of R&D, the impact of lags in the realization of the R&D outcome, and the choice between accelerated versus staged (sequential) R&D. These issues are also especially interesting since the history of decisions affects future decisions and the distributions of asset prices and induces path-dependency. We show that the existence of optional R&D efforts enhances the investment option value significantly. The impact of a dividend-like payout rate or of project volatility on optimal R&D decisions may be different with R&D timing flexibility than without. The attractiveness of sequential strategies is enhanced in the presence of learning-by-doing and decreasing marginal reversibility of capital effects.     

The resource allocation problem for demand forecasting and performance measurement

We consider a retailer investing in two monitoring functions for unobservable demand and salesperson’s effort. We show that improving effort monitoring is more effective. Moreover, demand monitoring may be less preferable when it becomes relatively cheaper and balancing these two may be worse when they become more complementary.     

Forest of stochastic meshes: A new method for valuing high-dimensional swing options

Swing options generalize American-style options as they allow the holder multiple exercise rights and control over the exercise amounts. In this work, we replace the standard (binomial) trees in the forest of trees algorithm with stochastic meshes, yielding the forest of stochastic meshes; a simulation-based method for valuing high-dimensional swing options. This new method handles general price processes and payoffs, produces high- and low-biased consistent estimators and a true option price confidence interval.     

Revisit of stochastic mesh method for pricing American options

From an importance sampling viewpoint, Broadie and Glasserman [M. Broadie, P. Glasserman, A stochastic mesh method for pricing high-dimensional American options, Journal of Computational Finance 7 (4) (2004) 35–72] proposed a stochastic mesh method to price American options. In this paper, we revisit the method from a conditioning viewpoint, and derive some new weights.     

The effect of forecasting and information sharing in SCM for multi-generation products

Proper selection of information sharing policy and forecasting method has a significant impact on supply chain performance, especially in the high-tech industry where the product life cycle is short and multiple generations of products coexist. This paper evaluates the value of information sharing with various forecasting methods where two generations of high-tech products compete with each other in the same market. We consider two market environmental factors and two supply chain factors for the Monte Carlo Simulation and find out the most ideal combination of information sharing policy and forecasting method producing the maximum profits and service level.     

The Valuation of American Options with Stochastic Stopping Time Constraints

Abstract

This paper concerns the pricing of American options with stochastic stopping time constraints expressed in terms of the states of a Markov process. Following the ideas of Menaldi et al., we transform the constrained into an unconstrained optimal stopping problem. The transformation replaces the original payoff by the value of a generalized barrier option. We also provide a Monte Carlo method to numerically calculate the option value for multidimensional Markov processes. We adapt the Longstaff–Schwartz algorithm to solve the stochastic Cauchy–Dirichlet problem related to the valuation problem of the barrier option along a set of simulated trajectories of the underlying Markov process.     

带补偿的随机二次规划的近似Lagrange-Newton方法

In this paper, two-stage stochastic quadratic programming problems with equality constraints are considered. By Monte Carlo simulation-based approximations of the objective function and its first (second)derivative,an inexact Lagrange-Newton type method is proposed.It is showed that this method is globally convergent with probability one. In particular, the convergence is local superlinear under an integral approximation error bound condition.Moreover, this method can be easily extended to solve stochastic quadratic programming problems with inequality constraints.     

Controlling for supplier switching in the presence of real options and asymmetric information

Applying a real option approach, this paper examines how asymmetric information alters key variables of a firm’s supplier switching process, such as the timing of contracting (hurried versus delayed contracting), transfer payments, set-up, switching, and abandonment decisions. In a symmetric information setting, delayed contracting is unambiguously beneficial. Abandoning the once established relation with the entrant supplier is never an issue. In contrast, under asymmetric information hurried contracting with potentially abandoning the relation can be beneficial. Consistent with adverse selection models, we find that under delayed contracting, in equilibrium, the firm switches less frequently to the entrant supplier (switching inertia). Surprisingly, we also find that under hurried contracting the firm switches more frequently to the entrant supplier (switching acceleration) and may abandon the relation. Finally, we study how these key variables of the supplier switching process change when also the incumbent supplier has private information (two-sided asymmetric information case).     

Model‐based quantification of the volatility of options at transaction level with extended count regression models

In this paper, we elaborate how Poisson regression models of different complexity can be used in order to model absolute transaction price changes of an exchange‐traded security. When combined with an adequate autoregressive conditional duration model, our modelling approach can be used to construct a complete modelling framework for a security's absolute returns at transaction level, and thus for a model‐based quantification of intraday volatility and risk. We apply our approach to absolute price changes of an option on the XETRA DAX index based on quote‐by‐quote data from the EUREX exchange and find that within our Bayesian framework a Poisson generalized linear model (GLM) with a latent AR(1) process in the mean is the best model for our data according to the deviance information criterion (DIC). While, according to our modelling results, the price development of the underlying, the intrinsic value of the option at the time of the trade, the number of new quotations between two price changes, the time between two price changes and the Bid–Ask spread have significant effects on the size of the price changes, this is not the case for the remaining time to maturity of the option. Copyright © 2006 John Wiley & Sons, Ltd.     

A Simple Derivation of and Improvements to Jamshidian's and Rogers' Upper Bound Methods for Bermudan Options

The additive method for upper bounds for Bermudan options is rephrased in terms of buyer's and seller's prices. It is shown how to deduce Jamshidian's upper bound result in a simple fashion from the additive method, including the case of possibly zero final pay‐off. Both methods are improved by ruling out exercise at sub‐optimal points. It is also shown that it is possible to use sub‐Monte Carlo simulations to estimate the value of the hedging portfolio at intermediate points in the Jamshidian method without jeopardizing its status as upper bound.     

Modelling stochastic skew of FX options using SLV models with stochastic spot/vol correlation and correlated jumps

It is known that the implied volatility skew of Forex (FX) options demonstrates a stochastic behaviour which is called stochastic skew. In this paper, we create stochastic skew by assuming the spot/instantaneous variance (InV) correlation to be stochastic. Accordingly, we consider a class of Stochastic Local Volatility (SLV) models with stochastic correlation where all drivers – the spot, InV and their correlation – are modelled by processes. We assume all diffusion components to be fully correlated, as well as all jump components. A new fully implicit splitting finite-difference scheme is proposed for solving forward PIDE which is used when calibrating the model to market prices of the FX options with different strikes and maturities. The scheme is unconditionally stable, of second order of approximation in time and space, and achieves a linear complexity in each spatial direction. The results of simulation obtained by using this model demonstrate the capacity of the presented approach in modelling stochastic skew.     

A stochastic model for operating room planning with elective and emergency demand for surgery

This paper describes a stochastic model for Operating Room (OR) planning with two types of demand for surgery: elective surgery and emergency surgery. Elective cases can be planned ahead and have a patient-related cost depending on the surgery date. Emergency cases arrive randomly and have to be performed on the day of arrival. The planning problem consists in assigning elective cases to different periods over a planning horizon in order to minimize the sum of elective patient related costs and overtime costs of operating rooms. A new stochastic mathematical programming model is first proposed. We then propose a Monte Carlo optimization method combining Monte Carlo simulation and Mixed Integer Programming. The solution of this method is proved to converge to a real optimum as the computation budget increases. Numerical results show that important gains can be realized by using a stochastic OR planning model.     

Global optimization by random perturbation of the gradient method with a fixed parameter

The paper deals with the global minimization of a differentiable cost function mapping a ball of a finite dimensional Euclidean space into an interval of real numbers. It is established that a suitable random perturbation of the gradient method with a fixed parameter generates a bounded minimizing sequence and leads to a global minimum: the perturbation avoids convergence to local minima. The stated results suggest an algorithm for the numerical approximation of global minima: experiments are performed for the problem of fitting a sum of exponentials to discrete data and to a nonlinear system involving about 5000 variables. The effect of the random perturbation is examined by comparison with the purely deterministic gradient method.     

Conditional Monte Carlo method for dynamic systems with random properties

A method is developed for calculating moments and other properties of states X(t) of dynamic systems with random coefficients depending on semi-Markov processes ξ(t) and subjected to Gaussian white noise. Random vibration theory is used to find probability laws of conditional processes X(t)∣ξ(·). Unconditional properties of X(t) are estimated by averaging conditional statistics of this process corresponding to samples of ξ(t). The method is particularly efficient for linear systems since X(t)∣ξ(·) is Gaussian during periods of constant values of ξ(t), so that and its probability law is completely defined by the process mean and covariance functions that can be obtained simply from equations of linear random vibration. The method is applied to find statistics of an Ornstein-Uhlenbeck process X(t) whose decay parameter is a semi-Markov process ξ(t). Numerical results show that X(t) is not Gaussian and that the law of this process depends essentially on features of ξ(t). A version of the method is used to calculate the failure probability for an oscillator with degrading stiffness subjected to Gaussian white noise.     

Dynamic pricing model and algorithm for perishable products with fuzzy demand

This paper studies the dynamic pricing problem of selling fixed stock of perishable items over a finite horizon, where the decision maker does not have the necessary historic data to estimate the distribution of uncertain demand, but has imprecise information about the quantity demand. We model this uncertainty using fuzzy variables. The dynamic pricing problem based on credibility theory is formulated using three fuzzy programming models, viz.: the fuzzy expected revenue maximization model, α‐optimistic revenue maximization model, and credibility maximization model. Fuzzy simulations for functions with fuzzy parameters are given and embedded into a genetic algorithm to design a hybrid intelligent algorithm to solve these three models. Finally, a real‐world example is presented to highlight the effectiveness of the developed model and algorithm. Copyright © 2009 John Wiley & Sons, Ltd.     

一类重尾风险因子的模拟及其投资高风险值和置信区间的估计

由于金融市场中的日周期或短周期对数回报率的样本数据多数呈现胖尾分布，于是现有的正态或对数正态分布模型都在不同程度上失效，为了准确模拟这种胖尾分布和提高投资风险估计及金融管理，本文引进了一种可根据实际金融市场数据作出调正的蒙特卡洛模拟方法．这个方法可以有效地复制金融产品价格的日周期对数回报率数据的胖尾分布．结合非参数估计方法，利用该模拟方法还得到投资高风险值以及高风险置信区间的准确估计。     

Bayesian modeling and forecasting of Value‐at‐Risk via threshold realized volatility

This study proposes a threshold realized generalized autoregressive conditional heteroscedastic (GARCH) model that jointly models daily returns and realized volatility, thereby taking into account the bias and asymmetry of realized volatility. We incorporate this threshold realized GARCH model with skew Student‐t innovations as the observation equation, view this model as a sharp transition model, and treat the realized volatility as a proxy for volatility under this nonlinear structure. Through the Bayesian Markov chain Monte Carlo method, the model can jointly estimate the parameters in the return equation, the volatility equation, and the measurement equation. As an illustration, we conduct a simulation study and apply the proposed method to the US and Japan stock markets. Based on quantile forecasting and volatility estimation, we find that the threshold heteroskedastic framework with realized volatility successfully models the asymmetric dynamic structure. We also investigate the predictive ability of volatility by comparing the proposed model with the traditional GARCH model as well as some popular asymmetric GARCH and realized GARCH models. This threshold realized GARCH model with skew Student‐t innovations outperforms the competing risk models in out‐of‐sample volatility and Value‐at‐Risk forecasting.     

Quantile forecasting based on a bivariate hysteretic autoregressive model with GARCH errors and time ‐varying correlations

To understand and predict chronological dependence in the second‐order moments of asset returns, this paper considers a multivariate hysteretic autoregressive (HAR) model with generalized autoregressive conditional heteroskedasticity (GARCH) specification and time‐varying correlations, by providing a new method to describe a nonlinear dynamic structure of the target time series. The hysteresis variable governs the nonlinear dynamics of the proposed model in which the regime switch can be delayed if the hysteresis variable lies in a hysteresis zone. The proposed setup combines three useful model components for modeling economic and financial data: (1) the multivariate HAR model, (2) the multivariate hysteretic volatility models, and (3) a dynamic conditional correlation structure. This research further incorporates an adapted multivariate Student t innovation based on a scale mixture normal presentation in the HAR model to tolerate for dependence and different shaped innovation components. This study carries out bivariate volatilities, Value at Risk, and marginal expected shortfall based on a Bayesian sampling scheme through adaptive Markov chain Monte Carlo (MCMC) methods, thus allowing to statistically estimate all unknown model parameters and forecasts simultaneously. Lastly, the proposed methods herein employ both simulated and real examples that help to jointly measure for industry downside tail risk.     

Optimal Partial Proxy Method for Computing Gammas of Financial Products with Discontinuous and Angular Payoffs

We extend the limit optimal partial proxy method to compute second-order sensitivities of financial products with discontinuous or angular payoffs by Monte Carlo simulation. The methodology is optimal in terms of minimizing the variance of the likelihood ratio weight. Applications are presented for both equity and interest-rate products with discontinuous payoff structures. The first-order optimal partial proxy method is also implemented to calculate the Hessians of insurance products with angular payoffs. Numerical results are presented which demonstrate the speed and efficacy of the method.