The pricing of derivatives on assets with quadratic volatility

The basic model of financial economics is the Samuelson model of geometric Brownian motion because of the celebrated Black-Scholes formula for pricing the call option. The asset's volatility is a linear function of the asset value and the model guarantees positive asset prices. In this paper, it is shown that the pricing partial differential equation can be solved for level-dependent volatility which is a quadratic polynomial. If zero is attainable, both absorption and negative asset values are possible. Explicit formulae are derived for the call option: a generalization of the Black-Scholes formula for an asset whose volatiliy is affine, the formula for the Bachelier model with constant volatility, and new formulae in the case of quadratic volatility. The implied Black-Scholes volatilities of the Bachelier and the affine model are frowns, the quadratic specifications imply smiles.     

Various passport options and their valuation

The passport option is a call option on the balance of a trading account. The option holder retains the gain from trading, while the writer is liable for the loss. Multi-asset passport options and passport options with discrete constraints are studied. For the first ones the pricing equations are Hamilton-Jacobi-Bellman equations. For those with discrete constraints, a linear complementary problem must be solved in order to price the option. The gain by selling passport options to utility maximizing investors and to investors who guess the market a certain percentage of the time is also examined.     

Option pricing and hedging in incomplete market driven by Normal Tempered Stable process with stochastic volatility

This paper develops a distribution class, termed Normal Tempered Stable, by subordinating a drifted Brownian motion through a strictly increasing Tempered Stable process that generalizes the Variance Gamma and the Normal Inverse Gaussian and is used to model the logarithm asset returns. The newly added parameter is to create subclasses for all the distributions discovered in financial market. The empirical test suggests that time series of Technology stock returns in US market reject both the Variance Gamma distribution and the Normal Inverse Gaussian distribution and admit instead another subclass of the Normal Tempered Stable distribution. Furthermore, we introduce stochastic volatilities into the Normal Tempered Stable process and derive explicit formulae for option pricing and hedging by means of the characteristic function based methods. To answer the question of how well different models work in practice, we investigate four models adopting data on daily equity option prices and obtain several findings from the numerical results. To sum up, the Normal Tempered Stable process with stochastic volatility is able to adequately capture implied volatility dynamics and seen as a superior model relative to the jump-diffusion stochastic volatility model, based on the construction methodology that incorporates more sophisticated and flexible jump structure and the systematic and realistic treatment of volatility dynamics. The Normal Tempered Stable model turns out to have the competitive performance in an efficient manner given that it only requires three parameters.     

Option pricing when the regime-switching risk is priced

We study the pricing of an option when the price dynamic of the underlying risky asset is governed by a Markov-modulated geometric Brownian motion. We suppose that the drift and volatility of the underlying risky asset are modulated by an observable continuous-time, finite-state Markov chain. We develop a two- stage pricing model which can price both the diffusion risk and the regime-switching risk based on the Esscher transform and the minimization of the maximum entropy between an equivalent martingale measure and the real-world probability measure over different states. Numerical experiments are conducted and their results reveal that the impact of pricing regime-switching risk on the option prices is significant.     

Can properly discounted projects follow geometric Brownian motion?

The geometric Brownian motion is routinely used as a dynamic model of underlying project value in real option analysis, perhaps for reasons of analytic tractability. By characterizing a stochastic state variable of future cash flows, this paper considers how transformations between a state variable and cash flows are related to project volatility and drift, and specifies necessary and sufficient conditions for project volatility and drift to be time-varying, a topic that is important for real option analysis because project value and its fluctuation can only seldom be estimated from data. This study also shows how fixed costs can cause project volatility to be mean-reverting. We conclude that the conditions of geometric Brownian motion can only rarely be met, and therefore real option analysis should be based on models of cash flow factors rather than a direct model of project value.     

Option pricing under joint dynamics of interest rates, dividends, and stock prices

This paper proposes a unified framework for option pricing, which integrates the stochastic dynamics of interest rates, dividends, and stock prices under the transversality condition. Using the Vasicek model for the spot rate dynamics, I compare the framework with two existing option pricing models. The main implication is that the stochastic spot rate affects options not only directly but also via an endogenously determined dividend yield and return volatility; consequently, call prices can be decreasing with respect to interest rates.     

An optimal portfolio model with stochastic volatility and stochastic interest rate

We consider a portfolio optimization problem under stochastic volatility as well as stochastic interest rate on an infinite time horizon. It is assumed that risky asset prices follow geometric Brownian motion and both volatility and interest rate vary according to ergodic Markov diffusion processes and are correlated with risky asset price. We use an asymptotic method to obtain an optimal consumption and investment policy and find some characteristics of the policy depending upon the correlation between the underlying risky asset price and the stochastic interest rate.     

Option pricing with mean reversion and stochastic volatility

Many underlying assets of option contracts, such as currencies, commodities, energy, temperature and even some stocks, exhibit both mean reversion and stochastic volatility. This paper investigates the valuation of options when the underlying asset follows a mean-reverting lognormal process with stochastic volatility. A closed-form solution is derived for European options by means of Fourier transform. The proposed model allows the option pricing formula to capture both the term structure of futures prices and the market implied volatility smile within a unified framework. A bivariate trinomial lattice approach is introduced to value path-dependent options with the proposed model. Numerical examples using European options, American options and barrier options demonstrate the use of the model and the quality of the numerical scheme.     

Heuristic option pricing with neural networks and the neuro-computer synapse 3

Today’s option and warrant pricing is based on models developed by Black, Scholes and Merton in 1973 and Cox, Ross and Rubinstein in 1979. The price movement of the underlying asset is modeled by continuous-time or discrete-time stochastic processes. Unfortunately these models are based on severely unrealistic assumptions. Permanently an unsatisfactory and quite artificial adaption to the true market conditions is necessary (future volatility of the underlying price). Here, an alternative heuristic approach with a highly accurate neural network approximation is presented. Market prices of options and warrants and the values of the influence variables form the usually very large output/ input data set. Thousands of multi-layer perceptrons with various topologies and with different weight initializations are trained with a fast sequential quadratic programming (SQP) method. The best networks are combined to an expert council network to synthesize market prices accurately. All options and warrants can be compared to single out overpriced and underpriced ones for each trading day. For each option and warrant overpriced and underpriced trading days can be used to ascertain a better buy and sell timing. Furthermore the neural model gains deep insight into the market price sen-sitivities (option Greeks), e.g., ?, Г, Θ and Ω. As an illustrative example we inves-tigate BASF stock call warrants. Time series from the beginning of 1996 to mid 1997 of 74 BASF call warrant prices at the Frankfurter Wertpapierborse (Frankfurt Stock Exchange) form the data basis. Finally a possible speed up of the training with the neuro-computer SYNAPSE 3 is briefly discussed     

A note on convergence of option prices and their Greeks for Lévy models

We study the robustness of option prices to model variation after a change of measure where the measure depends on the model choice. We consider geometric Lévy models in which the infinite activity of the small jumps is approximated by a scaled Brownian motion. For the Esscher transform, the minimal entropy martingale measure, the minimal martingale measure and the mean variance martingale measure, we show that the option prices and their corresponding deltas converge as the scaling of the Brownian motion part tends to zero. We give some examples illustrating our results.     

Stochastic Approximation Algorithms for Parameter Estimation in Option Pricing with Regime Switching

Abstract

This work is concerned with option pricing. Stochastic approximation/optimization algorithms are proposed and analyzed. The underlying stock price evolves according to two geometric Brownian motions coupled by a continuous-time finite state Markov chain. Recursive stochastic approximation algorithms are developed to estimate the implied volatility. Convergence of the algorithm is proved. Rate of convergence is also ascertained. Then real market data are used to compare our algorithms with other schemes.     

Mixing Monte-Carlo and Partial Differential Equations for Pricing Options

There is a need for very fast option pricers when the financial objects are modeled by complex systems of stochastic differential equations. Here the authors investigate option pricers based on mixed Monte-Carlo partial differential solvers for stochastic volatility models such as Heston’s. It is found that orders of magnitude in speed are gained on full Monte-Carlo algorithms by solving all equations but one by a Monte-Carlo method, and pricing the underlying asset by a partial differential equation with random coefficients, derived by Itô calculus. This strategy is investigated for vanilla options, barrier options and American options with stochastic volatilities and jumps optionally.     

High-order computational methods for option valuation under multifactor models

Many of the different numerical techniques in the partial differential equations framework for solving option pricing problems have employed only standard second-order discretization schemes. A higher-order discretization has the advantage of producing low size matrix systems for computing sufficiently accurate option prices and this paper proposes new computational schemes yielding high-order convergence rates for the solution of multi-factor option problems. These new schemes employ Galerkin finite element discretizations with quadratic basis functions for the approximation of the spatial derivatives in the pricing equations for stochastic volatility and two-asset option problems and time integration of the resulting semi-discrete systems requires the computation of a single matrix exponential. The computations indicate that this combination of high-order finite elements and exponential time integration leads to efficient algorithms for multi-factor problems. Highly accurate European prices are obtained with relatively coarse meshes and high-order convergence rates are also observed for options with the American early exercise feature. Various numerical examples are provided for illustrating the accuracy of the option prices for Heston’s and Bates stochastic volatility models and for two-asset problems under Merton’s jump-diffusion model.     

Optimal reinsurance–investment problem in a constant elasticity of variance stock market for jump‐diffusion risk model

In this paper, we consider the jump‐diffusion risk model with proportional reinsurance and stock price process following the constant elasticity of variance model. Compared with the geometric Brownian motion model, the advantage of the constant elasticity of variance model is that the volatility has correlation with the risky asset price, and thus, it can explain the empirical bias exhibited by the Black and Scholes model, such as volatility smile. Here, we study the optimal investment–reinsurance problem of maximizing the expected exponential utility of terminal wealth. By using techniques of stochastic control theory, we are able to derive the explicit expressions for the optimal strategy and value function. Numerical examples are presented to show the impact of model parameters on the optimal strategies. Copyright © 2011 John Wiley & Sons, Ltd.     

Indifference pricing and hedging in a multiple-priors model with trading constraints

This paper considers utility indifference valuation of derivatives under model uncertainty and trading constraints, where the utility is formulated as an additive stochastic differential utility of both intertemporal consumption and terminal wealth, and the uncertain prospects are ranked according to a multiple-priors model of Chen and Epstein (2002). The price is determined by two optimal stochastic control problems (mixed with optimal stopping time in the case of American option) of forward-backward stochastic differential equations. By means of backward stochastic differential equation and partial differential equation methods, we show that both bid and ask prices are closely related to the Black-Scholes risk-neutral price with modified dividend rates. The two prices will actually coincide with each other if there is no trading constraint or the model uncertainty disappears. Finally, two applications to European option and American option are discussed.     

Heston随机波动率模型下的资产负债管理问题

应用随机最优控制方法研究Heston随机波动率模型下带有负债过程的动态投资组合问题,其中假设股票价格服从Heston随机波动率模型,负债过程由带漂移的布朗运动所驱动.金融市场由一种无风险资产和一种风险资产组成.应用随机动态规划原理和变量替换法得出了上述问题在幂效用和指数效用函数下最优投资策略的显示解,并给出数值算例分别分析了市场参数在幂效用和指数效用函数下对最优投资策略的影响.     

Pricing Options on Defaultable Stocks*

? ?. This work was inspired by the SAMSI workshops on Financial Mathematics, Statistics and Econometrics (Fall 2005, Spring 2006 North Carolina). The author wishes to thank the organizers for the travel grant to participate in this stimulating event. I also would like to thank Bo Yang for his research assistance and the two anonymous referees and an anonymous associate editor for their valuable suggestions. Stock option price approximations are developed for a model which takes both the risk of default and the stochastic volatility into account. The intensity of defaults is assumed to be influenced by the volatility. It is shown that it might be possible to infer the risk neutral default intensity from the stock option prices. The proposed option price approximation has a rich implied volatility surface structure and fits the data implied volatility well. A calibration exercise shows that an effective hazard rate from bonds issued by a company can be used to explain the impliedvolatility skew of the option prices issued by the same company. It is also observed that the implied yield spread obtained from calibrating all the model parameters to the option prices matches the observed yield spread.     

HARA frontiers of optimal portfolios in stochastic markets

In this paper, we consider the optimal portfolio selection problem in continuous-time settings where the investor maximizes the expected utility of the terminal wealth in a stochastic market. The utility function has the structure of the HARA family and the market states change according to a Markov process. The states of the market describe the prevailing economic, financial, social and other conditions that affect the deterministic and probabilistic parameters of the model. This includes the distributions of the random asset returns as well as the utility function. We analyzed Black–Scholes type continuous-time models where the market parameters are driven by Markov processes. The Markov process that affects the state of the market is independent of the underlying Brownian motion that drives the stock prices. The problem of maximizing the expected utility of the terminal wealth is investigated and solved by stochastic optimal control methods for exponential, logarithmic and power utility functions. We found explicit solutions for optimal policy and the associated value functions. We also constructed the optimal wealth process explicitly and discussed some of its properties. In particular, it is shown that the optimal policy provides linear frontiers.     

Testing robustness in calibration of stochastic volatility models

Many numerical aspects are involved in parameter estimation of stochastic volatility models. We investigate a model for stochastic volatility suggested by Hobson and Rogers [Complete models with stochastic volatility, Mathematical Finance 8 (1998) 27] and we focus on its calibration performance with respect to numerical methodology.In recent financial literature there are many papers dealing with stochastic volatility models and their capability in capturing European option prices; in Figà-Talamanca and Guerra [Towards a coherent volatility pricing model: An empirical comparison, Financial Modelling, Phisyca-Verlag, 2000] a comparison between some of the most significant models is done. The model proposed by Hobson and Rogers seems to describe quite well the dynamics of volatility.In Figà-Talamanca and Guerra [Fitting the smile by a complete model, submitted] a deep investigation of the Hobson and Rogers model was put forward, introducing different ways of parameters' estimation. In this paper we test the robustness of the numerical procedures involved in calibration: the quadrature formula to compute the integral in the definition of some state variables, called offsets, that represent the weight of the historical log-returns, the discretization schemes adopted to solve the stochastic differential equation for volatility and the number of simulations in the Monte Carlo procedure introduced to obtain the option price.The main results can be summarized as follows. The choice of a high order of convergence scheme is not fully justified because the option prices computed via calibration method are not sensitive to the use of a scheme with 2.0 order of convergence or greater. The refining of the approximation rule for the integral, on the contrary, allows to compute option prices that are often closer to market prices. In conclusion, a number of 10 000 simulations seems to be sufficient to compute the option price and a higher number can only slow down the numerical procedure.     

Professor V. Lakshmikantham,Editor Stochastic Analysis and Applications (1983–2012)

We use the stochastic calculus of variations for the fractional Brownian motion to derive formulas for the replicating portfolios for a class of contingent claims in a Bachelier and a Black–Scholes markets modulated by fractional Brownian motion. An example of such a model is the Black–Scholes process whose volatility solves a stochastic differential equation driven by a fractional Brownian motion that may depend on the underlying Brownian motion.