Abstract
This article develops, a lattice-based approach for pricing contingent claims when parameters governing the logs of the underlying asset dynamics are modelled by generalized hyperbolic distribution and normal inverse Gaussian distribution. The pentanomial lattice is constructed using a moment matching procedure. Moment generating functions of generalized hyperbolic distribution and normal inverse Gaussian distribution are utilized to compute probabilities and jump parameters under historical measure . Minimal entropy martingale measure (MEMM) is used to value European call option with a view of comparing the results with some of the existing benchmark model such as Black Scholes model. Empirical data from S&P500 index, RUTSELL2000 index and RUI1000 index are used to demonstrate how the model works. There is a significant difference especially for long term maturity (six months and above) type of contracts, the proposed model outperform the benchmark model, while performing poorly at short term contracts. Pentanomial NIG models seems to outperform the other models especially for long dated maturities.
Public Interest Statement
This study develops a simplified approach to value financial contracts such as calls and puts. Features of daily relative changes of price are known to be non-normal in contrast to the common assumption. The proposed model accommodates such several aspects of reality observed in financial markets in recent times. Subjected to observed market data, the proposed valuation model seems to outperform existing benchmark models. The simplified approach may be used to value different types of contracts which may not have closed-form solutions.
1. Introduction
In the past decade option pricing has become one of the major areas of financial theory and practice. Since the introduction of the celebrated Black Scholes option pricing which assumes that the underlying stock prices follow a geometric Brownian motion; there is an explosive growth in trading derivatives in the world wide financial market. Due to its compact and computational simplicity; the Black and Scholes (Citation1973) model enjoys great popularity in the financial markets. Recently the formulae has been extended in various ways. Recent studies have shown that the normal distribution may not accurately describe observed properties of stock returns; see for example Barndorff-Nielsen (Citation1998), Cont (Citation2001), Carr, German, Madan, and Yor (Citation2002) for a well documented stylized facts about returns. The deviations from normality become more severe when high frequency financial data are used.
A lattice is a graphical representation of all possible paths that might be followed any stochastic process say for example stock price. We construct a pentanomial lattice that approximates the evolution stock price. Lattices are useful for valuing a wide variety of options such as path dependent contracts which may not have a closed form solution such as lookback options, American type options and barrier option. Since options cash flow are functions of the future value of the underlying asset, options can be valued in the lattice by taking the expectation of their payoff. The current option value equals the discounted expected option payoff.
Lattices for option pricing were first introduced in 1979 in the pioneering work of Cox, Ross, and Rubinstein (Citation1979). In particular, they used binomial lattice to model geometric Brownian motion and Rendleman and Bartter (Citation1979) used binomial lattice to model exponential Poisson process. An attractive property of their model is that the binomial lattice for geometric Brownian motion is consistent with the standard (Black & Scholes, Citation1973) formula for European options. Due to simplicity and versatility of lattice models, a number of extensions to the basic model have been proposed, see Derman and Kani (Citation1994), Ritchken and Trevor (Citation1999), Yamada and Primbs (Citation2001), Wu (Citation2006) for example. Florescu and Viens (Citation2008) use quadrinomial tree to model stochastic volatility in option pricing, while Primbs, Rathianam, and Yamada (Citation2007) price options with a pentanomial lattice. It is worthy noting that an efficient lattice method, may be significantly faster than a Monte Carlo method for valuing some types of path dependent options.
The objective of this paper is to develop an option pricing lattice model which combine skewness and the leptokurtic nature of daily log returns under an alternative distributional assumption, that is consistent with empirical stock returns. Minimal entropy martingale measure (MEMM) is used to change probability measure to a risk neutral economy within a pentanomial lattice framework. Parameters of the model are selected to match the first four central moments of the returns. Such a model, has the potential of estimating option prices that are more consistent with empirically observed stylized facts of returns.
This paper proceeds as follows. In Section 2, we establish the general dynamics of the asset price over a time interval . In Section 3, a brief review of binomial, and pentanomial lattice is outlined. In Section 4, option pricing formulaes are derived in pentanomial framework and minimal entropy martingale measure is applied to change measure to risk neutral world Section 4 introduces numerical procedures in relation to derived formulaes using real market data. European call option is priced and numerical results compared. Section draws conclusions.
2. Basic model setup
Consider the stochastic distribution of the price of non-divided paying stock in a risk-neutral economy. Let the stock price be S(t) at time t in a period . An option pricing model is generally based on assumed process of the stock price or return. The Black and Scholes (Citation1973), for example assume that the stock price (under risk neutral measure ), movement is governed by the following process
where r is the risk free rate and is the instantaneous volatility rate of the stock return distribution. This is equivalent to assuming daily log returns are normally distributed with mean and variance . The resulting price of a contingent claim is given by(2.1) (2.1)
where t refers to the current time. However such a process does not incorporate skewness and leptokurtic fat tails which are typically associated with empirical stock log returns.
We consider a discrete time economy for a period of [t, T] where the trading takes place at any of the trading nodes where . Suppose is a given probability space, where is the statistical or data generating probability measure. Here, the sample space represents the uncertainty in our financial model. Let be the time index set of our financial model such that all economic activities take place at each time point We equip our probability space with the information structure to get a filtered probability space denoted by . That is, for each represents the information set of all market information up to and including time j, where contains all -null sets in In general, asset price is sometimes assumed to follow the process(2.2) (2.2)
where , , and are sequences of independent Poisson random variables with parameter . Changes of daily log returns are known to be leptokurtic and assumed to be laced with Poisson mixture of normal distributions. For more detailed exposition on modeling the dynamics of the underlying risky asset see Hsieh (Citation1989), Nieuwland, Verschoor, and Wolff (Citation1994), Chan and Maheu (Citation2002), Duan, Ritchken, and Sun (Citation2006) etc. Without loss of generality, our focus is on developing lattice model for the underlying process, say(2.3) (2.3)
where is a well chosen constant. Let , and define , which are assumed to be random variable whose first four moments are known from the market data simply for the sake of developing a pentanomial lattice framework. For simplicity, we let be modelled by generalized hyperbolic distribution and normal inverse Gaussian which are as defined below.
2.1. The generalized hyperbolic distribution
In this subsection we provide definition of the generalized hyperbolic distribution along with normal inverse Gaussian as its special case.
Definition 2.1
The probability density function of the one-dimensional Generalized Hyperbolic distribution is given by the following:(2.4) (2.4)
where and is the modified Bessel function of third kind with index given by(2.5) (2.5) .
According to Barndorff-Nielsen, the parameters domain is given by
In all cases, is the location parameter and can take any real value, is a scale parameter; and determine the distribution shape and defines the subclasses of GH and is related to the tail flatness.
Characteristic function of the GH is given by(2.6) (2.6)
while mean and variance are given respectively by the following
We consider normal inverse Gaussian (hereafter NIG) which is a special case of Generalized hyperbolic distributions when .
Definition 2.2
The Normal inverse Gaussian (hereafer NIG) distribution is a flexible four parameter distribution that can describe a wide range of shapes. A random variable if(2.7) (2.7)
where is the modified Bessel function of third kind, with the index 1 given by
It is interesting to note that NIG distribution can take a variety of different shapes. Having a heavier tail than that of normal distribution is considered suitable for modeling data sets with many extremal observations. The moments of a random variable are(2.8) (2.8)
The characteristic function of NIG random variable say Y is given by(2.9) (2.9)
2.2. Parameterizations
Although the parametrization is mostly used in literature we have other parameterizations like which is invariant under the transformation of the scale and location parameters with and . McNeil, Frey, and Embrechts (Citation2005) used the following parameterizations where(2.10) (2.10)
The parametrization , is derived if we set
Similar parametrization is used in ghyp R package.
The main challenge is to construct branching probabilities in the lattice. Our approach would be using moment matching technique.
3. Multinomial lattices
One of the most important joint distributions is the multinomial distribution, which arises when a sequence of n independent and identical experiments are performed. Suppose that each experiment can result in any one of possible outcomes, with respective probabilities . If we let denote the number of the n experiments that result in outcome number i, then(3.1) (3.1)
In multinomial lattice model, we need to determine the up and down rates u and d, and the probabilities to fit the actual market data as closely as possible. This can be done by moment matching or directly from density function, see Kellezi and Webber (Citation2004) for different ways of constructing branching probabilities in the lattice. Note that u and d may be thought of as up and down factors at each step. Also it can be shown that the multinomial lattice still recombines even if u and d are time dependent when is satisfied for some constant where and are up and down factors in each time step, see Yamada and Primbs (Citation2001), Yamada and Primbs (Citation2004) for more details. Let the up and down rates, u and d, be given as(3.2) (3.2)
where is the number of branches, and and are real numbers.
We develop the basic theoretical set up to model the dynamics of the underlying with an objective to value options in discrete time. It is assumed that, trades occur only at discrete dates indexed by , and the stock price at date can take on values only in a discrete set specified exogenously by
where the variables , index time and state values respectively, while is the possible number of future states for from , i.e.(3.3) (3.3)
with probabilities , satisfying . In this case, the stock may achieve possible prices at time given by(3.4) (3.4)
Let , then its jth central moment,
For more information about multinomial approximating models see Kamrad and Ritchken (Citation1991), Kargin (Citation2005). We briefly illustrate moment matching methodology, by considering the binomial and pentanomial models for a two time steps in the following subsection.
3.1. Binomial lattices
The binomial option pricing model is an iterative solution that models the price evolution over the whole option validity period . Figure represents the price evolution of the underlying asset as the binomial lattices of all possible prices at equally spaced time steps from today () under the assumption that at each step, the price can move, either up or down at a fixed rate and with respective pseudo-probabilities and . A standard (Cox et al., Citation1979) binomial tree, consists of a set of nodes, representing possible future stock prices, with a constant logarithmic spacing between these nodes.
The necessary equations for the binomial lattice are (3.5) (3.5)
From these two equations, we obtain several possibilities of solutions e.g. or
where is the variance of Y. A European call option with exercise price K and date n will have payoff in state given by(3.6) (3.6)
3.2. Pentanomial lattice construction
We consider state space for risky stock price dynamics over two trading dates as shown in Figure . At each date , the stock price can take on values in an exogenously specified discrete set indexed by j. The price denotes the stock price in state j at date for and respectively.
To construct a pentanomial model of stock prices, we examine the behavior of the stock price in an interval . The discrete distribution of over the interval is approximated to be pentanomial,as illustrated in Figure . To model the stock price movement as a pentanomial lattice, the interval [t, T] is divided into n equal subintervals of length , where T is the maturity date of an option. For convenience, define . To develop an option pricing model dependent on skewness and kurtosis, the distribution of is assumed to be fully specified by parameters via central moments.
Proposition 3.1
Let , and where is the characteristic function of random variable Y, and is the central moment of Thus
It follows that skewness and kurtosis of is given by(3.7) (3.7)
The relation (3.8) is used to form system of linear equations.(3.8) (3.8)
where are as defined above.
To calibrate the pentanomial lattice, we need to solve the following five equations(3.9) (3.9)
which implies that,(3.10) (3.10)
Making the column of the probabilities the subject, we get(3.11) (3.11)
The third and fourth equations arise from matching the third and fourth central moments of the approximating distribution to the third and fourth central moments respectively of the empirical distribution. How these central moments are related to skewness and excess kurtosis is described in Equation (3.8). On solving these five equations, we get Equation (3.12),(3.12) (3.12)
where(3.13) (3.13)
Note that must be chosen in order to ensure positivity of probabilities and It so happens that if and then, there exists a range of values of (which includes ) which will ensure that all the probabilities are strictly positive (see Primbs et al., Citation2007; Yamada & Primbs, Citation2001; Yamada & Primbs, Citation2004). This translates to the following proposition.
Proposition 3.2
For the choice of , Equation (3.12) reduces to the following probabilities in (3.14) with guaranteed positivity, and the corresponding jump amplitudes for the pentanomial lattice in Equation (3.2) respectively.(3.14) (3.14) (3.15) (3.15)
The notion of change of measure from to in an incomplete market implies existence of an equivalent measure which is not unique, with absence of arbitrage. One such martingale measure is minimal entropy martingale measure.
3.3. Minimal entropy martingale measure
One of the most important economic insight underlying the preference free option pricing result, is the concept of perfect replication of contingent claims, by continuously adjusting a self-financing portfolio under the no-arbitrage principle. Cox et al. (Citation1979) provided further insight in the concept of perfect replication by introducing the notion of risk-neutral valuation and establishing its relationship with no-arbitrage principle in a transparent way under a discrete-time binomial setting.
Harrison and Kreps (Citation1979) and Harrison and Pliska (Citation1981) established a solid mathematical foundation for the relationship between no-arbitrage principal and the notion of risk-neutral valuation using the modern language of probability theory. They proposed the “Fundamental theorem for asset pricing” which states that the absence of arbitrage opportunities is equivalent to the existence of an equivalent martingale measure. If the securities market is complete, there is a unique martingale measure and hence the unique price of any contingent claim is given by its discounted payoff at expiry under the martingale measure. However, the assumption of market completeness is questionable in the real world securities market. Under an incomplete market, there is more than one equivalent martingale measure and hence a range of no-arbitrage prices for a contingent claim. One crucial issue is to identify an equivalent martingale measure which gives an economically consistent and justifiable price for the contingent claim.
Let be the cardinality of , (where r denotes single period interest rate) and be the price process of the risky asset. We assume that is known and the random variable takes five different positive values with the probability
The minimal entropy martingale measure (MEMM) for the pentanomial lattice is the solution to the objective function say (3.16) (3.16) (3.17) (3.17)
It can be shown quite easily that is given by(3.18) (3.18)
where is the unique real solution (that always exists under the assumption of no-arbitrage opportunities) of the following equation(3.19) (3.19)
This is part of lemma due to Frittelli (Citation2000) in which, he links existence and uniqueness to to no arbitrage assumption. See Ssebungenyi (Citation2008), Miyahara (Citation2001), Fujiwara and Miyahara (Citation2003), Esche and Schweizer (Citation2005), Choulli and Striker (Citation2006), Ssebugenyi, Mwaniki, and Konlack (Citation2013) for more application(s) of minimal entropy martingale measure.
Proposition 3.3
Let be the option value at the node where i refers to the time instant and j is the one of the nodes in period i. Let . The price of the underlying asset in pentanomial node is
where u and d parameters are given by
At maturity, we have and going backwards in time, entropy price of the contingent claim is given by
4. Empirical results
4.1. Data description
The data set consists of three daily adjusted closing price of three major indices, that is S&P500 January 2, 1990 up to April 16, 2016, RUT2000 index from 2 January 1990 up to 11 March 2016 and RUI1000 10 December 1992 up to 8 March 2016. Basic statistics of the resulting data set are computed as shown in Table . All the three data sets indicate that they are negatively skewed and highly leptokurtic. This implies that they are not normally distributed. Over the entire period, we have the daily closing (adjusted) values of the indexes which we use in estimating the volatility parameter.(4.1) (4.1)
Table provided maximum likelihood estimates of parameters of generalized hyperbolic distribution, normal inverse distribution for the three sets of log returns.
As discussed earlier, specifications of pentanomial lattices are developed using the numerical procedure outlined in the previous section, In Table are risk neutral probabilities for S&P500 index, similar computation can be done for other indices.
Once the parameters of discrete distributions are specified, pentanomial lattice building procedure is analogous to that of binomial lattices. Option values are obtained through a recursive procedure, and the corresponding graphical results presented in Figures , , .
4.2. European call option prices
A call option gives the owner the right, but not the obligation, to buy a particular security at a pre-specified price within a pre-specified time period. The value of such an option will be intimately related to the distribution of the price of the underlying instrument at the time of maturity. Specifically the more volatile the underlying price process, the more valuable the option. The standard approach for pricing options rely on risk neutral valuation methods. In this risk-neutralized probability measure, the price of a call option, that does not allow for early exercise and pays no dividends, will be qual to the discounted expected value of the payoffs at the maturity date. Our analysis is meant to illustrate a possibility of modeling volatility dependencies when calculating option prices.
To that end, we compare the performance of three lattice models for short time and long term maturity level at the money and out of the money European call options priced in Black and Scholes (Citation1973) world, i.e.
where we let t refer to the present time and the time to expiry date in days. An option is said to be at the money if the exercise price, K, equal the current value of the underlying security. Lattices are calibrated in data generating process and transformed to minimal entropy martingale measure. In all the models same parameters are used and results plotted against real market data and compared to Black Scholes model of 1973.
4.3. Empirical performance of the proposed model
The pricing performance of our model is tested relative to 200 European call options on the S&P500 index, RUT2000 index, and RUI1000 index at the close of the market on 11 April 2016. The data were taken from market watch website. On 11 April 2016,the closing price for the three indices (S&P500,RUT2000, RUI1000) were = $2021,$1073, and $1117 respectively. We assumed general annual risk free rate of with no divided yield. We took long term options with maturities, 431, 320 and 455 days respectively. The performance evaluated based on real option prices of data of the proposed model is measured with three indicators: (i) the dollar root mean squared absolute error (RMSE), (ii) the average relative pricing error (ARPE) and (iii) the average absolute error (APE) given below.(4.2) (4.2) (4.3) (4.3) (4.4) (4.4)
where N represents the total number of options and is the average option price. Table summarizes the overall pricing errors of the various models considered here. We notice option prices computed based on pentanomial NIG lattice and Pentanomial GH lattice outperform the ones calculated based on BSM73 model for long term maturity days considered here as presented.
5. Conclusion
In this paper we establish the asset dynamics under the physical probability measure in incomplete market, also apply minimal entropy martingale measure to change dynamics to risk neutral . We assumed log-Lévy model to calibrate dynamics of the underlying price process and MEMM to change the historical probability risk neutral probability measure.
The valuation of contingent claims whose value depend on multiple sources of uncertainty is an important problem in financial economics. Since numerical methods for valuing such claims can be computationally expensive, the need for an efficient algorithm is clear. We made simplifying assumptions in that direction, even though there is more to be refined.
Although the pentanomial lattice provided in this article are tractable as the standard binomial, the pentanomial lattices approach, and may be extended to the multinomial case. Pentanomial lattices can be considered useful for relatively long term contracts (200 days and above) which can be used to solve American type options problems incorporating skewness and kurtosis. In depth study and more data sets are required to fine tune this observation.
Since option prices may react sensitively to changes in volatility, a proper specification of the conditional means at each step may play a crucial role in the proposed pentanomial model. Under the proposed framework, the market is in general incomplete, which is challenging to handle for the implication is a multitude of equivalent martingale measures and thus, a variety of no-arbitrage prices.
We note that under the proposed underlying dynamics, the proposed pricing model outperform bench mark model such as Black scholes model for the long term contract, pentanomial NIG lattice outperforms the other two models. We leave model refinement and extensions for future research. It would be interesting to incorporate changing volatility in the pentanomial framework and compare the result.
Acknowledgements
Financial support from the International Science Programme (ISP) Sweden in collaboration with the Eastern African Universities Mathematics programme, is greatly acknowledged. This article was fine tuned while visiting University of Cape town, South Africa as part of my post doctoral study invited by the African Collaboration for Quantitative Finance and Risk Research (ACQuFRR) directed by Prof. David Taylor. Thanks to the anonymous reviewer for his helpful corrections and suggestions.
Additional information
Funding
Notes on contributors
Ivivi J. Mwaniki
Ivivi J. Mwaniki completed his PhD at the University of Nairobi Kenya in collaboration with Uppsala University Sweden under East African Universities Mathematics Programme/International Science Programme. Currently he works for the University of Nairobi School of Mathematics. This article is part of the ongoing research in relation to option pricing in general. His research interest include option pricing, optimization, time series analysis and mixtures of various distributions, insurance and risk management.
References
- Barndorff-Nielsen, O. (1998). Process of normal inverse Gaussian type. Finance Stochastics, 2, 41–68.
- Black, F., & Scholes, M. (1973). The pricing of options and corporate liabilities. Journal of Political Economy, 81, 659–683.
- Carr, P., German, H., Madan, D., & Yor, M. (2002). The fine structure of asset returns: An empirical investigation. Journal of Business, 75, 305–332.
- Chan, W., & Maheu, J. (2002). Conditional jump dynamics in stock market returns. Journal of Business and Economic Statistics, 20, 377–389.
- Choulli, T. & Striker, C. (2006). More on minimal entropy-hellinger martingale measure. Mathematical Finance, 16(1), 1–19.
- Cont, R. (2001). Empirical properties of asset returns: Stylized facts and statistical issues. Quantitative Finance, 1, 223–226.
- Cox, J., Ross, S., & Rubinstein, M. (1979). Option pricing: A simplified approach. Journal of Financial Economics, 7, 229–264.
- Derman, E., & Kani, I. (1994). Implied trinomial trinomial trees of the volatility smile. Goldman Sachs quantitative strategies research notes.
- Duan, J., Ritchken, P., & Sun, Z. (2006). Approximating garch-jump models, jump-diffusion processes and option pricing. Mathematical Finance, 16, 21–52.
- Esche, F., & Schweizer, M. (2005). Minimal entropy preserves the lévy property: how and why. Stochastics and their Applications, 115, 299–327.
- Florescu, I., & Viens, F. (2008). Stochastic volatility: Option pricing using a multinomial recombining tree. Applied Mathematical Finance, 15, 151–181.
- Frittelli, M. (2000). The minimal entropy martingale measure and valuation problem in incomplete markets. Mathematical Finance, 10, 39–52.
- Fujiwara, T., & Miyahara, Y. (2003). The minimal entropy martingale measures for geometric lévy processes. Finance and Stochastics, 7, 509–531.
- Harrison, M., & Kreps, M. (1979). Martingales and arbitrage in multiperiod securities markets. Economic Theory, 20, 381–408.
- Harrison, M., & Pliska, S. (1981). Martingales and stochastic integrals in the theory of continous trading. Stochastic Processes and their Applications, 11, 215–260.
- Hsieh, D. (1989). Testing for nonlinear dependance in daily foreign exchanges rates. Journal of Business, 62, 339–368.
- Kamrad, B., & Ritchken, P. (1991). Multinomial approximating models for options with k state variables. Management Science, 37, 1640–1652.
- Kargin, V. (2005). Lattice option pricing by multidimensional interpolation. Mathematical Finance, 15, 635–647.
- Kellezi, E., & Webber, N. (2004). Valuing Bermudan options when the asset returns are Lévy processes. Quantitative Finance, 4, 87–100.
- McNeil, A., Frey, R., & Embrechts, P. (2005). Embrechts. Quantitative risk management: Concepts techniques and tools. Priceton University Press.
- Miyahara, Y. (2001). Geometric lévy process & memm: Pricing model and related estimation problems. Asia-Pacific Financial Markets, 8, 45–60.
- Nieuwland, F., Verschoor, W., & Wolff, C. (1994). Stochastic trends and jumps in ems exchange rates. Journal of International Money and Finance, 13, 699–727.
- Primbs, J., Rathianam, M., & Yamada, Y. (2007). Option Pricing with a pentanomial lattice model that incorporates skewness and kurtosis. Applied Mathematical Finance, 14, 1–17.
- Rendleman, R., & Bartter, B. (1979). Two state option pricing pricing. Journal of Finance, 3, 1093–1110.
- Ritchken, P., & Trevor, R. (1999). Pricing options under generalized GARCH and stochastic volatility process. Journal of Finance, 54, 337–402.
- Ssebungenyi, C. (2008). Valuation of real options using the minimal entropy martingale measure. Applied Mathematical Sciences, 2, 2875–2800.
- Ssebugenyi, C. S., Mwaniki, I. J., & Konlack, V. S. (2013). On the minimal entropy martingale measure and multinomial lattices with cumulants. Applied Mathematical Finance, 20, 359–379.
- Wu, C. (2006). The GARCH option pricing model: A modification of lattice approach. Review of Qantitative Finance and Accounting, 26, 55–66.
- Yamada, Y. & Primbs, J. (2001). Properties of multinomial lattice random walks for optimal hedges. International Conference of Computational Science, 579–588.
- Yamada, Y., & Primbs, J. (2004). Properties of multinomial lattices with cumulants for option pricing and hedging. Asia-Pacific Fianacial Markets, 11, 335–365.