Next Article in Journal
A Longitudinal Analysis of the Impact of Distance Driven on the Probability of Car Accidents
Next Article in Special Issue
Determining Economic Security of a Business Based on Valuation of Intangible Assets according to the International Valuation Standards (IVS)
Previous Article in Journal / Special Issue
Modelling Australian Dollar Volatility at Multiple Horizons with High-Frequency Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Note on Simulation Pricing of π-Options

by
Zbigniew Palmowski
*,† and
Tomasz Serafin
Faculty of Pure and Applied Mathematics, Wrocław University of Science and Technology, ul. Wyb. Wyspiańskiego 27, 50-370 Wrocław, Poland
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Risks 2020, 8(3), 90; https://doi.org/10.3390/risks8030090
Submission received: 8 July 2020 / Revised: 12 August 2020 / Accepted: 24 August 2020 / Published: 28 August 2020
(This article belongs to the Special Issue Quantitative Methods in Economics and Finance)

Abstract

:
In this work, we adapt a Monte Carlo algorithm introduced by Broadie and Glasserman in 1997 to price a π -option. This method is based on the simulated price tree that comes from discretization and replication of possible trajectories of the underlying asset’s price. As a result, this algorithm produces the lower and the upper bounds that converge to the true price with the increasing depth of the tree. Under specific parametrization, this π -option is related to relative maximum drawdown and can be used in the real market environment to protect a portfolio against volatile and unexpected price drops. We also provide some numerical analysis.
MSC:
49L20; 60G40
JEL Classification:
G13; C61

1. Introduction

In this paper, we analyze π -options introduced by Guo and Zervos (2010) that depends on so-called relative drawdown and can be used in hedging against volatile and unexpected price drops or by speculators betting on falling prices. These options are the contracts with a payoff function:
g ( S T ) = ( M T a S T b K ) +
in case of the call option and
g ( S T ) = ( K M T a S T b ) +
in the case of put option, where
S t = S 0 exp r σ 2 2 t + σ B t
is an asset price in the Black-Scholes model under martingale measure, i.e., r is a risk-free interest rate, σ is an asset’s volatility and B t is a Brownian motion. Moreover,
M t = sup w t S w
is a running maximum of the asset price and T is its maturity. Finally, a and b are some chosen parameters.
A few very well-known options are particular cases of a π -option. In particular, taking a = 0 and b = 1 produces an American option and by choosing a = 1 and b = 0 we derive a lookback option. Another interesting case, related to the concept of drawdown (see Figure 1), is when a = b = 1 and K = 1 . Then the pay-out function ( K M T a S T b ) + = 1 S T b M T a = M T S t M T = D T R equals the relative drawdown D t R , defined as a quotient of the difference between maximum price and the present value of the asset and the past maximum price. In other words, D t R corresponds to the percentage drop in price from its maximum. We take a closer look at this specific parametrization of the π -option in the later sections, starting from Section 3.2.
Monte Carlo simulations are widely used in pricing in financial markets they have proved to be valuable and flexible computational tools to calculate the value of various options as witnessed by the contributions of Barraquand and Martineau (1995); Boyle (1977); Boyle et al. (1997); Broadie et al. (1997); Caflisch (1998); Clément et al. (2002); Dyer and Jacob (1991); Geske and Shastri (1985); Glasserman (2004); Jäckel (2002); Joy et al. (1996); Longstaff and Schwartz (2001); Niederreiter (1992); Raymar and Zwecher (1997); Rogers (2002); Tilley (1993); Tsitsiklis and van Roy (1999, 2001); Resenburg and Torrie (1993); Villani (2010). One of the first attempts of Monte Carlo simulation for American options is by Tsitsiklis and van Roy (1999) where the backward induction algorithm was introduced. However, as appears later, Tilley method suffers from exponentially increasing computational cost as the number of dimensions (assets) increases. Broadie et al. (1997) to overcome this problem offered a non-recombining binomial simulation approach instead combined with some pruning technique to reduce computation burden and other variance reduction techniques to increase precision. In the same year Broadie and Glasserman (1997) construct computationally cheap lower and upper bounds to the American option price. This method is used in this paper. An alternative way to formulate the American option pricing problem is in terms of optimal stopping times. This is done in Carriere (1996), where it was proved that finding the price of American option can be based on a backwards induction and calculating several conditional expectations. This observation gives another breakthrough in pricing early exercise derivatives by Monte Carlo done by Longstaff and Schwartz (2001). They propose least square Monte Carlo (LSM) method which has proved to be versatile and easy to implement. The idea is to estimate the conditional expectation of the payoff from continuing to keep the option alive at each possible exercise point from a cross-sectional least squares regression using the information in the simulated paths. To do so we have to then solve some minimization problem. Therefore, this method is still computationally expensive. Some improvements of this method have been also proposed; see also Stentoft (2004a, 2004b) who gave theoretical foundation of LSM and properties of its estimator.
There are other, various pricing methods in the case of American-type options; we refer Zhao (2018) for review. We must note though that not all of them are good for simulation of prices of general π -options as it is a path-dependent product. In particular, in pricing π -options one cannot use finite difference method introduced by Brennan and Schwartz (1978); Schwartz (1977) which uses a linear combination of the values of a function at three points to approximate a linear combination of the values of derivatives of the same function at another point. Similarly, the analytic method of lines of Carr and Faguet (1994) is not available for pricing general π -options. One can use though a binomial tree algorithm (or trinomial model) though which goes backwards in time by first discounting the price along each path and computing the continuation value. Then this algorithm compares the former with the latter values and decide for each path whether or not to exercise; see Broadie and Detemple (1996) for details and references therein. It is a common belief that Monte Carlo method is more efficient than binomial tree algorithm in case of path-dependent financial instruments. It has another known advantages as handling time-varying variants, asymmetry, abnormal distribution and extreme conditions.
In this paper, we adapt a Monte Carlo algorithm proposed in 1997 by Broadie and Glasserman (1997) to price π -options. This numerical method replicates possible trajectories of the underlying asset’s price by a simulated price tree. Then, the values of two estimators, based on the price tree, are obtained. They create an upper and a lower bound for the true price of the option and, under some additional conditions, converge to that price. The first estimator compares the early exercise payoff of the contract to its expected continuation value (based on the successor nodes) and decides if it is optimal to hold or to exercise the option. This estimation technique is one of the most popular ones used for pricing American-type derivatives. However, as shown by Broadie and Glasserman (1997), it overestimates the true price of the option. The second estimator also compares the expected continuation value and early exercise payoff, but in a slightly different way, which results in underestimation of the true price. Both Broadie–Glasserman Algorithms (BGAs) are explained and described precisely in Section 2. The price tree that we need to generate is parameterized by the number of nodes and by the number of branches in each node. Naturally, the bigger the numbers of nodes and branches, the more accurate price estimates we get. The obvious drawback of taking a bigger price tree is that the computation time increases significantly with the size of the tree. However, in this paper we show that one can take a relatively small price tree and still the results are satisfactory.
The Monte Carlo simulation presented in this paper can be used in corporate finance and especially in portfolio management and personal finance planning. Having American-type options in the portfolio, the analyst might use the Monte Carlo simulation to determine its expected value even though the allocated assets and options have varying degrees of risk, various correlations and many parameters. In fact determining a return profile is a key ingredient of building efficient portfolio. As we show in this paper portfolio with π -options out-performs typical portfolio with American put options in hedging investment portfolio losses since it allows investors to lock in profits whenever stock prices reaches its new maximum.
In this paper, we use BGA to price the π -option on relative drawdown for the Microsoft Corporation’s (MSFT) stock and for the West Texas Intermediate (WTI) crude oil futures. Input parameters for the algorithm are based on real market data. Moreover, we provide an exemplary situation in which we explain the possible application of the π -option on relative drawdown to the protection against volatile price movements. We also compare this type of option to an American put and outline the difference between these two contracts.
This paper is organized as follows. In the next section we present the Broadie–Glasserman Algorithm. In Section 3 we use this algorithm to numerically study π -options for the Microsoft Corporation’s stock and WTI futures. Finally, in the last section, we state our conclusions and recommendations for further research in this new and interesting topic.

2. Monte Carlo Algorithm

Formulas identifying the general price of π -option are known in some special cases and they are given in terms of so-called scale functions and hence in terms of the solution of some second order ordinary differential equations; see for example, (Christensen 2013, chp. 5) and Egami and Oryu (2017) for details and further references. Still, the formulas are complex, and a Monte Carlo method of pricing presented in this paper is very efficient and accurate alternative method. In this section we present a detailed description of the used algorithm. In particular, we give formulas for two estimators, one biased low and one biased high, that under certain conditions converge to the theoretical price of the option.

2.1. Preliminary Notations

We adapt the Monte Carlo method introduced by Broadie and Glasserman (1997) for pricing American options. In this algorithm, values of two estimators are calculated on the so-called price tree that represents the underlying’s behavior over time. This tree is parametrized by the number of nodes n and the number of branches in each node—denoted by l. For example, the tree with parameters n = 2 , l = 3 is depicted in Figure 2.
To apply the numerical algorithm, we must discretize the price process given in (3), by considering the time sequence t 0 = 0 < t 1 < < t n = T with t i = i T n for i = 0 , , n . By
S t i l 1 , , l i
we denote the asset’s price at the time t i = i T n . The upper index l 1 , , l i , associated with t i , describes the branch selection (see Figure 3) in each of the tree nodes and allows us to uniquely determine the path of the underlying’s price process up to time t i . Similarly, we define
M t i l 1 , , l i = max k i S t k l 1 , , l k .
We introduce the state variable S ˜ t i l 1 , , l i = ( S t i l 1 , , l i , M t i l 1 , , l i ) as well.
We relate with it the payoff of an immediate exercise (for π -put) at time t i in the state S ˜ t i l 1 , , l i given by
h t i ( S ˜ t i l 1 , , l i ) = ( K S t i l 1 , , l i a M t i l 1 , , l i b ) +
and the expected value of holding the option from t i to t i + 1 , given asset’s value S ˜ t i l 1 , , l i at time t i defined via
g t i ( S ˜ t i l 1 , , l i ) = E e r n f t i + 1 ( S ˜ t i + 1 l 1 , , l i + 1 ) | S ˜ t i l 1 , , l i ,
where
f t i ( S ˜ t i l 1 , , l i ) = max { h t i ( S ˜ t i l 1 , , l i ) , g t i ( S ˜ t i l 1 , , l i ) }
is the option value at time t i in state S ˜ t i l 1 , , l i . Please note that
f t n ( S ˜ t n l 1 , , l n ) = f T ( S ˜ T l 1 , , l n ) = h T ( S ˜ T l 1 , , l n ) = ( K S T l 1 , , l n a M T l 1 , , l n b ) + .

2.2. Estimators

We will now give the formulas for the estimators Θ and Φ which overestimate and underestimate the true price of the option, respectively. Then, we will state the main theorem showing that both estimators are asymptotically unbiased and that they converge to the theoretical price of the π -option. We also provide a detailed explanation of the estimation procedure based on the exemplary price tree. In all calculations we consider a π -put option with parameters a = 1 , b = 1 and K = 1 . Additionally, we assume that the risk-free rate used for discounting the payoffs equals 5 % .

2.2.1. The Θ Estimator

The formula for the estimator is recursive and given by:
Θ t i = max h t i ( S ˜ t i l 1 , , l i ) , e r n 1 l j = 1 l Θ t i + 1 l 1 , , l i , j , i = 0 , , n 1 .
At the option’s maturity, T, the value of the estimator is given by
Θ T = f T ( S ˜ T ) .
The Θ estimator, at each node of the price tree, chooses the maximum of the payoff of the option’s early exercise at time t i , h t i ( S ˜ t i l 1 , , l i ) , and the expected continuation value, i.e., the discounted average payoff of successor nodes. Figure 4 shows how the value of Θ estimator is obtained given the certain realization of a price tree. All calculations are also shown below:
  • a _ Holding value : 0 + 10 115 + 0 3 e 0.05 0.028 ̲ Early exercise : 0
  • b _ Holding value : 5 110 + 0 + 13 110 3 e 0.05 0.052 Early exercise : 10 110 0.091 ̲
  • c _ Holding value : 10 110 + 15 110 + 5 110 3 e 0.05 0.086 Early exercise : 20 110 0.182 ̲
  • d _ Holding value : 0.028 + 0.091 + 0.182 3 e 0.05 0.095 ̲ Early exercise : 10 110 0.091

2.2.2. The Φ Estimator

The Φ estimator is also defined recursively. Before we give the formula we need to introduce an auxiliary function ξ by
ξ t i l 1 , , l i j = { h t i ( S ˜ t i l 1 , , l i ) , if h t i ( S ˜ t i l 1 , , l i ) e r n 1 l 1 k = 1 k j l Φ t i + 1 l 1 , , l i , k e r n Φ t i + 1 l 1 , , l i , j , if h t i ( S ˜ t i l 1 , , l i ) < e r n 1 l 1 k = 1 k j l Φ t i + 1 l 1 , , l i , k
for j = 1 , , l . Now we can define the Φ estimator in the following way:
{ Φ t i l 1 , , l i = 1 l j = 1 l ξ t i l 1 , , l i j Φ T = f T ( S ˜ T ) .
The formula for this estimator is more complicated. Therefore, we provide a detailed explanation of the mechanism behind the algorithm in the following part of this section. In our explanation we refer to Figure 5. Please note that in the following example, underlined numbers correspond to the final values associated with the specific branches of the tree.
  • a _ Early exercise : 0 Holding value for branch j = 1 : 0.087 + 0 2 e 0.05 0.041 > 0 0 ̲ Holding value for branch j = 2 : 0 + 0 2 e 0.05 = 0 0 = h i ( S ˜ t i ) 0 ̲ Holding value for branch j = 3 : 0 + 0.087 2 e 0.05 0.041 > 0 0 ̲
    For the branch j = 1 we look at the two remaining ones to determine whether early exercising (payoff = 0) or holding the option (payoff = 0.087 + 0 2 e 0.05 ) is more profitable. Obviously, early exercise is not optimal, so we hold the option and thus, as the value of ξ t 1 1 1 we take the payoff of the branch j = 1 which is 0.
    For the branch j = 2 both early exercise value and holding value from two other branches equals 0. Thus, from (4) the value of ξ t 1 1 2 equals the payoff of early exercise, which is 0.
    For the third branch, again holding the option is a more profitable decision (based on the payoffs of the two remaining branches). Thus, ξ t 1 1 3 takes the value corresponding to the branch j = 3 and it is 0.
    Now the value of the estimator for node a _ is the sum of ξ t 1 1 j across all branches:
    Φ t 1 1 = 1 3 j = 1 3 ξ t 1 1 j = 0 .
    Similarly, we have the following values of our estimator.
  • b _ Early exercise : 10 110 0.091 ̲ Holding value for branch j = 1 : 0 + 0.118 2 e 0.05 0.056 < 0.091 0.091 ̲ Holding value for branch j = 2 : 0.045 + 0.118 2 e 0.05 0.078 < 0.091 0.091 ̲ Holding value for branch j = 3 : 0.045 + 0 2 e 0.05 0.021 < 0.091 0.091 ̲
    In this case, the value of the estimator for the b _ node equals 0.091.
  • c _ Early exercise : 20 110 0.182 ̲ Holding value for branch j = 1 : 0.045 + 0.136 2 e 0.05 0.086 < 0.182 0.182 ̲ Holding value for branch j = 2 : 0.091 + 0.045 2 e 0.05 0.065 < 0.182 0.182 ̲ Holding value for branch j = 3 : 0.091 + 0.136 2 e 0.05 0.108 < 0.182 0.182 ̲
    For node c _ the value of the estimator is 0.182.
  • d _ Early exercise : 10 110 0.091 ̲ Holding value for branch j = 1 : 0.091 + 0.182 2 e 0.05 0.13 > 0.091 0 ̲ Holding value for branch j = 2 : 0.0 + 0.182 2 e 0.05 0.087 < 0.091 0.091 ̲ Holding value for branch j = 3 : 0 + 0.091 2 e 0.05 0.043 < 0.091 0.091 ̲
    The value of the estimator for this node equals 0 + 0.091 + 0.091 3 = 0.061 . This is also the (under)estimated value of the option.
Following arguments of Broadie and Glasserman (1997), one can easily prove the following crucial fact.
Theorem 1.
Both Θ and Φ are consistent and asymptotically unbiased estimators of the option value. They both converge to the true price of the option as the number of price tree branches, l, increases to infinity. For a finite l:
  • The bias of the Θ estimator is always positive, i.e.,
    E [ Θ 0 ( l ) ] f 0 ( S ˜ 0 ) .
  • The bias of the Φ estimator is always negative, i.e.,
    E [ Φ 0 ( l ) ] f 0 ( S ˜ 0 ) .
On every realization of the price tree, the low estimator Φ is always less than or equal to the high estimator Θ, i.e.,
P ( Φ t i l 1 , , l i Θ t i l 1 , , l i ) = 1 .

3. Numerical Analysis

In this section, we will present results of the numerical analysis. First, we use the algorithm described above to price the American option with arbitrary parameters. This will allow us to confirm that our Monte Carlo algorithm produces precise estimates of options’ prices. We focus on options related to Microsoft Corporation stock. Next, we price π -options for several combinations of parameters. We also consider π -option on drawdown using the real market data and we compare it with an American put, which is one of the most popular tool for protecting our portfolio against price drops.

3.1. American Options

First of all, we decided to check the robustness of the Monte Carlo pricing algorithm. We estimate prices of the American call options with different strike prices. In the example, the uderlying asset price S 0 equals 100, σ = 20 % , risk-free rate r = 5 % and the maturity is 30 days. In Table 1 we present the results of the estimation. Please note that when using the Broadie–Glasserman algorithm, we obtain the upper and the lower boundaries of the option price. To obtain the American option price estimate we average both values.

3.2. π -Options

We will analyze put π -option for various combinations of parameters a and b. We assume that parameter a is varying from 1.1 to 0.9 and b parameter between 0.9 and 1.1 . The ranges of these parameters have been chosen arbitrarily for illustrative purposes. All input parameters for options pricing, S 0 , M 0 , volatility and interest rate are taken from the real market data for the Microsoft Corporation stock (MSFT) and are given in Table 2. The numerical results are presented in Figure 6.

3.3. π -Options on Relative Drawdown

Recall that for a = 1 and b = 1 the payoff of the π -option equals
K S t M t + ,
where S t / M t is the current value of the relative drawdown of the underlying asset. We believe that such contracts could be very efficiently used for hedging and managing portfolio risk against the volatile drops in underlying’s price (see Section 3.4). One can adjust the payoff function (6) by the appropriate choice of the strike K. The choice is arbitrary and solely dependent on the risk management goals of the option’s buyer. It allows the setting of the minimal size of drawdown we would like to protect against and let the buyer adjust and full control of the level of our exposure at risk associated with unexpected price drops. For example by setting K = 9 10 , the payoff of our option becomes greater than zero only if the drop in the price of the underlying from its maximum exceeds 10 % . Of course, the bigger the value of K, the more expensive the option is.
We take a closer look at the impact of M t and K on the price of this special case of π -option. Here, we assume that the maximum price M t is between 100 and 120 and K ranges between 0.8 and 1. This time, the remaining parameters, namely S 0 , r and σ , have been arbitrarily chosen for illustrative purposes and are given in Table 3. The results are shown in Figure 7.

3.4. π -Options on Relative Drawdown - Application

We now focus on the potential application of π -options and compare the prices of American put and π -option on relative drawdown. We compare these particular instruments due to the fact that their values increase with the decrease of the underlying asset’s price. As an exemplary environment for the options comparison we choose two time series containing daily closing prices of the Microsoft Corporation’s stock (see Figure 8) as well as daily closing prices of the West Texas Intermediate (WTI) crude oil futures (see Figure 9). Both datasets are taken from www.finance.yahoo.com and span approximately one year, from 6 November 2017 to 9 November 2018. We use the first 9 months (from 6 November 2017 to 3 August 2018) to calibrate the historical volatility for both assets, which is one of the input parameters in our pricing algorithm.
Then, using the historical volatility, we compute prices of π and American options (using assets’ prices from 3 August 2018), both expiring 3 months after the end of calibration period. Please note that the parameters for the π -option on a relative drawdown are a = 1 , b = 1 and K = 1 . Input parameters for calculation and estimated options prices for both assets are given in Table 4 and Table 5.
Since the payoff of π -option on relative drawdown with K = 1 is always less than 1, to compensate against the drop in underlying’s price, we need a certain number of these contracts per each unit of stock in our portfolio. This number must be equal to M 0 . Please note that in Table 4 and Table 5, the real price of the single π -option on relative drawdown contract should be 0.0735 for MSFT and 0.0949 for WTI. However, in order to be able to compare the results to the American put values, we initially need to make the instruments pay the same amount in case of a price drop, therefore we multiply the price of single π -option on drawdown by M 0 (110 and 74 for MSFT and WTI respectively). That is why in Table 4 and Table 5 the price of π -option equals 0.0735 · 110 = 8.09 for the stock and 0.094 · 74 = 6.95 for the oil futures contract.
It turns out that π -option is more expensive than vanilla put in case of both assets, which is not a surprise as it initially pays the amount equivalent to the present maximum drawdown. However, since the difference in price between these instruments is rather significant, a question emerges whether there exists a situation in which purchasing π -option on relative drawdown is more profitable than buying a simple vanilla put. To answer this question, let us focus on the dashed part of the Microsoft Corporation and WTI futures data from the beginning of this section. In Figure 10 and Figure 11 we show the amount each instrument would pay (on each day) throughout the whole 3-month period until options’ maturity.
To display the difference more clearly, we construct two portfolios V American and V π , both consisting of an underlying asset (a single Microsoft Corporation stock or a barrel of the WTI crude oil) and an option (American put and π -option on relative drawdown, respectively). We observe them at the end of the volatility calibration period. Assets’ prices and options prices are taken from Table 4 and Table 5. In Table 6 and Table 7 we show the initial net values of both V American and V π portfolios.
Then we analyze the behavior of the constructed portfolios, by calculating the net value of each portfolio for each day until options’ maturity; see Figure 12 and Figure 13.
Based in Figure 12 and Figure 13 we can observe that the maximum value of portfolio V American is greater than the one for V π . Thus, when focusing purely at the possible maximum profit over some period of time, then the portfolio containing American option performs better. However, we can notice that V American ’s value over time is much more volatile compared to V π and it directly follows the behavior of underlying asset (it increases when asset’s price rises and decreases in the opposite situation). The value V π of π -option portfolio is most of the time non-decreasing. Moreover, V π increases its value every time the asset’s price reaches a new maximum and essentially does not decrease in case of any price drop. In other words, combining the underlying asset and π -option on drawdown allow us to lock in our profit whenever the price reaches its new maximum.
This brings us to the conclusion that the purpose of using π -option on relative drawdown and an American put is completely different. Vanilla American option protects us from asset price drops and ensures us that the current worth of our portfolio will not be less than its initial value. Unfortunately, in this case our portfolio’s value is more volatile and reflects the volatility of the underlying asset. This may result in bigger gains when compared to the use of π -option on relative drawdown if the price of the underlying rises significantly and stays on that level until option’s maturity. However, in case of a drop in asset price after the upswing, we do not benefit from the fact that the new maximum has been reached and thus the value of our portfolio decreases together with the price of the underlying asset. When looking at the value of V π over time one can notice that combining stock or a commodity and π -option on relative drawdown protects us against price drops as well but the volatility of our portfolio is reduced significantly. Additionally, the contract allows us to benefit from the underlying’s price upswings and locks in the profit every time new maximum is reached.
We have analyzed two datasets, MSFT and WTI, and the above analysis shows that the behavior of a portfolio based on π -option is similar for various choices of underlying assets.

4. Conclusions

In this paper we focus on the numerical pricing of the new derivative instrument—a π -option. We adapted the Monte Carlo algorithm proposed by Broadie and Glasserman (1997) to price this new option. We focused on a specific parametrization of this option which we call the π -option on drawdown. We observed that this specific financial instrument is related to so-called relative maximum drawdown. We obtained prices of the π -option on relative drawdown for the Microsoft Corporation stock with different parameters to examine the influence of those parameters on option’s premium. Our next step involved the analysis of two portfolios: first one based on a π -option on relative drawdown and the second one based on an American put. We used the Microsoft Corporation data as well as the West Texas Intermediate crude oil futures dataset. It turned out that the portfolios behave in a completely different manner. The value of the portfolio containing the American put was highly correlated with the underlying’s price movements and thus had an unpredictable and volatile behavior. On the other hand, combining π -option on relative drawdown with the underlying asset not only ensures that the worth of the portfolio will not drop below the initial level, but it also allows us to take advantage of price upswings and to reduce the portfolio’s volatility at the same time. Similar analysis could be carried out for a geometric Lévy process of asset price. One can also consider the regime-switching market.

Author Contributions

Methodology, Z.P. and T.S.; Formal Analysis, Z.P., T.S.; Investigation, Z.P. and T.S.; Writing—Original Draft Preparation, T.S.; Writing—Review and Editing, Z.P.; Visualization, T.S.; Supervision, Z.P.; Project Administration, Z.P.; Funding Acquisition, Z.P. All authors have read and agreed to the published version of the manuscript.

Funding

This paper is supported by the National Science Centre under the grant 2016/23/B/HS4/00566 (2017–2020).

Conflicts 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.

References

  1. Barraquand, Jérôme, and Didier Martineau. 1995. Numerical valuation of high dimensional multivariate American securities. Journal of Financial and Quantitative Analysis 30: 383–405. [Google Scholar] [CrossRef] [Green Version]
  2. Brennan, Michael J., and Eduardo S. Schwartz. 1978. Finite Difference Method and Jump Processes Arising in the Pricing of Contingent Claims. Journal of Financial and Quantitative Analysis 13: 461–74. [Google Scholar] [CrossRef]
  3. Boyle, Phelim P. 1977. Options: A Monte Carlo approach. Journal of Financial Economics 4: 323–38. [Google Scholar] [CrossRef]
  4. Boyle, Phelim, Mark Broadie, and Paul Glasserman. 1997. Monte Carlo Methods for Security Pricing. Journal of Economic Dynamics and Control 21: 1267–321. [Google Scholar] [CrossRef]
  5. Broadie, Mark, and Jerome Detemple. 1996. American Option Valuation, New Bounds, Approximations, and a Comparison of Existing Methods. Review of Financial Studies 9: 1211–50. [Google Scholar] [CrossRef] [Green Version]
  6. Broadie, Mark, and Paul Glasserman. 1997. Pricing American-style securities using simulation. Journal of Economic Dynamics and Control 21: 1323–52. [Google Scholar] [CrossRef]
  7. Broadie, Mark, Paul Glasserman, and Gautam Jain. 1997. Enhanced Monte Carlo estimates for American option prices. Journal of Derivatives 5: 25–44. [Google Scholar] [CrossRef] [Green Version]
  8. Caflisch, Russel E. 1998. Monte Carlo, and Quasi-Monte Carlo Methods. Acta Numerica 7: 1–49. [Google Scholar] [CrossRef] [Green Version]
  9. Carr, Peter, and Dmitri Faguet. 1994. Fast Accurate Valuation of American Options. Working Paper. Ithaca: Cornell University. [Google Scholar]
  10. Carriere, Jacques F. 1996. Valuation of the early-exercise price for options using simulations and nonparametric regression. Insurance: Mathematics Economics 19: 19–30. [Google Scholar] [CrossRef]
  11. Christensen, Sören. 2013. Contributions to the Theory of Optimal Stopping. Ph.D. dissertation, Christian-Albrechts Universität Kiel, Kiel, Germany. [Google Scholar]
  12. Clément, Emmanuelle, Damien Lamberton, and Philip Protter. 2002. An Analysis of a Least Squares Regression Method for American Option Pricing. Finance and Stochastics 6: 449–71. [Google Scholar] [CrossRef]
  13. Dyer, Lawrence, and David Jacob. 1991. An Overview of Fixed Income Option Models. In The Handbook of Fixed Income Securities. Edited by Frank Fabozzi. Homewood: Business One Irwin, pp. 742–73. [Google Scholar]
  14. Egamim, Masahiko, and Tadao Oryu. 2017. A direct solution method for pricing options involving the maximum process. Finance and Stochastics 21: 967–93. [Google Scholar] [CrossRef]
  15. Geske, Robert, and Kuldeep Shastri. 1985. Valuation by Approximation: A Comparison of Alternative Option Valuation Techniques. Journal of Financial and Quantitative Analysis 20: 45–71. [Google Scholar] [CrossRef]
  16. Glasserman, Paul. 2004. Monte Carlo Methods in Financial Engineering. New York: Springer. [Google Scholar]
  17. Guo, Xin, and Mihail Zervos. 2010. π options. Stochastic Processes and their Applications 120: 1033–59. [Google Scholar] [CrossRef] [Green Version]
  18. Jäckel, Peter. 2002. Monte Carlo Methods in Finance. Chichester: John Wiley. [Google Scholar]
  19. Joy, Corwin, Phelim P. Boyle, and Ken Seng Tan. 1996. Quasi-Monte Carlo methods in numerical finance. Management Science 42: 926–38. [Google Scholar] [CrossRef] [Green Version]
  20. Longstaff, Francis A., and Eduardo S. Schwartz. 2001. Valuing American options by simulation: A simple least-squares approach. Review of Financial Studies 14: 113–47. [Google Scholar] [CrossRef] [Green Version]
  21. Niederreiter, Harald. 1992. Random Number Generation and Quasi-Monte Carlo Methods. Vienna: Society for Industrial and Applied Mathematics. [Google Scholar]
  22. Raymar, Steven B., and Michael J. Zwecher. 1997. A Monte Carlo valuation of American call options on the maximum of several stocks. Journal of Derivatives 5: 7–23. [Google Scholar] [CrossRef]
  23. Rogers, Leonard C.G. 2002. Monte Carlo Valuation of American Options. Mathematical Finance 12: 271–86. [Google Scholar] [CrossRef]
  24. Schwartz, Eduardo S. 1977. The Valuation of Warrants: Implementing a New Approach. Journal of Financial Economics 4: 79–93. [Google Scholar] [CrossRef]
  25. Stentoft, Lars. 2004a. Convergence of the Least Squares Monte Carlo Approach to American Option Valuation. Management Science 50: 1193–203. [Google Scholar] [CrossRef] [Green Version]
  26. Stentoft, Lars. 2004b. Assessing the Least Squares Monte Carlo Approach to American Option Valuation. Review of Derivatives Research 7: 129–68. [Google Scholar] [CrossRef]
  27. Tilley, James A. 1993. Valuing American Options in a Path simulation model. Transactions of the Society of Actuaries 15: 499–550. [Google Scholar]
  28. Tsitsiklis, John N., and Benjamin Van Roy. 1999. Optimal stopping of Markov processes: Hilbert space theory, approximation algorithms, and an application to pricing high-dimensional financial derivatives. IEEE Transactions on Automatic Control 44: 1840–51. [Google Scholar] [CrossRef] [Green Version]
  29. Tsitsiklis, John N., and Benjamin Van Roy. 2001. Regression methods for pricing complex American-style options. IEEE Transactions on Neural Networks 12: 694–703. [Google Scholar] [CrossRef] [Green Version]
  30. van Resenburg, Janse, and Glenn Torrie. 1993. Estimation of multidimensional integrals: Is Monte Carlo the best method? Journal of Physics A: Mathematics and General 26: 943. [Google Scholar]
  31. Villani, Giovanni. 2010. A Monte Carlo approach to value exchange options using a single stochastic factor. In Mathematical and Statistical Methods for Actuarial Sciences and Finance. Edited by Ballester Ferrer and Corazza Pizzi. Milan: Springer, pp. 305–31. [Google Scholar]
  32. Zhao, Jinsha. 2018. American Option Valuation Methods. International Journal of Economics and Finance 10: 1. [Google Scholar] [CrossRef] [Green Version]
Figure 1. A sample drawdown for the Microsoft Corporation stock is marked with black arrows and dashed lines. Data is taken from www.finance.yahoo.com.
Figure 1. A sample drawdown for the Microsoft Corporation stock is marked with black arrows and dashed lines. Data is taken from www.finance.yahoo.com.
Risks 08 00090 g001
Figure 2. An example of the price tree. Underlying’s price is marked with circles and corresponding maximums are marked with rectangles.
Figure 2. An example of the price tree. Underlying’s price is marked with circles and corresponding maximums are marked with rectangles.
Risks 08 00090 g002
Figure 3. Branch selecting.
Figure 3. Branch selecting.
Risks 08 00090 g003
Figure 4. Explanation of Θ estimator.
Figure 4. Explanation of Θ estimator.
Risks 08 00090 g004
Figure 5. Explanation of Φ estimator.
Figure 5. Explanation of Φ estimator.
Risks 08 00090 g005
Figure 6. π -option price estimations for varying a and b—Microsoft Corporation stock.
Figure 6. π -option price estimations for varying a and b—Microsoft Corporation stock.
Risks 08 00090 g006
Figure 7. π -option on relative drawdown price estimations for varying K and M 0 parameters.
Figure 7. π -option on relative drawdown price estimations for varying K and M 0 parameters.
Risks 08 00090 g007
Figure 8. Daily closing prices of the Microsoft Corporation’s stock. Data spans from 6.11.2017 to 9.11.2018. Vertical dashed line indicates the end of volatility calibration period. Option prices are calculated based on the volatility and the stock’s price on 3.08.2018.
Figure 8. Daily closing prices of the Microsoft Corporation’s stock. Data spans from 6.11.2017 to 9.11.2018. Vertical dashed line indicates the end of volatility calibration period. Option prices are calculated based on the volatility and the stock’s price on 3.08.2018.
Risks 08 00090 g008
Figure 9. Daily closing prices of the West Texas Intermediate crude oil futures contracts. Data spans from 6.11.2017 to 9.11.2018. Vertical dashed line indicates the end of volatility calibration period. Option prices are calculated based on the volatility and the asset’s price on 3.08.2018.
Figure 9. Daily closing prices of the West Texas Intermediate crude oil futures contracts. Data spans from 6.11.2017 to 9.11.2018. Vertical dashed line indicates the end of volatility calibration period. Option prices are calculated based on the volatility and the asset’s price on 3.08.2018.
Risks 08 00090 g009
Figure 10. Microsoft Corporation stock closing prices (top) and the corresponding payoffs of π -option on relative drawdown and American put (bottom) with the parameters from Table 4.
Figure 10. Microsoft Corporation stock closing prices (top) and the corresponding payoffs of π -option on relative drawdown and American put (bottom) with the parameters from Table 4.
Risks 08 00090 g010
Figure 11. WTI crude oil futures contract closing prices (top) and the corresponding payoffs of π -option on relative drawdown and American put (bottom) with the parameters from Table 5.
Figure 11. WTI crude oil futures contract closing prices (top) and the corresponding payoffs of π -option on relative drawdown and American put (bottom) with the parameters from Table 5.
Risks 08 00090 g011
Figure 12. Microsoft Corporation stock closing prices (top) and payoffs of portfolios with parameters from Table 4 (bottom).
Figure 12. Microsoft Corporation stock closing prices (top) and payoffs of portfolios with parameters from Table 4 (bottom).
Risks 08 00090 g012
Figure 13. WTI crude oil futures contract closing prices (top) and payoffs of portfolios with parameters from Table 5 (bottom).
Figure 13. WTI crude oil futures contract closing prices (top) and payoffs of portfolios with parameters from Table 5 (bottom).
Risks 08 00090 g013
Table 1. Comparison of the estimated and ’real’ American option prices with different strikes. Absolute percentage errors are also included.
Table 1. Comparison of the estimated and ’real’ American option prices with different strikes. Absolute percentage errors are also included.
StrikeLow Est.High Est.Estimated PriceReal PriceAbs. Perc. Err.
$80$20.16$20.55$20.36$20.330.14%
$85$15.10$15.54$15.32$15.350.19%
$90$10.28$10.62$10.45$10.430.19%
$95$5.84$6.02$5.93$5.890.68%
$100$2.54$2.60$2.57$2.512.36%
$105$0.76$0.77$0.77$0.735.33%
$110$0.16$0.16$0.16$0.1413.35%
Table 2. Input parameters for pricing π -put option on the Microsoft Corporation stock.
Table 2. Input parameters for pricing π -put option on the Microsoft Corporation stock.
ParameterValue
S 0 106.08
M 0 110.83
σ 17.03%
r1.5%
l65
K1
Table 3. Input parameters for pricing π -option on relative drawdown.
Table 3. Input parameters for pricing π -option on relative drawdown.
ParameterValue
S 0 100
σ 20%
r5%
l100
Table 4. Input parameters for computation and estimated options’ prices for the MSFT dataset.
Table 4. Input parameters for computation and estimated options’ prices for the MSFT dataset.
MSFT
American Put π on Drawdown
ParameterValueParameterValue
K108.13 M 0 110.83
S 0 108.13
σ 24.06%
r2.25%
l100
T3
Option price
$5.47$8.09
Table 5. Input parameters for computation and estimated options’ prices for the WTI dataset.
Table 5. Input parameters for computation and estimated options’ prices for the WTI dataset.
WTI
American Put π on Drawdown
ParameterValueParameterValue
K68.49 M 0 74.15
S 0 68.49
σ 23.97%
r2.25%
l100
T3
Option price
$3.07$6.95
Table 6. Portfolios and their initial net values for the MSFT dataset.
Table 6. Portfolios and their initial net values for the MSFT dataset.
Portfolio V American V π
Initial asset value108.13108.13
Option premium$5.47$8.09
Option initial payoff$0$2.68
Portfolio’s net value$102.66$102.72
Table 7. Portfolios and their initial net values for the WTI dataset.
Table 7. Portfolios and their initial net values for the WTI dataset.
Portfolio V American V π
Initial asset value68.4968.49
Option premium$3.07$6.95
Option initial payoff$0$5.66
Portfolio’s net value$65.42$67.20

Share and Cite

MDPI and ACS Style

Palmowski, Z.; Serafin, T. A Note on Simulation Pricing of π-Options. Risks 2020, 8, 90. https://doi.org/10.3390/risks8030090

AMA Style

Palmowski Z, Serafin T. A Note on Simulation Pricing of π-Options. Risks. 2020; 8(3):90. https://doi.org/10.3390/risks8030090

Chicago/Turabian Style

Palmowski, Zbigniew, and Tomasz Serafin. 2020. "A Note on Simulation Pricing of π-Options" Risks 8, no. 3: 90. https://doi.org/10.3390/risks8030090

APA Style

Palmowski, Z., & Serafin, T. (2020). A Note on Simulation Pricing of π-Options. Risks, 8(3), 90. https://doi.org/10.3390/risks8030090

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