Fractional Brownian motion with Hurst exponent H∈(1/2,1) is a good candidate for modeling financial time series with long-range dependence and self-similarity. The main purpose of this paper is to address the valuation of equity indexed annuity (EIA) designs under the market driven by fractional Brownian motion. As a result, this paper presents an explicit pricing expression for point-to-point EIA design and bounds for the pricing of high-water-marked EIA design. Some numerical examples are given to illustrate the impact of the parameters involved in the pricing problems.
1. Introduction
An 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.com's 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 Lévy process. All these processes are independent in increments and thus they are semimartingales, which enable the classical Itô 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 Aít-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 Itô 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={Xk,k∈ℕ} is said to exhibit the statistical long-range dependence, if its autocorrelation function ρ(k)=^Cov(Xn,Xk+n) satisfies
(1)limk→∞ρ(k)cρk-α=1,
for some cρ and α∈(0,1). This is to say that the dependence between Xn and Xk+n decays slowly as k→∞. In particular,
(2)∑k=0∞ρ(k)=∞.
In some literatures, (2) is also referred to as the definition of long-range dependence.
Definition 2 (self-similarity).
A centered stochastic process X={Xt,t∈(0,T)} is said to be statistically self-similar with Hurst exponent H, if
(3){Xt,t∈(0,T)}=d{α-HXαt,t∈(0,T)},
for all α>0. Here, =d 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={Xt,0≤t≤T} with stationary increments admits that
(4)
Cov(Xt,Xs)=12Var(X1)(t2H+s2H-|t-s|2H).
By (4), one can easily find that it is a necessary condition that H∈(0,1) for the process to be well defined.
Lemma 4.
Suppose that self-similar process X={Xk,k∈ℕ+} is centered square integrable process with stationary increments; then the increments
(5)Yk=^Xk-Xk-1
are stationary with autocorrelation function ρ(k)=(1/2)((k+1)2H-2k2H+(k-1)2H). Then, when k→∞, one has
(6)ρ(k)≈H(2H-1)k2H-2.
Thus, only if H∈(1/2,1), the increments {Yk,k∈ℕ+} exhibit long-range dependence.
Let (Ω,𝔉,(ℱt)0≤t≤T,ℙ) be a complete filtered probability space, a fBm BtH with Hurst exponent H∈(0,1) is a continuous, centered Gaussian process with covariance function
(7)Cov(BtH,BsH)=12(t2H+s2H-|t-s|2H),s,t∈(0,T),
where we assumed that 𝔉=ℱT, for some T∈(0,∞), and ℱt=σ(BsH,0≤s≤t≤T) satisfy the usual condition and ℙ is the real world probability.
If H=1/2, then the corresponding fBm is the usual standard Brownian motion. Furthermore, BtH 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 BtH, t∈(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, ω-ω 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 Øksendal [23]. Suppose that f∈C1([0,T]); the following change of variables formula for the path wise integration with respect to the fBm plays important role in what follows:
(8)f(BtH)-f(BsH)=∫stf′(BrH)dBrH.
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=MH, which is given by
(9)MtH=c1∫0tr1/2-H(t-r)1/2-HdBrH;
Now, MH is a martingale with quadratic process
(10)〈MH〉t=c2t2-2H,
where c1 and c2 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 μ(r) be a deterministic function and define a measure ℚ=ℚμ by
(11)dℚdℙ|ℱt=exp{∫0tμ(r)dMrH-12∫0tμ2(r)d〈M〉rH}=exp{∫0tμ(r)dMrH-12c2(2-2H)∫0tμ2(r)r1-2Hdr}.
Then the process
(12)BtH-∫0tμ(r)dr,0≤t≤T
is an fBm, if and only if under measure ℚ.
Lemma 6 (martingale representation for fractional Brownian motion).
Define B~tH=^𝔼[BTH|ℱtBH]; then
(13)B~tH=BtH+∫0tgT(t,s)dBsH,
where
(14)gT(t,s)=sin(π(H-1/2))πs1/2-H∫tTrH-1/2(r-t)H-1/2r-sdr.
3. Pricing EIA3.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 Wt by fBm BtH with index H∈(1/2,1) and then the dynamic of our risky asset is given by
(15)dSt=St(μdt+σdBtH).
The solution to (15) is called the geometric fractional Brownian motion. Similar to the assumption in classical model, we assume that the parameters r, μ, σ, 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
(16)St=S0exp{μt-σ22t2H+σBtH}
and the pricing model is free of arbitrage and complete. However, under the path wise integration framework, the solution to (15) is
(17)St=S0exp{μt+σBtH}
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.
(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.
(Practical needed) Price processes which are not semimartingales are appearing more regularly in the empirical literature estimating stock price processes (see Aít-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 Ø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 St be the value (or price) of an asset at time t that takes interest force r; that is, St are paid dividends rS(t)dt between time t and time t+dt, where r is nonnegative. Let α 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 eαlnST, α>0 or a fixed exercise price K (=elnK), 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 αlnSt, with the provision of a minimum guaranteed rate of return, lnK. Thus, we can express the valuation of this policy under fractal model by
(18)e-rT𝔼ℚ[max(St,K)].
In classical pricing formula, if the market is free of arbitrage, under the risk neutral measure ℚ, the discounted price process e-rtSt is a martingale with respect to the natural sigma filtration generated by Brownian motion {Bt,0≤t≤T}. Particularly, we have
(19)𝔼ℚ[e-rTST]=S0.
In factional setting, we cannot have the martingale property. However, the Girsanov theorem for fraction Brownian motion provides a unique probability measure ℚ, equivalent to ℙ such that (19) holds. By simple calculation, one can find that (19) holds, if process
(20)BtH+μ-rσt+σ2t2H
is the fBm under measure ℚ. By the Grisanov theorem, ℚ is defined as
(21)dℚdℙ=exp{-μ-rσMt+σH∫0ts2H-1dMsxxxxxxxxx-12c2(2-2H)xxxxxxxxx×((μ-rσ)2t2-2H+(μ-r)t+σ2t2H)∫0t}.
Under measure ℚ, we have
(22)𝔼ℚ[f(ST)]=𝔼ℚ[f(S0exp(μT+σBTH))]=𝔼ℚ[f(S0exp(σ22σ[BTH+μ-rσT+σ2T2H]xxxxxxxxxxxxxx+rT-σ22T2H))].
Note that, under measure ℚ, {BtH+((μ-r)/σ)t+(σ/2)t2H, 0≤t≤T} is a fractional Brownian motion; integrating its self-similarity property yields
(23)𝔼ℚ[f(ST)]=e-rT12π×∫-∞∞f(S0exp(σyTH+rT-σ22T2H))e-y2/2dy.
Given r, T, μ, and σ, denote by Pα the valuation of EIA with participation rate α. Without loss of generality, suppose that S0=1. Then, we obtain the valuation formula for EIA under point-to-point design directly by (23). Consider
(24)Pα=𝔼ℚ[max(eαlnST,K)]=𝔼ℚ[max(eα(μT+σBTH),K)]=e-rT12π×∫-∞∞max{exp(α(σyTH+rT-σ22T2H)),K}xxxxxxxx×e-y2/2dy.
Note that
(25)eα(σyTH+rT-(σ2/2)T2H)≤K⟺α(σyTH+rT-σ22T2H)≤lnKxxxxx⟺y≤lnK/ασTH-rσT1-H+σ2TH.
Thus, (24) can be reformulated as
(26)Ke-rTΦ(lnK/ασTH-rσT1-H+σ2TH)+e(α-1)rT+(σ2/2)T2Hα(α-1)×∫lnK/ασTH-(r/σ)T1-H+(σ/2)TH∞e-(y-ασTH)2/2dy=Ke-rTΦ(lnK/ασTH-rσT1-H+σ2TH)+e(α-1)rT+(σ2/2)T2Hα(α-1)×[1-Φ(lnK/ασTH-rσT1-H+σ2TH-ασTH)].
Remark 9.
If we take the dividend payment into account, we let δ be the nonnegative dividend yield rate, that is, the assets pay out dividend δS(t)dt between very small time interval [t,t+dt], then the “pure" riskless income of the asset is (r-δ)S(t)dt. One can easily find that, in this case, the pricing formula for point-to-point design is revived as
(27)Ke-rTΦ(lnK/ασTH-(r-δ)σT1-H+σ2TH)+e[(α-1)r-αδ]T+(σ2/2)T2Hα(α-1)×[1-Φ(lnK/ασTH-r-δσT1-H+σ2TH-ασTH)].
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 β, 0<β≤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
(28)K=βegT,
assuming an initial premium of 1. For most of the EIA designs sold by the insurance companies, β is set at around 90% and g at 3%. Now, we are going to make some observations about Pα. Generally speaking, the price of design Pα is increasing with respect to parameter K, σ, and α. As for r and T, the situations are not clear cut. Pα 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 α puts small value, for example 0.2 or 0.3, there is no obvious difference between two pricing models. While the value of α 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.
Comparison pricing Pα as a function of T given σ=20%, r=5%, δ=2%, and K=0.9e0.03T, for α=0.2,0.3, respectively.
Comparison pricing Pα as a function of T given σ=20%, r=5%, δ=2%, and K=0.9e0.03T, for α=0.6,0.9, respectively.
3.3. High-Water-Marked Design
Another popular design among the EIA products is the look-back or high-water-marked 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
(29)YT=max0≤t≤T(μt+σBtH)=max0≤t≤Tlog(StS0)
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 eαYT,0<α<1 or a fixed minimum guaranteed amount, K, whichever is larger. In this case, at time 0, denote by Ph the value of this policy; then
(30)Ph=𝔼ℚ[max(eαYT,K)]=𝔼ℚ[emax(αYT,lnK)]=𝔼ℚ[emax[αmax0≤t≤T(μt+σBtH),lnK]]=𝔼ℚ[emax[αmax0≤t≤T(σ(BtH+((μ-r)/σ)t+(σ/2)t2H)+rt-(σ2/2)t2H),lnK]]=𝔼ℚ[emax[ασmax0≤t≤T(B~tH+(r/σ)t-(σ/2)t2H),lnK]]=𝔼ℚ[eασmax[max0≤t≤T(B~tH+(r/σ)t-(σ/2)t2H),lnK/ασ]],
where B~tH=^BtH+((μ-r)/σ)t+(σ/2)t2H is a fBm under measure ℚ. 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-(σ2/2)t2H. 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 Z(T)=^max0≤t≤T(B~tH+(r/σ)t-(σ/2)t2H). Denote by F(·) the distribution of Z(T) under measure ℚ and by F¯(·) the tail distribution of F(·); that is, F¯(·)=1-F(·). Then (30) is reformulated as
(31)Kℚ(Z(T)≤lnKσα)+∫lnK/σα∞eσαydF(y)=Kℚ(Z(T)≤lnKσα)-∫lnK/σα∞eσαydF¯(y)=Kℚ(Z(T)≤lnKσα)-KF¯(lnKσα)+σα∫lnK/σα∞eσαyF¯(y)dy=2KF(lnKσα)-K+σα∫lnK/σα∞eσαyF¯(y)dy.
Note that
(32)ℚ(Z(T)≤lnKσα)=ℚ(max0≤t≤T(B~tH+rσt-σ2t2H)≤lnKσα)≤min0≤t≤Tℚ(B~tH+rσt-σ2t2H≤lnKσα)=min0≤t≤Tℚ(B~tH≤lnKσα-rσt+σ2t2H)=min0≤t≤T∫-∞lnK/σα+rt/σ-(σ/2)t2H12πtHe-y2/2t2Hdy=min0≤t≤T∫-∞(lnK/σα)t-H+rt1-H/σ-(σ/2)tH12πe-s2/2ds=min0≤t≤TΦ(lnKσαt-H+rt1-Hσ-σ2tH)=^min0≤t≤Th(lnKσα,t,T),
where Φ(·) is standard normal distribution. Denote by t0(y,T) the minimizer of h(y,t,T); then h(y,t0(y,T),T)|y=lnK/σα is an upper bound of ℚ(Z(T)≤lnK/σα). 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,t0(y,T),T)|y=lnK/σα, but it is easy to get their numerical solution.
On the other hand,
(33)ℚ(Z(T)≤lnKσα)=ℚ(max0≤t≤T(B~tH+rσt-σ2t2H)≤lnKσα)≥ℚ(max0≤t≤T(B~tH+rσt)≤lnKσα)=1-ℚ(max0≤t≤T(B~tH+rσt)>lnKσα).
By the Slepian inequality for fBm with drift (e.g., see Michna [28]), it follows that
(34)ℙ(max0≤t≤T(BtH+μt)>u)≤ℙ(max0≤t≤T(Bt2H1/2+μt)>u),
where Bs1/2 is standard Brownian motion. Denote ρ(u) as
(35)ρ(u)=inf{t>0:Bs1/2+μt1/2H>u};
then
(36)ℙ(max0≤t≤T(Bs2H1/2+μt)>u)=ℙ(max0≤t≤T2H(Bs1/2+μt1/2H)>u)=ℙ(ρ1/2H(u)≤T).
By using the integral equation for the first passage density for Brownian motion to the barrier f(s)=u-μs1/2H (see, e.g., Ferebee [29]), we are able to find a lower density of ρ(u), say p(s;u), which takes the form of
(37)p(s;u)=s-3/2(u-(1-1H)rσs1/2H)×ϕ(us1/2-rσs(1-H)/2H),
where ϕ(·) is the density function of standard normal distribution. By integrating (33), (35), and (37), we have
(38)ℚ(Z(T)≤lnKσα)≥1-ℚ(max0≤t≤T(B~tH+rσt)>lnKσα)≥1-ℙ(max0≤t≤T(Bt2H1/2+rσt)>lnKσα)=1-ℙ(ρ1/2H(lnKσα)<T)≥1-∫0Tp(s;lnKσα)ds.
Together with the definition of h(y,t0(y,T),T)|y=lnK/σα and (38), we have
(39)1-∫0Tp(s;lnKσα)ds≤ℚ(Z(T)≤lnKσα)≤h(y,t0(y,T),T)|y=lnK/σα,1-h(y,t0(y,T),T)≤F¯(y)≤∫0Tp(s;y)ds.
By integrating (30) and (39), we have the following bound estimation for the pricing of high-water-marked design:
(40)H1(K)≤Ph≤H2(K),
where
(41)H1(K)=K[1-2∫0Tp(s;lnKσα)ds]+σα∫lnK/σα∞∫0Tp(s;y)eσαydsdy,H2(K)=[2h(y,t0(y,T),T)|y=lnK/σα-1]K+σα∫lnK/σα∞(1-h(y,t0(y,T),T))dy.
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-water-marked 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 Vasicek model. 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.
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).
MarrionJ.4th quarter index annuity sales2007, http://www.indexannuity.orgTiongS.Valuing equity-indexed annuities20004414916310.1080/10920277.2000.10595945MR1989840ZBL1083.62545BoyleP.TianW.The design of equity-indexed annuities200843330331510.1016/j.insmatheco.2008.05.006MR2479577ZBL1152.91484GerberH. U.ShiuE. S. W.Pricing lookback options and dynamic guarantees2003714867With discussion by Griselda Deelstra10.1080/10920277.2003.10596076MR1986886ZBL1084.91507HardyM.2003Ontario, CanadaJohn Wiley & SonsJaimungalS.Pricing and hedging equity indexed annuities with variance gamma deviates2004, http://www.utstat.utoronto.ca/sjaimung/papers/eiaVG.pdfKijimaM.WongT.Pricing of ratchet equity-indexed annuities under stochastic interest rates200741331733810.1016/j.insmatheco.2006.11.005MR2364558ZBL1141.91457LeeH.Pricing equity-indexed annuities with path-dependent options20033336776902-s2.0-034537550510.1016/j.insmatheco.2003.09.006LinX. S.TanK. S.Valuation of equity-indexed annuities under stochastic interest rates200374729110.1080/10920277.2003.10596119MR2061238ZBL1084.60530MooreK. S.Optimal surrender strategies for equity-indexed annuity investors200944111810.1016/j.insmatheco.2008.08.006MR2492996ZBL1156.91379BiffisE.Affine processes for dynamic mortality and actuarial valuations200537344346810.1016/j.insmatheco.2005.05.003MR2185491ZBL1129.91024BiffisE.DenuitM.DevolderP.Stochastic mortality under measure changesCass Business School Research Paper, 2005HainautD.DevolderP.Mortality modelling with Lévy processes200842140941810.1016/j.insmatheco.2007.05.007MR2392097ZBL1141.91516JalenL.MamonR.Valuation of contingent claims with mortality and interest rate risks2009499-101893190410.1016/j.mcm.2008.10.014MR2532096ZBL1171.91349QianL.WangW.WangR.TangY.Valuation of equity-indexed annuity under stochastic mortality and interest rate201047212312910.1016/j.insmatheco.2010.06.005MR2732209ZBL1231.91446Aít-SahaliaY.Nonparametric pricing of interest rate derivative securities19966435275602-s2.0-0030369366AndrewL.Long-term memory in stock price19915912791313GrangerC. W. J.A typical spectral shape of an economic variable196634150161SottinenT.ValkeilaE.Fractional Brownain motion as a model in financehttp://www.mathstat.helsinki.fi/reports/Preprint302.psWillingerW.TaqquM.TeverovskyV.Stock market prices and long-range dependence19993113DudleyR.NorvaisaR.1998MaPhySto, Department of Mathematical Sciences, University of AarhusDuncanT. E.HuY.Pasik-DuncanB.Stochastic calculus for fractional Brownian motion. I. Theory200038258261210.1137/S036301299834171XMR1741154ZBL0947.60061HuY.ØksendalB.Fractional white noise calculus and applications to finance20036113210.1142/S0219025703001110MR1976868ZBL1045.60072NorrosI.ValkeilaE.VirtamoJ.An elementary approach to a Girsanov formula and other analytical results on fractional Brownian motions19995457158710.2307/3318691MR1704556ZBL0955.60034JarrowR. A.ProtterP.SayitH.No arbitrage without semimartingales200919259661610.1214/08-AAP554MR2521881ZBL1172.60027SonerH. M.ShreveS. E.CvitanićJ.There is no nontrivial hedging portfolio for option pricing with transaction costs199552327355MR1336872ZBL0837.90012NarayanO.Exact asymptotic queue length distribution for fractional Brownian traffic199813963MichnaZ.Self-similar processes in collective risk theory199811442944810.1155/S1048953398000367MR1663600ZBL0937.60029FerebeeB.An asymptotic expansion for one-sided Brownian exit densities198363111510.1007/BF00534172MR699781ZBL0488.60086AsmussenS.2000SingaporeWorld Scientific PublishingMR1794582