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 is traded. In our inflation problem, 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 , which will be denoted by . 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 (measured in ) 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 the real value function at time t, that is, the market price expressed in units of 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 the real value at time t of the future payoff , where is a sequence of real cash-flows to be received on a corresponding sequence of dates after t. The payoff 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 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
at the maturity date
. The time
t price of this ZCB with maturity
T is denoted by
If we consider a ZCB with deterministic payoff
to be received at time
T, we have, to prevent arbitrages,
that is,
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
, is given by
The
real instantaneous rate of interest at time
t is defined as the rate of return of a ZCB currently maturing, i.e.,
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
as the single risk factor. More specifically, we shall assume the following Markov property:
where the price
at time
t depends only on the current value
of
x. In this case,
is also referred to as the
state variable.
We assume that
is a diffusion process described by the following stochastic differential equation (SDE):
where
is a Wiener process. The choice of the functions
and
, together with the initial condition
, completely specifies the conditional probability distribution of
given
for
.
We shall assume for this process the “Vasicek specifications” (see
Vasicek (
1977)), i.e., a mean-reverting drift
and a constant diffusion coefficient
. Precisely,
With these
and
functions, the process (
2) is a special case of the Ornstein–Uhlenbeck (OU) process, which we characterize by the following SDE:
In
Appendix A, we provide some fundamental results on the
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
given
, for
, is normal, with mean
and variance
which is independent of
. 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
and
(and, obviously, interpreting the OU process
as
).
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
of any real security which depends only on
and
t must satisfy the no-arbitrage relation:
1
This expression prescribes that the instantaneous expected return of
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
expressing the
semi-elasticity of
with respect to
x, multiplied by the
risk premium coefficient . The form of this coefficient is to be specified, since, by Ito’s lemma,
expression (
7) can be written as
This deterministic (parabolic) partial differential equation is the GVE of our model. Once the coefficients
, and
are specified, the price
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
at time
, this condition is
.
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
can be expressed in integral (i.e., in expectation) form as
where
denotes the conditional expectation taken with respect to the so-called
risk-neutral probability measure
, that is, the probability determined by the
risk-adjusted drift
and the diffusion coefficient
. Since
and
have been chosen in the Vasicek form, to completely define the pricing problem, it remains to specify the form of the risk premium
.
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,
is given by the covariance of changes in the state variable
with percentage changes in the optimally allocated wealth
(the value of the so-called
market portfolio), i.e.,
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
:
With this choice, the risk-adjusted drift is
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,
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
and we find that the following closed-form expression holds, for
:
with
and
In order to derive this expression, let us denote by
, the moment-generating function of the random variable
X. As is well-known, in the normal case
, one has
. The ZCB price
is the exponential expectation
with
. Using Result A2 in
Appendix A.3, after the changes
we find that
is normal, so
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
Hence, the log-return at time
t, for the given
, is an affine function of the risk factor
. Moreover, it is immediately verified that
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
defined in (
15). By (
14),
Hence,
is equal to (minus) the semi-elasticity of the price
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
w.r.t.
x. For this reason,
is typically used as a measure of the real interest rate risk of the ZCB
.
is measured in time units, is zero for
, and is monotonically increasing with
. If
(genuine mean-reversion), then
for
and has a horizontal asymptote at level
; the mean-reversion produces a ceiling on the interest rate risk. If
(“mean repulsion”), then
for
and diverges for
. The function
for different values of
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
for
. For the corresponding log-returns, we have, by (
16),
where the conditional variance of
, given by (
6), depends on the parameters
. As illustrated in the right figure of
Table 1, for
, the factor
as a function of
takes values in
, is monotonically decreasing, and decreases more rapidly for increasing values of
. 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
at time
, can be simulated at time
t as:
The conditional mean and variance of the normal variable
are given by (
5) and (
6), and the random variable
is lognormal. Of course,
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 one unit of money (say, EUR 1). We denote by the price in units of money of one unit of the reference basket 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, is simply the number of needed for receiving one at time t. However, this is sufficient for modeling nominal bonds.
Let us denote by
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
because receiving one real unit,
, at time
t is equivalent to receiving
nominal units,
, at the same date. Let us consider the real value
of a bond which provides with certainty the unit nominal payoff
at maturity
. Since, by relation (
20),
, to prevent arbitrage, the real value of this bond can be represented by
Then, for the nominal unit ZCB price
, i.e., the money value at time
t of one money unit at time
T, we have, again by relation (
20),
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
—we have to extend the model for
by making appropriate assumptions on the stochastic dynamics of
.
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:
with
Since
, one has
Therefore,
is the
expected instantaneous rate of inflation at time
t. If
is deterministic, the CPI process is a geometric Brownian motion. However, in order to obtain a more realistic model, we assume that
is also stochastic.
Our assumption is that the expected instantaneous inflation rate is a diffusion process described by the following SDE:
with
With this Vasicek specification,
is also an Ornstein–Uhlembeck process.
As concerns the correlation structure, we assume
and
Assumptions (
26) state the independence between real and nominal quantities and have important consequences. First of all, the form of the market risk premium
is the same as in (
12) and the market risk premia for the new two state variables are zero:
since property (
11) also holds for the nominal quantities
and
and
is a real quantity.
2.4.2. The Inflation Discount Factor and the Fisher Relation
Since the real payoff
is independent of
p and
y, the price of real bonds
is still obtained by the univariate model for
. In particular,
still has the Vasicek form (
14).
As concerns
, we have, by extension of (
10):
where now the risk-neutral measure
is three-dimensional. However, by the independence assumptions (
26), we have
Hence,
For the expectation on the right-hand side, we observe that, since in the real economy, the risk premia
and
are set to zero, the relevant risk-neutral measure is equal to the natural probability measure. We define the folloing:
which will be referred to as the
inflation discount factor on the time interval
. Hence, we have
Expression (
28) can be referred to as the
Fisher relation in multiplicative form. As we shall see, since
is stochastic, future expectations, i.e.,
for
, 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,
expression (
28) provides
where
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
:
Of course,
Using Result A3 in
Appendix A.4, with
and
, for
, the random variable
is normal, with mean
and variance
where
is given by
with
, that is,
It is important to note that these expressions are independent of the current level
of the CPI. Furthermore, given the model parameters, the variance
only depends on
and the mean
depends on
t only through
.
Using (
31) and (
32), an expression for
is immediately obtained by computing the following moment-generating function:
We obtain
where
is given by (
33) and
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
, and future values
and
, for
, are, respectively, lognormally and normally distributed. Moreover, it is easily verified that
For the function
given by (
33), the same properties and the same interpretation hold as for
.
is a time measure of expected inflation risk, providing the sensitivity of the inflation discount factor
to the risk factor
. The only difference is that, being that
is the natural risk reversion parameter, it is positive by assumption, so only the “capped” behavior (asymptote at
) is possible for
; the difference
is strictly positive and monotonically increasing with
.
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
,
, of the CPI are also simulated as lognormal, that is,
with
standard normal. Taking account of (
34), expression (
38) can be also written as
This expression is important because if the discount factors 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 . In this case, since and do not enter into the expression of the variance , their observation or estimation is not required for the simulations of .
Remark 1. In the two-factor model used in DFM23, where is deterministic, the risk-neutral inflation discount factor is defined aswhere the expectation is taken w.r.t. the risk-neutral probability measure of the price process in the nominal economy setting and we use again the definition . However, in the two-factor model, one has , then , as in (30). Thus, the two models are formally consistent. Of course, in the three-factor model, the distribution of 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
, with drift parameters
the process
is described by the following SDE:
By Result A2 in
Appendix A.3, the stochastic integral
is normal with mean
and variance
Then, it is easily verified that expression (
35) can also be obtained as
The interesting point is that this is equivalent to saying
where
is as in (
37) and
is the expectation taken with respect to a probability measure
defined by the modified mean-reverting drift
. By the Feynman–Kac theorem, this expression of
is the solution, with terminal condition
, of the following partial differential equation:
which has exactly the same form of the general valuation Equation (
9) for the real ZCB price
. As is shown in
Appendix B, (
43) is actually the general valuation equation for the discount factor
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
and
factors are different from zero, and take exactly the form which determines the risk-adjusted drift
and the modified process
.
3 Then
is actually the risk-neutral measure in the nominal economy and
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:
The instantaneous nominal interest rate
is naturally defined as
with
given in (
29). Therefore, by (
17) and (
37), we obtain the Fisher relation for instantaneous returns:
Since
depends on the two risk factors
and
, two different risk measures must be used:
where
and
are the real interest rate measure and the expected inflation risk measure previously defined.
Given the normality properties of
and
,
is also normal, with mean equal to the sum of the means minus
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.