# Pricing of equity indexed annuity under fractional Brownian motion model.

1. IntroductionAn indexed annuity is a type of tax-deferred annuity whose credited interest is linked to an equity index--typically the S&P 500 or international index. It guarantees a minimum interest rate (typically between 1% and 3%) if held to the end of the surrender term and protects against a loss of principal. An equity indexed annuity (EIA for short) is a contract with an insurance or annuity company. The returns may be higher than fixed instruments such as CDs, money market accounts, and bonds but not as high as market returns. The guarantees in the contract are backed by the relative strength of the insurer. The contracts may be suitable for a portion of the asset portfolio for those who want to avoid risk and are in retirement or nearing retirement age. The objective of purchasing an equity index annuity is to realize greater gains than those provided by CDs, money markets, or bonds, while still protecting principal.

Because they provide minimum guarantee investment return and lock-in credit rate, EIA designs have been attractive products to customers. Sales have grown dramatically since their introduction in 1995. Indeed, EIA sales for 2007 were nearly 25.2 billion, a 380% increase over their 2000 level of 5.25 billion (see Marrion [1]). According to an industry survey, indexed annuities enjoyed record sales growth in 2010; thirty-nine indexed annuity carriers participated in the 54th edition of AnnuitySpecs.coms Indexed Sales & Market Report, representing 99% of indexed annuity production. Total fourth quarter sales were 8.3 billion, up to 19% from the same period of last year.

The mechanics of EIA are often complex and the returns can vary greatly depending on the month and year that the annuity is purchased. One main method for valuing (pricing) of EIA is the so-called structural method. The idea of the method is to introduce a stochastic process for modeling the behavior of equity involved in the EIA products. Original works on this aspect can be found in Tiong [2]. Stepped works can be found in Boyle and Tian [3], Gerber and Shiu [4], Hardy [5], Jaimungal [6], Kijima and Wong [7], Lee [8], Lin and Tan [9], and Moore [10]. These authors studied the pricing, hedging, and risk management of various features of EIA. There also exist a variety of literatures on EIA valuation by incorporating stochastic interest rate or stochastic mortality into model; for details, see the work of Biffis [11], Biffis et al. [12], Hainaut and Devolder [13], and Jalen and Mamon [14] and references therein.

In aforementioned papers, the dynamic of equity is assumed to be driven by drifted Brownian motion, Brownian motion with Poisson jumps, or Levy process. All these processes are independent in increments and thus they are semimartingales, which enable the classical Ito calculus for semimartingales to be applied for modeling the equity market. Another popular tool for the valuation of EIA is Esscher transform; see the work of Lin and Tan [9], Qian et al. [15], and Tiong [2] and references therein. The purpose of Esscher transform is to obtain a Radon-Nikodym derivative for constructing a new measure. Under the new measure, underlying risky asset is embedded into a risk-neutral world. This is also the idea of famous Black-Scholes-Merton option pricing framework. Esscher transform requires that the return process of risky asset (also named by the logarithm of asset price process) should be independent in increments. However, behavioral finance and econometrics as well as empirical studies found that not all financial data are consistent with this basic characteristic. Over the past decades, vast literatures show that many financial market time series display long-range dependence and momentum and self-similar properties; see the work of Ait-Sahalia [16], Andrew [17], and Granger [18], for example. In this case, the classical financial model with independent increments is invalid, consequently, the approach of deriving risk neutral measure by Esscher transform is not available. Thus, it is natural to propose new model for describing such new kinds of financial time series. By theoretical analysis and empirical test, it turns out that fractional Brownian motion (fBm for short) model is good candidate and consequently, the fractional Black-Scholes models are brought forward. The fractional Black-Scholes model is a generalization of the Black-Scholes model, which is based on replacing the standard Brownian motion by a fractional Brownian motion in financial model.

In this paper we introduce fBm into EIA valuing framework and obtain pricing expressions for two main kinds of EIA designs. This paper is organized as follows. In Section 2, we quickly cite some preliminary results which show why fBm captures the financial data with long-range dependence and self-similarity. Since fBm is not a semimartingale, the traditional Ito calculus method is not valid here. Thus, Section 2 also presents some existing basic calculus results on fBm which will be of great importance in our later discussion. In Section 3, an explicit pricing formula for point-to-point EIA design and bound estimation for high-water-marked EIA design are presented. Some numerical examples are presented for comparing our results with the classical model in this section. In Section 4, we present our conclusions and some remaining problems.

2. Preliminaries: Long-Range Dependence and fBm

This section just serves as a quick survey of the content presented in Sottinen and Valkeila [19] and references therein.

Definition 1 (long-range dependence). A stationary sequence X = {[X.sub.k], k [member of] N} is said to exhibit the statistical long-range dependence, if its autocorrelation function [rho](k) [??] Cov([X.sub.n], [X.sub.k+n]) satisfies

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (1)

for some [c.sub.[rho]] and [alpha] [member of] (0, 1). This is to say that the dependence between [X.sub.n] and [X.sub.k+n] decays slowly as k [right arrow] [infinity]. In particular,

[[infinity].summation over (k=0)] [rho](k) = [infinity]. (2)

In some literatures, (2) is also referred to as the definition of long-range dependence.

Definition 2 (self-similarity). A centered stochastic process X = {[X.sub.t], t [member of] (0, T)} is said to be statistically self-similar with Hurst exponent H, if

[[X.sub.t], t [member of] (0,T)} [??] {[[alpha].sup.-H] [X.sub.[alpha]t], t [member of] (0,T)}, (3)

for all [alpha] > 0. Here, [??] denotes the equivalence in distribution.

Following Lemmas 3 and 4 demonstrates the importance of fractional Brownian motion in modeling financial data with long-range dependence; for details we refer to Sottinen and Valkeila [19] and references therein.

Lemma 3. A square integrable self-similar process X = {[X.sub.t], 0 [less than or equal to] t [less than or equal to] T} with stationary increments admits that

Cov([X.sub.t], [X.sub.s]) = 1/2 Var ([X.sub.1] ([t.sup.2H] + [s.sup.2H] - [[absolute value of t - s].sup.2H]). (4)

By (4), one can easily find that it is a necessary condition that H [member of] (0,1) for the process to be well defined.

Lemma 4. Suppose that self-similar process X = {[X.sub.k], k [member of] [N.sup.+]} is centered square integrable process with stationary increments; then the increments

[Y.sub.k] [??] [X.sub.k] - [X.sub.k-1] (5)

are stationary with autocorrelation function [rho](k) = (1/2)((k + 1)[2.sup.H] - 2[k.sup.2H] + [(k - 1).sup.2H]). Then, when [right arrow] [infinity], one has

[rho](k) [approximately equal to] H (2H - 1) [k.sup.2H-2]. (6)

Thus, only if H [member of] (1/2,1), the increments {[Y.sub.k], k [member of] [N.sup.+]} exhibit long-range dependence.

Let ([OMEGA], T, [([F.sub.t]).sub.0[less than or equal to]t[less than or equal to]T], P) be a complete filtered probability space, a fBm [B.sup.H.sub.t] with Hurst exponent H [member of] (0, 1) is a continuous, centered Gaussian process with covariance function

Cov([B.sup.H.sub.t], [B.sup.H.sub.s]) = 1/2 ([t.sup.2H] + [s.sup.2H] - [[absolute value of t - s].sup.2H]), s, t [member of] (0, T), (7)

where we assumed that T = [F.sub.T], for some T [member of] (0, [infinity]), and [F.sub.t] = [sigma]([B.sup.H.sub.s], 0 [less than or equal to] s [less than or equal to] t [less than or equal to] T) satisfy the usual condition and P is the real world probability.

If H = 1/2, then the corresponding fBm is the usual standard Brownian motion. Furthermore, [B.sup.H.sub.t] has stationary increments and is H self-similar. Valuations of EIA under model driven by standard Brownian motion have been extensively studied, as we mentioned in introduction. If H > 1/2, the process [B.sup.H.sub.t], t [member of] (0, T) exhibits a long-range dependence. As it was mentioned in Willinger et al. [20], data with long-range dependence are widely spread in economics and finance and thus are an active research topic (e.g., see [10] for details). Long-range dependence seems also an important feature that explains the well-documented evidence of volatility persistence and momentum effects. Hereafter, we will only consider the case H > 1/2 in the rest of this paper, which is most frequently encountered in the real financial data.

To define a fractional analogue of the classical Black-Scholes Pricing model, we need to know how to integrate with respect to (w.r.t. for short) fBm as this is connected to hedging. There are two main ways in defining stochastic integration for fBm. One is path wise, that is, [omega]-[omega] stochastic integrals w.r.t. fBm as a refinement of Riemann-Stieltjes integrals by using _p-variation; for details, see Dudley and Norvaisa [21]. Another one is based on white noise analysis and Wick products; for details, see Duncan et al. [22] and, for applications, see Hu and 0ksendal [23]. Suppose that f [member of] [C.sup.1] ([0, T]); the following change of variables formula for the path wise integration with respect to the fBm plays important role in what follows:

f ([B.sup.H.sub.t]) - f([B.sup.H.sub.s]) = [[integral].sup.t.sub.s] f' ([B.sup.H.sub.r]). (8)

In option pricing or EIA valuation, the Girsanov theorem for semimartingales is of great importance for which transforms the underlying asset price process into an risk-neutral world; besides Girsanov theorem, Esscher transform is also one main method for obtaining the risk neutral measure. Since Esscher transform is invalid for fBm, to continue our discussion, analogues of Girsanov theorem for fBm is necessary. We start our discussion by defining the so-called fundamental martingale M = [M.sup.H], which is given by

[M.sup.H.sub.t] = [c.sub.1] [[integral].sup.t.sub.0] [r.sup.1/2-H] [(t - r).sup.1/2-H] d[B.sup.H.sub.r]; (9)

Now, [M.sup.H] is a martingale with quadratic process

[<[M.sup.H]>.sub.t] = [c.sub.2][t.sup.2-2H], (10)

where [c.sub.1] and [c.sub.2] are certain constants depending on Hurst Exponent H. The following two theorems present the Girsanov transform for fBM derived process and martingale representation for fBM; readers are referred to Norros et al. [24] for detailed discussion.

Lemma 5 (Girsanov theorem). Let [mu](r) be a deterministic function and define a measure Q = [Q.sup.[mu]] by

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (11)

Then the process

[B.sup.H.sub.t] - [[integral].sup.t.sub.0] [mu](r)dr, 0 [less than or equal to] t [less than or equal to] T (12)

is an fBm, if and only if under measure Q.

Lemma 6 (martingale representation for fractional Brownian motion). Define [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]; then

[[??].sup.H.sub.t] = [B.sup.H.sub.t] + [[integral].sup.t.sub.0] gT (t, s) d[B.sup.H.sub.s], (13)

where

gT (t, s) = sin ([pi](H-1/2))/[pi] [s.sup.1/2-H] [[integral].sup.T.sub.t] [r.sup.H-1/2][(r - t).sup.H-1/2]/r - s dr. (14)

3. Pricing EIA

3.1. EIA Designs and Model. There are variable EIAs products; in this paper, we discuss the pricing of some most commonly available EIA policies in the market. The simplest one is a point-to-point design, where the policy earns the realized return on the index (or some other risky asset) over a certain period of time at a prescribed participation rate, but with a minimum guarantee. For example, if the realized 5-year return on the S&P 500 is 100% and we assume a participation rate of 80%, then the actual interest credited to the policy will be 80% instead. Even when the market performs poorly, the policy still earns a return of at least the guaranteed rate. Another type of design that we will also touch on in this paper is the look-back or the high-water-marked. This is also referred to as the "no-regret" type of policy, which earns the highest return on the index attained during the life of the policy. The most favorable type of EIA product design seems to be the annual reset. However, due to technical difficulties, we cannot present satisfied framework and results for valuing the annual reset EIA products.

To price an EIA design, one has to bring forward the dynamic of indexed equity. To our knowledge, the risky assets are mostly driven by semimartingales. Due the properties of long-dependence and self-similarity of some financial time series, it is natural and urgent to introduce proper process for modeling such financial data. As it was shown in Section 2, fBm is the ideal continuous time candidate. Thus, we replace the source of randomness of Brownian motion [W.sub.t] by fBm [B.sup.H.sub.t] with index H [member of] (1/2, 1) and then the dynamic of our risky asset is given by

d[S.sub.t] = [S.sub.t] ([mu]dt + [sigma]d[B.sup.H.sub.t]). (15)

The solution to (15) is called the geometric fractional Brownian motion. Similar to the assumption in classical model, we assume that the parameters r, [mu], [sigma], and H are constants.

Remark 7. Since our idea is to apply the option pricing formula for valuing EIA designs, an inevitable problem is whether the pricing model is free of arbitrage and whether the pricing model is complete. Under the Wick product framework of Duncan et al. [22], the solution to (15) is

[S.sub.t] = [S.sub.0] exp {[mu]t - [[sigma].sup.2]/2 [t.sup.2H] + [sigma][B.sup.H.sub.t] (16)

and the pricing model is free of arbitrage and complete. However, under the path wise integration framework, the solution to (15) is

[S.sub.t] = [S.sub.0] exp {[mu]t + [sigma][B.sup.h.sub.t]} (17)

and there exists arbitrage for the pricing model. In 1997, Rogers [7] showed that the fractional Brownian motion could not be used as a price process for a risky security without introducing arbitrage opportunities.

Although under the framework of classical "fundamental asset pricing theorem," fBm is not good candidate for driving process of financial market; it is still meaningful to consider pricing problem under the market driven by fBm. Our reasons are listed as follows.

(i) (Nonarbitrage is possible in a fBm model.) Thanks to the work of Jarrow et al. [25], it is possible to find a general class of processes, which need not be semimartingales that do not permit arbitrage. The idea is to disallow continuous trading and, moreover, to require a minimal fixed time between successive trades. The fixed time can be as small as one likes, but once chosen, it cannot be changed. This disallows a clustering of trades around a fleeting arbitrage opportunity. One should note that continuous trading strategies generate infinite transaction costs under reasonable models of such costs, and those of unbounded variation generate infinite liquidity costs in any finite time interval (c.f. Soner et al. [26]). As such, these trading strategies could never be used in practice. Thus, the method of avoiding arbitrage proposed by Jarrow et al. [25] is naturally accepted in real world practice.

(ii) (Practical needed) Price processes which are not semimartingales are appearing more regularly in the empirical literature estimating stock price processes (see Ait-Sahalia [16], Andrew [17], Granger [18] Dudley and Norvaisa [21], and references therein); so it necessary to introduce fBm in modeling financial time series.

Remark 8. One can compute the price of EIA in path wise fractional model using a weak pricing principle, which coincides with the ones obtained in the generalized pricing model brought forward by Hu and [empty set]ksendal [23]. Although the path wise integration theory is not as satisfactory as the Wick product way, the former one takes its own advantages. Firstly, the functional analytic approach makes it impossible to consider the integrals as almost sure limits of the paths of the process under certain partitions of the integral. It should be noted that this interpretation is possible in Brownian motion. Secondly, and more subjectively, one wants to model the paths properties of the price of the indexed equity. The Wick products integration does not fit well to this aim since the path properties play no central role in the integrals.

3.2. Point-to-Point Design. Here, we shall discuss only the plain point-to-point design, where the index level at maturity is taken simply as the ending index. In all cases the starting index is the prevailing index level when the policy is issued. Let [S.sub.t] be the value (or price) of an asset at time t that takes interest force r; that is, [S.sub.t] are paid dividends rS(t)dt between time t and time t + dt, where r is nonnegative. Let [alpha] be the participation rate, which, in practice, is almost always less than or equal to 1. Suppose that at time T, T > 0, given an initial premium of 1, we have a policy that pays [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] or a fixed exercise price K (=[e.sup.ln K]), whichever is higher. Therefore, at maturity, the policy earns a percentage of the realized return on the asset over T periods (the term of the policy) which is [alpha] ln [S.sub.t], with the provision of a minimum guaranteed rate of return, ln K. Thus, we can express the valuation of this policy under fractal model by

[e.sup.-rT] [E.sup.Q] [max ([S.sub.t],K)]. (18)

In classical pricing formula, if the market is free of arbitrage, under the risk neutral measure Q, the discounted price process [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is a martingale with respect to the natural sigma filtration generated by Brownian motion {[B.sub.t], 0 [less than or equal to] t [less than or equal to] T}. Particularly, we have

[E.sup.Q] [[e.sup.-rT] [S.sub.T]] = [S.sub.0]. (19)

In factional setting, we cannot have the martingale property. However, the Girsanov theorem for fraction Brownian motion provides a unique probability measure Q, equivalent to P such that (19) holds. By simple calculation, one can find that (19) holds, if process

[B.sup.H.sub.t] + [mu] - r/[sigma] t + [sigma]/2 [t.sup.2H] (20)

is the fBm under measure Q. By the Grisanov theorem, Q is defined as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (21)

Under measure Q, we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (22)

Note that, under measure Q, {[B.sup.H.sub.t] + (([mu]-r)/[sigma])t+([sigma]/2)[t.sup.2H], 0 [less than or equal to] t [less than or equal to] T} is a fractional Brownian motion; integrating its self-similarity property yields

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (23)

Given r, T, [mu], and [sigma], denote by [P.sub.[alpha]] the valuation of EIA with participation rate [alpha]. Without loss of generality, suppose that [S.sub.0] = 1. Then, we obtain the valuation formula for EIA under point-to-point design directly by (23). Consider

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (24)

Note that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (25)

Thus, (24) can be reformulated as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (26)

Remark 9. If we take the dividend payment into account, we let [delta] be the nonnegative dividend yield rate, that is, the assets pay out dividend [delta]S(t)dt between very small time interval [t, t + dt], then the "pure" riskless income of the asset is (r - [delta])S(t)dt. One can easily find that, in this case, the pricing formula for point-to-point design is revived as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (27)

By putting H = 1/2, then our model is the famous Black-Scholes pricing model, which is also the indexed equity model in Tiong [2] and (3.4) of Tiong [2] is a special case of (27) by putting H = 1/2.

Let us now assume that the customer is guaranteed a minimum maturity value, which is a percentage [beta], 0 < [beta] [less than or equal to] 1, of the original premium compounded at a minimum guaranteed rate of return g, g > 0, for the duration of the policy. We can incorporate this feature into the policy by setting the exercise price to be

K = [[beta]e.sup.gT], (28)

assuming an initial premium of 1. For most of the EIA designs sold by the insurance companies, p is set at around 90% and g at 3%. Now, we are going to make some observations about [P.sub.[alpha]]. Generally speaking, the price of design [P.sub.[alpha]] is increasing with respect to parameter K, [sigma], and [alpha]. As for r and T, the situations are not clear cut. [P.sub.[alpha]] can be increasing or/and decreasing with respect to r and T, depending on other parameters. Although we can give some intuitive understanding of these relations by numerical examples, just as what has been done in Tiong [2], here, we want to emphasize on the difference of the EIA valuation under the model driven by Brownian motion and the one driven by fBm. Thus, we take the same parameters as the one taken in Figures 1 and 2 of Tiong [2], respectively.

In Figures 1 and 2, the red curves represent the pricing curves under fBm model and blue curves represent the ones under Brownian motion model. From these two figures, it follows easily that when the participation rate a puts small value, for example 0.2 or 0.3, there is no obvious difference between two pricing models. While the value of a is relatively large, for example 0.6 and 0.9, then, two pricing models are quite different. Generally speaking, at the same surrender term, the valuation of fBm model is higher than the one of Brownian motion.

3.3. High-Water-Marked Design. Another popular design among the EIA products is the look-back or high-watermarked method. The idea is that, at maturity, the interest earned on the policy will be based on the growth rate of the highest index value attained during the life of the policy over the index value at the start of the term, which we assume to be one in our calculations. In practice, the method usually looks at the index level at each policy anniversary, and the highest of these is then taken and figured as the index level on the maturity date. Here, we consider the continuous look-back case. Let

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (29)

be the maximum rate of return on the index attained over the time interval [0, T]. We assume that the payoff of the policy at time T is [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] or a fixed minimum guaranteed amount, K, whichever is larger. In this case, at time 0, denote by [P.sub.h] the value of this policy; then

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (30)

where [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] is a fBm under measure Q. Thus, to derive the pricing formula for high-water-marked design, it is sufficient to get the distribution of the maximum of an fBm with nonlinear drift rt - ([[sigma].sup.2]/2)[t.sup.2H]. Unfortunately, to our knowledge, there are no existing results on the topic and it is also not an easy job for us. Alternatively, we dedicate to derive a lower bound and upper bound for the pricing of EIA under high-water-marked design.

Let [MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. Denote by F(x) the distribution of Z(T) under measure Q and by [??](x) the tail distribution of F(x); that is, [??](x) = 1 - F(x). Then (30) is reformulated as

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (31)

Note that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (32)

where [PHI](x) is standard normal distribution. Denote by [t.sub.0](y, T) the minimizer of h(y, t, T); then h(y, [t.sub.0](y, T), T)[|.sub.y=ln K/[sigma][alpha]] is an upper bound of Q(Z(T) [less than or equal to] ln K/[sigma][alpha]). This upper bound seems to be quite rough; however, just as it was proved in Narayan [27] and Norros et al. [24], such bound is very tight. Since the power exponent of t is not an integer, it is not easy to get explicit expression of h(y, [t.sub.0](y, T), T)[|.sub.y=ln K/[sigma][alpha]], but it is easy to get their numerical solution.

On the other hand,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (33)

By the Slepian inequality for fBm with drift (e.g., see Michna [28]), it follows that

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (34)

where [B.sup.1/2.sub.s] is standard Brownian motion. Denote [rho](u) as

[rho](u) = inf {t > 0 : [B.sup.1/2.sub.s] + [[mu]t.sup.1/2H] > u}; (35)

then

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (36)

By using the integral equation for the first passage density for Brownian motion to the barrier f(s) = u - [[mu]s.sup.1/2H] (see, e.g., Ferebee [29]), we are able to find a lower density of [rho](u), say p(s; u), which takes the form of

p(s; u) = [s.sup.-3/2] (u-(1- 1/H)r/[sigma][s.sup.1/2H])

x[phi]([us.sup.1/2] - r/[sigma][s.sup.(1-H)/2H)], (37)

where [phi](x) is the density function of standard normal distribution. By integrating (33), (35), and (37), we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (38)

Together with the definition of h(y, [t.sub.0](y, T), T)|y=ln K/[sigma][alpha] and (38), we have

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (39)

By integrating (30) and (39), we have the following bound estimation for the pricing of high-water-marked design:

[H.sub.1] (K) [less than or equal to] [P.sub.h] [less than or equal to] [H.sub.2] (K), (40)

where

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (41)

4. Remarks and Conclusion

The valuing of high-water-marked EIA designs are highly depends on the distribution of the supremum of drifted (linear or nonlinear) fBm, although the readers can fall back on the results of Narayan [27] and references therein for the asymptotic distribution of the supremum of fBm and their applications. For example, see Willinger et al. [20] for their applications in performance evaluation and see Michna [28] for their applications in insurance. But, to our knowledge, the distribution of the supremum of fBm with nonlinear drift still remained a problem. Even a good estimation for the bounds of supremum distribution is unsolved. To some extent, our bounds estimation for the pricing of EIA under high-watermarked design is meaningful. In fact, bound estimation is also classical research topic in risk theory and risk management, for example, ruin probability; see Asmussen [30].

Usually, there are three kinds of EIA designs: point-to-point, high-water-marked, and annual reset. In this paper, we just focus on the former two designs. Up to now, we are not able to get an explicit expression for pricing formula or to get any proper bounds estimation for such design and thus it is a remaining problem. It is very important to consider stochastic interest rates when pricing EIA. Lin and Tan [9] considered the model of stochastic interest rates by postulating a Vasicek model which is correlated to the geometric Brownian motion of the risky asset. They argued that the effects of stochastic interest rates are crucial in EIA pricing by simulation results. Kijima and Wong [7] adopted the ordinary arbitrage-free pricing principle to price simple and compound annual reset EIA when the short rate follows the extended Vasicekmodel. In the above-mentioned literatures, mortality risk is considered to be deterministic or even not included. However, the expected life length has increased considerably in many countries during the past decades with the advances made in the health sciences and medicine; life insurance and annuities are exposed to unanticipated changes over time in the mortality rates of the appropriate reference population, which has forced life insurers to use a stochastic model to describe mortality laws. There exist a variety of literatures on this topic, for example, Biffis [11], Biffis et al. [12], Hainaut and Devolder [13], Jalen and Mamon [14], and Qian et al. [15]. Since up to now, we have no clear idea on the mechanics of both fBm and the stochastic mortality process or interest process, so we just focus on simple case. But these topics are interesting and worthy of further research.

http://dx.doi.org/10.1155/2014/380718

Conflict of Interests

The authors declare that there is no conflict of interests regarding the publication of this paper.

Acknowledgments

The authors are very grateful to editors and the anonymous referee's advice and help, which improved this paper greatly. Lin Xu was supported by Humanities and Social Science Projects of Ministry of Education of China (12YJC910012) and the National Natural Science Foundation of China (11201006). Guangjun Shen was supported by the National Natural Science Foundation of China (11271020) and the Anhui Provincial Natural Science Foundation (1208085MA11). Dingjun Yao was supported by the National Natural Science Foundation of China (11101205).

References

[1] J. Marrion, "4th quarter index annuity sales," 2007, http:// www.indexannuity.org.

[2] S. Tiong, "Valuing equity-indexed annuities," North American Actuarial Journal, vol. 4, no. 4, pp. 149-163, 2000.

[3] P. Boyle and W. Tian, "The design of equity-indexed annuities," Insurance: Mathematics & Economics, vol. 43, no. 3, pp. 303-315, 2008.

[4] H. U. Gerber and E. S. W. Shiu, "Pricing lookback options and dynamic guarantees," North American Actuarial Journal, vol. 7, no. 1, pp. 48-67, 2003, With discussion by Griselda Deelstra.

[5] M. Hardy, Investment Guarantees: Modeling and Risk Management for Equity-Linked Life Insurance, John Wiley & Sons, Ontario, Canada, 2003.

[6] S. Jaimungal, "Pricing and hedging equity indexed annuities with variance gamma deviates," 2004, http://www.utstat .utoronto.ca/sjaimung/papers/eiaVG.pdf.

[7] M. Kijima and T. Wong, "Pricing of ratchet equity-indexed annuities under stochastic interest rates," Insurance: Mathematics & Economics, vol. 41, no. 3, pp. 317-338, 2007.

[8] H. Lee, "Pricing equity-indexed annuities with path-dependent options," Insurance: Mathematics & Economics, vol. 33, no. 3, pp. 677-690, 2003.

[9] X. S. Lin and K. S. Tan, "Valuation of equity-indexed annuities under stochastic interest rates," North American Actuarial Journal, vol. 7, no. 4, pp. 72-91, 2003.

[10] K. S. Moore, "Optimal surrender strategies for equity-indexed annuity investors," Insurance: Mathematics & Economics, vol. 44, no. 1, pp. 1-18, 2009.

[11] E. Biffis, "Affine processes for dynamic mortality and actuarial valuations," Insurance: Mathematics & Economics, vol. 37, no. 3, pp. 443-468, 2005.

[12] E. Biffis, M. Denuit, and P. Devolder, "Stochastic mortality under measure changes," Cass Business School Research Paper, 2005.

[13] D. Hainaut and P. Devolder, "Mortality modelling with Levy processes," Insurance: Mathematics & Economics, vol. 42, no. 1, pp. 409-418, 2008.

[14] L. Jalen and R. Mamon, "Valuation of contingent claims with mortality and interest rate risks," Mathematical and Computer Modelling, vol. 49, no. 9-10, pp. 1893-1904, 2009.

[15] L. Qian, W. Wang, R. Wang, and Y. Tang, "Valuation of equity-indexed annuity under stochastic mortality and interest rate," Insurance: Mathematics & Economics, vol. 47, no. 2, pp. 123-129, 2010.

[16] Y. Ait-Sahalia, "Nonparametric pricing of interest rate derivative securities," Econometrica, vol. 64, no. 3, pp. 527-560, 1996.

[17] L. Andrew, "Long-term memory in stock price," Econometrica, vol. 59, pp. 1279-1313, 1991.

[18] C. W. J. Granger, "A typical spectral shape of an economic variable," Econometrica, vol. 34, pp. 150-161, 1966.

[19] T. Sottinen and E. Valkeila, "Fractional Brownain motion as a model in finance," Report, http://www.mathstat.helsinki.fi/ reports/Preprint302.ps.

[20] W. Willinger, M. Taqqu, and V Teverovsky, "Stock market prices and long-range dependence," Finance and Stochastics, vol. 3, pp. 1-13, 1999.

[21] R. Dudley and R. Norvaisa, An Introduction to P-Variation and Young Integrals with Emphasis on Sample Functions of Stochastic Processes, MaPhySto, Department of Mathematical Sciences, University of Aarhus, 1998.

[22] T. E. Duncan, Y. Hu, and B. Pasik-Duncan, "Stochastic calculus for fractional Brownian motion. I. Theory," SIAM Journal on Control and Optimization, vol. 38, no. 2, pp. 582-612, 2000.

[23] Y. Hu and B. 0ksendal, "Fractional white noise calculus and applications to finance," Infinite Dimensional Analysis, Quantum Probability and Related Topics, vol. 6, no. 1, pp. 1-32, 2003.

[24] I. Norros, E. Valkeila, and J. Virtamo, "An elementary approach to a Girsanov formula and other analytical results on fractional Brownian motions," Bernoulli, vol. 5, no. 4, pp. 571-587, 1999.

[25] R. A. Jarrow, P. Protter, and H. Sayit, "No arbitrage without semimartingales," The Annals of Applied Probability, vol. 19, no. 2, pp. 596-616, 2009.

[26] H. M. Soner, S. E. Shreve, and J. Cvitanic, "There is no nontrivial hedging portfolio for option pricing with transaction costs," The Annals of Applied Probability, vol. 5, no. 2, pp. 327-355, 1995.

[27] O. Narayan, "Exact asymptotic queue length distribution for fractional Brownian traffic," Advances in Performance Analysis, vol. 1, pp. 39-63, 1998.

[28] Z. Michna, "Self-similar processes in collective risk theory," Journal of Applied Mathematics and Stochastic Analysis, vol. 11, no. 4, pp. 429-448, 1998.

[29] B. Ferebee, "An asymptotic expansion for one-sided Brownian exit densities," Zeitschrift fur Wahrscheinlichkeitstheorie und Verwandte Gebiete, vol. 63, no. 1, pp. 1-15, 1983.

[30] S. Asmussen, Ruin Probabilities, World Scientific Publishing, Singapore, 2000.

Lin Xu, (1) Guangjun Shen, (1) and Dingjun Yao (2)

(1) School of Mathematics and Computer Science, Anhui Normal University, Wuhu, Anhui 241002, China

(2) School of Finance, Nanjing University of Finance and Economics, Nanjing, fiangsu 210046, China

Correspondence should be addressed to Guangjun Shen; gjshen@163.com

Received 13 December 2013; Revised 16 March 2014; Accepted 21 March 2014; Published 24 April 2014

Academic Editor: Yiming Ding

Printer friendly Cite/link Email Feedback | |

Title Annotation: | Research Article |
---|---|

Author: | Xu, Lin; Shen, Guangjun; Yao, Dingjun |

Publication: | Abstract and Applied Analysis |

Article Type: | Report |

Date: | Jan 1, 2014 |

Words: | 5801 |

Previous Article: | Three-step block method for solving nonlinear boundary value problems. |

Next Article: | On differential equations derived from the pseudospherical surfaces. |

Topics: |