Valuation and Risk Assessment of a Portfolio of Variable Annuities: A Vector Autoregression Approach

Affiliation(s)

^{1}
Istituto per le Applicazioni del Calcolo “Mauro Picone” CNR, Napoli, Italy.

^{2}
Department of Statistics and Actuarial Science, Simon Fraser University, Burnaby, Canada.

ABSTRACT

This paper focuses on assessing the financial position of an insurer issuing a portfolio of Variable Annuities (VAs). Two multivariate models for the underlying and the interest rate are considered. The first model uses a single total rate of return for the basket of assets. The second one, jointly models the rates of return on the n assets in the basket. For simplicity, the insurer is assumed to be able to implement a static hedging programme to manage the risk. As an example, a homogeneous portfolio of VAs with GMDB and GMMB guarantees offering different investment opportunities to the policyholders is studied. The insurer can choose to rebalance the basket of assets regularly or not. Results for these two cases are presented.

This paper focuses on assessing the financial position of an insurer issuing a portfolio of Variable Annuities (VAs). Two multivariate models for the underlying and the interest rate are considered. The first model uses a single total rate of return for the basket of assets. The second one, jointly models the rates of return on the n assets in the basket. For simplicity, the insurer is assumed to be able to implement a static hedging programme to manage the risk. As an example, a homogeneous portfolio of VAs with GMDB and GMMB guarantees offering different investment opportunities to the policyholders is studied. The insurer can choose to rebalance the basket of assets regularly or not. Results for these two cases are presented.

1. Introduction

Over the years, many practical and academic contributions have been offered describing VAs and their embedded guarantees.

Variable annuities (VAs) are life insurance products with investment guarantees providing substantial investment guarantees along with tax privileges and prudently managed assets. VAs have existed in USA since the 1950s. Recently they have been spreading across Europe. Some of the more significant and high profile launches have been AXA in France, Germany, Spain, Italy and Belgium as well as ING launches in Spain, Hungary and Poland, Generali launch (December 2007) in Italy and Ergo launch (February 2008) in Germany. This is in addition to the various launches by Aegon, Hartford, Metlife and Lincoln in the U.K. [1] .

VAs allow the allocation of premiums into a range of investment options which usually contain stocks, bonds, money market instruments or some combination of the three. The benefits to the policyholder will depend on the performance of the chosen investment option; typically, the benefit is the higher of the account value and the guaranteed amount.

The major types of VA guarantees are: Guaranteed Minimum Death Benefits (GMDBs) that guarantee a return of the principal invested upon the death of the policyholder; Guaranteed Minimum Accumulation Benefits (GMABs) where the guarantees typically bite on specified policy anniversaries or between specified dates if the policy is still in force. If guarantees are available at maturity they are called Guaranteed Minimum Maturity Benefits (GMMBs); Guaranteed Income Benefits (GMIBs) guarantee a minimum income stream (typically in the form of a life annuity) from a specified future point in time; Guaranteed Minimum Withdrawal Benefits (GMWBs) guarantee a minimum income stream through regular withdrawals from the account.

A detailed overview of variable annuities and a description of the VA international market are given in [1] and [2] . The main emphasis of the current literature is on pricing and hedging the options embedded in these life insurance contracts. A general framework in which any design of options and guarantees currently offered within VAs can be modeled, is presented in [3] . In the recent literature we find some interesting contributions on identification and assessment of risk. Among others, [4] deals with risk assessment, mostly via stochastic simulation. Their main contribution is to present an integrated approach to risks in VA products. Analytical methods for the calculation of risk measures for variable annuities guaranteed benefits and a quantitative model for market risk on the liability side and revenue risk in the asset side, are proposed in [5] .

The insurer will be holding the liability for the VAs guarantee until lapse, death or maturity. The impact of both mortality risks (due to random fluctuations in the number of deaths) and investment risks can be rather high in the presence of some types of guarantees and options. The resulting economic impact to the insurer selling the guarantee is based on how the future real world scenario plays out. Given these considerations, it may be useful for an insurer to assess real world risk associated with embedded guarantees of its VA portfolio, once pricing is done and a risk management strategy is adopted.

Most insurers manage the risk of their VA portfolios by adopting a sophisticated hedging programme. Some use static hedging to mitigate the risk, at least partially. Here, static hedging refers to strategies consisting of buying and selling long-term financial instruments, e.g. over the counter (OTC) options, at time of issue and holding them until maturity with little rebalancing, if any. Available static hedges often fail to closely match the liabilities of VA products and are expensive to set up. For these reasons, dynamic hedging strategies are more popular among insurers. Here, a portfolio of liquid assets is set up so that, over a short period, its change in value will approximately match the change in value of the liabilities. The portfolio of assets is then rebalanced at a frequency that takes into account trading costs and effectiveness of the hedge.

In this paper, it is assumed that the insurer is able to implement a static hedge. This allows us to compare the impact of the two multivariate processes used to model the rate of return of the basket of assets on the value of the portfolio. We can also study whether rebalancing the basket of assets held in the VA account or not affects the risk. Although using a dynamic hedging strategy is more realistic, we consider a static hedging program as a first step of our research which is not specifically focused on the kind of strategy chosen to hedge the portfolio. Dynamic hedging will be explored in future work.

In this paper we choose to focus on assessing the financial position at issue of an insurer selling a portfolio of VA products offering GMDB and GMMB guarantees. We measure the portfolio value at inception as the difference between the present values of the stochastic future inflows and outflows in a real world scenario.

Our main contribution consists of considering two multivariate VAR (1) processes to model the underlying and the interest rate and compare the portfolio values at inception.

Although variable annuities are written on more than one asset, in practice, single asset univariate stochastic investment models are mostly used for pricing and hedging purposes. This may lead to problems such as basis risk or a lack of flexibility with respect to asset allocations [6] . In light of these considerations, we compare a two-dimensional VAR (1) model, where the asset returns on each fund are combined in a single time series, to an (n + 1)-dimensional VAR (1) model where the asset returns on each fund are modelled separately. A special case of these models is proposed to better capture the mutual dependence of the variables involved.

Another issue for the insurer is the decision to rebalance or not the basket of assets constituting the fund itself. Rebalancing is the process of realigning the weightings of a portfolio of assets. It involves periodically buying and selling assets in a portfolio to maintain a specific asset allocation. In our work, we consider two different accumulation functions depending on whether the basket of assets is rebalanced or not. For sake of simplicity, no transaction costs are considered.

Finally, we assume that the VA writer offers to policyholders different premium allocations with different investment styles. These kind of Funds characterize some VA products sold in Europe (among others, the “Accumulator” sold by AXA).

The different cases described above are evaluated in a simulation framework. The results under static hedging are compared with those obtained when the insurer decides not to hedge.

The remainder of the paper is organized as follows. Section 2 introduces the product, Section 3 describes the VAR (1) models employed, Section 4 is devoted to the portfolio evaluation at inception. In Section 5, the results are presented and discussed. In Section 6, some concluding remarks are given.

2. The Contract

Let us consider a portfolio of VA contracts offering GMDB and GMMB guarantees issued to C independent lives aged x, each paying a unique premium P at issue.

The insurer invests the premiums, net of initial charges for general expenses, into several internal funds characterized by different investment styles.

Policyholders may choose in which funds to invest.

By virtue of the GMDB guarantee, if the insured dies in policy year t, the insurer pays a sum equal to the greater of the guaranteed amount or the fund value.

On the other hand, by virtue of the GMMB, at maturity T the insurer pays the greater of the guaranteed amount or the fund value, if the insured is still alive. The benefits are paid monthly. A monthly management charge is deducted and the money gathered by the insurer through this charge is used to cover the guarantees and other expenses linked to the Fund such as fund management and taxes.

The obligations the insurer has to front for the GMDB at time t are valued

$GMD{B}_{t}={N}_{D}\left(t\right)\cdot Max\left[{F}_{t},{G}_{t}\right]$ (1)

with $t\in \left[\mathrm{1,2,}\cdots \mathrm{,}T\right]$ .

For the GMMB we have:

$GMM{B}_{T}={N}_{S}\left(T\right)\cdot Max\left[{F}_{T},{G}_{T}\right]$ (2)

where ${N}_{D}\left(t\right)$ is the actual number of deaths in $\left[t-\mathrm{1,}t\right]$ and ${N}_{S}\left(T\right)$ is the actual number of survivors at time T.

According to the contractual design, the fund value used in (1) and (2) is:

${F}_{t}={P}^{\prime}\cdot {\left(1-m\right)}^{t}{e}^{f\left(t\right)}$ (3)

and the guaranteed value is:

${G}_{t}={P}^{\prime}\cdot {e}^{g\cdot t}$ (4)

with $t\in \left[\mathrm{1,2,}\cdots \mathrm{,}T\right]$ .

Here, ${F}_{t}$ is the fund value obtained by investing the net premiums ${P}^{\prime}=P-c$ into the chosen internal funds, c is the initial charge and P the single premium. ${e}^{f\left(t\right)}$ is the accumulation function and m the monthly management charge deducted. ${G}_{t}$ is the guaranteed amount obtained by accumulating the net premium at the guaranteed rate g (roll up).

Distinct funds are offered, each made up of a mix of several assets. The combination of the assets gives rise to various investment styles.

Moreover, with reference to the fund value ${F}_{t}$ in (3), we consider two cases: the rebalancing strategy of the basket of assets constituting the fund itself and the no rebalancing one (see Section 3).

By means of the well known put decomposition principle, we can rewrite the guarantees (1) and (2) as follows:

$GMD{B}_{t}={N}_{D}\left(t\right)\cdot \left({F}_{t}+{\left[{G}_{t}-{F}_{t},0\right]}_{+}\right)$ (5)

with $t\in \left[\mathrm{1,2,}\cdots \mathrm{,}T\right]$ and

$GMM{B}_{T}={N}_{S}\left(T\right)\cdot \left({F}_{T}+{\left[{G}_{T}-{F}_{T},0\right]}_{+}\right)$ . (6)

(5) and (6) imply that upon the death of the policyholder or expiration of the contract, the insurance Company pays the accumulated fund value ${F}_{t}$ to the investor plus an additional payoff of a put option with strike price equal to ${G}_{t}$ with $t\in \left[\mathrm{1,2,}\cdots \mathrm{,}T\right]$ and ${G}_{T}$ respectively.

Therefore if the fund performs so poorly that its value, when the contract matures or when death occurs, is below the corresponding guaranteed values, the insurer pays the difference.

For the GMDB, the insurer’s monthly payments are:

${L}_{t}^{GMDB}={N}_{D}\left(t\right)\cdot {\left[{G}_{t}-{F}_{t}\mathrm{,0}\right]}_{+}$ . (7)

with $t\in \left[\mathrm{1,2,}\cdots \mathrm{,}T\right]$ .

On the other hand for the GMMB we have:

${L}_{T}^{GMMB}={N}_{S}\left(T\right)\cdot {\left[{G}_{T}-{F}_{T}\mathrm{,0}\right]}_{+}$ (8)

The random liabilities evaluated at time zero are then given by:

${L}_{0}^{GMDB}={\displaystyle \underset{t=1}{\overset{T}{\sum}}}\text{\hspace{0.05em}}{N}_{D}\left(t\right)\cdot {\left({G}_{t}-{F}_{t},0\right)}_{+}\cdot {e}^{-y\left(t\right)}$ (9)

${L}_{0}^{GMMB}={N}_{S}\left(t\right)\cdot {\left({G}_{T}-{F}_{T},0\right)}_{+}\cdot {e}^{-y\left(T\right)}$ (10)

where ${e}^{-y\left(t\right)}$ $t\in \left[\mathrm{1,2,}\cdots \mathrm{,}T\right]$ is a random discounting factor.

3. Modeling the Financial Variables: The Vector Autoregressive Model

We choose a multivariate framework modelling jointly all financial variables, namely asset returns and interest rates, by means of a Vector Autoregressive model of order one, VAR (1).

The VAR model is one of the most flexible and easy to use models for analyzing multivariate time series. It is a natural extension of the univariate autoregressive model to study multivariate time series. The VAR model has proven to be especially useful for describing the dynamic behaviour of economic and financial time series and for forecasting [7] .

An n-dimensional time series random variable $\underset{\_}{{X}_{s}}$ is said to follow a VAR (1) model if

$\underset{\_}{{X}_{s}}-\underset{\_}{\mu}=\underset{\_}{\Phi}\cdot \left(\underset{\_}{{X}_{s-1}}-\underset{\_}{\mu}\right)+\underset{\_}{{a}_{s}}$ (11)

where: $s\in Z$ , $\underset{\_}{\mu}$ is an $\left(n\times 1\right)$ matrix of the long term mean values of $\underset{\_}{{X}_{s}}$ , $\underset{\_}{\Phi}$ is an $\left(n\times n\right)$ matrix of autoregressive coefficients, $\underset{\_}{{a}_{s}}$ is an n-dimensional white noise term which follows a multivariate normal distribution with zero mean and covariance matrix ${\underset{\_}{\Sigma}}_{a}$ .

A sufficient condition of stationarity is that all eigenvalues of $\underset{\_}{\Phi}$ are less than 1 in absolute value.

When conditioning on the initial value $\underset{\_}{{X}_{0}}$ , the mean of $\underset{\_}{{X}_{s}}$ is

$E\left[\underset{\_}{{X}_{s}}|\underset{\_}{{X}_{0}}\right]=\underset{\_}{{\Phi}^{s}}\cdot \left(\underset{\_}{{X}_{0}}-\underset{\_}{\mu}\right)+\underset{\_}{\mu}$ (12)

The covariance between $\underset{\_}{{X}_{s}}$ and $\underset{\_}{{X}_{s-k}}$ , where k is a non negative integer, $k<s$ , also conditional on $\underset{\_}{{X}_{0}}$ is:

$Cov\left[\underset{\_}{{X}_{s}},\underset{\_}{{X}_{s-k}}|\underset{\_}{{X}_{0}}\right]={\displaystyle \underset{i=0}{\overset{s-k-1}{\sum}}}\underset{\_}{{\Phi}^{s+i}}\cdot {\underset{\_}{\Sigma}}_{a}\cdot {\underset{\_}{\Phi}}^{i}$ (13)

We consider two models: a two-dimensional VAR (1) model for the composite rate of return of the n assets and the interest rate and a (n + 1)-dimensional VAR (1) for the rates of return and the interest rate.

3.1. Two-Dimensional VAR (1) Model

In this case we model the composite rate of return ${\delta}_{s}$ of the n assets constituting each fund and the interest rate ${r}_{s}$ by means of a two-dimensional VAR (1) model.

Referring to (11) we can write:

$\underset{\_}{{X}_{s}}=\left(\begin{array}{c}{\delta}_{s}\\ {r}_{s}\end{array}\right)$ (14)

$\underset{\_}{\Phi}=\left(\begin{array}{cc}{\Phi}_{1,1}& {\Phi}_{1,2}\\ {\Phi}_{2,1}& {\Phi}_{2,2}\end{array}\right)$ (15)

$\underset{\_}{{a}_{s}}=\left(\begin{array}{c}{a}_{{\delta}_{s}}\\ {a}_{{r}_{s}}\end{array}\right)$ (16)

The composite rate of return ${\delta}_{s}$ can be computed in two ways, depending on whether the basket of assets is rebalanced or not. When rebalancing the basket of assets we have:

${\delta}_{s}={\displaystyle \underset{i=1}{\overset{n}{\sum}}}\text{\hspace{0.05em}}{\omega}_{i}\cdot {\delta}_{s}^{i}$ (17)

with no rebalancing, the composite rate of return is:

${\delta}_{s}={\displaystyle \underset{i=1}{\overset{n}{\sum}}}\left({\displaystyle \underset{k=1}{\overset{s-1}{\prod}}}\left(1+{\delta}_{s-k}^{i}\right)*{\delta}_{s}^{i}\right)$ (18)

where ${\delta}_{1}={\displaystyle {\sum}_{i=1}^{n}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\omega}_{i}\cdot {\delta}_{1}^{i}$ and $t\in 1,2,\cdots ,T$ .

In (17) and (18) ${\delta}_{s}^{i}$ with $i\in 1,2,\cdots ,n$ are the monthly continuous rates of return of the n assets and ${\omega}_{i}$ $i\in 1,2,\cdots ,n$ are the proportion of the fund invested in each of the n assets.

If we model the financial variables by a two-dimensional VAR (1) model, ${F}_{t}$ in (3) can be rewritten as:

${F}_{t}={P}^{\prime}\cdot {\left(1-m\right)}^{t}\cdot {e}^{{\displaystyle {\sum}_{s=1}^{t}{\delta}_{s}}}$ (19)

with $t\in \left[1,2,\cdots ,T\right]$ .

In (19), ${\delta}_{s}$ is given by (17) or (18) depending on whether the fund is rebalanced or not. ${F}_{t}$ assumes different values depending on the selected investment strategy.

The discounting factor is given by:

${e}^{-y\left(t\right)}={e}^{-{\displaystyle {\sum}_{s=1}^{t}{r}_{s}}}$ (20)

with $t\in \left[1,2,\cdots ,T\right]$ .

3.2. (n + 1)-Dimensional VAR (1) Model

Here we consider a (n + 1)-dimensional VAR (1), modelling jointly the rates of return of the n assets in which the insurer invests ( ${\delta}_{s}^{i}$ with $i\in 1,2,\cdots ,n$ ) and the discounting interest rate ( ${r}_{s}$ ).

Referring to (11) we can write:

$\underset{\_}{{X}_{s}}={\left({\delta}_{s}^{1}\mathrm{,}{\delta}_{s}^{2}\mathrm{,}{\delta}_{s}^{3}\mathrm{,}\cdots ,{\delta}_{s}^{n}\mathrm{,}{r}_{s}\right)}^{\text{T}}$ (21)

$\underset{\_}{\Phi}=\left(\begin{array}{ccccc}{\Phi}_{1,1}& {\Phi}_{1,2}& \cdots & {\Phi}_{1,n}& {\Phi}_{1,n+1}\\ {\Phi}_{2,1}& {\Phi}_{2,2}& \cdots & {\Phi}_{2,n}& {\Phi}_{2,n+1}\\ \vdots & \vdots & \ddots & \vdots & \vdots \\ {\Phi}_{n,1}& {\Phi}_{n,2}& \cdots & {\Phi}_{n,n}& {\Phi}_{n,n+1}\\ {\Phi}_{n+1,1}& {\Phi}_{n+1,2}& \cdots & {\Phi}_{n+1,n}& {\Phi}_{n+1,n+1}\end{array}\right)$ (22)

$\underset{\_}{{a}_{s}}={\left({a}_{{\delta}_{s}^{1}},{a}_{{\delta}_{s}^{2}},{a}_{{\delta}_{s}^{3}},\cdots ,{a}_{{\delta}_{s}^{n}},{a}_{{r}_{s}}\right)}^{\text{T}}$ . (23)

Based on some observations on the estimated parameters, we also propose a modified version of the (n + 1)-VAR (1) which includes a Geometric Brownian Motion (GBM) to model the stock indexes (see Section 5).

The use of the (n + 1)-dimensional model implies that we have different values of ${F}_{t}$ in (3), in case of rebalancing or not the basket of assets

1) Rebalancing case:

${F}_{t}={P}^{\prime}\cdot {\left(1-m\right)}^{t}\cdot {e}^{{\displaystyle {\sum}_{s=1}^{t}{\displaystyle {\sum}_{i=1}^{n}{\omega}_{i}\cdot {\delta}_{s}^{i}}}}$ (24)

2) No rebalancing case:

${F}_{t}={\displaystyle \underset{i=1}{\overset{n}{\sum}}}\text{\hspace{0.05em}}\text{\hspace{0.05em}}{\omega}_{i}\cdot {P}^{\prime}\cdot {\left(1-m\right)}^{t}\cdot {e}^{{\displaystyle {\sum}_{s=1}^{t}{\delta}_{s}^{i}}}$ (25)

where ${\varpi}_{i}$ $i\in \left[\mathrm{1,2,}\cdots \mathrm{,}n\right]$ are the percentages of capital invested by the insurer in each of the n assets.

In both cases, according to the values of ${\omega}_{i}$ , we obtain different accumulation factors depending on the chosen investment style.

The discounting factor is given by:

${e}^{-y\left(t\right)}={e}^{-{\displaystyle {\sum}_{s=1}^{t}{r}_{s}}}$ (26)

with $t\in \left[1,2,\cdots ,T\right]$ .

4. Hedging the Guarantee

Bearing in mind that the VA contracts embed equity put options, a crucial problem is to hedge these options.

Hedging the put options means holding a hedge portfolio replicating the guarantee payoff. A possible solution is dynamic hedging which requires a continuous rebalancing of the hedge portfolio.

In this framework, we consider the case of static hedging; for example, the insurer buys equity puts on the market place.

4.1. The Cost of the Hedging Portfolio

In case of static hedging, the insurer buys a portfolio of put options, at time zero, to hedge the guarantee payoff. We know that, in this case, the option price may be too high because these options differ in nature and term from the typical over the counter options currently available. Moreover, the counter party risk may be a factor, given the long terms involved.

Static hedging is not always possible, neverthless this is a first step of our study. Future research will consider dynamic hedging, too.

We price the replicating portfolio by adapting the Black-Scholes results for the guarantees embedded in the VA contract like in [8] while considering the different assumptions we make about the financial variables.

The term to maturity of each option in the VA contract is random, depending on the remaining lifetime of the policyholders. For the GMDB guarantee it depends on the monthly death probabilities and for the GMMB it depends on the survival probability at the maturity of the contract.

Therefore, the expected cost of the hedge portfolio at time zero is:

$P{P}_{0}={\displaystyle \underset{t=1}{\overset{T}{\sum}}}\text{\hspace{0.05em}}{E}_{D}\left(t\right)\cdot BS{P}_{0}\left(t\right)+{E}_{S}\left(T\right)\cdot BS{P}_{0}\left(T\right)$ (27)

where ${E}_{D}\left(t\right)$ is the expected number of deaths at the end of the month t and ${E}_{S}\left(t\right)$ is the expected number of survivors at the maturity date of the contract T.

We assume that ${\left\{{N}_{D}\left(t\right)\right\}}_{t=1}^{T}\cup \left\{{N}_{S}\left(T\right)\right\}$ is multinomial with parameters

$\left(C\mathrm{;}{q}_{x}\mathrm{,}\cdots \mathrm{,}{}_{T-\mathrm{1|}}q{}_{x}{\mathrm{,}}_{n}{p}_{x}\right)$ where C is the number of policies issued at time zero, ${}_{t\mathrm{|}}{q}_{x}$ is the probability that a life aged x at issue dies in $\left[t-1,t\right]$ and ${}_{T}{p}_{x}$ is the probability that a life aged x at issue is alive at time T.

So we have that ${E}_{D}\left(t\right)=C\cdot {}_{t-\mathrm{1|}}q{}_{x}$ and ${E}_{S}\left(t\right)=C\cdot {}_{T}p{}_{x}$

As aforementioned, we calculate the quantity $BS{P}_{0}$ in (27) by an adaptation of the standard Black-Scholes formula.

We need to modify the standard formula for the put options to reflect the fact that the underlying asset is not the stock price itself, but the fund value, and this differs from the stock price through the deduction of the management charges.

Net of the monthly management charges the fund is:

${F}_{t}={F}_{0}\cdot \frac{{S}_{t}}{{S}_{0}}\cdot {\left(1-m\right)}^{t}$ (28)

where ${S}_{t}$ is the stock index.

Let ${F}_{0}={S}_{0}$ , then the option price with maturity t is:

$BS{P}_{0}\left(t\right)={e}^{-r\cdot t}\cdot {E}_{Q}\left[{\left({G}_{t}-{S}_{t}\cdot {\left(1-m\right)}^{t}\right)}_{+}\right]$ (29)

Using this approach, we replace ${S}_{0}$ by ${S}_{0}\cdot {\left(1-m\right)}^{t}$ in the standard Black-Scholes formula. Then the put option price is computed as follows:

$BS{P}_{0}\left(t\right)={G}_{t}\cdot {e}^{-r\cdot t}\cdot \Phi \left(-{d}_{2}\right)-{S}_{0}\cdot {\left(1-m\right)}^{t}\cdot \Phi \left(-{d}_{1}\right)$ (30)

where r is the risk-free (force) of interest, $\Phi (\u2022)$ denotes the cumulative normal density function, ${G}_{t}={P}^{\prime}\cdot {e}^{g\cdot t}$ is the increasing strike price, m is the monthly charge and

${d}_{1}=\frac{\mathrm{log}\left({S}_{0}\cdot {\left(1-m\right)}^{t}/{G}_{t}\right)+\left(r+{\sigma}_{S}^{2}/2\right)\cdot t}{\sigma \cdot \sqrt{t}},$ (31)

${d}_{2}={d}_{1}-\sigma \cdot \sqrt{t},$ (32)

${\sigma}_{S}^{2}$ is the estimated volatility of the asset return.

If we model the financial variables by a two-dimensional VAR model, ${\sigma}_{S}^{2}$ is the estimated volatility of the weighted continuously compounded rate of return ${\delta}_{t}$ . It assumes different values whether the basket of assets is rebalanced or not (Section 3.1)

Modelling the assets by a $\left(n+1\right)$ -dimensional VAR (1) model, ${\sigma}_{S}^{2}$ is the weighted volatility of the n returns on assets, that is ${\sigma}_{S}^{2}={\displaystyle {\sum}_{i=1}^{n}}{\displaystyle {\sum}_{j=1}^{n}}\text{\hspace{0.05em}}{\omega}_{i}{\omega}_{j}{\sigma}_{i,j}$ .

Let us now consider the case in which the initial allocation does not change and the financial variables are modelled by a $\left(n+1\right)$ -dimensional VAR (1) model.

In this case, at time zero, the insurer buys a portfolio of basket options.

As it is well known, a basket option is an option on a basket of assets, typically stocks.

In this case, the Fund value in (28) can be modified as follows:

${F}_{t}={\displaystyle \underset{i=1}{\overset{n}{\sum}}}\text{\hspace{0.05em}}{\omega}_{i}\cdot {F}_{0}^{i}\cdot \frac{{S}_{t}^{i}}{{S}_{0}^{i}}\cdot {\left(1-m\right)}^{t}$ (33)

where ${S}_{t}^{i}$ is the i-th stock index, ${\omega}_{i}$ is the proportion of the fund invested in asset i.

Letting $C\cdot {F}_{0}^{i}={S}_{0}^{i}$ , we can write:

${F}_{t}={\displaystyle \underset{i=1}{\overset{n}{\sum}}}\text{\hspace{0.05em}}{\omega}_{i}\cdot {S}_{t}^{i}\cdot {\left(1-m\right)}^{t}.$ (34)

Our interest is in the dynamics of the weighted arithmetic sum of the underlying ${S}_{t}^{i}\left(i=1,2,3,\cdots ,n\right)$ . Similar to a plain vanilla option, the value of a basket option is equal to its expected payoff discounted at the risk free rate. The expectation is computed with respect to the state-price density function (SPD) also known as the risk neutral probability density function. The SPD, in this case, is not known and closed form solutions to the problem are not available. A common approach is to estimate the state density function by a lognormal functional form, using different matching techniques. One way is to compute the first two moments of the basket’s payoff structure at maturity and then match them to those of the lognormal distribution.

Following [9] , we define the “pseudo-forward” time t of the basket:

${A}_{t}={\displaystyle \underset{i=1}{\overset{n}{\sum}}}\text{\hspace{0.05em}}{\omega}_{i}{A}_{t}^{i}$ (35)

where ${A}_{t}^{i}={S}_{0}^{i}\cdot {\left(1-m\right)}^{t}\cdot {e}^{r\cdot t}$ , being r the risk free rate.

Dividing ${F}_{t}$ by ${A}_{t}$ and denoting the ratio by ${F}_{t}^{\mathrm{*}}$ , we normalize it to have a mean of 1 (i.e. ${M}_{1}=1$ ).

The second moment of ${F}_{t}^{\mathrm{*}}$ is:

${M}_{2}=\frac{1}{{A}_{t}^{2}}\cdot {\displaystyle \underset{i=1}{\overset{n}{\sum}}}{\displaystyle \underset{j=1}{\overset{n}{\sum}}}\text{\hspace{0.05em}}{\omega}_{i}{\omega}_{j}{A}_{i}{A}_{j}{e}^{{\rho}_{i,j}{\sigma}_{i}{\sigma}_{j}t}$ (36)

where ${\rho}_{i\mathrm{,}j}$ is the correlation coefficient between ${A}_{i}$ and ${A}_{j}$ .

Assuming that the normalized basket ${F}_{t}^{\mathrm{*}}$ is lognormal with mean ${M}_{1}=1$ , it follows that the variance of the basket is

$v=\mathrm{ln}\left({M}_{2}\right)$ (37)

The moment-matched lognormal state-price density function gives the price of the put option, by means of the usual Black-Scholes formula:

$LNBS{P}_{0}\left(t\right)={G}_{t}\cdot {e}^{-r\cdot t}\cdot \Phi \left(-{d}_{2}\right)-{A}_{t}\cdot \Phi \left(-{d}_{1}\right)$ (38)

where $\Phi (\u2022)$ is the cumulative normal density function and

${d}_{1}=\frac{\mathrm{log}\left({A}_{t}\cdot {\left(1-m\right)}^{t}/{G}_{t}\right)+v/2}{\sqrt{v}}$ (39)

${d}_{2}={d}_{1}-sqrtv$ (40)

Equation (27) becomes

$P{P}_{0}={\displaystyle \underset{t=1}{\overset{T}{\sum}}}\text{\hspace{0.05em}}{E}_{D}\left(t\right)\cdot LNBS{P}_{0}\left(t\right)+{E}_{S}\left(T\right)\cdot LNBS{P}_{0}\left(T\right)$ (41)

4.2. The Portfolio Value at Inception

The portfolio value at time zero is given by the difference between the random inflows and outflows arising from the contract.

Let us now consider the inflows connected to the contract.

They are given by the expense charge paid at issue by each insured as a percentage c of the premium P plus the present value of the collected monthly charges m to the insured’s account. Denoting by $I{F}_{0}$ the inflows evaluated at time zero, in the case of a two-dimensional VAR (1) model and referring to (19) and (20) we have:

$\begin{array}{c}I{F}_{0}=C\cdot P\cdot c+{\displaystyle \underset{t=1}{\overset{T}{\sum}}}\left({N}_{D}\left(t\right)\cdot \left(1-{\left(1-m\right)}^{t}\right)\cdot {P}^{\prime}\cdot {e}^{{\displaystyle {\sum}_{s=1}^{t}{\delta}_{s}}}\cdot {e}^{-{\displaystyle {\sum}_{s=1}^{t}{r}_{s}}}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{N}_{S}\left(T\right)\cdot \left(1-{\left(1-m\right)}^{\text{T}}\right)\cdot {P}^{\prime}\cdot {e}^{{\displaystyle {\sum}_{s=1}^{t}{\delta}_{s}}}\cdot {e}^{-{\displaystyle {\sum}_{s=1}^{t}{r}_{s}}}\end{array}$ (42)

where the ${\delta}_{s}$ values depend on rebalancing or not and on the investment strategy (see Section 3.1).

In the case of the $\left(n+1\right)$ dimensional VAR (1) model, using Equations (24)-(26), we get:

1) Rebalancing case:

$\begin{array}{c}I{F}_{0}=C\cdot P\cdot c+{\displaystyle \underset{t=1}{\overset{T}{\sum}}}\left({N}_{D}\left(t\right)\cdot \left(1-{\left(1-m\right)}^{t}\right)\cdot {P}^{\prime}\cdot {e}^{{\displaystyle {\sum}_{s=1}^{t}{\displaystyle {\sum}_{i=1}^{n}{\omega}_{i}\cdot {\delta}_{s}^{i}}}}\cdot {e}^{-{\displaystyle {\sum}_{s=1}^{t}{r}_{s}}}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{N}_{S}\left(T\right)\cdot \left(1-{\left(1-m\right)}^{\text{T}}\right)\cdot {P}^{\prime}\cdot {e}^{{\displaystyle {\sum}_{s=1}^{T}{\displaystyle {\sum}_{i=1}^{n}{\omega}_{i}\cdot {\delta}_{s}^{i}}}}\cdot {e}^{-{\displaystyle {\sum}_{s=1}^{T}{r}_{s}}}\end{array}$ (43)

2) No Rebalancing case:

$\begin{array}{c}I{F}_{0}=C\cdot P\cdot c+{\displaystyle \underset{t=1}{\overset{T}{\sum}}}\text{\hspace{0.17em}}{\displaystyle \underset{i=1}{\overset{n}{\sum}}}\left({N}_{D}\left(t\right)\cdot \left(1-{\left(1-m\right)}^{t}\right)\cdot {\omega}_{i}\cdot {P}^{\prime}\cdot {e}^{{\displaystyle {\sum}_{s=1}^{t}{\delta}_{s}^{i}}}\cdot {e}^{-{\displaystyle {\sum}_{s=1}^{t}{r}_{s}}}\right)\\ \text{\hspace{0.17em}}\text{\hspace{0.17em}}+{N}_{S}\left(T\right)\cdot \left(1-{\left(1-m\right)}^{\text{T}}\right)\cdot {\omega}_{i}\cdot {P}^{\prime}\cdot {e}^{{\displaystyle {\sum}_{s=1}^{T}{\delta}_{s}^{i}}}\cdot {e}^{-{\displaystyle {\sum}_{s=1}^{T}{r}_{s}}}.\end{array}$ (44)

In (43) and (44) ${\omega}_{i}$ , $i=1,2,\cdots ,n$ , are the proportions invested in each of the n assets and their value changes according to the chosen investment strategy.

Let us now consider the cash flows for each t, connected to the contract. They are given by the difference between the obligations arising from the guarantees (linked to the actual number of deaths each month (GMDB) and to the actual number of survivors at time T (GMMB)) and the inflows arising from the puts payoff. At time zero we get:

$C{F}_{0}\left(GMDB\right)={\displaystyle \underset{t=1}{\overset{T}{\sum}}}\left({N}_{D}\left(t\right)-{E}_{D}\left(t\right)\right)\cdot {\left({G}_{t}-{F}_{t},0\right)}_{+}\cdot {e}^{-{\displaystyle {\sum}_{s=1}^{t}{r}_{s}}}$ (45)

$C{F}_{0}\left(GMMB\right)=\left({N}_{S}\left(t\right)-{E}_{S}\left(T\right)\right)\cdot {\left({G}_{T}-{F}_{T}\mathrm{,0}\right)}_{+}\cdot {e}^{-{\displaystyle {\sum}_{s=1}^{T}{r}_{s}}}$ . (46)

Therefore, the portfolio value at inception is:

$V{P}_{0}=I{F}_{0}-P{P}_{0}-C{F}_{0}\left(GMDB\right)-C{F}_{0}\left(GMMB\right)$ (47)

where the value of each quantity depends on the specific model and investment strategy, as previously explained.

We are interested in capturing the distribution of the potential future markets and mortality rates. To this aim we resort to a Monte Carlo simulation, assessing $V{P}_{0}$ across a real world scenario.

5. Illustrations

In this section we study an illustrative portfolio of VA contracts offering GMDB and GMMB guarantees (see Section 2) issued to 1000 independent lives aged 55 and paying a single premium P = 6000 at issue.

Table 1 summarizes the portfolio characteristics.

The premium P, net of the initial charge (3% of P), is invested into the fund. The benefits are paid monthly and a monthly charge m = 0.0025 is deducted. The guaranteed amount is obtained by accumulating the net premium at the guaranteed rate g. In our example we consider three different guarantee rates: 0.015, 0.02 and 0.025. We assume that the VA writer offers two asset allocations to the policyholders. The first one, called moderate, is characterized by a large proportion of assets invested in bonds. The second one, mixed, is characterized by equal proportions in bonds and stocks. These asset allocations characterize some funds available in VA products sold in Europe (among others the “Accumulator” sold by AXA)

According to the moderate allocation (say Allocation A) we assume that the Company invests 70% in the Euro Big index, 18% in the Eurostoxx index and 12% in the S & P index. For the mixed allocation (say Allocation B) the fund allocation is: 50% in the Euro Big index, the 30% in the Eurostoxx index and the 20% in the S & P index.

Moreover we use two different accumulation functions to calculate the expected returns on investment in each fund, depending on whether the basket of assets constituting the fund itself is rebalanced or not.

The Euro Big index (Broad investment Euro BIG Bond index) provides a well known benchmark for Euro Bonds fixed-income portfolios. It covers all sectors of the investment grade fixed income market that are accessible to institutional investors and accurately measures their performance and characteristics. The Eurostoxx index is a broad, yet liquid, subset of the Stoxx Europe 600 Index. With a variable number of components, the index represents large, mid and small capitalization companies of 12 Eurozone countries: Austria, Belgium, Finland, France, Germany, Greece, Ireland, Italy, Luxembourg, the Netherlands, Portugal and Spain. The S & P500 is a free-capitalization-weighted index of the

Table 1. The illustrative portfolio.

prices of 500 large-cap common stocks actively traded in the United States.

The discounting rate is modelled using the Italian “Rendistato” that is a weighted average yield on a basket of Italian government securities.

For each index, we use the data from June 2002 to April 2016. The data sources are: the yield book of Citi fixed income indices for the EuroBIG index, Yahoo finance for the Eurostoxx and the S & P500 indices and the historical series published by the Bank of Italy, for the Italian Rendistato.

Figure 1 illustrates the behavior of the three considered indexes as well as the Italian Rendistato series.

For the mortality rates, we refer to a non parametric mortality table: IPS55

5.1. Estimation of the VAR Model Parameters

The parameter estimation results are obtained using the VAR package of the software R. All the parameters are estimated on a monthly basis.

Before estimating the parameters, we subtract the long term mean vector $\underset{\_}{\mu}$ from the data. We consider the sample mean of each series as an estimator for the long term mean vector $\underset{\_}{\mu}$ .

Let us first consider the case of the two-dimensional VAR (1) model described in Section 3.1. The long term mean of the discounting rate is 0.0028. Table 2 shows the long term mean of the composite rate of return on assets in each case: “rebalancing” for allocation A and B (column 1) and “not rebalancing” for allocation A and B (column 2).

Table 3 and Table 4 show the parameter estimation results for Allocation A. We have a process that is stable since the eigenvalues of the $\underset{\_}{\Phi}$ matrix are less than one. They are(0.9929, 0.0977) for the rebalancing case and (0.9929, 0.0957)

Figure 1. EuroBIG, Eurostoxx, S & P 500 and Italian Rendistato series. June 2002-April 2016.

Table 2. μ estimation. Two-dimensional VAR (1).

in case of no rebalancing.

Table 5 and Table 6 show the parameter estimation results for Allocation B. We have a process that is stationary since the eigenvalues of the $\underset{\_}{\Phi}$ matrix are less than one. In this case they are (0.9929, 0.1009) for the rebalancing case and (0.9929, 0.1139) in case of no rebalancing.

Let us now consider the $\left(n+1\right)$ -dimensional VAR (1) model: in the specific example we are considering n = 3, therefore we will have a four-dimensional VAR (1) model.

We indicate by ${\delta}_{t}^{1}$ the Eurobig index, by ${\delta}_{t}^{2}$ the Eurostoxx50 index, by ${\delta}_{t}^{3}$ the S & P index and by ${r}_{t}$ the Italian Rendistato (Figure 1).

Table 3. $\underset{\_}{\Phi}$ matrix estimation. Allocation A (moderate). Two-dimensional VAR (1) model for the continuous weighted average rates of return in rebalancing and not rebalancing case and Italian Rendistato series. June 2002-April 2016.

Table 4. ${\underset{\_}{\Sigma}}_{a}$ matrix estimation. Allocation A (moderate). Two-dimensional VAR (1) model for the continuous weighted average rates of return in rebalancing and not rebalancing case and Italian Rendistato series. June 2002-April 2016.

Table 5. $\underset{\_}{\Phi}$ matrix estimation. Allocation B (mixed). Two-dimensional VAR (1) model for the continuous weighted average rates of return in rebalancing and not rebalancing case and Italian Rendistato series. June 2002-April 2016.

Table 6. ${\underset{\_}{\Sigma}}_{a}$ matrix estimation. Allocation B (mixed). Two-dimensional VAR (1) model for the continuous weighted average rates of return in rebalancing and not rebalancing case and Italian Rendistato series. June 2002-April 2016.

The long term means on the three assets are: ${\mu}^{1}=0.0024$ , ${\mu}^{2}=0.0022$ , ${\mu}^{3}=0.0044$ . The estimation results for the $\underset{\_}{\Phi}$ matrix and the covariance matrix of the residuals are shown in Table 7 and Table 8.

The process is stationary since the absolute values of the eigenvalues of the $\underset{\_}{\Phi}$ matrix, which are 0.9919, 0.9513, 0.1005 and 0.09064, are less than one.

Nevertheless let us look at Table 7. Each row shows the regression parameters of the equations describing the evolution in time of each variable.

Looking at the first and second rows, first of all we observe a strong autoregressive dependence of the two bond indexes, ${\delta}_{t}^{2}$ and ${r}_{t}$ . Indeed the parameters explaining the evolution based on their own lag are high: 0.98287 and 0.95810, respectively. The regression parameters explaining the dependence on the other variables show a weak dependence on stock indexes.

The third and fourth rows of Table 7 give the estimated parameters of stock indexes processes. The regression parameters explaining the dependence on the lag of bond indexes are: −10.70956 and 9.01693 for ${\delta}_{t}^{2}$ equation and −9.95162 and 9.77489 for the ${r}_{t}$ equation.

All these observations tell us that there is not a clear dependence of bond indexes ( ${\delta}^{2}$ and ${r}_{t}$ ) on stock indexes ( ${\delta}^{2}$ and ${\delta}^{3}$ ).

Let us now focus on stock indexes. The two indexes show a weak autoregressive dependence: −0.038769 and 0.05082, respectively. In order to better capture their mutual dependence, we modelled the two indexes by a two-dimensional VAR (1) model. The estimated parameters are shown in Table 9.

The autoregressive dependence of each index is still weak. Moreover the eigenvalues are low: 0.1322 and 0.06405.

Based on the aforementioned observations, we conclude that a VAR (1) is questionable when interested in describing the behavior of stock indexes.

We decide to model the stock indexes by a Geometric Brownian Motion model (GBM), the most widely used model for stock prices.

As it is well known, the stock indexes ${S}_{t}$ follows a GBM if:

$\text{d}{S}_{t}=\mu {S}_{t}\text{d}t+\sigma {S}_{t}\text{d}{W}_{t}$ (48)

where ${W}_{t}$ is a Brownian motion and μ and σ are positive constants.

This implies that the instantaneous rate of return ${\delta}_{t}$ is a Brownian Motion:

${\delta}_{t}=\frac{\text{d}{S}_{t}}{{S}_{t}}=\mu +\sigma \text{d}{W}_{t}.$ (49)

Table 7. $\underset{\_}{\Phi}$ matrix estimation. Four-dimensional VAR (1) model. Eurobig, Eurostoxx 50, SP and Italian Rendistato series. June 2002-April 2016.

Table 8. ${\underset{\_}{\Sigma}}_{a}$ matrix estimation. Four-dimensional VAR (1) model. Eurobig, Eurostoxx50, SP and Italian Rendistato series. June 2002-April 2016.

Table 9. $\underset{\_}{\Phi}$ matrix estimation. Two-dimensional VAR (1) model. Eurostoxx 50 and S & P series. June 2002-April 2016.

The solution we propose in order to include the GBM in the model, is a modified version of the four-dimensional VAR (1) model.

First of all, we estimate again the parameters of a two-dimensional VAR (1) model considering only the bond indexes, Eurobig and Rendistato.

The $\underset{\_}{\Phi}$ matrix estimation results are given in Table 10.

Again the process is stationary because the eigenvalues of the $\underset{\_}{\Phi}$ matrix are 0.992 and 0.9512.

Let us now come back to the four-dimensional VAR (1) model. We replace the estimated $\underset{\_}{\Phi}$ matrix in Table 7, with the ${\underset{\_}{\Phi}}^{\mathrm{*}}$ matrix of Table 11.

The ${\underset{\_}{\Phi}}^{\mathrm{*}}$ matrix is obtained by setting all the parameters connected to the Eurostoxx and the S & P equal to zero and by setting the other parameters equal to the ones given in Table 10.

With reference to (11) the residuals matrix is obtained as follows:

${\underset{\_}{{a}_{t}}}^{*}=\left(\underset{\_}{{X}_{t}}-\underset{\_}{\mu}\right)-{\underset{\_}{\Phi}}^{*}\cdot \left(\underset{\_}{{X}_{t-1}}-\underset{\_}{\mu}\right)$ (50)

where $\underset{\_}{{X}_{t}}$ is given by (21).

It is than possible to estimate the covariance matrix of the residuals ${\underset{\_}{\Sigma}}_{a}^{\mathrm{*}}$ (see Table 11).

In order to compare the VAR (1) model and the proposed modified version, we ran 10,000 simulations for values of ${\delta}_{t}^{1}$ , ${\delta}_{t}^{2}$ , ${\delta}_{t}^{3}$ and ${r}_{t}$ values with $t\in \left[1,2,\cdots ,120\right]$ assuming that their evolution in time is described by the ‘full’ four-dimensional VAR (1) model (Table 7 and Table 8) and by the modified version of the same model (Table 11 and Table 12).

Figure 2 and Figure 3 show the evolution over time of the expected values of the four variables in both cases and for different selected initial values. For each variable and for each t, the expected values are computed as the mean of the 10,000 simulated values.

We observe that the two processes have the same long term dynamics with different behavior only for a few months.

Table 10. $\underset{\_}{\Phi}$ matrix estimation. Two-dimensional VAR (1) model. Eurobig and Italian Rendistato series. June 2002-April 2016.

Table 11. $\underset{\_}{\Phi}$ matrix estimation. Four-dimensional VAR (1) model. Eurobig and Italian Rendistato series. June 2002-April 2016.

Table 12. ${\underset{\_}{\Sigma}}_{a}^{\mathrm{*}}$ matrix estimation. Modified four-dimensional VAR (1) model. Eurobig, Eurostoxx50, SP and Italian Rendistato series. June 2002-April 2016.

Figure 2. Four-dimensional VAR (1) model. Expected values with different initial values. Number of simulations = 10,000.

Figure 3. Modified four-dimensional VAR (1) model. Expected values with different initial values. Number of simulations = 10,000.

5.2. Some Results

In this section, we refer to the illustrative portfolio described in Section 5.1.

We are interested in the expected value and the conditional value at risk at a chosen confidence level of 95% of the portfolio at inception, for a static hedging strategy (Section 4.2). The software used for the following analysis is Matlab.

We resort to a Monte Carlo simulation, assessing the quantities of interest across a real-world scenario set for financial and mortality variables.

Here are the steps of the simulation procedure:

1) Randomly generate the monthly rates of return and the monthly discounting rates using the results of the estimated parameters;

2) Randomly generate the monthly death probabilities. For fractional ages a Uniform distribution of deaths assumption is made: for example the probability of dying in any month of a given year is the same.

3) Assess $V{P}_{0}$ value (Equation (47));

4) Repeat the procedure N = 10,000 times;

5) Estimate $E\left[V{P}_{0}\right]$ by taking the sample mean of the N simulated values of $V{P}_{0}$ ;

6) Estimate the conditional value at risk at a 95% confidence level from the empirical distribution of $V{P}_{0}$ . As is well known, the $CVaR\left(1-\alpha \right)\mathrm{\%}$ is the average of the $\left(1-\alpha \right)\mathrm{\%}$ worst scenarios. In our example we choose $\alpha =95\%$ and the CVaR 95% is the average of the 5% lowest values of $V{P}_{0}$ .

The simulation procedure is repeated for each assumption about the financial variables previously described, the two and four-dimensional VAR (1) models, both with rebalancing of the basket of assets and without rebalancing.

We then compare the results with the choice of no hedge, that is the case in which the insurer decides not to buy the hedge portfolio of put options.

Table 13 shows the cost of the hedging portfolio in each case. The portfolio price has been determined in a risk neutral framework on the basis of the results obtained in Sections 4.1 and 5.1.

First of all we observe that in each case the cost is lower than the intitial charge paid by the insurers, that is $C\ast P\ast c$ , that in our illustrative portfolio amounts to 180,000 euros.

As the guaranteed rate gets higher, the put options price increase and the cost of the entire portfolio increases, too.

In general, the cost of the hedging portfolio is higher for the rebalancing strategy than the no rebalancing one, as expected. Indeed, the estimated volatilities of the underlying are higher in the rebalancing case, for each considered option.

The same happens going from Allocation A to Allocation B. The second strategy is more risky because it is characterized by a higher proportion of assets invested in stocks and therefore a higher volatility.

Tables 14-17 compare the results for static hedging and those for when the Company decides not to hedge at all.

First of all, we observe that in each considered case the expected value of the portfolio $E\left[V{P}_{0}\right]$ is higher if we shift from no hedging to static hedging. This is expected because the no hedging case does not “adapt” to mortality experience over time, as static hedging does to some extent.

Moreover we observe that Allocation B gives lower values of $E\left[V{P}_{0}\right]$ and that

Table 13. Cost of the hedging Portfolio.

Table 14. Simulation results for $V{P}_{0}$ : 2DVAR (1), Rebalancing case, $N=10,000$ , $\alpha =5\%.$

Table 15. Simulation results for $V{P}_{0}$ : 2DVAR (1), No Rebalancing case, $N=10,000$ , $\alpha =5\%.$

Table 16. Simulation results for $V{P}_{0}$ : 4DVAR (1), Rebalancing case, $N=10,000$ , $\alpha =5\%.$

Table 17. Simulation results for $V{P}_{0}$ : 4DVAR (1), No Rebalancing case, $N=10,000$ , $\alpha =5\%.$

holds for each considered case. Indeed Allocation B is more risky, therefore we get higher expected liabilities.

In particular, in the no hedge case, Allocation B always leads to a loss. Regarding Allocation A, even if we get some positive values, they are lower than the initial sum of 180,000 euros.

Looking at static hedging, $E\left[V{P}_{0}\right]$ is always lower if allocation B is chosen, but we never get negative values. This worsening is due to the aforementioned higher cost of the hedging portfolio (Table 13).

Moving from the case of rebalancing the basket of assets to not rebalancing it, the results improve in each case. Referring to static hedging strategy, this is a consequence of the higher cost of the hedging portfolio (Table 13) due to a higher level of the underlying volatility in the rebalancing case. Moreover the higher volatility of the underlying allows for higher extreme returns and consequently higher liabilities. This effect is observable for the no hedging case as well.

Finally let us look at the results obtained by means of two-dimensional VAR (1) and four-dimensional VAR (1) models (2D VAR (1) and 4D VAR (1), respectively). For the static hedging case, 2D VAR (1) gives slightly higher values than 4D VAR both in case of rebalancing or not. Even if the cost of the hedging portfolio is lower in the 4D VAR (1) case, we can observe that in 2D the asset “volatilies” are averaged out in the estimators but in 4D they are estimated separately. This allows for more extreme returns in the 4D model and consequently higher liabilities. This effect is stronger when not hedging.

Looking at Conditional Value at Risk values, we observe the same behaviour as the expected values, of course. In particular, in the no hedge case the difference between the conditional values at risk and the expected values is very high. This is due to the higher risk associated with this choice. Moreover, for both the no hedge and the static hedge cases, we observe a worsening in terms of conditional values at risk, when changing from allocation A to the more risky allocation B.

6. Concluding Remarks

This paper studies, from the insurer point of view, a VA product offering GMDB and GMMB guarantees. In particular, the financial position of an insurer issuing a homogeneous portfolio of VAs offering those guarantees, is assessed. To this aim, the expected value and the conditional value at risk of this portfolio are quantified in a real world scenario. The results under a static hedging strategy are compared with those obtained when not hedging. Static hedging, as expected, gives better results.

Both in static hedging and no hedging cases, the portfolio value depends on many factors and we investigate some of them. Regarding the choices offered to the policyholders, we look at two different asset allocations and, as expected, the portfolio value decreases as the asset allocation becomes more risky.

A crucial aspect is the choice of a model for the asset returns and the discounting rate. This choice impacts the valuation of the portfolio for risk management purposes.

VA guarantees are linked to more that one asset and the use of single-asset (univariate) models for VA pricing and risk management purposes, may lead to results that are not reliable. We show the results obtained by modelling jointly the financial variables involved by means of two kinds of VAR (1) models. We observe that using a four-dimensional VAR (1) model gives lower values of the portfolio than a two-dimensional VAR (1) model. In the first case, the asset returns on each fund are combined in a single time series while in the four-dimensional VAR (1) model, the asset returns on each fund are modelled separately. This structure adds more information in modelling the financial variables and the results could be more reliable. This provides useful information to the insurer concerned about risk management.

We also consider two different accumulations functions depending on whether the basket of assets is rebalanced or not. The rebalancing case gives lower values of the portfolio.

Following the analysis perfomed in this paper, further work could be done. This may include, for instance, the implementation of a dynamic hedging strategy and the use of stochastic mortality models.

Cite this paper

Orlando, A. and Parker, G. (2018) Valuation and Risk Assessment of a Portfolio of Variable Annuities: A Vector Autoregression Approach.*Journal of Mathematical Finance*, **8**, 349-371. doi: 10.4236/jmf.2018.82023.

Orlando, A. and Parker, G. (2018) Valuation and Risk Assessment of a Portfolio of Variable Annuities: A Vector Autoregression Approach.

References

[1] Ledlie, M.C., Corry, D.P., Finkelstein, G.S., Ritchie, A.J., Su, K. and Wilson, D.C.E. (2008) Variable Annuities. British Actuarial Journal, 14, 327-389.

https://doi.org/10.1017/S1357321700001744

[2] Kalberer, T. and Ravindran, K. (2009) Variable Annuities. A Global Perspective. Risk Books.

[3] Bauer, D., Kling, A. and Russ, J. (2008) A Universal Pricing Framework for Guaranteed Minimum Benefits in Variable Annuities. Astin Bulletin, 38, 621-651.

https://doi.org/10.1017/S0515036100015312

[4] Bacinello, A.R., Millossovich, P., Olivieri, A. and Pitacco, E. (2010) Variable Annuities: Risk Identification and Risk Assessment. CAREFIN Research Paper 14/2010, 1-48.

https://ssrn.com/abstract=1799266

[5] Feng, R. and Volkmer, H.W. (2012) Analytical Calculation of Risk Measures for Variable Annuity Guaranteed Benefits. Insurance: Mathematics and Economics, 51, 636-648.

https://doi.org/10.1016/j.insmatheco.2012.09.007

[6] Li, J.S.-H. and Ng, A.C.-Y. (2013) Pricing and Hedging Variable annuity guarantees with multiasset stochastic investment models. North American Actuarial Journal, 17, 41-62.

https://doi.org/10.1080/10920277.2013.773240

[7] Zivot, E. and Wang, J. (2006) Vector Autoregressive Models for Multivariate Time Series. In: Modelling Financial Time Series with S-PLUS, Springer Science, 385-429.

[8] Hardy, M. (2003) Investment Guarantees: Modeling and Risk Management for Equity Linked Life Insurance. John Wiley & Sons, Hoboken, NJ.

[9] Milevsky, M.A. and Posner, S.E. (1998) A Closed-Form Approximation for Valuing Basket Options. The journal of Derivatives, 5, 54-61.

https://doi.org/10.3905/jod.1998.408005

[1] Ledlie, M.C., Corry, D.P., Finkelstein, G.S., Ritchie, A.J., Su, K. and Wilson, D.C.E. (2008) Variable Annuities. British Actuarial Journal, 14, 327-389.

https://doi.org/10.1017/S1357321700001744

[2] Kalberer, T. and Ravindran, K. (2009) Variable Annuities. A Global Perspective. Risk Books.

[3] Bauer, D., Kling, A. and Russ, J. (2008) A Universal Pricing Framework for Guaranteed Minimum Benefits in Variable Annuities. Astin Bulletin, 38, 621-651.

https://doi.org/10.1017/S0515036100015312

[4] Bacinello, A.R., Millossovich, P., Olivieri, A. and Pitacco, E. (2010) Variable Annuities: Risk Identification and Risk Assessment. CAREFIN Research Paper 14/2010, 1-48.

https://ssrn.com/abstract=1799266

[5] Feng, R. and Volkmer, H.W. (2012) Analytical Calculation of Risk Measures for Variable Annuity Guaranteed Benefits. Insurance: Mathematics and Economics, 51, 636-648.

https://doi.org/10.1016/j.insmatheco.2012.09.007

[6] Li, J.S.-H. and Ng, A.C.-Y. (2013) Pricing and Hedging Variable annuity guarantees with multiasset stochastic investment models. North American Actuarial Journal, 17, 41-62.

https://doi.org/10.1080/10920277.2013.773240

[7] Zivot, E. and Wang, J. (2006) Vector Autoregressive Models for Multivariate Time Series. In: Modelling Financial Time Series with S-PLUS, Springer Science, 385-429.

[8] Hardy, M. (2003) Investment Guarantees: Modeling and Risk Management for Equity Linked Life Insurance. John Wiley & Sons, Hoboken, NJ.

[9] Milevsky, M.A. and Posner, S.E. (1998) A Closed-Form Approximation for Valuing Basket Options. The journal of Derivatives, 5, 54-61.

https://doi.org/10.3905/jod.1998.408005