Deflation Risk and Implications for Life Insurers

Life insurers are exposed to deflation risk: falling prices could lead to insufficient investment returns, and inflation-indexed protections could make insurers vulnerable to deflation. In this spirit, this paper proposes a market-based methodology for measuring deflation risk based on a discrete framework: the latter accounts for the real interest rate, the inflation index level, its conditional variance, and the expected inflation rate. US inflation data are then used to estimate the model and show the importance of deflation risk. Specifically, the distribution of a fictitious life insurer’s future payments is investigated. We find that the proposed inflation model yields higher risk measures than the ones obtained using competing models, stressing the need for dynamic and market-consistent inflation modelling in the life insurance industry.


Introduction
Fluctuation in the general price level of goods and services is reflected through the inflation rate.Despite rather small changes in the inflation rate over the past two decades, the last recession created the potential for price instability.
Before the Great Recession, inflation risk was typically associated with high inflation rates.Events such as the Great Inflation of the 1970s are notable for steep inflation rates and could explain, to some extent, the common misconception that links inflation risk to high inflation rates.Yet, lately, inflation risk has become much more symmetric, and a fear that inflation may fall too low has emerged (so-called deflation).For instance, some recent actions undertaken by the Federal Reserve (Fed) in the US have been motivated by this fear of deflation, and an important number of European countries have had deflation issues in the past few years.Deflation is one of the most challenging risks faced by central banks nowadays.
In the long-run, insurance contracts, reserves, investment policies, and pension liabilities are sensitive to unexpected inflation. 1Long-term liabilities could be much worse if the inflation rate suddenly changes (e.g., see Bohnert et al. [1]).In addition, asset returns are tied up to the inflation, meaning that portfolio income could be quite dependent on the prevailing inflation rate.
Life insurers have an uncommonly large vulnerability to deflation risk: falling prices could lead to insufficient investment returns to support the rate credited to policyholders.Deflation is indeed "the worst possible world for a life insurance company" (see Mittelstaedt [2]).Japan is a notable example of how deflation can impact the life insurance business.The historical Japanese inflation rate remained below zero for most of the last twenty years, and this episode of deflation provoked the bankruptcy of at least eight Japanese insurers since 1996.
Life insurers are also vulnerable to embedded optionalities in the design of their products.For instance, inflation-indexed annuities often provide additional protection that increases the paid amount in the case of positive inflation, but does not reduce the amount in the case of deflation.This added protection is similar to an inflation floor, as the annuitant is covered against decreases in the general price level of goods and services.
On a similar note, inflation-indexed defined-benefit pension plans also provide the same kind of protection to their participants: in many instances, a decrease of the Consumer Price Index (CPI) would not reduce the retirement benefit amount. 2 As an increasing number of plan sponsors decide to de-risk their financial exposures by buying special annuities, they transfer their interest rate, longevity, and deflation risks to life insurance companies.These life insurers issue pension buy-in and buy-out annuities, similar to portfolios of inflation-indexed annuities, thus exposing themselves to additional deflation risk.
Oddly, inflation is only indirectly measured by insurance companies in general [3], and Solvency II does not directly account for inflation risk.As deflation risk could have a significant impact on life insurers' portfolios, it should be incorporated into the insurers' enterprise risk management (ERM) framework consistently.Market-consistency is important, as an ERM framework must be able to capture the changing economic conditions.Further, new solvency regulations force insurance companies to assess their liabilities using "mark-to-market" approaches.As a matter of fact, market-consistency requires the use of models with sufficient degrees of freedom to integrate all relevant market data [4].In this spirit, this paper presents a market-based methodology for measuring deflation risk based on a discrete framework inspired by Fisher's [5] equation.The model accounts for the real interest rate, the inflation index level, its conditional variance, and the expected inflation rate.The real interest rate and the expected inflation rate both follow a first-order autoregressive (AR) model.The conditional variance follows the generalized autoregressive conditional heteroskedasticity (GARCH) model of Heston and Nandi [6].This approach is somewhat different from common practice, as we model both the real interest rate and the expected inflation separately. 3The dynamics used in this paper are in line with the inflation modelling literature.Fleckenstein et al. [7] use a similar model to assess deflation risk, although their continuous-time model does not allow for heteroskedasticity in the index price level.A regime-switching first-order AR model for the expected inflation rate is proposed by Ahlgrim and D'Arcy [8].The authors' model is a generalization of Wilkie's [9] inflation model, as it allows for different parameters, depending on the prevailing regime.Singor et al. [10] consider a Heston [11] type framework to model the inflation index, thus considering a dynamic variance process as is the case in this study.Ang et al. [12] consider, among other things, autoregressive moving-average (ARMA) models and Phillips [13] curve to model the inflation rate. 4s is commonly the case with tail risks, deflation is hard to quantify using classic econometric methods.For instance, Ang et al. [12] use a large number of models based on time series of inflation to assess inflation forecasting, and conclude that these models perform weakly in forecasting the first moment of inflation.In addition, even though they conclude that survey data perform better, these provide little or no evidence of potential deflation risk.
In the light of this issue, assessing deflation risk would require more than just past inflation rate data.Fortunately, a considerable amount of inflation-dependent securities are available, and their 2 For instance, US Social Security benefits do not decrease if the cost-of-living adjustment is negative, as prescribed by the Social Security Act.It is the same in Canada: the Canada Pension Plan payments do not decrease in case of deflation.Note that other private and public pension plans may provide similar protections.3 Generally, the nominal interest rate is modelled directly, without special emphasis on the dynamics of underlying economic quantities.Modelling the real interest rate and the expected inflation rate separately allows for more flexibility, and has the advantage of capturing the risk associated specifically to these two factors.
prices shall inform us on the likelihood of deflation over a given time horizon.Inflation-indexed bonds are securities that protect the investors' purchasing-power by adjusting the principal based on changes in an inflation index.Most of the developed countries issue these bonds, yet US Treasury Inflation-Protected Securities (TIPS) are the most liquid ones, with about $500 billion in issuance.Inflation swaps are also available for most developed countries.In the US, they were introduced at the same time as TIPS, in 1997.The notional size of the US inflation swap market is estimated to be hundreds of billions [14].A parallel market for inflation option exists; it started in 2002 with the introduction of caps and floors on the realized inflation rate.According to Fleckenstein et al. [7], the inflation option market is sufficiently liquid: active quotations are available since 2009 for a broad range of maturities and strikes.
To estimate our model, multiple types of data source are used, since inflation index levels would only capture the average expected inflation rate and not the tails of the index level distribution.
Characterizing the entire index level distribution is of paramount importance in this study, because deflation risk is related to the left tail of this distribution.Thus, in addition to inflation index levels, we include nominal risk-free bonds, inflation swaps, and inflation options to the model estimation, as these securities contain relevant information about the tails of the distribution.The four types of data used in this paper are similar to what other inflation studies have utilized.For instance, Ang et al. [12] and Ahlgrim and D'Arcy [8] use time series of realized inflation rates, which is technically the same as using the index level. 5Inflation option market prices were used by Singor et al. [10], and Kitsul and Wright [15].Fleckenstein et al. [7] employed swap prices and option data, as in this study.
An estimation method based on the unscented Kalman filter (UKF) is considered.All physical and risk-neutral parameters are estimated using likelihood maximization in a unique stage.The proposed method contrasts with other methods used in the literature (i.e., calibration techniques and multi-stage estimations).On one hand, calibration methods such as the one used by Singor et al. [10] select the model parameters based on a given day's available prices.Even though this methodology can give a set of parameters that is consistent with market prices, it is impossible to know whether these are robust in time.In addition, only risk-neutral parameters can be inferred, since physical parameters are not directly involved in the pricing of derivatives.On the other hand, multi-stage methods, such as the one presented by Fleckenstein et al. [7], can be inadequate to capture the interaction among various model parameters, and also lack formal proofs of their effectiveness.
Another interesting consequence of using the UKF-based approach is that the expected inflation rate is forward-looking, as it is based on the current market prices.The model adapts readily to new market conditions, and the expected inflation rate is then consistent with investors' expectations.
To estimate the model, US inflation data are used.The unconditional annualized long-term trend of expected inflation is about 1.8%, which is consistent with the Fed desired target range for inflation-between 1.7% and 2.0%.The evolution of the expected inflation is dependent on the current economic conditions: for instance, during the pre-recession era, the average expected inflation was about 3.4%, whereas its average is about 0.8% in the post-recession period.Deflation risk changes over time, with average 1-year deflation probabilities of 4%, 18%, and 29% for the pre-recession, recession, and post-recession eras, respectively.
The effects of deflation risk are assessed in the life insurance industry.In this spirit, the distribution of a fictitious life insurer's future payments in present value is investigated.On average, the proposed inflation model yields risk measures that are higher than the ones obtained with competing models.The right tail of the future payments distribution is fatter, as low inflation scenarios are taken into consideration, and this is according to the current market conditions.The importance of modelling adequately embedded optionalities is also shown; an inflation-indexed annuity for which the paid amounts do not decrease when the inflation index level declines is used.The results show that 5 Time series of realized inflation rates is obtained by a simple transformation of inflation index levels.deflation could have a notable impact on life insurers when inflation protections are embedded in the product design.
The contributions of this paper are manifold.First, a modelling framework is constructed to capture the desired empirical facts, such as correlation between the real interest rate and the expected inflation, as well as conditional heteroskedasticity.This model allows for derivative prices in closed or semi-closed form solutions.Second, forward-looking and market-consistent assessment regarding deflation can be done using the framework.Option data are useful in capturing the tails of the expected inflation distribution, and the filter-based methodology allows for dynamic risk assessments that are consistent with market participants' expectations.This is of utmost importance, as new solvency regulations require insurers to assess their risk in a market-consistent way.Finally, our last contribution is to apply various models-full and nested frameworks-to assess, in the context of a life insurer, the importance of the two core components of the nominal interest rate: the real interest rate and the expected inflation.Based on two different exercises, we show that the full framework allows for more flexibility and a better understanding of deflation risk.
The remainder of the paper is organized as follows.Section 2 introduces the inflation framework.The estimation procedure is explained in Section 3. In Section 4, the estimation procedure is applied to US inflation data.Section 5 shows examples of how deflation risk can affect the life insurance industry.Finally, Section 6 concludes.

Framework
A discrete-time model for both inflation and real interest rate dynamics is presented in this section.Physical and risk-neutral dynamics are derived.The latter are of utmost importance for valuing derivatives such as bonds, inflation swaps, and inflation options.

Inflation and Real Interest Rate Dynamics
The Fisher [5] equation links nominal and real interest rates.Basically, the difference between the two interest rates is given by the expected inflation.Let r t be the nominal risk-free interest rate, defined as: where R t is the real risk-free interest rate at time t, and µ t is the expected inflation at time t.
Under the physical measure P, the real interest rate dynamics are given by an autoregressive model: where {ε R,t } ∞ t=0 is a sequence of standardized Gaussian random variables.This model is equivalent to a discretized version of Vasicek [16].Parameters θ and κ represent the unconditional level of real interest rate and the speed of mean reversion, respectively.The volatility of the real interest rate innovations is given by parameter σ.
There are many possible model formulations for modelling inflation, as explained in the introduction of this study. 6Even though the inflation index does not require an elaborate model per se, trying to fit various sources of information under a simple framework would be futile, as it cannot reflect the complexity of the inflation's distribution and its evolution in time.Therefore, a flexible model is proposed.The latter allows for heteroskedastic innovations in the inflation index level, as well as stochastic expected inflation rate.Since time series and prices are observed on Inflation modelling has been studied by Wilkie [9], Ang et al. [12], Ahlgrim and D'Arcy [8], Singor et al. [10], Fleckenstein et al. [7], and Bohnert et al. [1], to name a few. a discrete basis, the proposed framework is a discrete-time model.This would also come in handy in the estimation procedure explained in the subsequent section.
Let I t denote the inflation index level at time t.Under the physical measure P, the log-inflation index level dynamics are given by where h t is the conditional variance at time t.   The proposed framework could be interpreted as a generalization of the well-known Black and Scholes [17] dynamics: if µ t and h t were constant, the inflation index level would be lognormally distributed.As a first generalization, the expected inflation rate µ t is stochastic and follows a first-order AR process, meaning that the expected inflation rate depends linearly on its previous value.It is also mean-reverting and reverts to ν, the long-run trend in expected inflation.Parameter β controls for the speed of the mean reversion.This autoregressive expected inflation rate is similar to Wilkie's [9] inflation model.
The conditional variance h t depends on the time t innovation ε I,t .It is the GARCH process of Heston and Nandi [6].This type of GARCH is popular, since it allows semi-closed form solutions for option pricing.GARCH models were also used by Kilian and Manganelli [18] in the context of inflation modelling.
The proposed framework is related to the one of Fleckenstein et al. [7].However, their model does not consider heteroskedasticity in the inflation index level, and is a continuous-time model. 8On the other hand, they include a stochastic long-run trend of expected inflation. 9Bansal and Yaron [19] also propose a similar affine specification to model long-run risk.
Table 1 summarizes the various parameters of the model.The interpretation of the parameters is also available in the table.

Risk-Neutral Dynamics
To value inflation-dependent securities, it is convenient to rewrite the inflation model under some pricing measure.The market model is incomplete, which implies that there is an infinite number of pricing measures.In this paper, the Girsanov theorem is applied; the details of the measure change are given in Appendix A. 7 Over a short period of time, it is known from a first-order Taylor expansion that E P I t −I t−1 In an earlier version of their paper, Fleckenstein et al. [7] do have stochastic volatility.However, the variance could be negative in the model, which is inconsistent with the definition of this quantity.
Under the risk-neutral measure Q, the dynamics of R t -the real interest rate at time t-are given by where θ * = θ − ση R /(1 − κ), η R is the change of measure parameter associated to ε R,t , and {ε * R,t } ∞ t=0 is a sequence of standardized Gaussian random variables under the risk-neutral measure.
The dynamics of the log-inflation index level log(I t ) are given by where

Derivative Valuation
The proposed model-in addition to being flexible-allows for closed and semi-closed form solutions for bonds and inflation-dependent derivatives.These products are used to estimate the model, and should therefore be priced in a timely fashion.
To price a zero-coupon bond, the nominal risk-free interest rate r t is used.Its price can be expressed in a closed-form expression as In this study, we also use zero-coupon inflation swaps.These swaps are the most widely used type of inflation swap.It is executed between two counterparties at time zero; at expiry, a single cash flow of (I t+T /I t − F) is exchanged, where T is the maturity and F ≡ F(R t , µ t , T) is the inflation swap price.The price is set at initiation of the contract, and it is equal to F = (1 + f ) T , where f is the inflation swap rate. 10Note that the timing and index lag construction of the inflation index are chosen to precisely match the definitions applied to TIPS issues.The inflation swap price can be expressed in a closed-form expression as is a scalar constant.Again, one can compute it by recursion.Appendix B.2 shows the details on how to obtain the swap pricing formula.
As with inflation swaps, zero-coupon inflation options are used in this study.They are the most-actively traded inflation options (along with year-on-year inflation options). 11et C(R t , µ t , h t+1 , T) denote the value of a European inflation cap with strike K and maturity of T. 12The payoff of such an option at its expiration date is similar to the one of a European call option: max 0, I t+T /I t − (1 + K) T .The value of the cap option at time t is known in a semi-closed form solution.It is given by where and k = T log(1 + K).The integral equations are evaluated numerically using the trapezoidal rule. where are four scalar constants.Again, one can compute these by recursion, starting with s = t + T where The proof of this pricing formula is given in Appendix B.3.

Estimation Method
In this study, the real interest rate, the expected inflation rate, and the conditional variance of the inflation index level are latent variables. 13It is thus challenging to infer these unobservable variables.Moreover, in presence of latent variables, estimation can be somewhat more complicated.Filtering techniques are useful in these situations. 14 It is possible to link the observable quantities to the latent ones through a state space representation.The latter is useful to filter the real interest rate, the expected inflation rate and the conditional variance of the inflation index level by associating them to noisy security prices and the inflation index level.It is then possible to recover state variables' estimates based on investors' current expectations.Note that security prices contain forward-looking information, which allows for a better understanding of the three latent variables.Further, both P-and Q-parameters are estimated in a one-stage procedure.
The state propagation equations are given by Equations ( 2), (4), and ( 5).These equations explain how the latent variables evolve in time.The measurement equations show how the unobservable variables are linked to observable quantities.The first measurement equation is related to the inflation index level.It is given in Equation (3): as I t is observed.Then, the market security prices are incorporated in the measurements to capture the latent variables' true nature and the wedge between P-and Q-parameters.
The bond measurement equation is log where B (Market) t (T) is the market T-year zero-coupon bond price at time t, and υ B,t,T is a centred Gaussian random variable of standard deviation ζ B .To maintain positive prices, a logarithmic transformation is applied.Note that this measurement equation is linear in R t and µ t .
The swap measurement equation is log where F (Market) t (T) is the market T-year inflation swap price at time t and υ F,t,T is a centred Gaussian random variable of standard deviation ζ F .This measurement equation is also linear in R t and µ t .
Regarding option data, we would like to minimize the relative implied volatility root mean square error in the spirit of Christoffersen et al. [21] and Ornthanalai [22].However, since the real interest rate, the expected inflation rate, and the conditional variance are not predictable variables, we cannot use the method directly.Instead, we model the relative implied volatility error as a Gaussian random variable.The option measurement equation is therefore given by where σ (Market) t (K, T) is the market Black and Scholes implied volatility of a time t cap option of strike price K and maturity T. IV (C, K, T) is the model Black and Scholes implied volatility for a time t cap of price C, strike K, and maturity T.Moreover, υ IV,K,t,T is a centred Gaussian random variable of standard deviation σ 14 For a review of recent developments in filtering methods for economics and finance applications, see Creal [20].
and δ IV , δ T ,and δ K are three parameters to be estimated. 15Homoskedastic normally distributed errors could be restricting.Therefore, we allow for heteroskedasticity in option errors by letting the standard deviation be a function of the data; i.e., the absolute value of the strike K and the maturity of the option T. This formulation is similar to the standard deviation of the option errors proposed in Bardgett et al. [23].Note that for inflation floor, the rationale is the same: we simply replace cap prices by floor prices.
To estimate the model, we follow a simple filtering approach.The unscented Kalman filter of Julier and Uhlmann [24] is applied, since some measurement equations are nonlinear. 16The UKF handles the non-linearity and approximates the posterior state density using a set of deterministically chosen sample points.These points capture the mean and covariance of the Gaussian state variable, and when propagated through the measurement equations, it captures the posterior mean and covariance of the observations accurately up to the second order.According to Christoffersen et al. [26], the UKF may prove to be a good approach for a number of problems in fixed income pricing with nonlinear relationships between the state vector and the observations.In addition, many researchers used the UKF to infer the latent states in financial applications. 17 Note that the filter allows quasi-likelihood function computation in a somewhat direct manner. 18Hence, quasi-maximum likelihood estimation is done by numerical optimization.Overall, the parameters to be estimated are δ IV , δ T , and δ K .A single set of parameters is obtained by this estimation method.

Application to the US Inflation
The four major sources of information used in this study are introduced: inflation index levels, nominal risk-free bonds, inflation swaps, and inflation caps and floors.Then, the estimation results are shown for the full model and the nested specifications.Finally, the deflation probabilities are computed.

Data
We collect the Consumer Price Index for All Urban Consumers (CPI-U) levels on the Federal Reserve Bank of St. Louis website via Federal Reserve Economic Data (FRED) from April 2004 to March 2016. 19The CPI-U is computed using monthly average prices for each category-area combination.It is published on a monthly basis by the Bureau of Labour Statistics (BLS).
The inflation index level is constructed from the monthly CPI-U series.The underlying value used in the derivative pricing is a three-month lagged linearly interpolated value of the CPI-U.Therefore, the inflation index level used in the model respects the features of the inflation derivative market, and is defined as the three-month lagged linearly interpolated value of the CPI-U.We use Wednesday CPI-U levels.For more details on the different features of this market, see Kerkhof [29].
The top panel of Figure 1 shows the behaviour of the inflation index time series.It is increasing on average.The average increase is about 2% (annualized) over the 2004-2015 period, which is the target inflation rate set by the Fed.Note that the index decreased during the last financial crisis, going from 219.9 on 2 April 2008 to 212.1 on 29 October 2008. 15The standard deviation of the error is proportional to the size of the observed implied volatility.This is consistent with the assumption that relative error on implied volatility should be approximately Gaussian. 16The parameters of the UKF technique have been assumed to be κ UKF = 0.001, α UKF = 0, and β UKF = 2, as proposed in Wan and Van Der Merwe [25]. 17For an example using market data, see Carr and Wu [27] or Boudreault et al. [28]. 18Quasi-likelihood means here that the first two moments of the posterior distribution have a second order precision, and a posterior Gaussian distribution has to be assumed. 19 where F is the swap price and T is the maturity in years.The inflation swap rates are collected via the Bloomberg system from July 2004 to December 2015.The grey-shaded region highlights the National Bureau of Economic Research (NBER)-dated recession.
The nominal risk-free bond prices are computed from both London Interbank Offered Rate (LIBOR) rates (for maturities less than 1 year) and swap rates (for maturities greater or equal to 1 year).The LIBOR and interest swap rates are collected via FRED from July 2004 to December 2015. 20 Maturities of 1, 2, 3, and 6 months, as well as 1, 2, 3, 4, 5, 7, 10 and 30 years are used.
The middle panel of Figure 1 shows the 1-, 10-, and 30-year nominal risk-free interest rates.Interestingly, the term structure of the risk-free interest rates seems to invert between 2006 and 2008.Panel A of Table 2 shows descriptive statistics for risk-free interest rates.As the maturity increases, the dispersion of the rates decreases, as long-term rates are less uncertain than short-term rates.Additionally, the term structure of risk-free interest rates is upward sloping on average, because the average rates increase as a function of the maturity.Inflation swap rates are expressed in percentage and maturities in years.Zero-coupon inflation swap rates are extracted from the Bloomberg system from July 2004 to December 2015.The risk-free interest rate is proxied by both LIBOR rates (for maturities less than 1 year) and swap rates (for maturities more than or equal to 1 year).The LIBOR and interest swap rates are collected via FRED from July 2004 to December 2015.The mean, standard deviation (SD), the minimum (Min), the first quartile (Q 1 ), the median, the third quartile (Q 3 ), the maximum (Max), and the number of observations (N) are computed across time for each tenor.
According to Fleming and Sporn [30], the inflation swap market appears transparent and somewhat liquid, despite the market's modest activity and its over-the-counter nature.We again use Wednesday data. 22A total of 8895 observations are included in the final swap prices dataset.
The bottom panel of Figure 1 exhibits the evolution of the zero-coupon inflation swap rates for maturities of 1, 10, and 30 years.The inflation swap rates with shorter maturities are quite dependent on the inflation's current conditions.However, rates with longer maturities are less exposed to current conditions, and are more related to the long-term trends in expected inflation.Panel B of Table 2 shows descriptive statistics for inflation swap rates.As the maturity increases, the dispersion of the swap rates decreases; for instance, the 1-year swap rates have a standard deviation of 1.22, and the 30-year rates have a standard deviation of 0.34.One could argue that long-term swap rates are less affected by current economic conditions.The average swap rates exhibit a different pattern: they increase as a function of the maturity, going from 1.65% for 1-year swaps to 2.82% for 30-year swaps.For shorter maturities, the minimum value was reached during the last recession, when the US economy was in turmoil.
Zero-coupon inflation cap and floor prices are also extracted from the Bloomberg system.23 Maturities 1, 2, 3, 5, 7, 10, 12, 15, 20, and 30 are used for the period from October 2009 to December 2015.Unfortunately, option data before October 2009 are unavailable in the Bloomberg system.Strikes ranging from minus one percent to five percent in increments of 100 basis points are selected.Once more, Wednesday cap and floor prices are kept.
The quality of the dataset is checked by verifying that the selected caps and floors satisfy standard arbitrage conditions. 24Finally, Black and Scholes implied volatilities are computed for each option.
Table 3 exhibits the average implied volatilities for caps and floors.Means are computed for each maturity and strike across time.In general, the average implied volatility increases as a function of the maturity.Moreover, the well-known volatility smile effect is present in the inflation option data.For instance, the average 20-year implied volatilities decrease from K = −1% to K = 3%, and then increase, starting at K = 4%.This is related to deficiencies in the Black and Scholes option pricing model, which assumes constant volatility.Finally, throughout this study, the sample is divided in three subperiods: pre-recession (from 2004 to December 2007), recession (December 2007 to June 2009), and post-recession (June 2009 to 2015). 25

Estimation Results
The model is estimated using CPI-U levels, nominal risk-free bond prices, inflation swap prices, and option implied volatilities. 26The parameter estimates and their standard errors are given in Table 4.The robust standard errors (in parentheses) are calculated from the outer product of the gradient at the optimum parameter value.
Parameter κ is close to one (i.e., 0.9989) meaning that the effect of the current real interest rate is persistent.The unconditional real interest rate average level θ is 3.54 × 10 −4 under the physical measure, leading to an annualized average of 1.84%.The Q-measure unconditional real interest rate average level θ * is slightly lower-i.e., 3.24 × 10 −4 -which translates into an annualized average of 1.68%.These numbers are consistent with the low nominal interest and inflation rates observed during the 2004-2015 period.The volatility parameter σ is about 3.64 × 10 −5 .The real rate is not expected to change drastically from one week to another, to some extent explaining the low σ.
Parameters are estimated using the UKF with CPI-U levels, nominal risk-free bond prices and inflation swap prices only from July 2004 to December 2015, as well as option implied volatilities from October 2009 to December 2015.The log-likelihood function is maximized using multiple simplex search method optimizations (fminsearch in Matlab).Robust standard errors (in parentheses) are calculated from the outer product of the gradient at the optimum parameter value.For more details on this technique, see Appendix B.5 of Rémillard [33].
Parameter β is close to one (i.e., 0.9925), meaning that the current expected inflation rate is also persistent.This is consistent with the results of Ang et al. [12], who find that inflation is persistent over time.The long-term trend of expected inflation ν is 3.48 × 10 −4 , which translates into an annualized figure of 1.8%.This long-term trend is within the Fed target range for inflation-between 1.7% and 2.0%.Parameter ν * is slightly higher than its physical version ν; i.e., 7.05 × 10 −4 , or 3.66% when annualized.Volatility parameter γ is about 7.55 × 10 −5 .Again, we do not expect the expected inflation rate to vary severely from one week to another.In a similar framework, but using a different estimation technique, Fleckenstein et al. [7] find a slightly larger value for the expected inflation volatility parameter.The varying conditional volatility of the index level used in this study could explain this difference.The correlation coefficient between real interest rate innovations and inflation rate innovations is rather high: 0.708.Real interest and inflation rates tend to evolve in the same direction, as these two quantities are intimately related to the economic conditions of a given country.
Estimated GARCH parameter b is 0.9871 implying a high persistence in the conditional variance.The unconditional level of variance w is about 9.13 × 10 −7 , which translates into an annualized unconditional volatility of 0.69%.Even if parameter a is small in magnitude, it is somewhat substantial when compared to the unconditional variance average level.The positive asymmetry parameter c (i.e., 111.83) implies a negative correlation between the time t + 1 conditional variance and time t log-CPI level.In other words, when the level of the inflation index decreases, the conditional variance tends to increase.The parameter related to the error on bond log-prices ζ B is about 0.8%.Fitting a large number of bond prices in both time series and cross section could be ambitious, although the model is accommodating different bond prices, leaving little space for mistakes.The standard deviation on swap log-prices' errors ζ F is also small: about 0.2%.Parameter δ IV is about 48%, implying some aberrations in the model's implied volatilities.However, we feel the need to stress that is multiplied by observed option implied volatilities, leading to small standard deviations for option measurement equation errors.Parameter δ T is rather small, but positive.Therefore, the error on implied volatility is a function of the option's maturity: for instance, the standard deviation is increased by a factor of 1.2 for a maturity of 30 years.Finally, δ K is about 13.5, implying a positive relationship between ζ IV,K,t,T and the absolute value of K: the standard deviation of the relative error is raised by a factor of 1.4 when considering a strike of 5%.Mistakes in the prices, the over-the-counter nature of the market, and the long maturities used in this study could explain the implied volatility error terms.
The top and middle panels of Figure 2 show the behaviour of the filtered real risk-free interest rate and the expected inflation rate, respectively.The real interest rate R t is somehow stable at the beginning of the sample; then, during the last recession, the real interest rate decreased to reach −1.2% when annualized.In the post-recession era, the real interest rate increased briefly, and decreased again in 2012.The behaviour of the expected inflation rate µ t is rather different.It increased from about 2.0% to 4.5% during the pre-recession era.Then, starting in 2008, the expected inflation rate decreased and reached a low of −0.5% in 2010.Since then, the filtered value of µ t is about 0.8% on average in the last four years of the sample.The expected inflation µ t overall average over 2004-2015 is 1.6%, which is close to the Fed target regarding inflation.Yet, a lower average is consistent with the turmoil in the economy and the low realized inflation rates during the period considered.
The bottom panel of Figure 2 exhibits the annualized conditional volatility.As stated earlier, b is close to one, and there is therefore some persistence in the volatility process.The first half of the sample is somehow more volatile, with an average conditional volatility of 0.82%; the second half has an average of 0.57%.The last recession could explain this result, to some extent.Indeed, the conditional volatility jumped from about 0.5% before the recession to 1.54% in 2008.4.
To test the distribution fitting, the filtered values of ε I,t are investigated.From a visual inspection of the histogram of the ε I,t (left panel of Figure 3), we can conclude that the noise terms seem to be normally distributed.A QQ-plot (right panel of Figure 3) reveals the same conclusion.In addition, two statistical tests are used to verify that the ε I,t are indeed normally distributed: the chi-square and the Jarque-Bera goodness of fit tests fail to reject the null hypothesis, at a significance level of 5%: there is no statistical evidence of non-normality in the inflation index noise terms. 27

Estimation of Nested Models
To compare the performance of the model presented in Section 2, two competing models are now presented.Throughout the rest of the paper, these nested models are used to compare the full model: a discretized version of the Vasicek [16] interest rate model and a generalization of the Wilkie [9] inflation framework.

Real Interest Rate Expected Inflation Conditional Variance Error Standard Deviations
Parameters of Panel A (Nominal interest rate model) are estimated using the UKF with nominal risk-free bond prices from July 2004 to December 2015.Parameters of Panel B (Generalized Wilkie) are estimated using the UKF with CPI-U levels and inflation swap prices from July 2004 to December 2015.The log-likelihood function is maximized using multiple simplex search method optimizations (fminsearch in Matlab).Robust standard errors (in parentheses) are calculated from the outer product of the gradient at the optimum parameter value.For more details on this technique, see Appendix B.5 of Rémillard [33].
First, as it is common practice to model the nominal interest rate by using a unique process, we compare our approach to a simple autoregressive model:28 This specification is actually a discretized version of the Vasicek [16] model.It is estimated by using the same UKF-based approach as before, except that we discard data that require inflation modelling (i.e., inflation swaps and options).The estimated parameters are presented in Panel A of Table 5.
The second nested specification is a generalized version of Wilkie's [9] model.We generalize the latter to account for the fact that the expected and the realized inflation might not be the same: where w is the variance associated with realized inflation.This specification is nested in our framework: by setting R t = 0 for each t, and by assuming that θ = θ * = κ = σ = a = b = c = 0, we can recover the generalized Wilkie framework for inflation.The UKF-based approach is once again used to estimate the parameters.The time series inflation index level and inflation swaps are used in the estimation. 29anel B of Table 5 shows the estimated parameters for the generalized Wilkie (GW) model.The parameters of the two nested specifications seem reasonable.Yet, these models are misspecified: the Jarque-Bera test rejects normality for the two nested cases on r,t and I,t , respectively.The chi-square goodness of fit test rejects the null hypothesis for the first model.

Deflation Risk through Time
Once the models are estimated, it is possible to compute the probability density functions and the cumulative distributions for the inflation index level I t , given the filtered values of the real interest rate, the expected inflation rate, and the conditional variance.The latent variables evolve over time (as shown in Section 4.2), and imply changes in the distribution of the inflation.Details on the derivation of the inflation index level distribution are given in Appendix C.
A way of reporting the changing perceptions of deflation risk is by plotting the time series of deflation probabilities for a given time horizon.Deflation is simply defined as an annualized average inflation rate being below zero for a given time horizon τ; i.e., P (I t+τ ≤ I t | R t , µ t , h t ).Probabilities must be computed under the physical measure, since the evolution of the "real world" inflation is of interest here.The model and the estimation procedure employed in this study allow for the computation of such probabilities in a somewhat direct manner.The probabilities are forward-looking and based on current market's expectations.
The deflation probabilities are shown in Figure 4 for various time horizons and for two specifications of the model: the full version and the generalized Wilkie model.For the full model, the pattern of the different deflation probability time series of Figure 4 are similar across time horizons, even though their magnitudes tend to be larger for shorter maturities.Additionally, the deflation probabilities skyrocketed after the last NBER recession.For instance, on June 9, 2010, the deflation probabilities reached 54%, 47%, 42%, and 35% for 1-, 2-, 3-, and 5-year time horizons, respectively.For the same day, the GW's deflation probabilities were 1%, 1%, 2%, and 4% for 1-, 2-, 3-, and 5-year time horizons, respectively.
During the pre-recession period, the average 1-year deflation probability for the full model is 4%; it is about 0% for the generalized Wilkie model.The averages were about 18% and 2% during the recession for the full and the GW model, respectively.Lastly, the average 1-year probability for the full model increased to 29% during the post-recession era, while the average for the GW model was about 3%, on average. 30he values obtained from both models are different, and this is not a surprise: the GW model is not able to capture the left tail of the index inflation level distribution.To adequately model the tails, two things are needed: a model that accounts for changes in the level of uncertainty, and data that contains information about the tails.Indeed, option data contain information about the likelihood of extreme events, and are thus of utmost importance in assessing the deflation probabilities.   4 and Panel C of Table 5 were used to filter the real inflation rate, the expected inflation rate, and the conditional variance.Then, the latent variables were utilized to compute time series of probabilities.Deflation corresponds to an annualized average inflation that is less than 0% for a given time horizon, under the physical measure P. The grey-shaded region highlights the NBER-dated recession.
Other studies computed deflation probabilities using US inflation data.Kitsul and Wright [15] extract high inflation and deflation probabilities over the 2009-2013 period using option data.They find somewhat higher values than what is obtained in this study.However, the authors use an assumption of risk neutrality, and thus derive the probabilities under the Q-measure.The different measures could explain the differences in the results.In fact, Fleckenstein et al. [7] show that the ratio of risk-neutral to physical probabilities of deflation is often on the order of 1.30. 31The previously mentioned authors also find deflation P-probabilities similar to the ones presented in this study for most weeks, even though their model and estimation procedure are somewhat different.Christensen et al. [34] construct probability forecasts using yields on nominal and real US bonds.Their 1-year deflation probabilities are stuck around zero during the 2004-2008 period; this result is indeed somewhat similar to what we have obtained.They estimate the 1-year deflation probability to be close to one at the end of 2008, which is slightly higher than the probabilities of the full model presented in Figure 4. Discrepancies in the datasets and in the models used could explain these differences.

Impacts of Deflation on Life Insurance
The risk of deflation is real and could be rather important, as shown in Section 4. In this section, we focus on the implications of deflation risk on a fictitious life insurer.Indeed, deflation risk can have a material impact on the insurance industry.Life insurance is threatened by deflation risk, as negative inflation could have a significant impact on the prevailing nominal interest rates, and could thus increase the present value of the payments to be made in a distant future.Generally, the guaranteed interest rate is fixed for the whole term of the insurance contract.Thus, an inflation rate decline would make the guaranteed rate too high, as it cannot be reduced in general.Life insurers would therefore guarantee more interest than what they can earn.Additionally, inflation-indexed insurance policies and annuities might include implicit floor options.As a matter of fact, this kind of inflation protection would hurt life insurers in the eventuality of deflation.
A good example of the impact of deflation risk on life insurers is Japan: the nominal 10-year Japanese government bond yield dropped from a level of 8% in the early 1990s to less than 2% over the past decade.In addition, its historical inflation rate remained below zero for most of the last two decades.As a consequence of this extreme economic environment, eight Japanese insurers have gone bankrupt since 1996.For an analysis of deflation in the context of Japanese life insurers, see Hoshi and Kashyap [35].
Few papers discuss the role of deflation on life insurers.Ahlgrim and D'Arcy [8] acknowledge that the life insurance industry might be more affected by uninterrupted deflationary pressures, and that deflation makes it difficult to earn promised rates.Additionally, Dorfman et al. [36] find that deflation could have strong negative effects on life insurers while using a simple model for deflation.In their framework, the deflation risk is captured through a flat interest rate term structure that can decrease to lower levels.

Impacts on the Discount Rate
In this subsection (and the next), the focus is put on annuities.To this end, we consider a simple structure to assess the role of deflation risk on the present value of a life insurer's future payments.The insurer portfolio is a collection of immediate annuities for which the payments are made at the end of each year until the moment of death. 32Thus, the present value at time t of the life insurer's future payments is given by where N is the number of annuities sold, T i is the time to death of individual i, X i is the annual amount paid to the ith individual, and r * s is the time s insurer discount rate for this line of business.To simulate the life insurer portfolio, some additional assumptions are made.We assume that this portfolio contains only females aged 65.The mortality in the portfolio is random, and follows a Gompertz model fitted to US data.The parameters estimated by Pflaumer [37] are used in this study.33 The annual payment is set to X i = 1000, and the number of life insurance sold is N = 1000.The annualized discount rate r * s used by the life insurance company is assumed to be 2% over the nominal risk-free interest rate r s .
The purpose of this simple exercise is to assess the implications of deflation on life insurers without having to rely on complex assumptions.Complicated assumptions regarding mortality or the insured population could be easily taken into account, although it is not the purpose of the current study: deflation risk is our main concern.
The nominal interest rate is given by the full model: the nominal interest rate is thus defined as the sum of the real interest rate R s and the expected inflation rate µ s , as in Equation (1).In addition to the full model, three other models are used to compare the impact of the modelling assumptions on the risk measures of the fictitious life insurer.First, a model with stochastic real interest rate and deterministic inflation rate is used: the so-called anticipated inflation rate (AI) is calculated from the full model.This inflation rate is defined as the conditional expected value of the future inflation given today's economic conditions.It is therefore deterministic, but consistent with market participants' expectations.Second, the nominal interest rate (NIR) model presented in Section 4.3 is employed.In this framework, the sum of the real interest rate and the expected inflation rate is modelled directly through a univariate process (instead of a bivariate process, as in the full model).Third, the generalized Wilkie (GW) model described in Section 4.3 is used.As the latter does not model the real interest rate, we assume that the annualized real risk-free interest rate is constant and equal to 0.65%. 34,35 At this point, it is relevant to stress that the application does not focus on pricing.Instead, it focuses on risk assessment and management. 36We focus our attention at the impact of deflation upon two risk measures: the value-at-risk (VaR) and the conditional tail expectation (CTE).These risk measures are important to set up reserves and understand the risk profile of the insurer's portfolio.Remember that, in our setting, the nominal interest rate is not known.It is rather a random variable that has an impact on the tails of the future payments distribution.Specifically, decreases in the nominal interest-or in its constituents, the real interest rate and the expected inflation-create scenarios that would increase the present value of the portfolio, thus rising the risk measures.Figure 5 shows an example of four histograms of the life insurer's portfolio, in present value (5 February 2014).The four inflation assumptions are used: the full model (top left panel), the anticipated inflation assumption (top right panel), the nominal interest rate model (bottom left panel), and the generalized Wilkie model (bottom right panel).The full model allows for more dispersion than the other cases-especially the AI model and the GW model-leading to a larger VaR and CTE at a confidence level of 95%.The VaR(95%) is 21.9 when the full model is used.It is 17.2, 19.9, and 15.6 for the AI, the NIR, and the GW models, respectively.The right tail of the AI's distribution is very thin, yielding a CTE(95%) of only 19.3.As AI considers a stochastic real interest rate and a deterministic inflation rate, we can assume that a (market-consistent) deterministic inflation rate is not flexible enough to capture the tail risk adequately in the insurer's portfolio.The right tail of the GW's distribution is also thin: a CTE(95%) of 17.2 is associated to this model.Under this framework, the expected inflation is stochastic, but the real interest rate is deterministic.Therefore, to capture the tail behaviour of the insurer's portfolio, we also need a stochastic real interest rate.The expected inflation and the real interest rate may be modelled jointly via the nominal interest rate model, as is commonly done in practice.The NIR model yields a higher CTE(95%) than AI and GW.The full model and the nominal interest rate model histograms are also somewhat similar, even though the risk measures of the full model are slightly larger.This shows that modelling the sum of the expected inflation and the real interest rate through a unidimensional process could most probably be enough if one is only interested in the average behaviour of the future changes in the discount rate.Yet, the risk measures are not exactly the same, and separating the real interest rate and the expected inflation can have an impact on the specific behaviour of the risk measures, which depend explicitly on the levels of R t and µ t .Figure 6 displays risk measures for different levels of real interest rate and expected inflation, and for the four different inflation assumptions.These surfaces show the marginal impact of changes in the real interest rate and in the expected inflation.The AI surfaces' level is lower than the one of the full model, implying thinner tails for AI (27% and 38% lower on average for the VaR and the CTE, respectively).This observation is consistent with the results of Figure 5.For the generalized Wilkie model, the real interest rate has no marginal impact, as the GW model does not account for this dimension.Overall, the level of the GW surfaces is lower than the one of the full model (20% and 33% lower on average, respectively).On average, the full model VaR and CTE surfaces' level is about 5% and 6% higher than the one of the NIR model, respectively.4 and 5. Real interest rate (R t ) and expected inflation rates (µ t ) displayed in the figure are in percentages (annualized).The conditional variance is set to its unconditional average w (i.e., 9.1318 × 10 −7 ).Monte Carlo simulations are used to compute the risk measures.A sample of 10 5 observations are generated to compute each measure.
Even though the surfaces associated with the full model and NIR are tilted in the same direction-the risk measures increase when either the real interest rate or the expected inflation decrease-a marginal decrease in each of these two dimensions is not equivalent.Interestingly, for a given level of R t + µ t , the average VaR (CTE) of the full model is similar to the average VaR (CTE) of the NIR model.This implies that decomposing the nominal rate into its core components is important to capture market participants' expectations with respect to each source of risk.According to Figure 6, a marginal decrease of 1% in the real interest rate has the same impact as a marginal decrease of 1% in the expected inflation rate for the NIR model by construction.For the full model, this is not the case: a marginal decrease of 1% in R t is different when compared to a marginal decrease of 1% in µ t .4 and 5. Monte Carlo simulations were used to compute the risk measures.A sample of 10 5 observations were generated to compute each measure.The grey-shaded region highlights the NBER-dated recession.
Figure 7 shows the time series evolution of the VaR(95%) and the CTE(95%), computed using the four modelling assumptions.The four models are able to capture the dynamic nature of the risk, especially during the last NBER recession: there was an increase in both the VaR and the CTE during 2009.The full model yields risk measures that are larger than the other three models.Indeed, for AI and GW, this is not a surprise: these frameworks do not adequately capture the future evolution of the nominal interest rate, as one core component-either the real inflation rate or the expected inflation-is not modelled through a stochastic process.Even if the values used in the computation are market-consistent, a deterministic assumption is unable to capture the risk associated with low nominal rates.The NIR model does a better job; yet, the risk measure estimates using NIR are lower than the ones of the full model.A couple of reasons could explain these differences.First, the tails of the expected inflation distribution are better captured by the full model, as inflation caps and floors are used in the estimation of this model.Second, the separation of the real interest rate and the expected inflation brings a better understanding of the core constituents of the nominal rate, as these two components do not have the same marginal impact on the risk of the insurer's portfolio.
In summary, the nominal interest rate needs to be fully stochastic to adequately capture the right tail of the insurer's future payments distribution.For instance, AI and the GW model are not adequate, as an essential component of the discount rate is not modelled explicitly.Modelling the sum of the real interest rate and the expected inflation-instead of each component individually-captures the average behaviour of the risk in the right tail.Yet, the individual impact of each component does not seem to be symmetric, as shown in the top panels of Figure 6.Accounting for the real interest rate and the expected inflation separately allows for a better understanding of the portfolio's risk.

Impacts of Embedded Optionalities
Life insurers are exposed to deflationary pressures, as falling prices could lead to insufficient investment returns to support the rate credited to policyholders.Yet, this is not the only channel of transmission for deflation risk.Inflation-indexed protections could also make insurers vulnerable to deflation.This kind of insurance policy is constructed to cope with the risk of purchasing-power erosion.In practice, the paid amounts of inflation-indexed products follow some inflation index, such as the CPI.
Inflation-indexed annuities often provide additional protection that increases the paid amount in the case of positive inflation, but does not reduce the amount in the case of deflation.This added protection-similar to an inflation floor-covers the annuitant against decreases in the CPI level. 37This kind of protection-which is common in practice-is risky for life insurers, as deflationary pressures would increase the value of the embedded inflation floors.
Even though these inflation-indexed products are available, they are far less popular than equity-indexed products in the general population [38].Yet, there is still demand for similar inflation-indexed products through reinsurance channels.For example, an increasing number of defined-benefit pension plan sponsors decide to de-risk their financial exposures by buying pension buy-in and buy-out annuities.These special annuities transfer the risks borne by pension funds to life insurers.In the case of inflation-indexed pension benefits, the life insurer is exposed to the dependence between the inflation index level and the nominal interest rate.As most of these inflation-indexed pension plans do not reduce the retirement benefit amount in case of deflation, the life insurer thus bears additional deflation risk in a similar fashion to the indexed annuities described above-a consequence of the embedded inflation floor option.
The focus in this exercise is still on annuities.The portfolio is a collection of annuities for which the indexed paid amount is made at the end of each year until the moment of death.The insurer increases the paid amount using the inflation index when it is favourable for the annuitant.Thus, the present value at time t of the life insurer's future payments is given by The term max(1; I t+τ /I t ) is a combination of the protection, I t+τ /I t , and the embedded inflation floor paying the maximum between 1 − I t+τ /I t and zero.
The only model that can be used to compare the full model is the generalized Wilkie model, as it is the only one modelling the inflation index as a full-fledged variable. 38Figure 8 shows the time series evolution of the VaR and the CTE for both models. 39To some extent, the full model as well as the GW model are able to capture the changes in the risk as economic conditions unfold.Yet, the risk measures of the full model capture the increased deflation probabilities in the post-recession era, whereas the GW's risk measures are much more stable in time.The GW model is not able to capture deflation risk in the post-recession era, as the risk measures are back to their pre-recession levels.

CTE(95%)
Full model Generalized Wilkie Values-at-risk and conditional tail expectations were computed under two inflation assumptions: the full model and the generalized Wilkie model.The confidence level is set at 95%.Values are reported in millions.Computations were done using the parameters of Tables 4 and 5. Monte Carlo simulations were used to compute the risk measures.A sample of 10 5 observations were generated to compute each measure.The grey-shaded region highlights the NBER-dated recession.
The risk measures computed with the full model are higher than the ones of the GW model.As seen in Section 4.4, the deflation probabilities are not well captured by the GW model.Three main reasons might explain this: (1) the use of option in the model estimation; (2) the conditional heteroskedasticity of the inflation index noise terms; and (3) a fully stochastic nominal interest rate. 40 By using information coming from nominal bonds and inflation derivatives, we are thus able to break down the contributions of each component of the nominal interest rate and understand their specific impact.This is of utmost importance, as the future payments of life insurers depend on both the inflation index level and the nominal inflation rate.
The approach used in this study is different from the common practice that combines both the real interest rate and the expected inflation by modelling the nominal rate directly.Dividing these two effects allows for more flexibility and permits a better understanding of the underlying economic risk factors.It is also the only alternative to account for embedded inflation options such as the floors depicted in this exercise.

Concluding Remarks
For the past 20 years, inflation remained virtually constant.However, the recent recession unveiled the potential for unexpected changes in the price level.
The longstanding nature of insurance activities exposes insurers to deflation risk.Dynamic and market-consistent estimates of insurers' risk profile is of paramount importance to capture the changing economic conditions.To this end, a market-based methodology for assessing deflation risk is presented in this study.The proposed model accounts for the real interest rate, the inflation index level, its conditional variance, and the expected inflation rate.Pricing formulas for nominal bonds, inflation swaps, and inflation options are available in closed or semi-closed form expressions.
An estimation method based on the UKF was adopted.It captures the dynamics of the three latent variables.Multiple types of data source were utilized in the estimation: the use of inflation options-along with nominal risk-free bonds, inflation index levels, and inflation swaps-allows us to characterize the tails of the inflation index level distribution more appropriately.Additionally, the methodology incorporates the information as it becomes available in the market, leading to a market-consistent forward-looking assessment of deflation risk.This is important for proper estimation of deflation risk and its time-varying nature.
Then, we used US inflation data to estimate the framework.The evolution of the three latent variables is quite dependent on the market conditions.For instance, the filtered expected inflation rate decreased during and after the last recession, implying increased probabilities of deflation.Overall, the average 1-year deflation probabilities for the full model are estimated to be 18% and 29% for recession and post-recession eras, respectively.
The effects of deflation risk are analyzed in the context of life insurance.In general, the full model yields larger risk measures than the ones of the three nested models.We show that the nominal interest rate needs to be fully stochastic to capture the right tail of the insurer's future payments distribution.Additionally, modelling the core components of the nominal interest rate separately-in opposition to the common practice of modelling the nominal interest rate directly-allows for a better understanding of the portfolio's risk profile.By directly modelling the nominal interest rate, one cannot capture the asymmetric marginal impact of changes in both the real interest rate and the expected inflation, as shown in Figure 6, for instance.
Some life insurance products rely on the complex dependence structure between the inflation index level and the nominal interest rate.In these cases, modelling the core components of the nominal interest rate separately is the only alternative.For instance, the inflation-indexed annuities are impacted by both the inflation index level and the nominal interest rate, as investigated in this study. 40As explained in previous sections, options contain critical information about the tails of the inflation index distribution.
Additionally, flexible modelling assumptions are of paramount importance to capture the evolution of market-consistent risk measures.
Market-consistency links the likelihood of deflation scenarios to current market expectations.This study shows how deflation impacts the risk measures of a fictitious life insurer: the increases shown in the last section stress the need to include deflation in risk modelling and management-especially nowadays as deflation risk is more uncertain than during the last two decades.
which is the moment generating function of a Gaussian random variable of mean −Λ x,t and variance 1.
To get a noise term centred at zero, we define We can also find that First, we conjecture that the solution of E Q t exp − ∑ t+T s=t+1 R s has the following exponential affine form: where A R (t, t + T) and B R (t, t + T) are scalar coefficients.We know that One way to solve this expectation is by inverting the moment generating function.The inversion used in this study is closely related to the one of Gil-Pelaez [40] (as in Heston and Nandi [6] for instance).
The only thing missing is the moment generating function of log

Figure 3 .
Figure 3. Histogram and QQ-plot of inflation index noise terms ε I,t .The inflation index noise terms ε I,t are extracted from the UKF using the parameters of Table4.

Figure 4 .
Figure 4. Probability of deflation for 1-, 2-, 3-, and 5-year horizons.Parameters of Table4and Panel C of Table5were used to filter the real inflation rate, the expected inflation rate, and the conditional variance.Then, the latent variables were utilized to compute time series of probabilities.Deflation corresponds to an annualized average inflation that is less than 0% for a given time horizon, under the physical measure P. The grey-shaded region highlights the NBER-dated recession.

Figure 5 .
Figure 5. Histograms of life insurer's future payments on 5 February 2014, in present value.Portfolio present values are simulated using the method presented in Section 5 for four different inflation assumptions: full inflation (top left panel), anticipated inflation (top right panel), nominal interest rate (bottom left panel), and generalized Wilkie model (bottom right panel).Then, histograms are plotted and risk-measures are calculated.Values are reported in millions.CTE refers to the conditional tail expectation and VaR means value-at-risk.

Figure 6 .
Figure 6.Risk-measures for the life insurer's future payments, in present value.Values-at-risk and conditional tail expectations are computed under four inflation assumptions: full model, anticipated inflation, nominal interest rate, and the generalized Wilkie model.The confidence level is set at 95%.Values are reported in millions.Computations are done using the parameters of Tables4 and 5. Real interest rate (R t ) and expected inflation rates (µ t ) displayed in the figure are in percentages (annualized).The conditional variance is set to its unconditional average w (i.e., 9.1318 × 10 −7 ).Monte Carlo simulations are used to compute the risk measures.A sample of 10 5 observations are generated to compute each measure.

Figure 8 .
Figure 8.Time series of risk-measures for the life insurer's future payments, in present value.Values-at-risk and conditional tail expectations were computed under two inflation assumptions: the full model and the generalized Wilkie model.The confidence level is set at 95%.Values are reported in millions.Computations were done using the parameters of Tables4 and 5. Monte Carlo simulations were used to compute the risk measures.A sample of 10 5 observations were generated to compute each measure.The grey-shaded region highlights the NBER-dated recession.
Moreover, {ε I,t } ∞ t=0 and {ε µ,t } ∞ t=0 are two sequences of standardized Gaussian random variables.The real interest rate noises are correlated with the ones of the expected inflation rate; i.e., ρ = Corr P ε R,t , ε µ,t .The other sequences are uncorrelated.

Table 1 .
and η µ is the change of measure parameter associated to ε µ,t .Summary of the model parameters and their interpretation.
1, t + T) − 1)β are three scalar constants.One can compute these by recursion, starting with s = t + T, where A B (t + T, t + T) = B B (t + T, t + T) = C B (t + T, t + T) = 0. Appendix B.1 shows how to obtain the bond pricing formula.
Series ID: CPIAUCNS.Evolution of the inflation index level I, the risk-free interest rate for maturities of 1, 10, and 30 years in percentages, and the zero-coupon inflation swap rate for maturities of 1, 10, and 30 years in percentages.The inflation index level is defined as the three-month lagged linearly interpolated value of the Consumer Price Index for All Urban Consumers (CPI-U).The CPI-U level is collected via Federal Reserve Economic Data (FRED) from April 2004 to March 2016.The risk-free interest rate is proxied by both London Interbank Offered Rate (LIBOR) rates (for maturities less than 1 year) and swap rates (for maturities greater or equal to 1 year).The LIBOR and interest swap rates are collected via FRED from July 2004 to December 2015.The zero-coupon inflation swap rate is defined as f

Table 2 .
Summary statistics for risk-free interest rates and inflation swap rates.

Table 3 .
Summary statistics for inflation option implied volatilities.
Inflation option implied volatilities and strikes are expressed in percentage, and maturities in years.The zero-coupon inflation cap and floor prices are collected via the Bloomberg system from October 2009 to December 2015.The inflation cap and floor implied volatility are computed by inverting Black and Scholes' pricing formula.The average implied volatilities and the number of observations (N) are computed across time and strikes for each maturity.

Table 4 .
Parameter estimates and standard errors for the full model when using CPI-U levels, nominal risk-free bond prices, inflation swap prices, and option implied volatilities.

Interest Rate Expected Inflation Conditional Variance Error Standard Deviations
Evolution of the latent variables (annualized, in percentage) using CPI-U levels, nominal bond prices, inflation swap prices, and option implied volatilities.Real risk-free short rate R t , expected inflation rate µ t , and conditional variance h t are filtered using the UKF with CPI-U levels, nominal bond prices, and inflation swap prices from July 2004 to December 2015, as well as option implied volatilities from October 2009 to December 2015.Parameters of Table4are used.Values are annualized (i.e., R t and µ t is multiplied by 52; conditional volatility is √ 52h t ).The grey-shaded region highlights the NBER-dated recession.

Table 5 .
Parameters estimates and standard errors for nested versions of the model.

Nominal Interest Rate Model Nominal Interest Rate Expected Inflation Conditional Variance Error Standard Deviations
Time series of risk measures for the life insurer's future payments, in present value.Values-at-risk and conditional tail expectations were computed under four inflation assumptions: full model, anticipated inflation, nominal interest rate, and the generalized Wilkie model.The confidence level is set at 95%.Values are reported in millions.Computations were done using the parameters of Tables , t + T) = (B R (t + 1, t + T) − 1)κ, Hence,F(R t , µ t , T) = exp (A R (t, t + T) − A B (t, t + T) + (B R (t, t + T) − B B (t, t + T)) R t − C B (t, t + T)µ t ) .However, since B R (t, t + T) = B B (t, t + T), F(R t , µ t , T) = exp (A R (t, t + T) − A B (t, t + T) − C B (t, t + T)µ t )and the inflation swap price does not depend explicitly on R t .Appendix B.3.Cap and Floor OptionsThe payoff of a European inflation cap or call option with expiration T and strike price K is + µ s ) max 0,I t+T I t − (1 + K) T if C(R t , µ t , h t+1 , K, T)is the price of the inflation cap and r s = R s + µ s .The T-forward measure Q T is defined through the following Radon-Nikodym derivative: − ∑ t+T s=t+1 (R s + µ s ) B(R t , µ t , T) .