Abstract

Background

This article investigates the Least-Squares Monte Carlo Method by using different polynomial basis in American Asian Options pricing. The standard approach in the option pricing literature is to choose the basis arbitrarily. By comparing four different polynomial basis we show that the choice of basis interferes in the option's price.

Methods

We assess Least-Squares Method performance in pricing four different American Asian Options by using four polynomial basis: Power, Laguerre, Legendre and Hermite A. To every American Asian Option priced, three sets of parameters are used in order to evaluate it properly.

Results

We show that the choice of the basis interferes in the option's price by showing that one of them converges to the option's value faster than any other by using fewer simulated paths. In the case of an Amerasian call option, for example, we find that the preferable polynomial basis is Hermite A. For an Amerasian put option, the Power polynomial basis is recommended. Such empirical outcome is theoretically unpredictable, since in principle all basis can be indistinctly used when pricing the derivative.

Conclusion

In this article The Least-Squares Monte Carlo Method performance is assessed in pricing four different types of American Asian Options by using four different polynomial basis through three different sets of parameters. Our results suggest that one polynomial basis is best suited to perform the method when pricing an American Asian option. Theoretically all basis can be indistinctly used when pricing the derivative. However, our results does not confirm these. We find that when pricing an American Asian put option, Power A is better than the other basis we have studied here whereas when pricing an American Asian call, Hermite A is better.

Keywords

Background

Asian options are often used for cash flow hedges in companies whose purchase programming is set to mitigate the fluctuation of raw materials’ prices. Its versatility is confirmed by its presence in markets like commodities, electric power, interest rates and currency rates (McDonald 2006). Because they are complex (or exotic), the Asian options are usually traded over the counter. The characteristics of the contract (subject, premium, strike price, deadlines and maturity) are freely agreed between the parties, emphasizing their non-standardization.

Traditional techniques such as the finite-differences method and lattice become less attractive when dealing with pricing derivatives with multiple stochastic variables, problems with many dimensions, or even path-dependent American options, as it is the case of American Asian (Amerasian) options. The most flexible technique for pricing exotic options, such as American options, is the use of stochastic simulation with optimization algorithm. This technique includes different methods, such as the Least-Squares Monte Carlo method (LSM), first introduced by Longstaff and Schwartz (2001). Besides being faster and more precise to compute than other methodologies, the LSM methodology helps assess path-dependent American options with multiple dimensions and multiple state variables, being also applied to Markovian and non-Markovian problems.

The Least-Squares Monte Carlo Method has been used to price American Asian options. Longstaff and Schwartz (2001) exemplified the use of their technique in pricing an American Asian Arithmetic Average Fixed Strike call option with a specific polynomial basis, i.e., power. Moreno and Navas (2003) access the performance of Least Squares Monte Carlo Method numerically, by using two different polynomial basis, i.e., Laguerre and Hermite B, to value the same American Asian option priced by Longstaff and Schwartz (2001). Moreno and Navas (2003) found that the choice of the polynomial basis and the degree of these polynomials influence the estimated prices. Chaudhary (2005) used quasi-random sequences to improve the performance of this technique by pricing an American Asian Arithmetic Average call option without varying the polynomial basis used. In their experiment, they used power polynomial basis. Also to access the performance of this techniques, Cerrato and Cheung (2007) priced an American Asian Arithmetic Average call option using three sets of simulated trajectories. They applied two different polynomial basis, i.e., exponential and Laguerre. Although not dedicated to the task of pricing an American Asian option, other authors, such as Glasserman (2004) states that “accuracy depends on the choice of basis functions, which may require experimentation or good information about the structure of the problem”.

In this paper, we apply the Least-Squares Monte Carlo method intensively to price four different American Asian (Amerasian) options with four different polynomial basis. We do this for both put and call option. To access the performance of this algorithm, we use three sets of parameters.

The Least-Squares Monte Carlo Method

Every minute prior to the expiration date of an American option, the option holder must choose the optimal outcome between exercising the option or keeping it for another term. The choice is based on the immediate exercise payoff and the future payoff expected. If the former is greater than the latter, then the option is exercised.

Consistent with the no-arbitrage valuation theory, the continuation value is the risk-neutral expectation of the discounted future cash flow, i.e.,:

The idea of the Least-Squares Monte Carlo Method is to approximate the continuation value by using least-squares regression at every moment in which it is possible to exercise the option. The regression is done based on the data obtained for the state variables via Monte Carlo simulation and by choosing the trajectories where the option is in the money.

Within this context, whenever it is possible to exercise the option, the continuation value can be expressed as a linear combination of orthogonal basis functions, such as Power, Legendre and Laguerre polynomials. This follows from the Finance Literature that considers the payoff functions that belong to the function space of finite variance, represented by L2(Ω, ℑ, Q) 1. Since this is a Hilbert space2, any function F that belongs to this space can be written as a linear combination of orthogonal basis functions (Lima 2007); (Stentoft 2004). Thus, the same function F can be rewritten as:

where the polynomial basis, represented by Φj, is a function of one of the state variables. Note that the coefficients of the basis, αj, are not previously known, but they can be estimated by linear regression. The properties of the estimator and the convergence of the method are addressed algebraically by Clement et al. (2002) and Stentoft (2004).

When considering the number of basis to be used, it is important to point out a problem mentioned by Moreno and Navas (2003). They analyze the use of a large number of basis and conclude that several degrees of the chosen polynomial basis and its respective crossed products can increase the accuracy of the estimation. In some cases, however, an excessive increase in the number of basis can reduce the precision of the method, making it computationally expensive.

Methods

The Least-Squares Monte Carlo Method is applied in four different cases of Fixed-Start Time Window3 American Asian options: Case 1 –Arithmetic Average Floating Strike American Asian option; Case 2 –Geometric Average Floating Strike American Asian Option; Case 3 –Arithmetic Average Fixed Strike American Asian Option; and Case 4 – Geometric Average Fixed Strike American Asian Option.

Following Moreno and Navas (2003), we use the Bermuda-Asian option as a discrete approximation of the Amerasian option, whereby the price of the asset-object is monthly noted. All the priced options have a term of one year, and can be exercised from the third month of their issuance. It is assumed that the price of the asset-object follows a risk-neutral Geometric Brownian Motion.

Based on matrix notation, in which every column represents a specific time step and every line represents a simulated path, their respective payoff functions are formally expressed below, and the variables are set according to the following definitions:

According to Longstaff and Schwartz (2001), the structure used in regressions, expressed below by equation (1), comprises a constant, the two first degrees of the chosen polynomial basis and their crossed product up to the third degree, totaling eight basis. The polynomials used are Power, Legendre, Laguerre and Hermite A.

The matrix of the underlying asset prices, S, is defined recursively. The ith row vector corresponds to the ith simulated path and its elements are defined by:

Note that each column vector of the matrix corresponds to a specific time step.

Once the matrix of the underlying asset prices is set, the matrix of averages, \( \overline{\mathbf{S}} \), is determined by using relations introduced latter, in which the form depends solely on what kind of average is considered, i.e., whether it is an arithmetic average or it is a geometric one.

where, Φn(.) is the polynomial basis of nth degree and it is a function of one of the state variables. The dimension of matrices is set by j time steps and m simulated paths. So, at jth time step, the regression is performed by using the column vector related to the same time step. Therefore,

In this work, the Monte Carlo simulation was performed with pseudorandom number sequences and the variance reduction technique used is antithetic variables4. We used three sets of parameters varying them solely on strike prices used in experiments as follows:

Polynomial basis used

As mentioned above, the polynomials used in this study are Power, Legendre, Laguerre and Hermite A. All of them can be alternatively expressed by Rodrigues’ formula, explicit form or by the recurrence law (Abramowitz and Stegun 1972).

Following Abramowitz and Stegun (1972), Rodrigues’ formula is expressed by:

The format used to write a polynomial is chosen in such a way that the pricing procedure becomes operationally more practical. This method proved to be effective throughout the whole process, excluding the possibility that the rank of the coefficient matrix of the least square regression is insufficient. Thus, the explicit form of the Power and Laguerre polynomials is applied, the recurrence law is used to express the Hermite A polynomial, and Rodrigues’ formula is used to write the Legendre polynomial.

Results and discussion

Tables 4, 5, 6 and 7 show the average prices and their respective standard deviations of call and put options for the four types of Amerasian options studied. Tables 4, 5, 6 and 7 show that when different basis are used in all sets of simulated trajectories, and taking into account all sets of parameters used, there is a slight difference among the estimates of the option price. These results follow what Moreno and Navas (2003) stated i.e., the choice of the polynomial basis affects the estimates of the option prices. This difference becomes insignificant when a greater number of simulated trajectories is used.

Such pattern is the same for both put options and for call options of all the American Asian options considered. This result suggests two basic conclusions. First, it ratifies one of the main features of the method: if all of the polynomial basis are orthogonal, then theoretically all of them could be indistinctly used in derivative pricing. Second, it shows that the gradual reduction of differences among the price estimates for each of the basis due to the incremental increase in the number of simulated paths is a direct result of the law of large numbers.

The increase in the number of simulated trajectories markedly decreases the standard deviation of the price estimates for all of the basis used. Therefore, the estimates converge to their real price. This is true for both call and put American Asian options. To exemplify it graphically, Figs. 1 and 2 show the convergence of the estimated prices for the set of parameters in which K = 100 to all simulated trajectories in each studied cases considering the four polynomial basis.

The results in Table 8 suggest a certain degree of homogeneity regarding the most accurate basis considering not only the cases studied in a specific strike, but also when analyzing one case varying solely its strike price, whose basis presented in experiments made the lowest standard deviation at 50.000 simulated paths. This happens for most of cases and throughout strikes when the very same case is considered.

Table 8

The polynomial basis whose standard deviation is the lowest at 50.000 simulated paths

PUT

Strike

Case 1

Case 2

Case 3

Case 4

95

Power

Power

Legendre

Laguerre

100

Hermite A

Hermite A

Hermite A

Legendre

105

Hermite A

Hermite A

Legendre

Hermite A

CALL

Strike

Case 1

Case 2

Case 3

Case 4

95

Power

Hermite A

Power

Legendre

100

Power

Power

Power

Legendre

105

Legendre

Hermite A

Legendre

Hermite A

When analyzing one of the cases in study throughout strike prices, we observe that a specific basis delivers the lowest standard deviation in most of experiments. For example, in case 1, if it is a put option, then Hermite A delivers a price whose standard deviation is the lowest in two of the three set of parameters used. Whereas, for example, in case 3 for a call option across all strike prices considered, the polynomial basis that delivers a price with the lowest standard deviation is power. By repeating this process of column analysis of Table 8, in most of strikes considered, when pricing an American Asian put option, Hermite A provides a better performance to the algorithm. On the other hand, when pricing an American Asian call option, Power provides a better performance through different strikes.

Considering a specific strike price throughout cases in Table 8, we figure out that the same feature happens, i.e., there is, in most cases, one polynomial basis whose price delivered has the lowest standard deviation. For example, for a put option whose strike price is 100 in three out of the four cases studied, Hermite A delivers the lowest standard deviation for the price estimates. Whereas, for a call option with the same strike price, the polynomial basis power provides a better performance, i.e., the lowest standard deviation in most cases. By proceeding this process of row analysis in Table 8, we got in most cases when pricing an American Asian put option that Hermite A had the best performance. Whereas when pricing an American Asian call option, Power delivers a better performance in most of cases.

These results suggest that pricing American Asian options using the Least-Squares Monte Carlo Method enables the selection of one polynomial basis, regardless of the kind of Amerasian. In fact, the main concern is whether we are pricing put or call option.

Conclusions

Our experiments using the four chosen polynomial basis ratified one of the method’s elements, i.e., that any orthogonal polynomial basis can be used for pricing American Asian options, since in all of the simulated paths the estimated prices of the options are virtually the same.

For all of the Amerasian options priced, and for each of the basis used, the convergence of the price estimate results from the decrease of the standard deviation while the number of simulated paths increases.

Considering 50,000 simulated trajectories, our study suggests that it is possible to empirically choose one specific polynomial basis for pricing an American Asian option. The results show that one polynomial basis is marginally more accurate than others, by providing a lower standard deviation it delivers a better performance to the algorithm when pricing a complex option as American Asian options. Such empirical outcome is theoretically unpredictable, since in principle all basis can be indistinctly used when pricing the derivative.

Therefore, for practical purposes, if a trader is looking for a faster and more accurate way to valuing an American Asian call option, the polynomial basis suggested to perform this procedure is Hermite A. On the other hand, if their concern is valuing an Amerasian put option, the Power polynomial basis is recommended. This speediness means being capable of measuring the worthiness of any derivative in order to screen accurately business opportunities not only for hedge purposes but also for arbitrage purposes.

Declarations

Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

This article is part of Ms. USMdL dissetation of Master’s degree got in Pontifical Catholic University of Rio de Janeiro, whose advisor was Professor CPS Ph.D. Both authors read and approved the final manuscript.

Authors’ information

Ursula Silveira Monteiro de Lima holds a master’s degree in Industrial Engeneering at Pontifical Catholic University of Rio de Janeiro - DEI/PUC-Rio - research field: Finance and Investment Analysis. Besides, Ms. Lima holds a B.A. in Economics at Getulio Vargas Foundation – FGV/EPGE and a Graduate Specialization in Finance and the Capital Markets at Getulio Vargas Foundation of São Paulo – FGV-SP.

Carlos Patricio Samanez is graduated in Industrial Engineering, M.Sc. in Production Engineering and Ph.D. in Financial Economics. Associate professor at Faculty of Economics –State University of Rio de Janeiro - UERJ, and Department of Industrial Engineering - Pontifical Catholic University of Rio de Janeiro - PUC-Rio.

Authors’ Affiliations

(1)

Department of Industrial Engineering, Pontifical Catholic University of Rio de Janeiro, Gávea, Brazil