Next Article in Journal
Microinsurance and Economic Growth in Africa
Next Article in Special Issue
New Classes of Distortion Risk Measures and Their Estimation
Previous Article in Journal
Multifactor Risk Attribution Applied to Systemic, Climate and Geopolitical Tail Risks for the Eurozone Banking Sector
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Three-Factor Market Model for Incorporating Explicit General Inflation in Non-Life Claims Reserving

Department of Economics, University of Perugia, 06123 Perugia, Italy
Risks 2023, 11(10), 174; https://doi.org/10.3390/risks11100174
Submission received: 1 September 2023 / Revised: 27 September 2023 / Accepted: 3 October 2023 / Published: 7 October 2023
(This article belongs to the Special Issue Advancements in Actuarial Mathematics and Insurance Risk Management)

Abstract

:
In a recent paper “Stochastic Chain-Ladder Reserving with Modeled General Inflation”, the effects of modeled general inflation on non-life claims reserving were studied using, along with the so called “market approach”, a stochastic two-factor market model, characterized by deterministic expected inflation. In the present paper, we repeat the same study, again with the market approach, using a three-factor market model which extends the two-factor model by including stochastic expected inflation. After detailing the theoretical model and estimating the relevant parameters on the same market data used in “Stochastic Chain-Ladder Reserving with Modeled General Inflation”, we repeat the application to claims reserving presented in that paper and compare the results obtained with the two models. With these data, it is found that the inclusion of stochastic expected inflation produces a non-negligible increase in the reserve solvency capital requirement under the one-year view.

1. Introduction

In a recent paper (De Felice and Moriconi (2023), DFM23 for short), the problem is considered of explicitly including general inflation in non-life stochastic claims reserving using market-consistent approaches. An important issue discussed in DFM23 is that, given the presence of inflation-linked securities regularly quoted in the market, inflation risk should not be considered unhedgeable, in the sense defined in the Solvency II framework. Therefore, the technical provisions should not be obtained by adding to the best estimate an inflation risk margin calculated by the insurer. For hedgeable liabilities, the risk margin is actually a market risk premium, which is to be considered incorporated in the market prices, and which therefore must be derived using a pricing model estimated on market data. This problem is well-known in the market-consistent valuation of life insurance products, where the interest rate risk plays an important role, and many insurance benefits, such as those provided by profit-sharing policies, must be valued—and hedged—using complex market pricing models. These methods, however, are not part of the toolkit of a typical non-life actuary, essentially because in non-life insurance applications, interest rate risk is considered immaterial. The problem of including inflation risk in non-life claims reserving now seems to reopen the question.
A correct assessment of the inflation effects requires the use of appropriate models for real and nominal interest rates. In particular, only by using these methodologies can the effect on the capital requirements, as defined in a Solvency II internal model, be properly measured and controlled; only these methods can show the possible immateriality of these risks. Obviously, there remains the need to consider an explicit risk margin for technical risks (i.e., pure claims development uncertainty), since those risks remain intrinsically unhedgeable.
All of these issues have been addressed and analyzed in DFM23, but, for simplicity’s sake, a model with deterministic (future) expected inflation was used in that paper. In the present paper, an alternative model is proposed that also includes uncertainty on the expected inflation rates. In this valuation framework, a correct assessment of the market risk premia is clearly a central issue. As is well-known, in pricing models based only on the no-arbitrage principle (also referred to as partial equilibrium models), which are typically used for asset trading purposes, the functional form of the risk premia is undetermined, since it cannot be derived by the model assumptions, and for this reason, it is usually given exogenously. This could imply some degree of arbitrariness in the valuations. To avoid this problem, we propose in this paper a general equilibrium model, where the market premia for the relevant risk factors are consistently derived by the model assumptions. Three risk factors are required to correctly model our problem, and some basic concepts in financial economics need to be used. However, we try to maintain our approach at the simplest possible level, referring to an elementary real and nominal economic framework. Moreover, the real interest rate and the expected inflation rate are modeled as Ornstein–Uhlembeck processes, which provides an essentially Gaussian model that is easy to understand and apply.
A three-factor model with an identical structure but based on “mean-reverting square-root” processes was proposed in (Moriconi 1995). As is well-known, these Cox–Ingersoll–Ross (CIR)-type processes have the property of precluding negative nominal interest rates, a feature considered necessary at the time but which became unacceptable after the period of negative nominal rates that began in 2016 in the market. Instead of correcting the model à la CIR using mathematically complex devices, we choose here the path of simplicity, using Gaussian, Vasicek-type processes, notwithstanding some of their empirical weaknesses.
The proposed model can be easily used for pricing and hedging interest rate derivatives and inflation-linked securities. An example of a not very different model that pursues these objectives is the one proposed by Jarrow and Yildirim (2003). However, since our aim is to provide a market-consistent stochastic version of the actuarial approach typically used in non-life insurance (and illustrated in DFM23), we will limit ourselves to developing pricing formulas for real, nominal, and inflation pure discount bonds, without dealing with problems, for example, of financial option pricing. On the other hand, since our main purpose is to derive the probability distribution of the year-end obligations of a non-life insurer, we are particularly interested, in addition to the risk-neutral probabilities required for pricing, in the natural (“real world”) probabilities, which are the basis for risk management.
The present paper is organized as follows. The structure of the three-factor model is presented in Section 2. We first describe the real interest rate component of the model and derive the explicit expression of real zero-coupon bonds. We then introduce money into this real economy framework and derive the two-factor model for the CPI, which provides the inflation component of the model. The explicit expression for the inflation discount factors is derived and the nominal interest rates are obtained from the Fisher relation. The estimation of the model is presented in Section 3, where the parameters of the inflation component and the real interest rate component are estimated in two separated steps and the real interest rate term structure is then calibrated on market data using the Hull–White model. In Section 4, the model is applied to non-life claims reserving, using the same data and actuarial assumptions as in DFM23. The differences in the numerical results from the two-factor model are highlighted and discussed. In order to make the paper more self-contained, some fundamental results on Ornstein–Uhlembeck processes are provided in Appendix A. In Appendix B, the general valuation equation in the nominal economy setting is derived, which is useful for clarifying the “economic” derivation of the inflation discount factor presented in the text of the paper. The essential details for the application of the Hull–White model are presented in Appendix C.

2. The Three-Factor Market Model

In order to model the inflation effects, measured by the percentage changes of a specified Consumer Price Index (CPI), we consider a simple general equilibrium model under uncertainty.

2.1. The Real Economy Setting

We start considering an elementary real economy, where a single consumption good B is traded. In our inflation problem, B will be interpreted as the basket of goods and services which defines the reference CPI.
  • Assumptions on the market. In this economy, all values (and all securities) are measured in terms of units of B , which will be denoted by 1 ¯ . We shall consider a market defined in continuous time which is perfect and frictionless; essentially, all agents are price taker, securities are infinitely divisible, there are no taxes or transaction costs, and riskless arbitrage opportunities are precluded. The representative agent in this market is non-satiated (he prefers more to less) and risk-averse, with logarithmic individual utility function. This agent allocates his wealth W ¯ ( t ) (measured in 1 ¯ ) seeking to maximize the expected utility of consumption over a fixed time horizon. Moreover, all the securities to be priced in this market are not explicitly dependent on wealth.
  • Real bonds and interest rates. Let us denote by V ¯ ( t ) the real value function at time t, that is, the market price expressed in units of B at time t. The pricing problem of a real bond traded in this market can be represented with some degree of generality by denoting by V ¯ ( t ; Y ¯ τ ) the real value at time t of the future payoff Y ¯ τ , where Y ¯ is a sequence of real cash-flows to be received on a corresponding sequence τ of dates after t. The payoff Y ¯ τ can also be stochastic, i.e., not known at time t. However, we shall consider only (default-)risk-free bonds, i.e., bonds for which Y ¯ τ is known on dates τ and paid for certainty on those dates.
The simplest kinds of real bonds are the unit real zero-coupon bonds (ZCBs), i.e., the bonds which provide the single payoff 1 ¯ at the maturity date T t . The time t price of this ZCB with maturity T is denoted by
v ¯ ( t , T ) : = V ¯ ( t ; 1 ¯ T ) .
If we consider a ZCB with deterministic payoff Y ¯ T to be received at time T, we have, to prevent arbitrages,
V ( t ; Y ¯ T ) = v ¯ ( t , T ) Y ¯ T ,
that is, v ¯ ( t , T ) is the discount factor to be used to obtain the market price (or present value) of this bond (for this reason, the unit ZCB is also referred to as a discount bond).
The (logarithmic) rate of return, the “log-return”, corresponding to the discount factor v ¯ ( t , T ) , is given by
h x ( t , T ) : = log v ¯ ( t , T ) T t .
The real instantaneous rate of interest  x ( t ) at time t is defined as the rate of return of a ZCB currently maturing, i.e.,
x ( t ) : = lim T t h ¯ ( t , T ) .

2.2. The Single-Factor Model for Real Interest Rates

Under uncertainty, at time t, the future prices and returns of real bonds are not known. In order to model these prices and rates, they are represented as known functions of one or more sources of uncertainty, or risk factors (also referred to as base variables). We shall consider for the real interest rates a univariate stochastic model having x ( t ) as the single risk factor. More specifically, we shall assume the following Markov property:
V ¯ ( t ; Y ¯ τ ) = V ( x ( t ) , t ; Y ¯ τ ) ,
where the price V ¯ ( t ) at time t depends only on the current value x ( t ) of x. In this case, x ( t ) is also referred to as the state variable.
We assume that x ( t ) is a diffusion process described by the following stochastic differential equation (SDE):
d x ( t ) = f x ( x ( t ) , t ) d t + g x ( x ( t ) , t ) d W x ( t ) ,
where W x ( t ) is a Wiener process. The choice of the functions f x and g x , together with the initial condition x ( t ) = x t , completely specifies the conditional probability distribution of x ( T ) given x t for T > t .
We shall assume for this process the “Vasicek specifications” (see Vasicek (1977)), i.e., a mean-reverting drift f x and a constant diffusion coefficient g x 2 . Precisely,
f x ( x ( t ) , t ) = α x ( γ x x ( t ) ) , α x , γ x constant , α x 0 , g x ( x ( t ) , t ) = σ x 0 constant .
With these f x and g x functions, the process (2) is a special case of the Ornstein–Uhlenbeck (OU) process, which we characterize by the following SDE:
d z ( t ) = ( a z ( t ) + b ) d t + σ d W z ( t ) , a , b R , σ > 0 .
In Appendix A, we provide some fundamental results on the z ( t ) process that will be useful for deriving many important properties of our model.
A first property is that in the Vasicek case, the conditional probability distribution of x ( T ) given x ( t ) , for T = t + τ t , is normal, with mean
E x ( T ) | x ( t ) = γ x γ x x ( t ) e α x τ
and variance
Var x ( T ) | x ( t ) = σ x 2 2 α x 1 e 2 α x τ ,
which is independent of x ( t ) . This normality property is Result A1 in Appendix A.1, and expressions (5) and (6) are obtained by (A4) and (A5), respectively, after the changes a = α x and b = α x γ x (and, obviously, interpreting the OU process z ( t ) as x ( t ) ).
A fundamental tool for no-arbitrage pricing in a perfect market is the so-called hedging argument. It essentially consists in composing a portfolio of securities that is instantaneously riskless and then requiring that, to avoid arbitrages, the instantaneous return on this portfolio is equal to the risk-free rate. In a diffusion model, this condition provides the general valuation equation (GVE), a deterministic differential equation that must be satisfied by the price of all securities traded in the market.
In our model, the hedging argument requires that the price V ¯ ( t ) of any real security which depends only on x ( t ) and t must satisfy the no-arbitrage relation:1
E t d V ¯ ( t , x t ) V ¯ ( t , x t ) d t = x ( t ) + ϕ x ( t , x t ) V ¯ ( t , x t ) V ¯ ( t , x t ) x .
This expression prescribes that the instantaneous expected return of V ¯ must be equal to the instantaneous real interest rate x plus an additional return that reflects the market risk premium for the risk factor x. This excess return is given by a factor ( V ¯ / x ) / V ¯ expressing the semi-elasticity of V ¯ with respect to x, multiplied by the risk premium coefficient  ϕ x ( t , x t ) . The form of this coefficient is to be specified, since, by Ito’s lemma,
E t d V ¯ ( t , x t ) d t = 1 2 g x 2 2 V ¯ x 2 + f x V ¯ x + V ¯ t ,
expression (7) can be written as
1 2 g x 2 2 V ¯ x 2 + ( f x ϕ x ) V ¯ x + V ¯ t = x V ¯ .
This deterministic (parabolic) partial differential equation is the GVE of our model. Once the coefficients f x , g x , and ϕ x are specified, the price V ¯ ( t ) of any security will be obtained as the solution of (9) under the appropriate boundary conditions (essentially determined by the contractual characteristics of the security). In particular, for the zero-coupon bond providing the payoff Y ¯ T at time T t , this condition is V ¯ ( T ) = Y ¯ T .
Another fundamental result for our pricing problem is obtained by the well-known Feynman–Kac theorem. According to this theorem, the solution of the GVE (9) with terminal condition V ¯ ( T ) = Y ¯ T can be expressed in integral (i.e., in expectation) form as
V ¯ ( t ; Y ¯ T ) = E ˇ t R e t T x ( s ) d s Y ¯ T ,
where E ˇ t R ( · ) denotes the conditional expectation taken with respect to the so-called risk-neutral probability measure P ˇ R , that is, the probability determined by the risk-adjusted drift f ˇ x : = f x ϕ x and the diffusion coefficient g x 2 . Since f x and g x have been chosen in the Vasicek form, to completely define the pricing problem, it remains to specify the form of the risk premium ϕ x .

2.3. The Market Risk Premium and the ZCB Pricing Formula

2.3.1. Equilibrium and the Market Risk Premium

In a partial equilibrium model, where prices are determined solely by no-arbitrage arguments, the form of the risk premia must be given exogenously, which introduces some degree of arbitrariness into the pricing problem. In a general equilibrium model, instead, the market risk premium for a given state variable is endogenously determined, consistent with the model assumptions. This consistency property motivated our choice of a general equilibrium approach, given that risk premia can have a significant impact on the level of prices. It can be shown (see Cox et al. (1985)) that in our model, ϕ x is given by the covariance of changes in the state variable x ( t ) with percentage changes in the optimally allocated wealth W ¯ * ( t ) (the value of the so-called market portfolio), i.e.,
ϕ x ( x ( t ) , t ) = Cov d W ¯ * ( t ) W ¯ * ( t ) d t , d x ( t ) .
With arguments similar to those in (Cox et al. 1985), we find that in our model, the real interest rate risk premium must be an affine function of x ( t ) :
ϕ x = η π x ( t ) , η , π constant .
With this choice, the risk-adjusted drift is
f ˇ x = ( α x γ x η ) ( α x π ) x ( t ) ,
and the pricing problem (9), or (10), is completely defined. It is convenient to write the risk-adjusted drift as a function of risk-adjusted “mean-reverting” parameters,2 that is,
f ˇ x = α ˇ x γ ˇ x α ˇ x x ( t ) , with α ˇ x : = α x π , γ ˇ x : = α x γ x η α x π .

2.3.2. The ZCB Pricing Formula

For our purposes, it is sufficient to consider the pricing of the real unit ZCB. In this case, we have
v ¯ ( t , T ) = E ˇ t R e t T x ( s ) d s
and we find that the following closed-form expression holds, for T = t + τ 0 :
v ¯ ( t , t + τ ) = A x ( τ ) e x ( t ) B x ( τ ) ,
with
B x ( τ ) : = 1 e α ˇ x τ α ˇ x
and
A x ( τ ) : = exp γ ˇ x σ x 2 2 α ˇ x 2 B x ( τ ) τ σ x 2 4 α ˇ x B x 2 ( τ ) .
In order to derive this expression, let us denote by φ X ( s ) = E e s X , s R , the moment-generating function of the random variable X. As is well-known, in the normal case X N ( μ , σ 2 ) , one has φ X ( s ) = e s μ + s 2 σ 2 / 2 . The ZCB price v ¯ ( t , T ) is the exponential expectation E ˇ t R e X ( t , T ) with X ( t , T ) : = t T x ( s ) d s . Using Result A2 in Appendix A.3, after the changes
z ( t ) = x ( t ) , Z ( t , T ) = X ( t , T ) , a = α ˇ x , b = α ˇ x γ ˇ x ,
we find that X ( t , T ) is normal, so
v ¯ ( t , T ) = φ X ( t , T ) ( 1 ) = exp E t X ( t , T ) + 1 2 Var t X ( t , T ) ,
which provides (14) if the expressions (A14) and (A15) of Result A2 are used, respectively, for the mean and the variance.
Expression (14) plays the role of the fundamental pricing formula in our three-factor model since, as we will see, the expression of the inflation discount factor and that of nominal ZCB price can also be considered variants of this formula.
Expression (14) implies
h x ( t , t + τ ) = log A x ( τ ) τ + x ( t ) B ( τ ) τ .
Hence, the log-return at time t, for the given τ , is an affine function of the risk factor x ( t ) . Moreover, it is immediately verified that
lim τ 0 h x ( t , t + τ ) = x ( t ) .

2.3.3. The Sensitivity Function and the Effects of Mean-Reversion

For a better interpretation even of some of the results of our actuarial applications, it is useful to recall some properties of the function B x ( τ ) defined in (15). By (14),
B x ( τ ) = v ¯ ( t , t + τ ) v ¯ ( t , t + τ ) x .
Hence, B x ( τ ) is equal to (minus) the semi-elasticity of the price v ¯ with respect to x. It provides the percentage changes in the bond price attributable to an unexpected shift in the risk factor x and is then also referred to as the sensitivity of v ¯ w.r.t. x. For this reason, B x ( τ ) is typically used as a measure of the real interest rate risk of the ZCB v ¯ ( t , t + τ ) . B x ( τ ) is measured in time units, is zero for τ = 0 , and is monotonically increasing with τ . If α ˇ x > 0 (genuine mean-reversion), then B x ( τ ) < τ for τ > 0 and has a horizontal asymptote at level 1 / α ˇ x ; the mean-reversion produces a ceiling on the interest rate risk. If α ˇ x < 0 (“mean repulsion”), then B x ( τ ) > τ for τ > 0 and diverges for τ . The function B x ( τ ) for different values of α ˇ x is illustrated in the left figure of Table 1.
In our applications, it will be important to consider the variance of future values of the discount factors v ¯ ( T , T + τ ) for T = t + τ 0 . For the corresponding log-returns, we have, by (16),
Var h ( T , T + τ ) | x ( t ) = Var x ( T ) | x ( t ) B x ( τ ) τ 2 ,
where the conditional variance of x ( T ) , given by (6), depends on the parameters σ x , α x . As illustrated in the right figure of Table 1, for α ˇ x > 0 , the factor B x ( τ ) / τ as a function of τ takes values in ( 0 , 1 ] , is monotonically decreasing, and decreases more rapidly for increasing values of α ˇ x . This illustrates the important role that the risk-neutral mean-reversion plays in reducing the variability of returns and discount factors when the maturity changes.

2.3.4. Simulation of Future Prices and Rates

Using expression (14), at time t, future values of ZCB prices, i.e., future discount factors v ¯ ( t , t + τ ) at time t > t , can be simulated at time t as:
v ¯ ˜ ( t , t + τ ) = A x ( τ ) e x ˜ ( t ) B x ( τ ) .
The conditional mean and variance of the normal variable x ( t ) are given by (5) and (6), and the random variable v ¯ ( t , t + τ ) is lognormal. Of course, h x ( t , t + τ ) is normal.

2.4. Introducing Money: The Two-Factor Model for the Price Index

Obviously, inflation effects can only occur when money is present in the economy. Therefore, to model inflation, we introduce money in this real market and denote by 1 one unit of money (say, EUR 1). We denote by p ( t ) the price in units of money of one unit of the reference basket B which defines the CPI. In this simple model, we assume that the introduction of money has no effects on the underlying equilibrium and money is only a measure of nominal value, that is, p ( t ) is simply the number of 1 needed for receiving one 1 ¯ at time t. However, this is sufficient for modeling nominal bonds.
Let us denote by V ( t ) the nominal value function at time t, that is, the market price expressed in EUR at time t. The general relation between nominal and real price is
V ( t ) = p ( t ) V ¯ ( t ) ,
because receiving one real unit, 1 ¯ , at time t is equivalent to receiving p ( t ) nominal units, 1 , at the same date. Let us consider the real value V ¯ ( t ; 1 T ) of a bond which provides with certainty the unit nominal payoff 1 T at maturity T t . Since, by relation (20), 1 ¯ T = p ( T ) 1 T , to prevent arbitrage, the real value of this bond can be represented by
V ¯ ( t ; 1 T ) = V ¯ t ; 1 / p ( T ) .
Then, for the nominal unit ZCB price v ( t , T ) , i.e., the money value at time t of one money unit at time T, we have, again by relation (20),
v ( t , T ) = p ( t ) V ¯ t ; 1 / p ( T ) .
The same argument applies to more general nominal bonds. Therefore, in order to model nominal bond prices—or, equivalently, to model prices of the form V ¯ t ; 1 / p ( T ) —we have to extend the model for x ( t ) by making appropriate assumptions on the stochastic dynamics of p ( t ) .

2.4.1. The Two-Factor Model for the Price Index

We assume that the CPI is a diffusion process described by the following SDE:
d p ( t ) = f p ( p ( t ) , t ) d t + g p ( p ( t ) , t ) d W p ( t ) ,
with
f p ( p ( t ) , t ) = y ( t ) p ( t ) , g p ( p ( t ) , t ) = σ p p ( t ) , σ p > 0 .
Since E t d W p ( t ) = 0 , one has
y ( t ) = E t d p ( t ) p ( t ) d t .
Therefore, y ( t ) is the expected instantaneous rate of inflation at time t. If y ( t ) is deterministic, the CPI process is a geometric Brownian motion. However, in order to obtain a more realistic model, we assume that y ( t ) is also stochastic.
Our assumption is that the expected instantaneous inflation rate is a diffusion process described by the following SDE:
d y ( t ) = f y ( y ( t ) , t ) d t + g y ( y ( t ) , t ) d W y ( t ) ,
with
f y ( y ( t ) , t ) = α y ( γ y y ( t ) ) , α y , γ y constant , α y > 0 , g y ( y ( t ) , t ) = σ y , σ y > 0 constant .
With this Vasicek specification, y ( t ) is also an Ornstein–Uhlembeck process.
As concerns the correlation structure, we assume
Cov t [ d W y ( t ) , d W p ( t ) ] = ρ d t , ρ [ 1 , 1 ] constant ,
and
Cov t [ d W x ( t ) , d W y ( t ) ] = Cov t [ d W x ( t ) , d W p ( t ) ] = 0 .
Assumptions (26) state the independence between real and nominal quantities and have important consequences. First of all, the form of the market risk premium ϕ x is the same as in (12) and the market risk premia for the new two state variables are zero:
ϕ y ( t ) = ϕ p ( t ) = 0 ,
since property (11) also holds for the nominal quantities ϕ y and ϕ p and W ¯ * ( t ) is a real quantity.

2.4.2. The Inflation Discount Factor and the Fisher Relation

Since the real payoff Y ¯ τ is independent of p and y, the price of real bonds V ¯ ( t ; Y ¯ τ ) is still obtained by the univariate model for x ( t ) . In particular, v ¯ ( t , T ) still has the Vasicek form (14).
As concerns V ¯ ( t ; 1 / p ( T ) ) , we have, by extension of (10):
V ¯ t ; 1 p ( T ) = E ˇ t R e t T x ( s ) d s 1 p ( T ) ,
where now the risk-neutral measure P ˇ R is three-dimensional. However, by the independence assumptions (26), we have
V ¯ t ; 1 p ( T ) = E ˇ t R e t T x ( s ) d s E ˇ t R 1 p ( T ) = v ¯ ( t , T ) E ˇ t R 1 p ( T ) .
Hence,
v ( t , T ) = v ¯ ( t , T ) E ˇ t R p ( t ) p ( T ) .
For the expectation on the right-hand side, we observe that, since in the real economy, the risk premia ϕ y and ϕ p are set to zero, the relevant risk-neutral measure is equal to the natural probability measure. We define the folloing:
u ( t , T ) : = E ˇ t R p ( t ) p ( T ) = E t p ( t ) p ( T ) ,
which will be referred to as the inflation discount factor on the time interval [ t , T ] . Hence, we have
v ( t , T ) = v ¯ ( t , T ) u ( t , T ) .
Expression (28) can be referred to as the Fisher relation in multiplicative form. As we shall see, since y ( t ) is stochastic, future expectations, i.e., u ( t , T ) for t < t < T , are not known at time t. That is, like future real discount factors, future inflation discount factors too are stochastic.
If we consider the inflation and the nominal log-returns,
h y ( t , T ) : = log u ( t , T ) T t , h r ( t , T ) : = log v ( t , T ) T t ,
expression (28) provides
h r ( t , T ) = h x ( t , T ) + h y ( t , T ) ,
where h x ( t , T ) is given by (1). This is the celebrated Fisher relation in its classical form, representing the nominal interest rate as the sum of the real interest rate and the expected inflation rate.

2.4.3. The Expression for the Inflation Discount Factor

Let us define the CPI logarithmic price ratio on the time period [ t , T ] :
L ( t , T ) = log p ( T ) p ( t ) .
Of course,
u ( t , T ) : = E t p ( t ) p ( T ) = E t e L ( t , T ) .
Using Result A3 in Appendix A.4, with a = α y and b = α y γ y , for T = t + τ t , the random variable L ( t , T ) is normal, with mean
μ L ( t , T ) : = E t L ( t , T ) = y ( t ) B y ( τ ) + γ y τ B y ( τ ) 1 2 σ p 2 τ
and variance
σ L 2 ( t , T ) : = Var t L ( t , T ) = σ p 2 τ σ y 2 α y 2 + 2 ρ σ y σ p α y ( B y ( τ ) τ ) σ y 2 2 α y B y 2 ( τ ) ,
where B y ( τ ) is given by β ( τ ) with a = α y , that is,
B y ( τ ) : = 1 e α y τ α y .
It is important to note that these expressions are independent of the current level p ( t ) of the CPI. Furthermore, given the model parameters, the variance σ L 2 ( t , T ) only depends on τ = T t and the mean μ L ( t , T ) depends on t only through y ( t ) .
Using (31) and (32), an expression for u ( t , T ) is immediately obtained by computing the following moment-generating function:
φ L ( t , T ) ( 1 ) = e μ L ( t , T ) + 1 2 σ L 2 ( t , T ) .
We obtain
u ( t , t + τ ) = A y ( τ ) e y ( t ) σ p 2 B y ( τ ) ,
where B y ( τ ) is given by (33) and
A y ( τ ) : = exp γ y σ p 2 ρ σ y σ p α y σ y 2 2 α y 2 B y ( τ ) τ σ y 2 4 α y B y 2 ( τ ) .
For the inflation discount factor, the same properties already illustrated for the real discount factor hold; at time t, the inflation log-returns are affine functions of y ( t ) , and future values u ( t , T ) and h y ( t , T ) , for t < t < T , are, respectively, lognormally and normally distributed. Moreover, it is easily verified that
lim τ 0 h y ( t , t + τ ) = y * ( t ) : = y ( t ) σ p 2 .
For the function B y ( τ ) given by (33), the same properties and the same interpretation hold as for B x ( τ ) . B y ( τ ) is a time measure of expected inflation risk, providing the sensitivity of the inflation discount factor u ( t , t + τ ) to the risk factor y ( t ) . The only difference is that, being that α y is the natural risk reversion parameter, it is positive by assumption, so only the “capped” behavior (asymptote at 1 / α y ) is possible for B y ( τ ) ; the difference τ B y ( τ ) is strictly positive and monotonically increasing with τ > 0 .

2.4.4. Simulation of Future Rates and Prices

Also in this case, expression (35) allows us to easily simulate at time t future inflation discount factors as lognormal variables. More interesting for our applications, future values p ( T ) , T > t , of the CPI are also simulated as lognormal, that is,
p ˜ ( T ) = p ( t ) e L ˜ ( t , T ) ) = p ( t ) e μ l ( t , T ) + σ L ( t , T ) ε ˜ ,
with ε standard normal. Taking account of (34), expression (38) can be also written as
p ˜ ( T ) = p ( t ) 1 u ( t , T ) e 1 2 σ L 2 ( t , T ) + σ L ( t , T ) ε ˜ .
This expression is important because if the discount factors u ( t , T ) can be observed in the market (as will be the case in our applications), future values of the CPI can be simulated in a market-consistent manner, i.e., consistent with the current market expectations on future inflation. This is equivalent to posing μ L ( t , T ) = log u ( t , T ) + σ L 2 ( t , T ) / 2 . In this case, since y ( t ) and γ y do not enter into the expression of the variance σ L 2 ( t , T ) , their observation or estimation is not required for the simulations of p ( T ) .
Remark 1. 
In the two-factor model used in DFM23, where y ( t ) is deterministic, the risk-neutral inflation discount factor is defined as
u ˇ ( t , T ) : = 1 / E ˇ t p p ( T ) p ( t ) = 1 / E ˇ t p e L ( t , T ) ,
where the expectation is taken w.r.t. the risk-neutral probability measure P ˇ p of the price process p ( t ) in the nominal economy setting and we use again the definition L ( t , T ) = log p ( T ) / p ( t ) . However, in the two-factor model, one has E ˇ t p e L ( t , T ) = e t T ( y ( s ) σ p 2 ) d s , then u ˇ ( t , T ) = e t T ( y ( s ) σ p 2 ) d s = E ˇ t p e L ( t , T ) , as in (30). Thus, the two models are formally consistent. Of course, in the three-factor model, the distribution of L ( t , T ) is more complex than in the case with deterministic expected inflation.

2.4.5. An “Economic” Derivation of the Inflation Discount Factor

Alternative derivations of expression (35) are available in the literature. We illustrate one which is interesting for its economic meaning. This derivation originates by the following observation. If we define a “new” Vasicek-type stochastic process y ˇ ( t ) , with drift parameters
α ˇ y = α y , γ ˇ y = γ y σ p 2 ρ σ y σ p α y ,
the process y ˇ * ( t ) : = y ˇ ( t ) σ p 2 is described by the following SDE:
d y ˇ * ( t ) = α ˇ y γ ˇ y y ˇ * ( t ) d t + σ y d W y ˇ * ( t ) .
By Result A2 in Appendix A.3, the stochastic integral Y ˇ * ( t , T ) : = t T y ˇ * ( s ) d s is normal with mean
μ Y * ( t , T ) : = E t Y ˇ * ( t , T ) = y ˇ * ( t ) B y ( τ ) + γ ˇ y τ B y ( τ )
and variance
σ Y * 2 ( t , T ) : = Var t Y ˇ * ( t , T ) = σ y 2 α y 2 ( B y ( τ ) τ ) σ y 2 2 α y B y 2 ( τ ) .
Then, it is easily verified that expression (35) can also be obtained as
u ( t , T ) = E t e t T y ˇ * ( s ) d s .
The interesting point is that this is equivalent to saying
u ( t , T ) = E ˇ t N e t T y * ( s ) d s ,
where y * ( t ) is as in (37) and E ˇ N ( · ) is the expectation taken with respect to a probability measure P ˇ N defined by the modified mean-reverting drift f ˇ y ( y * ( t ) , t ) = α ˇ y ( γ ˇ y y * ( t ) ) . By the Feynman–Kac theorem, this expression of u ( t , T ) is the solution, with terminal condition u ( T , T ) = 1 , of the following partial differential equation:
1 2 g y 2 2 u y 2 + f ˇ y u y + u t = y * u ,
which has exactly the same form of the general valuation Equation (9) for the real ZCB price v ¯ . As is shown in Appendix B, (43) is actually the general valuation equation for the discount factor u ( t , T ) in the nominal economy setting, i.e., when the security prices are measured in money units. In this setting, the market price of risk for the p ( t ) and y ( t ) factors are different from zero, and take exactly the form which determines the risk-adjusted drift f ˇ y and the modified process y * ( t ) .3 Then P ˇ N is actually the risk-neutral measure in the nominal economy and E ˇ N ( · ) is the corresponding expectation. These findings derive from the more general priciple that our pricing model can be reformulated in a nominal economy framework provided that the appropriate change of probability measure is made.

2.5. Nominal Interest Rates

By the Fisher relation (28) and expressions (14) and (35), we obtain the following closed-form expression for the nominal discount factor:
v ( t , t + τ ) = A x ( τ ) A y ( τ ) e x ( t ) B x ( τ ) y ( t ) σ p 2 B y ( τ ) .
The instantaneous nominal interest rate r ( t ) is naturally defined as
r ( t ) = lim τ 0 h r ( t , t + τ ) ,
with h r ( t , t + τ ) given in (29). Therefore, by (17) and (37), we obtain the Fisher relation for instantaneous returns:
r ( t ) = x ( t ) + y * ( t ) .
Since v ( t , t + τ ) depends on the two risk factors x ( t ) and y ( t ) , two different risk measures must be used:
v ( t , t + τ ) v ( t , t + τ ) x = B x ( τ ) , v ( t , t + τ ) v ( t , t + τ ) y = B y ( τ ) ,
where B x ( τ ) and B y ( τ ) are the real interest rate measure and the expected inflation risk measure previously defined.
Given the normality properties of x ( T ) | x ( t ) and y ( T ) | y ( t ) , r ( T ) | r ( t ) is also normal, with mean equal to the sum of the means minus σ p 2 and, by the independence assumption, with variance equal to the sum of the variances. Given the normality of the exponent in (44), future nominal log-returns are normally distributed and future nominal ZCB prices are lognormally distributed.

3. Estimation of Model Parameters

3.1. Market Data

In order to estimate the parameters of the three-factor model, we used the same market data, by type and observation period, as in DFM23. This allows the results obtained with the three-factor model to be compared with those obtained with the two-factor model used in that paper.
  • Price index process. As the reference CPI we choose the Harmonised Index of Consumer Prices (HICP), which is compiled by Eurostat and the euro-area national statistical institutes and is published at the end of each month.
    We use the monthly time series { p ( θ ) } θ Θ m , where Θ m is the set of month-end dates in [January 2005, December 2022].
  • Inflation discount factors. To estimate the parameters of the inflation component of the model, we use the Zero-Coupon Inflation-Indexed Swaps (ZCIIS). On each trading day t the ZCIIS rate i M ( t , t + τ ) is quoted on the market for time-to-maturity τ T M : = { 1 , , 10 , 12 , 15 , 20 , 25 } years. As shown in DFM23 (see also Brigo and Mercurio (2007)), the following no-arbitrage relations hold:
    v ( t , t + τ ) 1 + i M ( t , t + τ ) τ = v ¯ ( t , t + τ ) , τ T M .
    Hence, 1 + i M ( t , t + τ ) τ provides the inflation discount factor u M ( t , t + τ ) prevailing on the market at time t for time-to-maturity τ .
    We use the daily time series { i M ( θ , θ + τ ) } τ T M , θ Θ , where Θ is the set of trading days in [1 January 2005, 30 December 2022]. As a usual market practice, swap rates for the missing maturities in T M are computed by linear interpolation of quoted swap rates. Hence, at each date θ Θ the swap rates i M ( θ , θ + τ ) are available for τ T : = { 1 , 2 , , 25 } years. These maturities are usually sufficient for non-life actuarial applications.
  • Nominal discount factors. Since the real discount factors are not directly observable on the market, we take data on the nominal risk-free discount factors, or rates, and then use the Fisher relation.
    In order to calibrate the model on the market at the valuation date t = 30 December 2022, we consider the risk-free nominal discount factors v M ( t , t + τ ) published by EIOPA. We choose the version including volatility adjustment. However, in order to obtain a complete estimate of the model, including the parameters of the natural probability measure, we also take the triple A yield curves computed by the European Central Bank (ECB) on each trading day (for details, see the Technical Notes downloadable from the ECB website European Central Bank (2023)). These curves are obtained by interpolating the observed returns on triple A bonds with different maturities using the Svensson model, a six-parameter smooth function. For the log-returns the Svensson function has the following form:
    h r ( t , t + τ ; π ) = β 0 + β 1 ϑ 1 1 e τ / ϑ 1 τ + β 2 ϑ 1 1 e τ / ϑ 1 τ e τ / ϑ 1 + β 3 ϑ 2 1 e τ / ϑ 2 τ e τ / ϑ 2 ,
    with the parameter set π = ( β 0 , β 1 , β 2 , β 3 , ϑ 1 , ϑ 2 ) (to avoid confusion, we have renamed the ECB parameters τ 1 and τ 2 as ϑ 1 and ϑ 2 , respectively). Observe that lim τ 0 h r ( t , t + τ ; π ) = r ( t ) = β 0 + β 1 . Thus, with this interpolating function there is no problem in determining the initial point of the yield curve. Differently from DFM23, where only the parameters β 0 , β 1 are required, we use the daily time series of the entire parameter sets { π ^ ( θ ) } θ Θ estimated by ECB. The corresponding yield curves are denoted by h r M ( θ , θ + τ ) and the instantaneous interest rates as r M ( θ ) .4
For further details and illustrations, see De Felice and Moriconi (2023). In particular, Figure 1 in that paper shows the yield curves h r M ( t , t + τ ) , h y M ( t , t + τ ) , h x M ( t , t + τ ) as well as the risk-neutral market expectations of the future annual inflation rates:
μ ˇ τ : = log u M ( t , t + τ ) u M ( t , t + τ 1 ) ,
for τ up to 25 years. In our three-factor model, however, the risk-neutral market expectations μ ˇ τ are equal to the natural expectations μ τ of future inflation rates.

3.2. Parameters of the Inflation Component

Our aim in this section is to derive an estimate of the parameters of the CPI distribution, including the instantaneous expected inflation rate, which is not directly observable. In a first step we estimate α y and σ y , then we estimate σ p and ρ , and finally, we provide an estimate of γ y and y ( θ ) for all θ Θ .

3.2.1. Time Series Estimate of α y and σ y

We consider the time series of market inflation log-returns h y M ( θ , θ + τ ) for θ Θ and τ T . As a first step of the estimation procedure, we estimate α y and σ y on the time series of the one-year log-returns h y M ( θ , θ + 1 ) . In principle, a time series of instantaneous inflation rates y M ( θ ) : = lim τ 0 h y M ( θ , θ + τ ) would be needed. However, y M ( θ ) is not directly observable and should be extrapolated by observed returns over finite time periods. Since an appropriate extrapolation function is not available at this stage, we defer to Section 3.2.3 the estimate of a time series of instantaneous inflation rates y M ( θ ) and consider here directly the market log-returns for one-year time-to-maturity. The rationale is that h y M ( θ , θ + 1 ) should be the observed inflation rate for which the correlation with the instantaneous rate y M ( θ ) is highest.5
We use the affinity property of the log-returns. We rewrite expression (35) as
u ( t , t + τ ) = A y ( τ ) e y ( t ) B y ( τ ) ,
with
A y ( τ ) : = e σ p 2 B y ( τ ) A y ( τ ) ,
where A y ( τ ) is given by (36). Then we can write
h y ( t , t + τ ) = C y ( τ ) + y ( t ) B y ( τ ) τ ,
with
C y ( τ ) : = log A y ( τ ) τ .
By Ito’s lemma, we have the following SDE for h y ( t , t + τ ) (where we omit the dependence on t + τ ):
d h y ( t ) = f y B y ( τ ) τ d t + σ y B y ( τ ) τ d W y ( t ) = α y γ y B y ( τ ) τ α y B y ( τ ) τ y ( t ) d t + σ y B y ( τ ) τ d W y ( t ) .
Since, by (48), y ( t ) B y ( τ ) / τ = h y ( t ) C y ( τ ) we obtain
d h y ( t ) = α y γ y B y ( τ ) τ + C y ( τ ) α y h y ( t ) d t + σ y B y ( τ ) τ d W y ( t ) ,
or
d h y ( t ) = a h y ( t ) + b d t + σ h d W y ( t ) ,
with
a : = α y , b : = α y γ y B y ( τ ) τ + C y ( τ ) , σ h : = σ y B y ( τ ) τ .
Observe that, since α y > 0 , then B y ( τ ) / τ < 1 ; hence, σ h < σ y .
As shown in Appendix A.2, for a chosen time-to-maturity τ and a chosen time step Δ t , (51) corresponds to the following first-order autoregressive model:
h y ( θ + Δ t , θ + Δ t + τ ) = β 0 + β 1 h y ( θ , θ + τ ) + e ( θ ) ,
where the error terms e ( θ ) are independent with zero mean and variance Ω 2 . Using the time series of market values h y M ( θ , θ + τ ) , the parameter estimates β ^ 0 , β ^ 1 , and Ω ^ 2 are obtained by linear regression and the corresponding estimate of the mean-reversion coefficient is given by
α ^ y = a ^ = log β ^ 1 Δ t .
Then we obtain
B ^ y ( τ ) = 1 e α ^ y τ α ^ y , σ ^ h 2 = Ω ^ 2 2 α ^ y 1 e 2 α ^ y Δ t , σ ^ y 2 = σ ^ h 2 B ^ y 2 ( τ ) / τ 2 .
This estimation procedure can be applied for any value of τ for which market observations are available. However, as previously noticed, the most reliable results should be that obtained for τ = 1 year. The corresponding estimation procedure (assuming 260 trading days per year, that is, Δ t = 1 / 260 ) provided the following results:
α y = 0.53878 , σ y = 0.01802 .
Remark 2. 
The illustrated estimation procedure is viable because in our model, the inflation market risk premium ϕ y N is independent of y ( t ) , and hence the mean-reversion coefficient α y is the same under the natural and the risk-neutral measure. This approach could not be directly applied to x ( t ) , since the time series estimate allows us to obtain the real-world parameter α x , while B x ( τ ) depends on its risk-adjusted version α ˇ x = α x π .

3.2.2. Time Series Estimate of σ p and ρ

From the monthly time series { p ( θ ) } θ Θ m of HICP values observed at each month-end in Θ m = [January 2005, December 2023], we derived the time series of monthly log-price-ratios ( Δ t = 1 / 12 ):
L ( θ , θ + Δ t ) = log p ( θ + Δ t ) p ( θ ) .
We show how σ p and ρ can be estimated in our model considering this time series and a corresponding time series of stochastic integrals.
Let us consider the two random variables L ( t , t + Δ t ) and Y ( t , t + Δ t ) : = t t + Δ t y ( s ) d s . For chosen Δ t , these variables are normal and, using expression (A22) in Appendix A.4 (with the usual specifications Z = Y and a = α y ), we have
Var t L ( t , t + Δ t ) = σ p 2 Δ t + Cov t L ( t , t + Δ t ) , Y ( t , t + Δ t ) ,
where, by expression (A21),
Cov t L ( t , t + Δ t ) , Y ( t , t + Δ t ) = Var t Y ( t , t + Δ t ) 2 ρ σ y σ p α y B y ( Δ t ) Δ t .
Therefore, σ p 2 can be obtained as
σ p 2 = 1 Δ t Var t L ( t , t + Δ t ) Cov t L ( t , t + Δ t ) , Y ( t , t + Δ t )
and, in turn, ρ is given by
ρ = Cov t L ( t , t + Δ t ) , Y ( t , t + Δ t ) Var t Y ( t , t + Δ t ) 2 ( σ y σ p / α y ) Δ t B y ( Δ t ) .
Since y ( t ) is not observable, we use
H ( t , t + Δ t ) : = t t + Δ t h y ( s , s + τ ) d s
for a given τ (which, for the same reasons illustrated in Section 3.2.1, should be 1 year). By (48), Cov t ( L , H ) = Cov t ( L , Y ) B y ( τ ) / τ and Var t ( H ) = Var t ( Y ) B y 2 ( τ ) / τ 2 ; hence,
σ p 2 = 1 Δ t Var t L ( t , t + Δ t ) τ / B y ( τ ) Cov t ( t , t + Δ t ) , H ( t , t + Δ t )
and
ρ = Cov t L ( t , t + Δ t ) , H ( t , t + Δ t ) τ / B y ( τ ) Var t H ( t , t + Δ t ) τ / B y ( τ ) 2 2 ( σ y σ p / α y ) Δ t B y ( Δ t ) .
The variance Var t ( H ) and the covariance Cov t ( L , H ) can be estimated on a monthly ( Δ t = 1 / 12 ) time series L ( θ , θ + Δ t ) θ Θ m of CPI log-price-ratios and a corresponding monthly time series H ( θ , θ + Δ t ) θ Θ m , where
H ( θ , θ + Δ t ) : = δ k = 0 n 1 h y M ( θ + k δ , θ + k δ + τ ) ,
with δ = Δ t / n , and n represents the number of trading days in each month. The estimates obtained in the previous section were used for α y and σ y .
This estimation procedure provided the following results:
σ ^ p = 0.01731 , ρ ^ = 0.61428 .

3.2.3. Estimate of y ( t ) and γ y on a Time Series of Cross-Sections

Having estimated the parameters α y , σ y , σ p , and ρ , our model provides an appropriate extrapolation procedure that allows us to estimate the instantaneous inflation rate y ( t ) , together with the parameter γ y , by a regression on an observed cross-section of log-returns { h y M ( t , t + τ ) } τ T . Since these returns are affine functions of y ( t ) , by some manipulations, this estimation procedure takes the form of a linear regression.
Using again expression (35), let us consider the following:
log u ( t , t + τ ) = τ h y ( t , t + τ ) = log A y ( τ ; γ y ) σ p 2 B y ( τ ) + y ( t ) B y ( τ ) .
We have
log A y ( τ ) + σ p 2 B y ( τ ) = γ y σ y 2 2 α y 2 B y ( τ ) τ + σ y 2 4 α y B y 2 ( τ ) + σ p 2 B y ( τ ) = B y ( τ ) τ x γ ( τ ) γ y + σ y 2 2 α y 2 B y ( τ ) τ + σ y 2 4 α y B y 2 ( τ ) + σ p 2 B y ( τ ) c 0 ( τ ) .
Then, defining u ( τ ) : = τ h y ( t , t + τ ) c 0 ( τ ) as the dependent variable and
x y ( τ ) : = B y ( τ ) , x γ ( τ ) : = τ B y ( τ )
as the independent variables, we are led to the following two-dimensional linear model (without intercept):
u ( τ ) = y x y ( τ ) + γ y x γ ( τ ) + e ( τ ) ,
with the usual assumptions on the error terms e ( τ ) .
For the estimation at the valuation date t, the data for the dependent variable are the cross-section { u M ( t , t + τ ) : = τ h y M ( t , t + τ ) c ^ 0 ( τ ) } τ T , where c ^ 0 ( τ ) is computed using the parameter estimates previously derived. However, we need not only an estimate at time t = 30 December 2022, but also an estimate for each θ Θ , i.e., at each date of the time series. So, we repeated the linear regression estimate (56) for each cross-section h y M ( θ , θ + τ ) τ T , θ Θ , thus obtaining the time series of estimated values γ ^ y ( θ ) , y ^ ( θ ) θ Θ . These estimates are illustrated in Figure 1.
Remark 3. 
Of course, a γ y parameter different for each time point is clearly inconsistent with the model assumptions. However, this forcing of the assumptions is only used here to obtain an extrapolation function not too inconsistent with the model. On the other hand, our model for the inflation discount factors is too simple to capture the market dynamics well enough to produce constant parameter estimates across any observation date. We point out, however, that, since in Section 3.3.3 we will calibrate the real discount factors v ¯ M ( t , t + τ ) = v M ( t , t + τ ) / u M ( t , t + τ ) on the observed term structure of real interest rates using the Hull and White model, the estimates of γ y ( θ ) will not be used.

3.3. Parameters of the Real Interest Rate Component

In this section, we estimate both the natural and the risk-neutral parameters of real interest rate component of the model, introduced in Section 2.2 and Section 2.3.

3.3.1. Derivation of Real Interest Rate Data by Nominal Interest Rate Data

By the estimated Svensson parameters π ^ published each trading day by ECB, we computed the risk-free nominal log-returns h r M ( θ , θ + τ ) according to (46) and then derived the corresponding real log-returns using the Fisher relation:
h x M ( θ , θ + τ ) = h r M ( θ , θ + τ ) h y M ( θ , θ + τ ) , τ T , θ Θ .
The Fisher relation was also used to derive the corresponding time series of estimated real instantaneous interest rates:
x ^ ( θ ) = r M ( θ ) y ^ ( θ ) σ ^ p 2 , θ Θ ,
where y ^ ( θ ) is estimated in the previous section and σ ^ p in Section 3.2.2. The time series of these estimated instantaneous rates are illustrated in Figure 2.

3.3.2. Time Series Estimate of Natural Real Interest Rate Parameters

In order to estimate the real-world parameters α x , γ x , σ x of the real interest rate process introduced in Section 2.2, we applied to the time series { x ^ ( θ ) } the classical OU estimators just used in Section 3.2.1. This estimation procedure provided the following results:
α ^ x = 1.19764 , γ ^ x = 0.00159 , σ ^ x = 0.02184 .
By the independence assumption, σ ^ r = σ ^ x 2 + σ ^ y 2 = 0.02832 .

3.3.3. Calibration of the Real Interest Rate Model

To simulate at the valuation date t the future real discount factors v ¯ ( t , t + τ ) , t t , we have to calibrate the real interest rate model on time t data. In order to obtain a perfect fit to the term structure derived from the market data, we used the extension of Vasicek’s model provided by Hull and White (Hull and White 1990). The details of this model are given in Appendix C. As shown in the appendix, the Hull–White estimates of future discount factors v ¯ ( t , t + τ ) , for t t , τ 0 , are computed as
v ¯ H W ( t , t + τ ) = v ¯ M ( t , t + τ ) v ¯ M ( t , t ) A x H W ( t , t + τ ) e ω x ( t ) B x ( τ ) ,
where ω x ( t ) is normal with mean zero and variance
Var ω x ( t ) = σ x 2 1 e 2 α x ( t t ) 2 α x ,
and
A x H W ( t , t + τ ) : = e 1 2 V ( τ ) V ( t + τ ) + V ( t ) ,
with
V ( τ ) : = σ x 2 α x ˇ 2 τ + 2 α ˇ x e α ˇ x τ 1 2 α ˇ x e 2 α ˇ x τ 3 2 α ˇ x .
The sensitivity function in (59) has the usual form B x ( τ ) = 1 e α ˇ x τ / α ˇ x .
The parameters σ x and α x have already been estimated (as σ ^ x = 0.02184 and α ^ x = 1.19764 ). It therefore remains to estimate the risk-neutral mean-reversion parameter α ˇ x . In DFM23, where the same problem arises for the calibration of nominal interest rates, this was accomplished by a best-fitting procedure of the model yield function h r ( t , t + τ ) to the current market term structure h r M ( t , t + τ ) , τ T . This procedure applied to our real interest rate data has proven to be scarcely reliable, producing rather unstable estimates of α ˇ x . We therefore found it convenient to perform the estimation using the entire time series of cross-sections h r M ( θ , θ + τ ) τ T for θ Θ . Using the affinity property of v ¯ ( t , t + τ ) as we did in Section 3.2.1 for the inflation discount factor u ( t , t + τ ) , we are led to the following SDE for the real log-returns, which is the analogue of (51):
d h x ( t , t + τ ) = a h x ( t , t + τ ) + b d t + σ x B x ( τ ) τ d W x ( t ) , τ 0 .
The coefficients a and b, which have the same form as in (52), are known at time t. Thus, we have
Var t d h x ( t , t + τ ) = σ x 2 B x ( τ ) τ 2 d t , τ 0 .
These variances are independent of t. For a sufficiently small finite time interval Δ t , we can consider the following approximations:
σ Δ h x 2 ( τ ) : = Var t Δ h x ( t , t + τ ) σ x 2 B x ( τ ) τ 2 Δ t , τ 0 .
We derived a time series estimate σ ^ Δ h x 2 ( τ ) of these variables for τ T by taking the sample variance of the daily increments:
Δ h r M ( θ , θ + τ ) : = h r M ( θ + Δ t , θ + Δ t + τ ) h r M ( θ , θ + τ ) θ Θ , Δ t = 1 / 260 ,
and we thus obtained an estimate of α ˇ x by solving the ordinary-least-squares problem:
α ˇ ^ x = argmin α ˇ x τ T σ x B x ( τ ; α ˇ x ) τ Δ t σ ^ Δ h x ( τ ) 2 .
We obtained α ˇ ^ x = 0.33109 .

4. Applying the Market Model to Claims Reserving

4.1. The Essential Elements of the Claims Reserving Application

For the claims reserving applications of the model, we adopt the definitions, theoretical framework, and simulation schemes proposed in DFM23. However, we do not consider the “actuarial approach”, which is characterized by the exclusion of interest rate risk. As also noted in DFM3, if we accept this approximation, it would not make sense to use a stochastic model for interest rates, much less one with stochastic expected inflation. Therefore, here we will only consider the so-called “market approach” and remove the superscript M used in the DFM23 notation when it is necessary to distinguish this approach from the actuarial approach (denoted there by the superscript A). We recall here for convenience the essential definitions introduced in DFM23.
The chain-ladder reserve at the current date t = I , the end of year 2022, is given by
R I c c = k = 1 n 1 C I D I + k c c v ¯ M ( I , I + k ) ,
where
  • the superscript cc indicates that the chain-ladder algorithm is applied to a triangle of claim payments expressed at current costs, i.e., adjusted for the past inflation as measured by the HICP index (the CPI to which the ZCIIS are linked);
  • C I D I + k c c denotes the chain-ladder estimate at time I of the sum D I + k c c of the incremental payments (at time I costs) to be performed in the k-th future calendar year.
Of course, v ¯ M ( I , I + k ) are the current market real discount factors, computed as v M ( I , I + k ) / u M ( I , I + k ) , where v M ( I , I + k ) are the risk-free discount factors provided at time I by EIOPA (including volatility adjustment) and u M ( I , I + k ) are the inflation discount factors derived by the ZCIIS observed at time I.
As regards reserve risk under the ultimate point of view, the predictive distribution of the discounted ultimate obligations (DUOs) provided by S simulations is given by
DUO = s U ˜ I c c = k = 1 n 1 s D ˜ I + k c c v ¯ M ( I , I + k ) ; s = 1 , , S .
It is assumed that the claim payments s D ˜ I + k c c are provided by the bootstrap simulation of a chosen stochastic chain-ladder model. Under the short-term view of reserve risk, as prescribed by Solvency II, the predictive distribution of the year-end obligations (YEOs) is given by
YEO = s Y ˜ I + 1 c c : = s D ˜ I + 1 c c s p ˜ ( I + 1 ) p ( I ) + s R ˜ I + 1 c c ; s = 1 , , S ,
where the residual reserve at time I + 1 in the s-th simulation is given by
s R ˜ I + 1 c c : = k = 2 n 1 s C ˜ I + 1 D I + k c c s p ˜ ( I + 1 ) p ( I ) s v ¯ ˜ M ( I + 1 , I + k ) .
Here, s v ¯ ˜ M ( I + 1 , I + k ) are the real discount factors at time I + 1 simulated consistently with the current term structure v ¯ M ( I , I + k ) using the Hull and White Formula (59) and the price index s p ˜ ( I + 1 ) at time I + 1 is simulated as suggested in Section 2.4.4 using expression (39). The inflation discount factors u M ( t , t + τ ) derived by the ZCIIS quoted in the market represent, by definition, risk-neutral expectations. However, in the three-factor model, we have u M ( t , t + τ ) = u ( t , t + τ ) , since the risk premia ϕ p and ϕ y are zero. So, in effect, expression (39) is written as
p ˜ ( T ) p ( t ) = 1 u M ( t , T ) e 1 2 σ L 2 ( t , T ) + σ L ( t , T ) ε ˜ .
Remark 4. 
In the two-factor model (see Section 4.2 in DFM23), the expression corresponding to (39) is
p ˜ ( T ) p ( t ) = 1 u ( t , T ) e 1 2 σ L 2 ( t , T ) + σ L ( t , T ) ε ˜ .
But in that model, the natural and the risk-neutral measures are different and the relation holds as u ( t , T ) = u M ( t , T ) e σ L 2 ( t , T ) . Thus, we obtain a relation formally identical to (62). This means that in the two models, the market data u M ( t , T ) are used in the same way. Of course, as already pointed out, the expression of σ L 2 ( t , T ) is different in the two models. However, it is immediately verified that if in the three-factor model we set σ y = 0 , i.e., if we make expected inflation deterministic, we obtain σ L 2 ( t , T ) = σ p 2 ( T t ) , as in the two-factor model.

4.2. Numerical Results with the Market Approach

For the numerical calculations, we used the same triangle of paid losses used in DFM23, which contains the MTPL claim payments observed at the end of 2022 (time I) over the previous 22 accident years (so the data also include the period of the COVID-19 pandemic). As in that paper, the stochastic chain-ladder model used is the ODP reserving model introduced by Renshaw and Verral (1998), in the bootstrap version proposed by England and Verral (2002). Since the definitions of the quantities to be estimated, the market data and the paid losses data, are the same in the two articles, the differences in the results obtained can only be attributed to the difference in the market model used and the estimates of the relevant parameters.

4.2.1. Recall of Results with the Two-Factor Model

In DFM23, the results of the market approach are as follows (see the paper for further details):
Results with deterministic computations
-  Undiscounted reserve with implicit inflation100,000
-  Discounted reserve with implicit inflation: R I h c 91,091
-  Undiscounted reserve with triangle at current costs:100,087
-  Discounted reserve with triangle at current costs: R I c c | n o r e p 91,391
-  Undiscounted reserve with modeled inflation109,445
-  Discounted reserve with modeled inflation (BE): R I c c 99,402
Results with the two-factor model (100,000 simulations)
-  DUO sample mean: M U I c c 99,304
-  DUO Std5101
-  DUO CV5.14%
-  YEO sample mean: M Y I + 1 c c 102,809
-  YEO Std4773
-  YEO CV4.64%
-   99.5 % percentile of YEO 115,880
-  Solvency Capital Requirement13,071
-  Present value of YEO sample mean: M Y I + 1 c c v ( I , I + 1 ) 99,461

4.2.2. Using the Three-Factor Model

Using the three-factor model, the deterministic results are obviously the same. The results on reserve risk under the ultimate point of view also do not change, as the discount factors used in expression (60) are those observed on the market. For the reserve risk under the one-year view, the results are different, however.
Results with the three-factor model (100,000 simulations)
-  YEO sample mean: M Y I + 1 c c 102,762
-  YEO Std4950
-  YEO CV4.82%
-   99.5 % percentile of YEO 116,314
-  Solvency Capital Requirement 13,551
-  Present value of YEO sample mean: M Y I + 1 c c v ( I , I + 1 ) 99,416
While the YEO sample mean M Y I + 1 c c and its discounted value are the same, within the Monte Carlo error, as in the two-factor model, the introduction of stochastic expected inflation in the valuation model produces a non-negligible increase in the variability of the predictive loss distribution. The coefficient of variation increases from 4.64 % to 4.82 % , and the SCR, as a percentage of the BE, increases by about 50 basis points, from 13.14 % to 13.63 % . These increases are roughly similar to those found in DFM23 when passing from the actuarial approach to the market approach with the two-factor model.
With regard to the kurtosis of the predictive distribution and the relationships with the corresponding moment-fitted lognormal distribution, there are no relevant differences to what was found with the two-factor model.

4.2.3. The Determinants of the Increase in Variability

It may be of some interest to analyze the determinants of the increase in variability produced by the model with stochastic expected inflation. The variance of the empirical distribution s R ˜ I + 1 c c defined in (61) is determined by the variance of the one-year repricing factor p ( I + 1 ) / p ( I ) and the year-end real discount factors v ¯ M ( I + 1 , I + k ) , k = 2 , 3 , , n 1 . Expression (61) applies to both the tree-factor model and the two-factor model in DFM23, the only difference being in the way these factors are calculated. Let us introduce the shorthand notation p 1 for p ( I + 1 ) / p ( I ) , v ¯ 1 k for v ¯ ( I + 1 , I + k ) , and p 1 v ¯ 1 k for the product p ( I + 1 ) / p ( I ) v ¯ ( I + 1 , I + k ) , and let us denote by σ p 1 , σ v ¯ 1 k , and σ p 1 v ¯ 1 k the corresponding standard deviations. In the three-factor model, σ p 1 depends on the parameters α y , σ y , σ p , ρ through the relation (39) and the expression (32) of σ L 2 ( I , I + 1 ) , and σ v ¯ 1 k depends on the parameters α x , α ˇ x , σ x entering the Hull–White expression (59) of v ¯ ( I + 1 , I + k ) . In particular, α x and σ x determine the variance of ω ( I + 1 ) . In the two-factor model, σ p 1 depends on the parameter σ L 2 ( I , I + 1 ) = σ p and σ v ¯ 1 k depends on the parameters α r , α ˇ r , σ r in the Hull–White expression of the nominal discount factors v ( I + 1 , I + k ) analogous to (59) (expression (32) in DFM23). For these parameters, the estimates derived in DFM23 are
σ ^ L 2 ( I , I + 1 ) = σ ^ p = 0.01808 , σ ^ r = 0.01552 , α ^ r = 0.61346 , α ˇ ^ r = 0.17433 ,
while with the three-factor model, we derived in the previous sections
σ ^ L ( I , I + 1 ) = 0.02313 , σ ^ x = 0.02184 , α ^ x = 1.19764 , α ˇ ^ x = 0.33109 .
Thus, volatilities are higher in the three-factor model than in the two-factor model, but so are the mean-reversion parameters, both natural and risk-neutral. Since, as we have seen, mean-reversion produces a decrease in variability for increasing maturities, it is useful to compare the combined effect of these parameter values. This can be easily carried out, since the standard deviations of the factors involved in (61) can be computed in closed form, using the relation Std ( e X ) = E ( e X ) e σ X 2 1 , valid for X N ( μ X , σ X 2 ) . The results of the computations are reported in Figure 3, where the standard deviations for different maturities are illustrated in blue color for the three-factor model and in red for the two-factor model. It can be seen that, although the standard deviation of ω ( I + 1 ) is larger in the three-factor than in the two-factor model, the standard deviations of the discount factors in the three-factor model are smaller even before the 2-year maturity. The change in dominance is even more evident when looking at the overall standard deviations, indicated by continuous lines, where the crossing of the curves occurs for a maturity between 3 and 4 years. Since by far the largest portion of the year-end obligations have maturity of less than 3 years, the variability of the overall present value is greater for the blue curve than for the red curve. That is, we have an increase in SCR with the three-factor model.

5. Conclusions

When simulating the distribution of year-end obligations, we find with our data that the three-factor model produces higher volatility than the two-factor model for both the repricing factor of the CPI and the next-year real discount factors to be applied to the claim payments of the remaining calendar years. However, the mean-reversion, both natural and risk-neutral, of real interest rates in the three-factor model is also higher than that derived from nominal interest rates in the two-factor model. This may cause the overall volatility of the three-factor model to be lower over long maturities. In our case, the overall effect is that, with the market approach of the three-factor model, the reserve SCR (under the one-year view) has a non-negligible increase over that derived with the market approach of the two-factor model. This increase is grosso modo similar to that found in DFM23 when switching from the actuarial to the market approach in the two-factor model.

Funding

This research received no external funding.

Acknowledgments

The author would like to thank Luca Passalacqua for the useful discussions on the development and the estimation of the three-factor model presented in this paper.

Conflicts of Interest

The author declares no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
BEBest Estimate
DFM23De Felice and Moriconi (2023) paper
HICPHarmonised Index of Consumer Prices
SCRSolvency Capital Requirement
SDEStochastic Differential Equation
ZCBZero-Coupon Bond
ZCIISZero-Coupon Inflation-Indexed Swap

Appendix A. Some Fundamental Results on Ornstein–Uhlembeck Process

Appendix A.1. Integral Representation of the Ornstein–Uhlembeck Process

Let z ( t ) be the OU process described by the SDE:
d z ( t ) = ( a z ( t ) + b ) d t + σ z d W z ( t ) , a , b R , σ > 0 .
As is well-known, this is a shorthand expression for the following integral representation:
z ( T ) = z ( t ) + t T ( a z ( s ) + b ) d s + σ z t T d W z ( s ) , T t
(where the stochastic integral is defined in Ito’s sense). We prove the following fundamental result concerning the conditional random variable z ( T ) | z ( t ) .
Result A1. 
For T = t + τ > t , we have
z ( T ) | z ( t ) = μ z ( t , T ) + σ z ε z ( t , T ) ,
with
ε z ( t , T ) : = e a T t T e a s d W z ( s ) N ( 0 , v ε ( τ ) ) .
Therefore,
z ( T ) | z ( t ) N μ z ( t , T ) , σ z 2 v ε ( τ ) ,
where the conditional mean is
μ z ( t , T ) : = E ( z ( T ) | z ( t ) ) = z ( t ) e a τ + b β ( τ ) ,
with
β ( τ ) : = e a τ 1 a ,
hence depending only on z ( t ) and τ , and the conditional variance is
v z ( t , T ) : = Var z ( T ) | z ( t ) = σ z 2 e 2 a τ 1 2 a ,
which depends only on τ .
Proof. 
We first apply a transformation providing an SDE with constant drift. Let u ( z ( t ) , t ) = z ( t ) e a t . Differentiating:
d u ( z , t ) = d z ( t ) e a t z ( t ) a e a t d t = ( a z ( t ) + b ) e a t d t + σ z e a t d W z ( t ) z ( t ) a e a t d t = b e a t d t + σ z d W z ( t ) e a t .
Then u ( t ) is described by the following SDE:
d u ( t ) = b e a t d t + σ z e a t d W z ( t ) .
For T = t + τ t ,
u ( T ) = u ( t ) + t T d u ( s ) d s = u ( t ) + b e a t e a T a + σ z t T e a s d W z ( s ) .
Replacing u ( t ) with z ( t ) e a t , we then have
z ( T ) = z ( t ) e a τ + b e a τ 1 a + σ z t T e a ( T s ) d W z ( s ) .
This gives expression (A2), with ε z ( t , T ) given by (A3) and μ z ( t , T ) given by (A4). Expression (A5) is immediately obtained by the following property:
Var t t T f ( x ) d W ( x ) = t T f 2 ( x ) d x ,
which gives
v ε ( τ ) : = Var t ε z ( t , T ) = e 2 a T t T e 2 a s d s = e 2 a τ 1 2 a .

Appendix A.2. Parameter Estimate for the Ornstein–Uhlembeck Process

Taking T t = Δ t and interpreting expression (A6) as an iterative relation, we see that the discrete equivalent of the SDE (A1) is the first-order autoregressive equation:
z t + Δ t = β 0 + β 1 z t + ε t ,
with
β 0 : = b e a Δ t 1 a ,
β 1 : = e a Δ t ,
and where the error terms ε t are independent normal with zero mean and constant variance
Var ( ε t ) : = Ω 2 = σ 2 e 2 a Δ t 1 2 a .
If a time series of observations of z ( t ) with time step Δ t is available, the standard linear regression defined by (A8) provides the corresponding consistent estimates β ^ 0 , β ^ 1 and Ω ^ 2 . By (A9)–(A11), we obtain
a ^ = log β ^ 1 Δ t , b ^ = β ^ 0 a ^ e a ^ Δ t 1 , σ ^ 2 = Ω ^ 2 2 a ^ e 2 a ^ Δ t 1 .
If Equation (4) is expressed in the form
d z ( t ) = α ( γ z ( t ) ) d t + σ d W ( t ) ,
the estimates for α and γ are obtained as
α ^ = a ^ , γ ^ = b ^ / α ^ .
If time is measured in years and we have a daily (weekly, monthly) time series, then we pose Δ t = 1 / 260 ( Δ t = 1 / 52 , 1 / 12 , resp.). Consequently, all the estimates will be consistently expressed on annual basis.

Appendix A.3. Integral of the Ornstein–Uhlembeck Process

For T = t + τ > t , let us consider the stochastic integral:
Z ( t , T ) = t T z ( s ) d s ,
with z ( t ) described by (A1). The following fundamental result holds for the conditional random variable Z ( T ) | z ( t ) .
Result A2. 
For T = t + τ > t we have
Z ( t , T ) | z ( t ) = μ Z ( t , T ) + σ z a ε z ( t , T ) δ z ( t , T ) ,
with ε z ( t , T ) given by (A3) and
δ z ( t , T ) : = t T d W z ( s ) N ( 0 , τ ) .
Hence,
Z ( t , T ) | z ( t ) N μ Z ( t , T ) , v Z ( t , T ) ,
with conditional mean
μ Z ( t , T ) : = E Z ( t , T ) | z ( t ) = z ( t ) β ( τ ) + b a β ( τ ) τ ,
which depends only on z ( t ) and τ , and conditional variance
v Z ( t , T ) : = Var Z ( t , T ) | z ( t ) = σ 2 2 a β 2 ( τ ) σ 2 a 2 β ( τ ) τ ,
which depends only on τ .
Proof. 
The normality of Z ( T ) is an obvious consequence of the normality of z ( T ) . Expression (A14) is obtained by integrating μ z ( t , T ) in (A4). The stochastic term in (A13) is obtained by integrating ε z ( t , s ) in (A3). In fact, we have (using Fubini’s theorem)
σ z t T t s e a ( s u ) d W z ( u ) d s = σ z t T u T e a ( s u ) d s d W z ( u ) = σ z a t T e a ( T u ) 1 d W z ( u ) = σ z a ε z ( t , T ) δ z ( t , T ) .
Expression (A15) for the variance can be obtained by computing the following integral:
Var t Z ( t , T ) = σ z 2 a 2 t T e a ( T u ) 1 2 d u .
More interestingly, the same expression is obtained by observing that
Cov t ε z ( t , T ) , δ z ( t , T ) = β ( τ ) .
In fact, since W z ( t ) is a zero-mean i.i.d. process,
E t t T e a ( T s ) d W z ( s ) · t T d W z ( s ) = E t t T e a ( T s ) d W z 2 ( s ) = t T e a ( T s ) d s = e a τ 1 a = β ( τ ) .
Then we have
Var t Z ( t , T ) = σ z 2 a 2 Var t ε z ( t , T ) δ z ( t , T ) = σ z 2 a 2 Var t ε z ( τ ) + Var t δ z ( t , T ) 2 Cov t ε z ( t , T ) , δ z ( t , T ) = σ z 2 a 2 v ε ( τ ) + τ 2 β ( τ ) ,
which gives (A15), observing that (A7) can be written as v ε ( τ ) = 1 2 a β 2 ( τ ) + β ( τ ) . □

Appendix A.4. Log-Returns of the Price Index Process

Let us consider the process p ( t ) defined by the following SDE:
d p ( t ) = z ( t ) p ( t ) d t + σ p p ( t ) d W p ( t ) ,
with z ( t ) given by (A1) and Cov d W z , d W p = ρ d t . Let us define the “log-return”, i.e., the logarithm of the price ratio, on the time period [ t , T ] :
L ( t , T ) = log p ( T ) p ( t ) .
We prove the following fundamental result for the random variable L ( t , T ) | z ( t ) , p ( t ) .
Result A3. 
For T = t + τ > t , we have
L ( t , T ) | z ( t ) , p ( t ) = μ Z ( t , T ) 1 2 σ p 2 τ + σ z a ε z ( t , T ) δ z ( t , T ) + 1 ρ 2 σ p δ p ( t , T ) + ρ σ p δ z ( t , T ) ,
with
ε z ( t , T ) : = e a T t T e a s d W z ( s ) , δ z ( t , T ) : = t T d W z ( s ) , δ p ( t , T ) : = t T d W p ( s ) ,
and W p ( t ) , W z ( t ) independent. Then
L ( t , T ) | z ( t ) , p ( t ) N μ L ( t , T ) , v L ( t , T ) ,
with conditional mean
μ L ( t , T ) : = E L ( t , T ) | z ( t ) , p ( t ) = z ( t ) β ( τ ) + b a β ( τ ) τ 1 2 σ p 2 τ ,
which depends only on z ( t ) and τ , and conditional variance
v L ( t , T ) : = Var L ( t , T ) | z ( t ) , p ( t ) = σ p 2 τ + σ z 2 2 a β 2 ( τ ) σ z 2 a 2 β ( τ ) τ + 2 ρ σ z σ p a β z ( τ ) τ ,
which depends only on τ .
Proof. 
We rewrite the stochastic dynamics of z ( t ) , p ( t ) using a Cholesky decomposition of the bivariate Wiener process W z ( t ) , W p ( t ) . This gives
d z ( t ) = ( a y ( t ) + b ) d t + σ z d W z ( t ) , d p ( t ) = z ( t ) p ( t ) d t + 1 ρ 2 σ p p ( t ) d W p ( t ) + ρ σ z d W z ( t ) ,
where W z ( t ) and W p ( t ) are independent Wiener processes. By the second SDE in (A20), we have, applying Ito’s lemma,
d log p ( t ) = z ( t ) 1 2 σ p 2 d t + 1 ρ 2 σ p d W p ( t ) + ρ σ p d W z ( t ) .
The solution of this SDE for T = t + τ t , given p ( t ) , is
log p ( T ) = log p ( t ) + t T z ( s ) d s 1 2 σ p 2 τ + 1 ρ 2 σ p t T d W p ( s ) + ρ σ p t T d W z ( s ) .
For the log-price-ratio L ( t , T ) , we then have
L ( t , T ) = Z ( t , T ) 1 2 σ p 2 τ + 1 ρ 2 σ p t T d W p ( s ) + ρ σ p t T d W z ( s ) ,
which is obviously normal. Expression (A18) is immediately obtained, since E ( L ( t , T ) | z ( t ) , p ( t ) ) = μ Z ( t , T ) 1 2 σ p 2 τ , with μ Z ( t , T ) given by (A14). Expression (A19) is easily proven considering that, by (A13),
Var t Z ( t , T ) = Var t σ z a ε z ( t , T ) δ z ( t , T ) .
Then, by the independence of W z and W p :
v L ( t , T ) : = Var t σ z a ε z + ρ σ p σ z a δ z + ( 1 ρ 2 ) σ p 2 Var t δ p .
That is,
v L ( t , T ) = σ z 2 a 2 Var t ( ε z ) + ρ σ p σ z a 2 Var t ( δ z ) + 2 σ z a ρ σ p σ z a Cov ( ε z , δ z ) + ( 1 ρ 2 ) σ p 2 Var t δ p = σ z 2 a 2 v ε ( τ ) + ρ σ p σ z a 2 τ + 2 σ z a ρ σ p σ z a β ( τ ) + ( 1 ρ 2 ) σ p 2 τ ,
where we used (A16). Since v ε ( τ ) = 1 2 a β 2 ( τ ) + β ( τ ) , as in (A7), expression (A19) is obtained. □
A useful alternative expression for v L ( t , T ) is obtained by computing the covariance between L ( t , T ) and Z ( t , T ) . By the representation (A17), we have
Cov t ( L , Z ) = Var t ( Z ) + 2 ρ σ z σ p a Cov t ε z δ z , δ z = Var t ( Z ) + 2 ρ σ z σ p a Cov t ε z , δ z Var t ( δ z ) = v Z ( t , T ) + 2 ρ σ z σ p a β z ( τ ) τ .
Using (A15) and comparing with (A19), we then obtain
Var t L ( t , T ) = σ p 2 τ + Cov t L ( t , T ) , Z ( t , T ) .

Appendix B. The General Valuation Equation in the Nominal Economy

Let us consider at time t a real ZCB which provides at time T t a real payoff X ¯ T which is defined as a specified function contractually specified of x ( T ) , y ( T ) and p ( T ) . Then the price V ¯ ( t ) of this bond is a function V ¯ ( t , x ( t ) , y ( t ) , p ( t ) ) of all three risk factors. By the hedging argument, V ¯ ( t ) is obtained as the solution of the following general valuation equation:
1 2 g x 2 2 V ¯ x 2 + 1 2 g y 2 2 V ¯ y 2 + 1 2 g p 2 2 V ¯ p 2 + ρ g y g p 2 V ¯ y p + ( f x ϕ x ) V ¯ x + ( f y ϕ y ) V ¯ y + ( f p ϕ p ) V ¯ p + V ¯ t = x V ¯ ,
with terminal condition V ¯ ( T ) = X ¯ T . By the Feynman–Kac theorem, this solution is given by
V ¯ ( t ) = E ˇ t R e t T x ( s ) d s X ¯ T ,
where E ˇ t R is taken with respect to the three-dimensional risk-neutral distribution P ˇ R . The valuation Equation (9) and the corresponding solution for V ¯ = v ¯ , given by (14) after ϕ x has been specified as in (12), is immediately obtained as a particular case of (A23) if one requires X ¯ T = 1 ¯ , since in this case, V ¯ ( t ) only depends on x ( t ) .
If we skip to the corresponding nominal economy setting, the general valuation equation for the nominal price V ( t ) is derived by Equation (A23) by posing V ¯ = V / p and computing the corresponding derivatives. Making the calculations and rearranging, we obtain
1 2 g x 2 2 V x 2 + 1 2 g y 2 2 V y 2 + 1 2 g p 2 2 V p 2 + ρ g y g p 2 V y p + ( f x ϕ x ) V x + ( f y ϕ y ρ g y g p p ) V y + ( f p ϕ p g p 2 p ) V p + V t = x + f p ϕ p p g p 2 p 2 V ,
or
1 2 g x 2 2 V x 2 + 1 2 g y 2 2 V y 2 + 1 2 g p 2 2 V p 2 + ρ g y g p 2 V y p + ( f x ϕ x ) V x + ( f y ϕ y N ) V y + ( f p ϕ p N ) V p + V t = x + f p ϕ p p g p 2 p 2 V ,
with
ϕ y N : = ϕ y + ρ g y g p p = ρ g y g p p , ϕ p N = ϕ p + g p 2 p = g p 2 p .
These expressions for the nominal risk premia for y and p can also be obtained by computing the corresponding covariance w.r.t. the nominal wealth W * : = p W ¯ * . This is carried out in a footnote in Section 2.4.5.
Let us specify the price V ( t ) as the price v ( t , T ) of the nominal unit ZCB maturing at time T. By the Fisher relation we know that v ( t , T ) = v ¯ ( t , T ) u ( t , T ) , where v ¯ ( t ) depends only on x. By the independence assumption, u ( t ) is independent of x. Moreover, as shown in Section 2.4.3, u ( t ) is also independent of p. Hence, we have v t , x ( t ) , y ( t ) = v ¯ t , x ( t ) u t , y ( t ) . As a consequence of this separation property, the nominal valuation Equation (A24) provides two separated equations, one for v ¯ , corresponding again to Equation (9), and the other for u, having the following form:
1 2 g y 2 2 u y 2 + ( f y ϕ y N ) u y + u t = f p ϕ p N p g p 2 p 2 u .
Using our specifications for f y , g y and ϕ y N , we have
1 2 σ y 2 2 u y 2 + ( α y γ y ρ σ y σ p ) α y y u y + u t = y * u ,
with y * : = y σ p 2 . Since deriving with respect to y * is the same as deriving with respect to y, this equation can be written as
1 2 σ y 2 2 u y 2 + ( α y γ y α y σ p 2 ρ σ y σ p ) α y y * u y + u t = y * u ,
or
1 2 σ y 2 2 u y 2 + α y γ ˇ y α y y * u y + u t = y * u ,
with γ ˇ y = γ y σ p 2 ρ σ y σ p / α y . The solution of this equation, under terminal condition u ( T , T ) = 1 , is given by u ( t , T ) = E ˇ N e t T y * ( s ) d s , which gives (42).

Appendix C. Applying the Hull and White Model

This section largely mimics Appendix B in DFM23. In order to have a perfect fit to the market discount factors at the valuation date, we consider the Hull and White extension of the Vasicek model for the z ( t ) process (Hull and White 1990). As usual, z ( t ) can be interpreted as x ( t ) or y ( t ) in the main text of the paper. Correspondingly, the discount factor v which will be used in this appendix can be interpreted as v ¯ or u, respectively. In Section 3.3.3, the first interpretation is used.
In the Hull and White model, the long-term rate γ ˇ is assumed to be time-dependent and z ( t ) is expressed by
z ( t ) = ω ( t ) + φ ( t ) , z ( 0 ) = z 0 ,
where ω ( t ) is the following OU process:
d ω ( t ) = α ω ( t ) d t + σ d W ( t ) , ω ( 0 ) = 0
and φ ( t ) is a deterministic function with φ ( 0 ) = z 0 . Time 0 is the reference date when z and φ are equal and it is also the date where we observe on the market the discount factors v M ( 0 , τ ) for τ 0 . For 0 t T = t + τ our aim is to derive an expression v 0 H W ( t , t + τ ) which for t = 0 perfectly fits the observed term structure v M ( 0 , τ ) .
As a trivial extension of Result A1 in Appendix A.1, z ( T ) | z ( t ) and ω ( T ) | ω ( t ) are normal, with conditional moments:
E ( ω ( T ) | ω ( t ) ) = ω ( t ) e α τ , E ( z ( T ) | z ( t ) ) = z ( t ) e α τ + φ ( t ) ,
and
Var ω ( T ) | ω ( t ) = Var z ( T ) | z ( t ) = σ 2 1 e 2 α τ 2 α .
In this model, we have
v 0 H W ( t , T ) = E ˇ t e t T z ( s ) d s = e t T φ ( s ) d s E ˇ t e t T ω ( s ) d s ,
where the expectation is taken w.r.t. the appropriate risk-neutral measure. By the usual calculations, we obtain
t T ω ( s ) d s N ω ( t ) B ( τ ) , V ( τ ) ,
with
V ( τ ) = σ 2 α ˇ 2 τ + 2 α ˇ e α ˇ τ 1 2 α ˇ e 2 α ˇ τ 3 2 α ˇ
and
B ( τ ) = 1 e α ˇ τ α ˇ .
Hence,
E ˇ t e t T ω ( s ) d s = e 1 2 V ( τ ) e ω ( t ) B ( τ ) .
We also obtain
exp t T φ ( s ) d s = v M ( 0 , T ) v M ( 0 , t ) exp 1 2 V ( T ) V ( t ) .
Therefore,
v 0 H W ( t , T ) = v M ( 0 , T ) v M ( 0 , t ) A H W ( t , T ) e ω ( t ) B ( τ ) ,
with
A H W ( t , T ) = e 1 2 V ( τ ) V ( T ) + V ( t ) .
For t = 0 we have V ( 0 ) = 0 , V ( T ) = V ( τ ) and v M ( 0 , 0 ) = 1 ; therefore, since ω ( 0 ) = 0 , (A25) provides v 0 H W ( 0 , T ) = v M ( 0 , T ) , i.e., perfect fitting, as desired.
At time 0, the future discount factor v 0 H W ( t , T ) , with T = t + τ t 0 , is lognormal, with moments (under the natural distribution):
E 0 log v 0 H W ( t , t + τ ) = log v M ( 0 , t + τ ) v M ( 0 , t ) + 1 2 V ( τ ) V ( t + τ ) + V ( t ) ,
since E 0 ω ( t ) = 0 , and:
Var log v 0 H W ( t , t + τ ) = σ 2 1 e 2 α t 2 α B 2 ( τ ) .
Therefore computation of expectation only requires the risk-adjusted parameter α ˇ , while for the computation of variance both α and α ˇ are required.
It can be shown that if these moments are computed under the risk-neutral measure the corresponding expressions fulfil the theoretical property E ˇ 0 [ v ( t , t + τ ) ] < v ( 0 , t , t + τ ) .

Notes

1
To simplify notation we shall denote by E t G ( t , T ) , where G ( t , T ) is a random variable depending on the sample path in ( t , T ] of a given stochastic process x ( t ) , the conditional expectation of G ( t , T ) given F t , where F t is the σ -field generated by x ( t ) up to time t.
2
Actually, in this risk-neutral setting we may not have real mean-reversion, because α x π could be negative.
3
This is immediately seen by observing that in the nominal economy the optimally allocated wealth is measured in money unit and is then given by W * ( t ) = p ( t ) W ¯ * ( t ) . Since the market risk premia are given by the covariance of changes in the relevant risk factor with percentage changes in W * ( t ) , we obtain:
ϕ p N ( p ( t ) , t ) = Cov d p ( t ) W ¯ * ( t ) p ( t ) W ¯ * ( t ) d t , d p ( t ) = σ p 2 p ,
ϕ y N ( y ( t ) , t ) = Cov d p ( t ) W ¯ * ( t ) p ( t ) W ¯ * ( t ) d t , d y ( t ) = ρ σ y σ p ,
where ϕ p N and ϕ y N are the market risk premia for p and y, respectively, in the nominal economy setting.
4
The ECB risk-free yield curves are slightly different from the yield curves provided by EIOPA, since EIOPA yield curves are primarily derived from observed interest rate swaps (after appropriate adjustment for credit risk), inter/extrapolated by the Smith-Wilson method. However, these differences have little relevance to our estimation problem.
5
As it is well-known, in a one-factor model rates of return of different maturities are perfectly correlated. This property is clearly unrealistic, but could be considered acceptable for maturities that are very close to each other.

References

  1. Brigo, Damiano, and Fabio Mercurio. 2007. Interest Rate Models: Theory and Practice. With Smile, Inflation and Credit, 2nd ed. New York: Springer. [Google Scholar]
  2. Cox, John C., Jonathan E. Ingersoll, and Stephen A. Ross. 1985. A theory of the term structure of interest rates. Econometrica 53: 385–406. [Google Scholar] [CrossRef]
  3. De Felice, Massimo, and Franco Moriconi. 2023. Stochastic Chain-Ladder Reserving with Modeled General Inflation. July 21. Available online: https://ssrn.com/abstract=4517501 (accessed on 5 October 2023). [CrossRef]
  4. England, Peter D., and Richard J. Verrall. 2002. Stochastic Claims Reserving in General Insurance (with discussion). British Actuarial Journal 8: 443–544. [Google Scholar] [CrossRef]
  5. European Central Bank. 2023. Technical Notes on ECB Yield Curve Methodology. Available online: https://www.ecb.europa.eu/stats/financial_markets_and_interest_rates/euro_area_yield_curves/shared/pdf/technical_notes.pdf (accessed on 5 October 2023).
  6. Hull, John, and Alan White. 1990. Pricing Interest Rate Derivative Securities. Review of Financial Studies 3: 573–92. [Google Scholar] [CrossRef]
  7. Jarrow, Robert, and Yildiray Yildirim. 2003. Pricing Treasury Inflation Protected Securities and Related Derivatives using an HJM Model. Journal of Financial and Quantitative Analysis 38: 337–58. [Google Scholar] [CrossRef]
  8. Moriconi, Franco. 1995. Analyzing default-free bond markets by diffusion models. In Financial Risk in Insurance. Edited by Giuseppe Ottaviani. Berlin: Springer. [Google Scholar]
  9. Renshaw, Arthur E., and Richard J. Verrall. 1998. A stochastic model underlying the chain ladder technique. British Actuarial Journal 4: 903–23. [Google Scholar] [CrossRef]
  10. Vasicek, Oldrich. 1977. An equilibrium characterization of the term structure. Journal of Financial Economics 5: 177–88. [Google Scholar] [CrossRef]
Figure 1. Daily cross-sectional estimates of y ( θ ) (blue line) and γ y ( θ ) (red line).
Figure 1. Daily cross-sectional estimates of y ( θ ) (blue line) and γ y ( θ ) (red line).
Risks 11 00174 g001
Figure 2. Daily estimates of y ( θ ) (blue), r ( θ ) (red), and x ( θ ) (green line).
Figure 2. Daily estimates of y ( θ ) (blue), r ( θ ) (red), and x ( θ ) (green line).
Risks 11 00174 g002
Figure 3. Standard deviations of the stochastic factors v ¯ 1 k and p 1 v ¯ 1 k in the three-factor model (blue lines) and the two-factor model (red lines), as a function of the time-to-maturity k. The dashed lines represent σ v ¯ 1 k and the continuous lines represent the overall volatility σ p 1 v ¯ 1 k in the two models. The dotted lines indicate the level of the standard deviation of ω ( I + 1 ) .
Figure 3. Standard deviations of the stochastic factors v ¯ 1 k and p 1 v ¯ 1 k in the three-factor model (blue lines) and the two-factor model (red lines), as a function of the time-to-maturity k. The dashed lines represent σ v ¯ 1 k and the continuous lines represent the overall volatility σ p 1 v ¯ 1 k in the two models. The dotted lines indicate the level of the standard deviation of ω ( I + 1 ) .
Risks 11 00174 g003
Table 1. Functions B x ( τ ) and their horizontal asymptotes (left) and functions B x ( τ ) / τ (right) for α ˇ x = 0.05 (blue lines), α ˇ x = 0.1 (green lines), and α ˇ x = 0.5 (red lines).
Table 1. Functions B x ( τ ) and their horizontal asymptotes (left) and functions B x ( τ ) / τ (right) for α ˇ x = 0.05 (blue lines), α ˇ x = 0.1 (green lines), and α ˇ x = 0.5 (red lines).
Risks 11 00174 i001Risks 11 00174 i002
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Moriconi, F. A Three-Factor Market Model for Incorporating Explicit General Inflation in Non-Life Claims Reserving. Risks 2023, 11, 174. https://doi.org/10.3390/risks11100174

AMA Style

Moriconi F. A Three-Factor Market Model for Incorporating Explicit General Inflation in Non-Life Claims Reserving. Risks. 2023; 11(10):174. https://doi.org/10.3390/risks11100174

Chicago/Turabian Style

Moriconi, Franco. 2023. "A Three-Factor Market Model for Incorporating Explicit General Inflation in Non-Life Claims Reserving" Risks 11, no. 10: 174. https://doi.org/10.3390/risks11100174

APA Style

Moriconi, F. (2023). A Three-Factor Market Model for Incorporating Explicit General Inflation in Non-Life Claims Reserving. Risks, 11(10), 174. https://doi.org/10.3390/risks11100174

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop