`Journal of Applied Mathematics and Stochastic AnalysisVolume 2008, Article ID 275217, 20 pageshttp://dx.doi.org/10.1155/2008/275217`
Research Article

## A Time-Series Approach to Non-Self-Financing Hedging in a Discrete-Time Incomplete Market

Department of Mathematical Sciences, Bentley College, 175 Forest Street, Waltham, MA 02452-4705, USA

Received 16 May 2008; Accepted 30 July 2008

Copyright © 2008 N. Josephy et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

We present an algorithm producing a dynamic non-self-financing hedging strategy in an incomplete market corresponding to investor-relevant risk criterion. The optimization is a two-stage process that first determines market calibrated model parameters that correspond to the market price of the option being hedged. In the second stage, an optimal set of model parameters is chosen from the market calibrated set. This choice is based on stock price simulations using a time-series model for stock price jump evolution. Results are presented for options traded on the New York Stock Exchange.

#### 1. Introduction

This paper is a continuation of the work originally presented in [1] where we developed an algorithm producing a dynamic non-self-financing hedging strategy in an incomplete market optimizing a suitable investor-relevant risk criterion. The algorithm expanded on theoretical investigations of A. V. Nagaev et al. (see [25]) where the authors studied asymptotic behavior of the residual value of a minimum cost super-hedge. The residual value occurs as a result of a non-self-financing dynamic hedging strategy of an option seller in an incomplete market introduced and discussed in [2].

The financial market model considered by Nagaev et al. consists of a bond, a stock, and a European type derivative security with a convex payoff function. The stock price jumps are assumed to be distributed over a bounded interval. This model is a natural extension of the classical binomial model of Cox-Ross-Rubinstein as well as its further multinomial extensions (see, e.g., [68], and references therein).

The model considered in this paper as well as its multinomial predecessors produces an incomplete market. A significant proportion of research on option pricing and hedging in incomplete markets constructs self-financing trading strategies that satisfy both a primary no-arbitrage condition and secondary conditions on portfolio risk and return. A comprehensive survey of modern methodologies can be found in [9]. A number of articles that deal with frictions in markets, shortfall risks, and quadratic hedging (all producing incomplete markets) can be found in the recent compendium [10].

Less prevalent is the study of non-self-financing trading strategies in similar economic environments. The encyclopedic reference [11] and the more modest [12] both illuminate option pricing with consumption, the model which is similar to the work presented here. The application of constrained quadratic programming to the pricing of options by constructing non-self-financing portfolios in an incomplete market can be found in [13]. We note that the alternative approach of specifying a deterministic model of incomplete markets is an active area of research, as illustrated by [1417].

In [1] as well as in the present work, we study the short-term behavior of the residual value of a dynamic non-self-financing hedging strategy, whose long-term behavior was studied by Nagaev in [5]. Our goal is to develop and evaluate an algorithm that will determine a dynamic non-self-financing hedging strategy consisting of a portfolio of our stock and bond assets in an incomplete market. The portfolio will approximately hedge the derivative security and will satisfy additional criteria that are meaningful to the investor, based on the deviation of the portfolio value from the required hedging value.

Our algorithm is a two-stage process that first determines a set of market calibrated model parameters that correspond to the market price of the option being hedged. In the second stage, an optimal set of model parameters is chosen from this set. This choice is based on stock price simulations. In our initial work [1], we used a bootstrap method to simulate future stock price paths assuming that the stock price jumps are independent identical distributed random variables. In the present work, we relax the i.i.d. assumption on the stock price jumps and use a time-series model for stock price jump evolution.

We fit an ARMA model to historical stock price jump series and extract the correlation structure in the stock price jump process. Further, we simulate ARMA model driving noise to construct future values of the price jump process by using the nonparametric bootstrap method of sampling with replacement from the ARMA model residual sequence. The above method (unlike common parametric methods based on an assumed Gaussian distribution) captures long tails observed in typical stock price jump distributions.

The remainder of the paper is organized as follows. We develop the discrete time financial model in Section 2. The notion of the residual value of a minimum cost super-hedge is described in Section 3. Non-self-financing dynamic hedging strategies and residual values associated with them are discussed in Section 4. Section 5 contains a brief overview of our hedging algorithm. The first stage of the algorithm, where we find the set of market calibrated model parameters, is presented in Section 6. Section 7 describes the second stage of the algorithm: the numerical optimization over the set of market calibrated model parameters in order to determine an optimal hedging strategy based on a risk criterion chosen by an investor. Section 8 is devoted to the ARMA model based stock price simulation. The results of applying our algorithm to a variety of stock options are presented in Section 9. Section 10 presents the comparison of our optimal hedging strategy with the Black-Scholes hedging strategy. We conclude with some remarks in Section 11.

#### 2. Discrete-Time Market Model and No-Arbitrage Option Prices

Our discrete-time financial market model consists of two fundamental assets and a derivative security.(1) A risk-free bond with fixed periodic interest rate , evolving from an initial value at time to at time as(2) A risky stock evolving from an initial value at time to at time aswhere the stock price jumps are assumed to be random variables distributed over a bounded interval , at every time step No further assumptions are made on the distribution for the .(3)A European type option with expiry and convex payoff function . We will require the conditionto be fulfilled to guarantee no-arbitrage in this market.

This discrete-time market model is incomplete. Indeed it is well known that the famous binomial model is complete while its trinomial (as well as more general multinomial) extensions are not (see, e.g., [6]). Our model generalizes the multinomial models allowing the stock price jumps at every time step to be distributed over a bounded interval , The incompleteness of the model is manifested in an open interval () of no-arbitrage option prices (see, e.g., [2, 18]). It is shown in [18] that the end points of the interval are given by the following formulas:where and define the support of the stock price jump distribution (see above), is a convex payoff function of the option, is the number of periods to the option expiration, is the stock price at time , and is defined as follows:(Here, is the binomial coefficient.)

For the option seller, the upper bound is the demarcation between risk sharing with the option buyer (if the option sale price is below ) and the potential for arbitrage profit (if the option sale price is at or above ).

#### 3. Minimum Cost Super Hedge

For the rest of the paper, we take the position of an option seller who wishes to hedge the potential liability of the sold option being exercised. In this section, we will consider an extreme case: suppose the option with the convex payoff function is sold at time for the upper bound price The option seller uses the amount to finance the setup cost of a portfolio consisting of stocks and bonds with the intention of hedging the short position in the option: The seller rebalances the portfolio at each time instant () creating a dynamic trading strategy ().

We will choose the quantities and according to the formulaswhere is defined in (2.6). S. A. Nagaev and A. V. Nagaev [2] showed (based on a convexity argument) that the dynamic trading strategy () represents a minimum cost super-hedging strategy in the following sense: the associated portfolio value at every time instant is greater than or equal to the value of the option. It is worth explaining the above statement in more detail.

Suppose at each time instant , the option seller liquidates the portfolio constructed in the prior period and uses the proceeds to construct a new portfolio for the current period . Let us denote by the liquidation value of the prior period portfolio:Combining (3.4) with (3.2) and (3.3), one getsOn the other hand, the funds required to finance the current period portfolio, or setup cost, are given by the upper bound of the no-arbitrage price interval corresponding to the time moment : (see (2.4)). The difference between the liquidation value (3.5) and the setup cost (3.6) is a residual amount : Since we assume that , it follows from the convexity of the payoff function that the residual is nonnegative:

In this fashion, at every time step , the option seller obtains a nonnegative residual , which is withdrawn after each portfolio liquidation prior to the construction of the next time period super-hedge. The above constructed super-hedging strategy is in general non-self-financing.

The accumulated value of the withdrawn residuals at maturity, which we will refer to as the minimum cost super-hedge residual, is given by

#### 4. Non-Self-Financing Hedging Strategies and their Residuals

In the previous section, we considered a hypothetical situation where the option has been sold at time for the price equal to the upper bound of the no-arbitrage option price interval. We saw how the option seller could use this option premium to construct a minimum cost super hedge based on a non-self-financing strategy.

In this section, we consider a more realistic situation where the initial (time zero) option price is lower than the upper bound , but still falls within the open interval of no-arbitrage option prices In this case, the option seller cannot in general construct a super hedge, but it is possible to construct a non-self-financing trading strategy that will produce a possibly nonnegative residual amount.

In order to explain how such trading strategy can be constructed, we need to present the following short mathematical introduction. Let us consider the quantity given as follows:where is the stock price at time , and is defined by (2.6) with boundary parameters , replaced with the values , such thatThe next proposition follows from the convexity arguments.Proposition 4.1. Let be a convex function. The function maps the set of pairs satisfying (4.2) onto the option price interval defined in (2.4) and (2.5). When and when then

The above proposition infers that for any choice of and satisfyingthe quantity (4.1) falls within the no-arbitrage option price intervalConversely, every no-arbitrage market option price (every point ) can be put into correspondence with at least one pair satisfying (4.3):

Suppose the option with the convex payoff function was sold at time for the price If we assume that there is no arbitrage on the market, then there is at least one pair satisfying (4.3) that allows the amount to be identified with the no-arbitrage option price satisfying (4.4) with In other words,We remark here that the exact values of the boundary parameters and are not important for practical purposes of option hedging, since they play purely theoretical role in our setting.

Let us return to (4.5). There are an infinite number of pairs satisfying (4.5). For a fixed observed market option price we will call the set of pairs satisfying (4.5) the market calibrated set of model parameters. Our goal will be to choose the best market calibrated pair based on an optimization criterion explained in Section 7.

The option seller uses the amount to setup the hedging portfolio :which will be rebalanced at every time instant following the dynamic portfolio strategywhere and are defined in (3.2) and (3.3), respectively, with the boundary parameters , replaced with the values , :

We remark that in fact there is an infinite number of dynamic portfolio strategies (4.7) defined by the formulas (4.8). These strategies are distinguished by the values of the market calibrated parameters We will call them the market calibrated dynamic portfolio strategies, or simply market calibrated hedging strategies.

By analogy with the arguments of Section 3, for every one can calculate the difference between the liquidation value of the portfolio constructed in the period and the setup cost of the portfolio constructed for the period This difference constitutes a residual valueThe latter formula (4.9) was obtained by analogy with the formula (3.7), which represents the residual amount corresponding to the boundary case of the minimum cost super hedge.

It is straightforward to show that(i) if (ii) if or (iii) if or In order to maintain the dynamic portfolio strategy defined by (4.7), at each time step the investor will either withdraw the residual (4.9) from the liquidated proceeds when or add the amount when . The local residuals produce an accumulated residual at option expiration:Note that the accumulated residual defined in (4.10) differs from the minimum cost super-hedge residual defined in (3.9) as illustrated by the characterizations of given above. S. A. Nagaev and A. V. Nagaev studied asymptotic properties of the minimum cost super-hedge residual (3.9) (or riskless profit of the investor, in his terminology) extensively in [24]. That work was extended to the asymptotic properties of the accumulated residual (4.10) in [5]. Our work investigates the short-term behavior of the accumulated residual and its usefulness in constructing practical hedging strategies for this market model.Remark 4.2. We would like to stress here that the dynamic portfolio strategy constructed in (4.7) is in general non-self-financing. Also, although the strategy provides an approximate hedging of the short position in an option, for the sake of simplicity we will still call it a hedging strategy.

Given the market option price theoretically one has an infinite choice of market calibrated model parameters ( pairs) each determining a market calibrated hedging strategy. Each market calibrated hedging strategy produces a residual sequence. An investor (an option seller) will want to choose values for and that determine a residual sequence with additional desirable risk/return characteristics. It is the choice of the model parameter values and based on the risk/return characteristics of the residual sequence that constitutes our algorithm design.

For the remainder of this paper, we will assume a European call option payoff function : where is the option strike price. Also, for all our computations we will use the value for the risk-free interest rate. The latter is justified by the results of our computational experiments confirming that a realistic variation in values of has a minimal effect on the results for the case of short-term option contracts studied in this paper.

#### 5. Brief Overview of the Hedging Algorithm

We will start with the brief summary of the previous sections. As was explained in Section 2, the incompleteness of our market model results in an infinite number of possible no-arbitrage option prices located within the open interval Assuming that there is no arbitrage on the market, we can associate a market option price with one of the points within the theoretical interval of no-arbitrage option prices: . The infinite set of pairs for which the equality holds is called a market calibrated set of model parameters.

Each no-arbitrage option price gives rise to a non-self-financing dynamic hedging strategy , , that is distinguished by the value of the parameters If the no-arbitrage option price under consideration is identified with the actual market option price so that we will call the associated hedging strategies market calibrated.

Each market calibrated hedging strategy produces a residual sequence and an accumulated residual amount when applied to a given stock price path. By imposing a risk criterion on the residual sequence, we can choose a unique optimal pair of parameters by optimizing the risk criterion over the set of market calibrated parameters. This results in the choice of the corresponding optimal market calibrated hedging strategy

This optimization is accomplished computationally by selecting a representative finite subset from the infinite set of market calibrated model parameters ( pairs) and simulating a large number of stock price paths for each selected pair. Each stock price path determines a residual sequence, from which we can estimate a criterion value using a suitable sample statistic.

Figure 1 gives a brief scheme of our proposed two-stage algorithm.

Figure 1: Two-stage algorithm.

We will now proceed with a more detailed explanation of each stage of the algorithm.

#### 6. Market Calibrated Set of the Model Parameters: Stage 1 of the Algorithm

Let us place ourselves in the setting described in Section 4. At time the option seller receives a premium which he/she is willing to use in order to build a dynamic hedging strategy that depends on the unknown parameters and Our goal at the first stage of the algorithm is to choose a set of market calibrated parameters and In order to do so we will identify the market option premium with a no-arbitrage option price : (see (4.1)). The market calibrated parameters must satisfy (6.1), where is given by It is more convenient to deal with the normalized quantity which we will denote by whereWe will be interested in the level curves of the surface (we will call them contours):The set of market calibrated model parameters ( pairs) is represented by the contour withWe recall that is the market option price at time and is the time zero stock price.

Several contours for a representative option are shown in Figure 2. The option BAC+JV is a Bank of America call option expiring on October 15, 2004 with a strike price of \$42.50 and a current stock price of \$44.16. The contours shown each correspond to a value of (indicated on the curve scaled by a factor of 100), and the contour shown in bold is the set of market calibrated model parameters described by (6.5) and (6.6) with a value of . The value of (being less than 1) indicates that the call option is currently in the money. The shape of the market calibrated contour is similar for all data sets examined.

Figure 2: Contours for BAC+JV option.

For the sake of comparison, Figure 3 depicts the set of contours for an Intel call option (INTC+JE) that is out of the money (). The option expired October 15, 2004 with a strike price of \$25. The market calibrated contour (indicated in bold) corresponds to the value The shape of the contour differs slightly from the market calibrated contour shown in Figure 2. The shape of this contour is similar for other out of the money options.

Figure 3: Contours for INTC+JE option.

The first stage of our algorithm is accomplished computationally by utilizing contour construction software to compute a finite number of market calibrated pairs satisfying (6.5) and (6.6). There are typically between 90 and 100 pairs identified on the market calibrated contour. It is this set of pairs that is used by the second stage of our algorithm.

#### 7. Choosing the Optimal Hedging Strategy: Stage 2 of the Algorithm

Let us recall that each pair on the market calibrated contour defined in (6.5) and (6.6) determines a market calibrated dynamic hedging strategy , . For a given stock price path , this hedging strategy produces a sequence of residuals , each representing a residual profit/loss for an investor. This sequence is the economic consequence of choosing model parameters and the associated dynamic hedging portfolio.

The second stage of the two-stage algorithm selects a unique pair on the market calibrated contour This chosen pair of parameters uniquely defines a hedging strategy , that will numerically optimize a risk criterion chosen by the investor.

We consider three possible criteria to choose from. These criteria convert the residual sequence into a scalar measure of investor risk, each reflecting some aspect of the option seller attitude towards risk.(i) Maximize the likelihood of a positive accumulated residual:(ii) Minimize expected shortfall:(iii) Maximize the expected accumulated profit:

The first criterion (7.1) interprets a positive residual as a profit and chooses a pair that has the highest probability of a net profit. In the absence of arbitrage, a large accumulated profit is not attainable with high probability. There is the possibility, however, of an investor achieving a small positive profit. The optimization problem presented here produces a market calibrated hedging strategy that maximizes the likelihood of a positive accumulated profit.

The second criterion (7.2) reflects an investor's desire to minimize the amount of single period additional funding needed to rebalance the portfolio over the life of the option. A negative residual represents the cash shortfall in rebalancing the portfolio at time The largest negative is the largest shortfall value. Optimizing this criterion produces a hedging portfolio with minimal expected single period additional funding.

Our final criterion (7.3) maximizes the expected accumulated residual, which reflects total net profit from using the dynamic portfolio strategy based on the chosen . It was shown in [2] that the expected accumulated profit is asymptotically constant on contours of constant no-arbitrage price. We thus anticipate minimal differences in the expected accumulated profit at each pair on our market calibrated contours when is large. For small empirical results show that it is possible to have a market contour with nonconstant expected accumulated profit.

Let us choose one of the listed criteria. The process of selecting a unique optimal pair of parameters consists of the following. For each pair in the market calibrated subset, we simulate a number of stock price time series , . We then apply the hedging strategy associated with the given pair to each of the simulated paths in order to determine a sequence of residuals . A particular choice of a risk criterion reduces the sequence to a single scalar value of risk (described in (7.1) through (7.3)). For example, the appropriate scalar value of risk for the criterion (7.3) is the accumulated residual For a given pair, we collect the corresponding sample of scalar risk values (the sample size equals the number of the simulated stock price paths). An appropriate sample statistic (mean value or probability of a desirable event) is then computed from the sample as the utility value of the pair. The best value of the sample statistic, as the pair is varied over the market calibrated contour, is chosen as the optimal . For example, the optimal pair for criterion (7.3) will correspond to the largest averaged accumulated residual

A description of the modeling process for the stock price time series is described in the following section.

#### 8. Stock Price Process Modeling

In our initial study (see [1]), we examined hedging portfolios produced by a risk minimization algorithm that simulated a stock price process under the assumption of independence of stock price jumps. A bootstrap procedure was used to sample stock price jumps with replacement from historical data to generate stock price paths. This approach assumed the independence of the sample stock price jumps. In this paper, we consider a more sophisticated model of the stock price process. We use time-series model to extract the correlation structure in the stock price jump process, producing a more structured model of the stock price process that accounts for autocorrelation in the stock price jump data.

Time-series models are developed for stock price jumps . We fit an ARMA model of the formto historical price jump series. Here, and are the parameters of the model, while are the error terms. Examining the data for the stock price jumps, it is apparent that the jump series is stationary having approximately constant mean and variance. It is also important to note that we are dealing with short-term behavior, so long-term trends and seasonality are not an issue.

The appropriate model order is determined individually for each data set based on the standard analysis of several factors (see [19]). A preliminary estimate of is made by examining the autocorrelation (ACF) and partial autocorrelation (PACF) functions. The model fit for several values of is evaluated using diagnostic statistics. The evaluation criteria include the statistical significance level of individual ARMA terms, Akaike information criterion (AIC) value associated with the model, and validation of the noise assumptions. The lack of correlation in the residual noise sequence is established using the Box-Ljung statistics. Although many of the data sets considered indicate a model, each data set exhibits unique patterns, and higher order models are appropriate in some cases. The model parameters are evaluated based on the model order. The time-series model for the stock price jump process and an appropriate residual noise process can be used to simulate future values of the stock price jumps (and then stock prices).

Common parametric methods for generating residual noise values include random number generation using an assumed underlying distribution. Examining the residuals in our time-series jump data, we see that the distribution has a long tail indicative of a nonnormal distribution (see Figure 4). To guarantee the appropriate inclusion of such tail valued price jumps in our simulation, the method chosen for simulating ARMA model driving noise for future values of the price jump process is the nonparametric bootstrap method of sampling with replacement from the residual sequence .

Figure 4: ARMA residual density function.

Using the ARMA residual sequence as noise, we simulate future values of the stock jump process. The jump process is accumulated to form a stock price process, which can be used to calculate a residual sequence or an accumulated residual using our risk minimization algorithm. This process is summarized in the flow chart shown in Figure 5.

Figure 5: Simulation of stock price paths.

We conclude this section with a brief overview of the existing literature on autoregressive processes with given marginal distributions. In our model, we assume that the stock price jump distribution at every time instant has a bounded support. A natural question arises: is it possible to construct an ARMA (or AR) type process with marginal distributions of bounded support? This issue is addressed in a number of recent works (see, e.g., [2023], and references therein). Namely, in [20] the author uses a special random coefficient autoregression to model a first-order AR process with beta marginal distributions. In [21], the authors extend on the earlier works [22, 23] and present a first-order autoregressive time-series model with the uniform marginal distributions. Autoregressive time-series processes with marginal distributions of other types (not necessarily with bounded supports) have been studies extensively in the literature. Namely, [24] discusses non-Gaussian ARMA processes with marginal distributions of the Laplace and l-Laplace types. The first-order autoregressive models producing time series with logistic, hyperbolic secant, exponential, Laplace, and Gamma marginal distributions are considered in [25]. In [26], by means of special choice of noise, autoregressive processes with Student type marginal distributions are constructed. An extensive survey of non-Gaussian conditional linear autoregressive models can be found in [27].

#### 9. Numerical Results

##### 9.1. Market Environment and Option Characteristics

In order to test the algorithm under varying market conditions and with a range of option characteristics, data were collected for 55 call options traded on the New York Stock Exchange. The volatility index, or VIX, was used as an indicator of the market environment. The volatility index for a period from September 2002 through April 2007 is shown in Figure 6. Two time periods were selected with differing VIX characteristics. The selected time periods are indicated by the darker sections of the graph between the dashed vertical lines. The first period encompasses dates ranging from October 2002 through October 2004. As shown in the figure, the volatility index was high and widely varying at the beginning of the time period, with lower values towards the end of the period. The second chosen time period was July 2006 through March 2007, with smaller and less variable volatility index values.

Figure 6: Volatility index.

Option data was collected for options in the money, at the money and out of the money, from varying industrial sectors and time to expiration. Option characteristics are detailed in Table 4. The data collected for each option was strike price and expiration date, historical stock price data, and historical option prices. Over 500 daily stock prices were recorded. To capture the market behavior corresponding to the volatility index fluctuations occurring at the beginning of 2003, all 500 data values were used in constructing time-series models for the stock price jump process for the options expiring on October 15, 2004. We used 150 historical stock prices in constructing time-series models for the stock price jump process for the options expiring on March 17, 2007 since the volatility index was lower and more stable in this regime. Tests using larger amounts of historical data produced similar results, indicating that a large amount of historical data was not required.

For all options in this study, 100 stock price paths were simulated. The market option price data for a period of 40 days prior to expiration were collected, and the market option price days prior to expiration was used in identifying the appropriate contour for each of the reported values of (as described in Section 6).

##### 9.2. Discussion of Algorithm Results

The risk criteria values given in (7.1) through (7.3) were evaluated for each pair on the market calibrated contour for each option listed in Table 4. Table 1 presents the optimal risk criteria values for all options with days to expiration. In addition, the minimum value for each risk criterion is presented in order to gauge the magnitude of the investor's gain associated with choosing the optimal hedging strategy produced by the algorithm. Note that the number at the end of each option ticker indicates the number of days to expiration.

Table 1: Algorithmic results for options with 30 days to expiration.

Considering the first option in Table 1 (BAC+CJ30), we see that the optimal value of or the expected accumulated residual is 0.3446 as reported in column 3 (Resid-max), which is over twice as large as the smallest possible value of 0.1701 reported in column 2 (Resid-min). Similarly, we can examine the improvement in shortfall associated with the algorithmic results. The minimum expected shortfall is (Shortfall-min), approximately 10% of the worst case scenario with a shortfall of (Shortfall-max). The probability of a positive accumulated residual increases from a possible low of 0.71 (PosProb-min) presented in column 6 to 1 (PosProb-max) using the optimal portfolio identified by the algorithm.

In most cases, the expected accumulated residual value is improved by a factor between 1.5 and 3 when the optimal hedging portfolio is chosen. Several cases produced much more dramatic improvements in the objective value. The most significant improvement in expected accumulated residual is found in the results for the XOM+CN30 option. The value associated with the optimal hedging portfolio is 0.5323, almost 100 times as large as the smallest possible value of .0058. Other notable cases include the Microsoft options MSQ+CK30 and MSQ+CY30 where the optimal accumulated residual is approximately 20 times as large as the smallest possible value and the Cisco option CYQ+CY30 where the optimal is more than 13 times as large as the smallest possible value.

On average, the shortfall is reduced to approximately 13% of the largest possible shortfall by following the hedging strategy associated with the optimal hedging portfolio. While the shortfall is reduced to just over 1% of the worst case for the GE option GE+CF30, the optimal is approximately 46% of the worst case for the Cisco option CYQ+CY30. It is important to note that the hedging portfolio associated with minimizing shortfall is not the same portfolio choice for maximizing expected accumulated residual.

The Intel option INTC+JE30 provides an interesting example for the benefit of choosing the hedging strategy produced by the algorithm to maximize the probability of a positive accumulated residual. The objective value associated with the optimal strategy is 0.85 as seen in column 7 (ProbPos-max) as compared to values as low as 0.06 (PosProb-min) associated with other strategies.

To investigate the implications of a shorter time horizon, results are summarized for two options (BAC+JV and GE+CF) in Table 2 for days to expiration. The number of days to expiration is indicated by the value at the end of the option ticker in the first column. Examining the first row of the table, we see that for the option BAC+JV with ten days to expiration, the expected accumulated residual has optimal value associated with algorithmic hedging portfolio of 0.3202 (Resid-max), minimum shortfall of (Shortfall-min), and the probability of a positive accumulated residual of 1.0 (PosProb-max). Following the hedging strategy associated with the optimal choice of chosen by the algorithm still provides some advantage as compared to other possible hedging portfolios, although not as much as with a longer time to expiration, as would be expected. The largest improvement in expected accumulated residual for less than 30 days to expiration occurs with the BAC+JV option with 15 days to expiration. Choosing the optimal hedging strategy provides a value 1.4 times as large as the smallest possible expected accumulated residual. For all cases presented, the minimum shortfall produced by the algorithm is very close to zero and the probability of a positive accumulated residual is very close to 1.

Table 2: Algorithmic results for options with varying time to expiration.

#### 10. Evaluation of Numerical Results

The results given in Section 9 document the advantages of using the optimal hedging strategy in comparison to choosing from the range of other market calibrated hedging strategies. To further evaluate the results of our risk minimization algorithm, we take the position of an investor possessing a hedging portfolio based on the optimal strategy identified by the algorithm and compute the accumulated residual value using a representative actual stock price path. To obtain representative stock price values, we divide the collected stock price data into two sets: a large set of historical values used to compute the optimal strategy and a set of values, where is the number of days to expiration of the option, to be interpreted as the actual values of the stock (we use days to expiration for this study).

Figure 7 illustrates this data splitting for a set of Bank of America stock prices from October 15, 2002 through October 15, 2004. To evaluate the optimal hedging portfolio produced by the risk minimization algorithm for an option days to expiration, the stock prices for October 15, 2002 through September 1, 2004 are used in the algorithm to compute the optimal hedging strategy, and stock prices from September 2, 2004 through October 15, 2004 (highlighted in gray) are used as the actual stock price values. The stock price and option price on the first day of the actual stock price path are used as the current stock and option price. The data used as the actual stock path is shown in more detail in the lower graph in Figure 7.

Figure 7: Stock price data.

Having computed the accumulated residual for the single actual path, we compare these to the hedging cost of maintaining the classic Black-Scholes hedge when rebalancing is done daily. Thus, in parallel with the computation of our residuals, we compute the accumulated rebalancing cost of the Black-Scholes hedge, as described in [28]. To be consistent with the method described in [28], we account for the future value (value at expiration) of the option premium received when we take the position of the option seller. It consists of adding the future value of the option premium to the accumulated residuals at expiration. We call the resulting value the accumulated hedging residual. Positive accumulated hedging residuals indicate a gain for the portfolio holder, while a negative accumulated hedging residual indicates a loss for the portfolio holder.

Accumulated residual for both the optimal hedging strategy produced by the algorithm and the Black-Scholes hedge is presented in Table 3. Of the twenty-seven options listed in the table, there are sixteen for which our accumulated hedge residual exceeds that of the Black-Scholes hedge. The results illustrate that following the hedging portfolio strategy produced by the algorithm is beneficial to an investor, producing residual values that exceed the classic Black-Scholes value in 59% of the test cases. The number of favorable cases provides support for the use of our algorithm.

Table 3: Model evaluation results.
Table 4: Options used in data analysis.

#### 11. Conclusions

We have developed an algorithm based on an ARMA time-series model for a stock price jump process that produces a non-self-financing hedging strategy in an incomplete market corresponding to one of several investor risk criteria. The two-stage algorithm optimizes an investor chosen statistical property of the portfolio residual profit or shortfall. The algorithm was tested on a number of options traded on the New York Stock Exchange.

#### Acknowledgments

The authors would like to thank N. N. Leonenko for his helpful discussions and insight as well as for providing us with important references. They would also like to thank the anonymous referees for their suggestions that have improved the quality of the presentation. This work is partially supported by the Bentley Fund for Strategic Research.

#### References

1. N. Josephy, L. Kimball, V. Steblovskaya, A. V. Nagaev, and M. Pasniewski, “An algorithmic approach to non-self-financing hedging in a discrete-time incomplete market,” Discrete Mathematics and Applications, vol. 17, no. 2, pp. 189–207, 2007.
2. A. V. Nagaev and S. A. Nagaev, “Asymptotics of riskless profit under selling of discrete time call options,” Applicationes Mathematicae, vol. 30, no. 2, pp. 173–191, 2003.
3. S. A. Nagaev, A. V. Nagaev, and R. M. Kunst, “A diffusion approximation to the Markov chains model of the financial market and the expected riskless profit under selling of call and put options,” Economics Series 165, Institute for Advanced Studies, Vienna, Austria, 2005.
4. S. A. Nagaev, A. V. Nagaev, and R. M. Kunst, “A diffusion approximation for the riskless profit under selling of discrete time call options: non-identically distributed jumps,” Economics Series 164, Institute for Advanced Studies, Vienna, Austria, 2005.
5. A. V. Nagaev, “A diffusion approximation of the expected risky profit on an investor under selling of discrete time options” (Russian), working paper.
6. G. Tessitore and J. Zabczyk, “Pricing options for multinomial models,” Bulletin of the Polish Academy of Sciences, Mathematics, vol. 44, no. 3, pp. 363–380, 1996.
7. G. Wolczynska, “An explicit formula for option pricing in discrete incomplete markets,” International Journal of Theoretical and Applied Finance, vol. 1, no. 2, pp. 283–288, 1998.
8. O. Hammarlid, “On minimizing risk in incomplete markets option pricing models,” International Journal of Theoretical and Applied Finance, vol. 1, no. 2, pp. 227–233, 1998.
9. J. Staum, “Incomplete markets,” in Handbooks in Operations Research and Management Science, vol. 15, pp. 511–564, Elsevier Science, Amsterdam, The Netherlands, 2006.
10. E. Jouini, J. Cvitanić, and M. Musiela, Eds., Option Pricing, Interest Rates and Risk Management, Handbooks in Mathematical Finance, Cambridge University Press, Cambridge, UK, 2001.
11. A. N. Shiryaev, Essentials of Stochastic Finance: Facts, Models, Theory, vol. 3 of Advanced Series on Statistical Science & Applied Probability, World Scientific, River Edge, NJ, USA, 1999.
12. A. V. Mel'nikov, S. N. Volkov, and M. L. Nechaev, Mathematics of Financial Obligations, vol. 212 of Translations of Mathematical Monographs, American Mathematical Society, Providence, RI, USA, 2002.
13. V. Ryabchenko, S. Sarykalin, and S. Uryasev, “Pricing European options by numerical replication: quadratic programming with constraints,” Asia-Paciffc Financial Markets, vol. 11, no. 3, pp. 301–333, 2004.
14. W. B. Arthur, J. H. Holland, B. LeBaron, R. Palmer, and P. Tayler, “Asset pricing under endogenous expectations in an artificial stock market,” in The Economy as an Evolving, Complex System II, Addison-Wesley, Reading, Mass, USA, 1997.
15. J. Benhabib, Cycles and Chaos in Economic Equilibrium, Princeton University Press, Princeton, NJ, USA, 1992.
16. W. A. Brock and C. Sayers, “Is the business cycle characterized by deterministic chaos?” Journal of Monetary Economics, vol. 22, no. 1, pp. 71–90, 1988.
17. C. W. Granger, “Is chaotic economic theory relevant for economics? A review essay,” Journal of International and Comparative Economics, vol. 3, pp. 139–145, 1994.
18. L. Rüschendorf, “On upper and lower prices in discrete-time models,” Proceedings of the Steklov Institute of Mathematics, vol. 237, pp. 134–139, 2002.
19. R. Carmona, Statistical Analysis of Financial Data in S-Plus, Springer Texts in Statistics, Springer, New York, NY, USA, 2004.
20. E. McKenzie, “An autoregressive process for beta random variables,” Management Science, vol. 31, no. 8, pp. 988–997, 1985.
21. M. Ristić and B. Popović, “A new uniform AR(1) time series model (NUAR(1)),” Publications de l'Institut Mathématique, vol. 68(82), pp. 145–152, 2000.
22. M. R. Chernick, “A limit theorem for the maximum of autoregressive processes with uniform marginal distributions,” The Annals of Probability, vol. 9, no. 1, pp. 145–149, 1981.
23. A. J. Lawrance, “Uniformly distributed first-order autoregressive time series models and multiplicative congruential random number generators,” Journal of Applied Probability, vol. 29, no. 4, pp. 896–903, 1992.
24. L. S. Dewald Sr., P. A. W. Lewis, and E. McKenzie, “Marginally specific alternatives to normal ARMA processes,” in Proceedings of the 19th Winter Simulation Conference (WSC '87), A. Thesen, H. Grant, and W. D. Kelton, Eds., pp. 300–301, Atlanta, Ga, USA, December 1987.
25. C. H. Sim, “Modelling non-normal first-order autoregressive time series,” Journal of Forcasting, vol. 13, no. 4, pp. 369–381, 1994.
26. C. C. Heyde and N. N. Leonenko, “Student processes,” Advances in Applied Probability, vol. 37, no. 2, pp. 342–365, 2005.
27. G. K. Grunwald, R. J. Hyndman, L. Tedesco, and R. L. Tweedie, “Non-Gaussian conditional linear AR(1) models,” Australian & New Zealand Journal of Statistics, vol. 42, no. 4, pp. 479–495, 2000.
28. J. C. Hull, Options, Futures and Other Derivatives, Prentice Hall, Upper Saddle River, NJ, USA, 6th edition, 2005.