Startseite The Optimal Strategy of Reinsurance-Investment Problem for an Insurer with Dynamic Income Under Stochastic Interest Rate
Artikel Öffentlich zugänglich

The Optimal Strategy of Reinsurance-Investment Problem for an Insurer with Dynamic Income Under Stochastic Interest Rate

  • Delei Sheng EMAIL logo
Veröffentlicht/Copyright: 25. Juni 2016
Veröffentlichen auch Sie bei De Gruyter Brill

Abstract

This paper considers the reinsurance-investment problem for an insurer with dynamic income to balance the profit of insurance company and policy-holders. The insurer’s dynamic income is given by a net premium minus a dynamic reward budget item and the net premium is obtained according to the expected premium principle. Applying the stochastic control technique, a Hamilton-Jacobi-Bellman equation is established under stochastic interest rate model and the explicit solution is obtained by maximizing the insurer’s power utility of terminal wealth. In addition, the comparison with corresponding results under constant interest rate helps us to understand the role and influence of stochastic interest rates more in-depth.

1 Introduction

Existing literatures investigate the optimal reinsurance-investment problem only from the view of insurer, but lack of meticulous care to balance the profit of both insurance company and policy-holders. An effective way to balance both profits of insurer and policy-holders is reward budget. The reward budget is a dynamic item which can be subtracted from the actuarial premium income. As a consequence, insurer’s income is dynamic, which given by an actuarial income under classic premium principle minus a dynamic reward budget item.

Portfolio has been a popular research field for more than sixty years after Markowitz’s pioneering work. It is still a very interesting problem which attracts a lot of research work in recent years. For example, Bjork, Murgoci, Zhou[1] gave a more realistic model for the case when the risk aversion depends dynamically on current wealth and studied the mean–variance portfolio optimization to find the optimal time consistent strategy. Zhu[2] studied the uncertain optimal control problem by applying Bellman’s principle of optimality and used the equation of optimality to solve a portfolio selection model. Hu, Øksendal, and Sulem[3] presented a mathematical model for a Black–Scholes market driven by fractional Brownian motion and found explicitly the optimal consumption rate and the optimal portfolio in such a market for an agent with utility of power functions. Cao and Xu[4] derived the forms of proportional and excess-of-loss reinsurance contracts under the assumption that investment fund follows the logarithm-normal distribution. Countless earlier literature need not too much listed here.

However, very few studies on the portfolio optimization under the stochastic interest rates before and small amount of literature emerges in nearly more than a decade. Korn and Kraft[5] considered investment problems under a stochastic interest rate to maximize her utility from terminal wealth without the usual Lipschitz assumptions. Grasselli[6] studied an investment problem for the HARA utility functions assuming the interest rates following the Cox-Ingersoll-Ross dynamics. Pang[7] considered a portfolio optimization problem on an infinite time horizon by maximizing the infinite horizon expected discounted log utility of consumption when the interest rate described by an ergodic Markov diffusion process. Li and Wu[8] obtained a closed-form expression of the optimal policy that maximizes a power utility for the optimal investment problem under the stochastic interest rate formulated as CIR model. Guan and Liang[9] investigated an optimal investment strategy of DC pension plan in the stochastic interest rate and stochastic volatility framework. They also investigated the optimal reinsurance and investment problem[10] when interest rate risk and inflation risk are considered in the same year.

In this paper, we study the optimal reinsurance strategy and investment strategy for an insurer with dynamic income under Ho-Lee stochastic interest rate model. A Hamilton-Jacobi-Bellman equation is established and the explicit solution for the HJB equation is obtained by maximizing the insurer’s power utility of terminal wealth. At last, some numerical simulations and graphics are given to show the impacts of different parameters. By contrast, the results under Ho-Lee model are always higher than results under constant interest rate.

The remainder of this paper is organized as follows. In Section 2, we formulate a general framework and give some models will be used hereinafter. In Section 3, a Hamilton-Jacobi-Bellman equation is established and the explicit solution has been obtained. At last, numerical simulations and graphics are provided in Section 4.

2 Formulation of the Model

Throughout this paper, (Ω, ℱ, P, {ℱt}0 ≤ tT) denotes a complete probability space satisfying the usual condition, where T > 0 is a finite constant representing the investment time horizon; ℱt stands for the information available until time t.

2.1 The Financial Market

We consider a financial market in which transaction amounts are so small that we can consider they have no influence on the prices, and the market with no arbitrage, frictionless, and continuously open as well.

Without loss of generality, we assume that the financial market is composed of two kinds of assets: Cash and equities. For the sake of simplicity, we will only consider one equity asset which can indeed represent the index of a stock market.

The instantaneous risk-free rate r(t) followed by the stochastic interest rate models of constant volatility, which described by the Ho and Lee[11] model

dr(t)=a(t)dt+bdW^(t),r(0)=r0>0,

where W^(t) is a standard Brownian motion, a(t) > 0 is a real-valued function and b > 0 is a constant.

The price S0(t) of the cash asset is given by

dS0(t)S0(t)=r(t)dt,S0(0)=1.(1)

Here r(t) is the stochastic interest rate.

The stock market index price S(t) at time t ≥ 0 can be described by the model

dS(t)S(t)=[r(t)+u(t)]dt+vdW~(t),S(0)=s0,(2)

where u(t) is another positive real-valued function, and the constant v > 0 denotes the volatility rate of risk asset. Two Brownian motions satisfying E[ W^(t), W^(t)] = ρ2t, ρ2 ∈ [−1, 1] is the correlation coefficient of W^(t) and W^(t).

In the financial model (2), u(t) > 0 is a natural hypothesis, which means that the investment on stock can make more money than investing on cash asset, which is consistent with the economic principle of high-risk high-yield. In fact, only when u(t) > 0, more people are willing to invest on equities.

2.2 The Surplus Process of Insurer

We use the classical Cramer–Lunderburg model to describe the surplus process of insurer.

dR(t)=cdtd(i=1N(t)Yi),

where N(t) is a homogeneous Poisson process with intensity λ, the claim sizes {Yi,i ≥ 1} are independent and identically distributed positive random variables with E[Yi] = μ and E[Yi2]=σ2. According to the expected value premium principle, the pure premium rate of insurer is c = (1 + θμ and θ > 0 is the safety loading of insurer.

The diffusion approximation of the compound poisson process d(i=1N(t)Yi) is given by

d(i=1N(t)Yi)=λμdtλσdW(t),(3)

where W(t) is a standard Brownian motion and independent with W^(t). But E[W(t), W^(t)] = ρ1t, ρ1 ∈ [−1, 1] is the correlation coefficient of W(t) and W^(t).

In order to avoid unbearable risk, the insurance company can buy proportional reinsurance to transfer part of the potential risk to reinsurance company. The insurer’s retention proportion is p ∈ [0, 1], the remaining 1 − p shall be borne by the reinsurer. 0 ≤ θη ≤ 1 is the safety loading of reinsurer.

dR(t)=[(ηθ)λμ+ηpλμ]dt+pλσdW(t).(4)

In the above equation, ηpλμ − (ηθμ = c − (1 + η)(1 − pμpλμ is a difference after reinsurance premium and average claim size being subtracted from the net premium.

A new notation d(ηθ)λμηpλμ introduced here for that it’s interesting and really existing for every combination of an insurance contract with an proportional reinsurance transaction. It’s known that 0 < θη < 1, (ηθ) ≥ 0, and then d > 0 which indicates the equilibrium game of both benefits of insurer and reinsurer. If η = θ, then d = 0, such proportional reinsurance transaction being more favourable for insurer but not favourable for reinsurer. Else if η > θ, then d > 0, however, actual reinsurance transaction always balance both benefits so that 0 < d < 1 is hold invariably.

2.3 The Insurer’s Wealth Process with Dynamic Income

The proportion of wealth invested on stock asset denoted by π(t) which is a progressively measurable real function, X(t) with X(0) = x0 represents the wealth of the insurer, then we have

dX(t)=[(1d)ηpλμ+(π(t)u(t)+r(t))X(t)]dt+pλσdW(t)+vπ(t)X(t)dW~(t).

The dynamic income of an insurer is equal to the net premium according to premium principle minus reward budget. Here, the dynamic income of insurer is (1 + θμkX(t), where k is the reward budget coefficient and (1 + θμ is the insurer’s net premium according to expected premium principle. The central element which makes the insurer’s income be dynamic is the reward budget which change with the dynamic process of insurer’s wealth. To determine the reward budget coefficient is a prudent event for any insurer, the reward budget deducted from premium income at a fixed share depend on both bygone experience and prospective prognosis.

dX(t)=[(1d)ηpλμ+(1k)X(t)r(t)kX(t)+π(t)(1k)X(t)u(t)]dt+pλσdW(t)+π(t)(1k)X(t)vdW~(t).(5)

A strategy α = (p(t), π (t)) is said to be admissible, if ∀ t∈ [0, T], α is ℱ(t) progressively measurable and E[0π2(t)X2(t)]<,p(t)[0,1],π(t)[0,+).Λ denotes the set of all admissible strategies and (5) has a unique (strong) solution. Suppose the insurer has a utility function U(x) which is strictly concave and continuously differentiable on (−∞, +∞) and aims to maximize the expected utility of his/her terminal wealth,

3 Main Results

In this subsection, we consider the reinsurance-investment problem with stochastic interest rate. Assuming the insurer’s objective is to maximize the expected utility from terminal wealth, where the utility U(x) is a power utility (constant relative risk aversion) function U(x) = xδ, the constant δ ∈ (0, 1) is the relative risk aversion coefficient of a risk averse insurer. This utility function plays a vital role in actuarial mathematics and insurance practice. For an admissible strategy α = (p(t), π (t)), the value function Jα(t,x,r) from state (x, r) at time t is defined as

Jα(t,x,r)=E[U(X(T))X(t)=x,r(t)=r],

and the optimal value function is given by

H(t,x,r)=max(π,p)Jα(t,x,r),(6)

with the boundary condition

H(T,x,r)=U(x).(7)

The goal of the insurer is to find an optimal strategy α* = (p*(t), π*(t)) such that Jα*(t, x, r) = H(t, x, r), where p*(t) is called the optimal reinsurance strategy and π*(t) is called the optimal investment strategy.

3.1 An Insurer with Reward Budget Under the Ho-Lee Model

If the stochastic interest rate followed by the Ho-Lee model

dr(t)=a(t)dt+bdW^(t),r(0)=r0>0.

For any H(t, x, r) ∈ C1,2,2([0, T] × R+ × R+), according to the classical results of Fleming and Soner[12], the generator can be defined as a variational operator

AHα(t,x,r)=Ht+Hx[r(1k)xkx]+Hra(t)+Hrr12b2+{Hx(1d)pηλμ+Hxrρ1bpλσ+Hxx12p2λσ2}+{Hxπu(t)(1k)x+Hxrρ2bvπ(1k)x+Hxx12v2π2(1k)2x2},

where Ht, Hx, Hr, Hxx, Hrr and Hxr denote the corresponding first and second-order partial derivatives with respect to (w.r.t.) the corresponding variables t, x, r, respectively. The Hamilton-Jacobi-Bellman (HJB) equation for problem (6) is given by

supαΛ{AHα(t,x,r)}=0(8)

with the boundary condition (7).

In the remainder of this section, we try to find the explicit solutions to problem (6) with stochastic interest by stochastic control technique.

Theorem 1

For the optimal control problem(6), assume that the objective is to maximize utility of terminal wealth, at the fixed terminal time T, and the power utility function is given by

U(x)=xδ,0<δ<1,
then the optimal value function is given by
H(t,x,r)=e16δ6k(tT)b2(1+k)2(tT)3δeb2(1+k)2(tT)3δ3ρ226(1+δ)e(1+d)3η3λ3/2μ3(1+d)ηλμ+b(1+k)(t+T)δρ1σ36b(1+k)(1+δ)ρ1σ3etT2b(1k)vδ2ρ2(Tτ)u[τ]+δu[τ]22v2(1+δ)dτetT(1k)δ(Tτ)a[τ]dτe(1k)δ(Tt)rxδ.

The optimal investment strategy is

π=u(t)v2(1k)(δ1)ρ2bδ(Tt)v(δ1),(9)
and the optimal reinsurance strategy is
p=(1d)ηλμx+ρ1bλσx(1k)δ(Tt)λσ2(δ1),(10)
where ρ1 ∈ [−1, 1], ρ2 ∈ [−1, 1], b > 0, v > 0, μ > 0, λ > 0, σ > 0, d > 0 and 0 < ηθ < 1 are constants, and u(t) > 0, a(t) > 0 are deterministic real functions.

Proof

The HJB equation (8) rewritten by

Ht+Hx[r(1k)xkx]+Hra(t)+Hrr12b2+maxp{Hx(1d)pηλμ+Hxrρ1bpλσ+Hxx12p2λσ2}+maxπ{Hxπu(t)(1k)x+Hxrρ2bvπ(1k)x+Hxx12v2π2(1k)2x2}=0.(11)

According to the first-order necessary conditions of the extremum, the optimal strategies

p=(1d)ηλμHx+ρ1bλσHxrλσ2Hxx,(12)
π=u(t)Hx+ρ2bvHxrv2(1k)xHxx.(13)

Inserting (12) and (13) into (11), the HJB equation becomes

Ht+Hra(t)+Hx[r(1k)xkx]+Hrr12b2[(1d)ηλμHx+ρ1bλσHxr]22λσ2Hxx[u(t)Hx+ρ2bvHxr]22v2Hxx=0.(14)

Conjecture of the solution has the following form

H(t,x,r)=f(t)eh(t)rxδ,H(T,x,r)=xδ,0<δ<1,f(T)=1,h(T)=0.(15)

The derivatives of the conjecture with respect to t, x, r, respectively.

Ht=(f(t)+f(t)h(t)r)eh(t)rxδ,Hx=f(t)eh(t)rδxδ1,Hr=f(t)eh(t)rxδh(t),Hxx=f(t)eh(t)rδ(δ1)xδ2,Hrr=f(t)eh(t)rxδh(t)2,Hxr=f(t)eh(t)rδxδ1h(t).

Plugging these derivatives into (14)

f(t)+f(t)h(t)r+f(t)δ[r(1k)k]+f(t)h(t)a(t)+b22f(t)h(t)2[(1d)ηλμ+ρ1bλσh(t)]2δf(t)2λσ2(δ1)[u(t)+ρ2bvh(t)]2δf(t)2v2(δ1)=0.(16)

Splitting (16) into two equations

[h(t)+δ(1k)]f(t)r=0,

i.e.,

h(t)+δ(1k)=0,h(T)=0.(17)
f(t)f(t)δk+f(t)h(t)a(t)+b22f(t)h(t)2[(1d)ηλμ+ρ1bλσh(t)]2δf(t)2λσ2(δ1)[u(t)+ρ2bvh(t)]2δf(t)2v2(δ1)=0,f(T)=1.(18)

The solution of (17) is

h(t)=(1k)δ(Tt).(19)

Inserting (19) into (18),

f(t)f(t)=δka(t)(1k)δ(Tt)b22[(1k)δ(Tt)]2+[(1d)ηλμδ+ρ1bλσδ(1k)δ(Tt)]22λσ2δ(δ1)+[δu(t)+ρ2bvδ(1k)δ(Tt)]22v2δ(δ1),f(T)=1.(20)

The solution of (20)

f(t)=e16δ6k(tT)b2(1+k)2(tT)3δeb2(1+k)2(tT)3δ3ρ226(1+δ)e(1+d)3η3λ3/2μ3(1+d)ηλμ+b(1+k)(t+T)δρ1σ36b(1+k)(1+δ)ρ1σ3etT2b(1k)vδ2ρ2(Tτ)u[τ]+δu[τ]22v2(1+δ)dτetT(1k)δ(Tτ)a[τ]dτ.(21)

Plugging (21) and (19) into (15), the optimal value function is obtained.

H(t,x,r)=e16δ6k(tT)b2(1+k)2(tT)3δeb2(1+k)2(tT)3δ3ρ226(1+δ)e(1+d)3η3λ3/2μ3(1+d)ηλμ+b(1+k)(t+T)δρ1σ36b(1+k)(1+δ)ρ1σ3etT2b(1k)vδ2ρ2(Tτ)u[τ]+δu[τ]22v2(1+δ)dτetT(1k)δ(Tτ)a[τ]dτe(1k)δ(Tt)rxδ.

Plugging (19) into (12) and (13), respectively. The optimal strategy (p*, π*)

p=(1d)ηλμHx+ρ1bλσHxrλσ2Hxx=(1d)ηλμf(t)eh(t)rδxδ1+ρ1bλσf(t)eh(t)rδxδ1h(t)λσ2f(t)eh(t)rδ(δ1)xδ2=(1d)ηλμx+ρ1bλσxh(t)λσ2(δ1)=(1d)ηλμx+ρ1bλσx(1k)δ(Tt)λσ2(δ1),
π=u(t)Hx+ρ2bvHxrv2(1k)xHxx=u(t)f(t)eh(t)rδxδ1+ρ2bvf(t)eh(t)rδxδ1h(t)v2(1k)xf(t)eh(t)rδ(δ1)xδ2=u(t)+ρ2bvh(t)v2(1k)(δ1)=u(t)+ρ2bv(1k)δ(Tt)v2(1k)(δ1)=u(t)v2(1k)(δ1)ρ2bδ(Tt)v(δ1).

The proof is finished.

Remark 1

For the optimal investment strategy (9)

π=u(t)v2(1k)(δ1)ρ2bδ(Tt)v(δ1)=1v(1δ)[1(1k)u(t)v+bρ2δ(Tt)],

where u(t) > 0 is a positive real-valued function and the constant b > 0 is the volatility of interest rate. The proportion u(t)v is excess-appreciate-rate of stock divided by the volatility rate of stock.

Observed carefully from the model (2) and (1), we conclude that u(t) is the excess-appreciate-rate of stock relative to cash. The principal parts of the optimal investment strategy π* are two components: one is the relative proportion u(t)v, the other is the volatility of interest rate b. Because the remainder parts 1(1k) and ρ2δ (Tt) are usually fixed in a specific insurer’s implementing process. In addition, both of them are so small actually.

In other words, the optimal investment strategy π* is mainly determined by two components: the relative ratio of stock’s excess-appreciate-rate divided by the stock’s volatility rate and the volatility rate of stochastic interest.

Remark 2

For the optimal reinsurance strategy (10), if we let m = λμ which represents mean-value of insurer’s claim size and n=λσ which denotes variance of insurer’s claim size, then

p=(1d)ηλμx+ρ1bλσx(1k)δ(Tt)λσ2(δ1)=ηmx+ρ1bnx(1k)δ(Tt)n2(1δ)=x[η(1δ)mn2+ρ1bn(1k)δ(1δ)(Tt)].

It’s nature that the optimal reinsurance strategy p* increases with respect to the wealth level at time t. The most interesting discovery is that we find the optimal reinsurance strategy mainly determined by mn2 and ρ1bn. Because η(1δ) and (1k)δ(1δ)(Tt) are usually fixed in a specific insurer’s implementing process. In fact, ρ1bn shows us the correlation of insurer’s claim volatility and the volatility of interest rate. mn2 is the proportion for mean and variance of insurer’s claim size. After analyzing carefully, we conclude that the volatility of insurer’s claim size is the origin force to change of the optimal reinsurance strategy p*.

Corollary 1

Assuming that the objective is to maximize utility of terminal wealth, at the fixed terminal time T, and the power utility function is given by

U(x)=xδ,0<δ<1,

If a(t) = a > 0 andu(t) = u > 0, where a and u are constants, the control problem(6)has a simple version, then the optimal value function is given by in simple form

H(t,x,r)=e16δ6k(tT)b2(1+k)2(tT)3δeb2(1+k)2(tT)3δ3ρ226(1+δ)e(1+d)3η3λ3/2μ3(1+d)ηλμ+b(1+k)(t+T)δρ1σ36b(1+k)(1+δ)ρ1σ3e(tT)uδ(u+b(1+k)(tT)vδρ2)2v2(1+δ)e12a(1+k)(tT)2δe(1k)δ(Tt)rxδ,(22)
and the optimal strategy (p*, π*) in simple form given by
p=(1d)ημxσ2(1δ)+ρ1bx(1k)δ(Tt)λσ(1δ),π=uv2(1k)(1δ)+ρ2bδ(Tt)v(1δ),
where ρ1 ∈ [−1, 1], ρ2 ∈ [−1, 1], b > 0, v > 0, μ > 0, λ > 0, σ > 0, d > 0 and 0 < ηθ < 1 are constants.

Proof

In fact, if the deterministic function a(t) = a > 0 and u(t) = u > 0, where a and u are constants, similar to the above proof completely, we obtain that f(t) can be given in more simple form.

f(t)=e16δ6k(tT)b2(1+k)2(tT)3δeb2(1+k)2(tT)3δ3ρ226(1+δ)e(1+d)3η3λ3/2μ3(1+d)ηλμ+b(1+k)(t+T)δρ1σ36b(1+k)(1+δ)ρ1σ3e(tT)uδ(u+b(1+k)(tT)vδρ2)2v2(1+δ)e12a(1+k)(tT)2δ,

and h(t) = (1 − k) δ (Tt). So the desired value function and optimal strategy are obtained.

Theorem 2

(verification theorem) Let H(t,x,r) be a convex, twice differential solution to(8)with the boundary condition(7)such that |H(0, x0, r0)| < ∞. Then, for all t ∈ [0, T], x ∈ (0, ∞),

  1. H(t, x, r) ≤ Jα(t, x, r) for every admissible strategy α.

  2. If there exists an admissible control α* = (p*, π*) such that

    αargmax[AαH(t,x(t),r(t))],
    then the value function H(t,x,r) = Jα*(t, x, r) and the policy α* = (p*, π*) is the optimal strategy corresponding to Xα*(t) which is the solution of(5).

Proof

Details you can refer to Fleming and Soner[12], Gu, Guo, Li and Zeng[13], and so on. The proof is omitted here.

3.3 A Special Case of an Insurer with Reward Budget Under a Constant Interest Rate

In this subsection, we investigate a special case that the interest rate r(t) is a constant r. All assumptions are the same as above except that r(t) ≡ r throughout this part, thus we have

dS0(t)S0(t)=rdt,S0(0)=1,

and

dS(t)S(t)=[r+u(t)]dt+vdW~(t),S(0)=s0,

where r, v are positive constants, u(t) is a positive real-valued function. The surplus processes with claim and reinsurance are the same as (3) and (4), W(t) and W^(t) are independent standard Brownian motions. So we get the wealth process with reward budget and investing in the financial market as follows:

dX(t)=[(1k)X(t)rkX(t)+(1d)ηpλμ+π(t)(1k)X(t)u(t)]dt+π(t)(1k)X(t)vdW~(t)+pλσdW(t).

The utility U(x) is also a power utility (constant relative risk aversion) function U(x) = xδ to maximize the expected utility of terminal wealth, where the constant δ ∈ (0, 1) is the relative risk aversion coefficient of a risk averse insurer. For the optimal control problem

H(t,x)=max(π,p)Jα(t,x)

with the boundary condition H(T,x,r) = xδ.

Applying the classical tools of stochastic optimal control, if the optimal value function H(t,x) ∈ C1,2([0, TR), then H satisfies the following Hamilton-Jacobi-Bellman (HJB) equation:

Ht+Hx[r(1+r)k]x+maxp{Hx(1d)pηλμ+Hxx12p2λσ2}+maxπ{Hxπu(t)(1k)x+Hxx12v2π2(1k)2x2}=0.

According to the first order necessary condition of extreme, we obtain that

p=(1d)ημHxσ2Hxx,π=u(t)Hxv2(1k)xHxx.

Conjecture of the solution has the following form

H(t,x)=f(t)xδ,H(T,x)=xδ,0<δ<1,f(T)=1.

Through the similar calculation as before, we obtain the value function

H(t,x)=xδexp{[r+(1+r)k(1d)2η2λμ22σ2(δ1)]δ(Tt)δtT[u(τ)22v2(δ1)]dτ},(23)

and the optimal strategy α* = (p*, π*) as follows:

p=1(1δ)(1d)ημxσ2andπ=u(t)v2(1k)(1δ).

Remark 3

Comparing the results of constant interest rate with the results of stochastic interest rate, it is easy to find that the optimal strategy under the stochastic interest rate exists one term more than the optimal strategy under a constant interest rate.

pstochasticpconstant=ρ1bλσδ(Tt)(1δ)(1k)x.

If ρ1 > 0, then p* stochastic > p* constant, else if ρ1 < 0, then p* stochastic < p* constant.

Let’s pay attention to bλσ which is a ratio of interest rate volatility rate divided by the volatility rate of claim size. So in other words, it’s a share of risk-free asset volatility relative to the claim size volatility at time t.

Similarly,

πstochasticπconstant=ρ2bvδ(Tt)(1δ).

If ρ2 > 0, then π* stochastic > π* constant, else if ρ2 < 0, then π* stochastic < π* constant.

It is noteworthy that bv is a ratio of interest rate volatility rate divided by the volatility rate of risky assets, in other words, it’s a share of risk-free asset volatility relative to the risk assets volatility.

In one word, whether the expression ρ1bλσδ(Tt)(1δ)(1k)x or the formula ρ2bvδ(Tt)(1δ) reflects the interest rate volatility’s influence on the optimal strategy. However, this influence is a synergy effect based on the three major factors: Interest rate volatility, risky asset volatility and claim size volatility.

4 Numerical Simulations

In this section we illustrate some experiments on simulated data for the case of u(t) = u and a(t) = a where u and a are positive constants. Unless otherwise stated, the related model parameters are as follows:

Table 1

Numerical values of parameters

xηλμbvukTδσρ1ρ2a
30.120.20.10.10.150.05100.520.30.10.06

First, we propose two experiments to study the optimal value function, the one to analyze the formula (22) under the stochastic interest rate and the other to analyze the formulas of (23) under a constant interest rate. The simulations mainly focus on the changes of the value function with respect to x, t, r, u, respectively. See Figure 1, Figure 2, Figure 3 and Figure 4.

Figure 1
Figure 1
Figure 2
Figure 2
Figure 3
Figure 3
Figure 4
Figure 4

The analysis below will always be a comparison of the results under the stochastic interest rate and the corresponding results under a constant interest rate. Each numerical graphic consists of two curves, dashed line representing the result under the stochastic interest rates and solid line representing the result under the constant interest rate. With the curve under constant interest rate as reference, it can help us to understand the role and influence of stochastic interest rates more clearly and more in-depth.

Figure 1 shows that the value function increases as the wealth level increasing, but the speed becoming slower and slower. This is the basic characteristics of utility function. The results obtained in line with the features of utility function, which confirms the correctness of our result.

Figure 2 shows the trend that the value function changes with time. When other parameters determined, you can see that the utility under the stochastic interest rate decreases with time but there is little change of the utility under a constant interest rate.

Figure 3 and Figure 4 shows the utility increases with the return level. Both the increase of interest rate r and the increase in risk premium u will lead to a significant increase in utility.

But the utility increasing with risk premia is more significant.

In the following, let’s review the optimal strategy α* = (p*, π*) with a constant interest rate

p=(1d)ημx(1δ)σ2,π=uv2(1k)(1δ),

and the optimal strategy α* = (p*, π*) with stochastic interest rate

p=ημxσ2(1δ)+ρ1bλσδ(Tt)(1δ)(1k)x,π=uv2(1k)(1δ)+ρ2bvδ(Tt)(1δ).

The rest of this part mainly focus on the analysis and comparisons of the optimal strategies, which will be illustrated through several graphics.

For the optimal reinsurance proportion p*, the simulations mainly focus on the changes with respect to x, δ, μ, σ, respectively. See Figure 5, Figure 6, Figure 7 and Figure 8.

Figure 5
Figure 5
Figure 6
Figure 6
Figure 7
Figure 7
Figure 8
Figure 8

Figure 5 shows that whether a constant interest rate or a stochastic interest rate, the optimal reinsurance proportion always increases with the wealth level x. But the trend under the stochastic interest rate is more significant, which in turn suggests that the growth of insurer’s wealth under stochastic interest rates need to bear even greater risk.

Figure 6 reflects the sensitivity of the reinsurance proportion on the degree of risk aversion. Figure 7 and Figure 8 show the trend of reinsurance proportion changing with the claim size and the claim volatility, respectively.

For the optimal investment strategy π*, the simulations mainly focus on the changes with respect to u, v, k, δ, respectively. See Figure 9, Figure 10, Figure 11 and Figure 12.

Figure 9
Figure 9
Figure 10
Figure 10
Figure 11
Figure 11
Figure 12
Figure 12

Figure 9 and Figure 10 show the trend of investment strategy changing with the risk premium u and the risk volatility v, respectively. Figure 11 shows that the reward budget increasing inevitably leads to the increase of the risk investment share. The increase of investment income can make up for the increase of budget share. Figure 12 shows the impacts of risk preference on the investment share.

Remark 4

Overall, the curve under the stochastic interest rate is always above the corresponding curve under a constant interest rate. Comparing with constant interest rate, the stochastic interest rate brings more risk into the wealth process, so that the insurer has to transfer more of the risk to the reinsurer and the reinsurance proportion p will be higher. However, transferring more risks to reinsurer must increase the insurance company’s spending, which requires insurer to seek much greater investment return. This eventually leads to a further increase of the investment share π on risk asset.

5 Conclusion

Based on the viewpoint of interaction between insurance company and policy-holder, we introduce reward budget into the insurance risk management process in this paper, which makes the insurer’s income to be dynamic. Some results are obtained under the stochastic interest rate model, but the follow-up work remains to be further perfect.


Supported by the National Natural Science Foundation of China (11301376)


Acknowledgements

We thank the referees for their time and comments.

References

[1] Bjork T, Murgoci A, Zhou X Y. mean-varicance portfolio optimization with state-dependent risk aversion. Mathematical Finance, 2014, 24(1): 1–24.10.1111/j.1467-9965.2011.00515.xSuche in Google Scholar

[2] Zhu Y. Uncertain optimal control with application to a portfolio selection model. Cybernetics and Systems, 2010, 41(7): 535–547.10.1080/01969722.2010.511552Suche in Google Scholar

[3] Hu Y, Øksendal B, Sulem A. Optimal consumption and portfolio in a Black-Scholes market driven by fractional Brownian motion. Infinite Dimensional Analysis Quantum Probability and Related Topics, 2011, 6(4): 519–536.10.1142/S0219025703001432Suche in Google Scholar

[4] Cao Y, Xu J. Proportional and excess-of-loss reinsurance under investment gains. Applied Mathematics and Computation, 2010, 217(6): 2546–2550.10.1016/j.amc.2010.07.067Suche in Google Scholar

[5] Korn R, Kraft H. A stochastic control approach to portfolio problems with stochastic interest rates. SIAM Journal on Control and Optimization, 2002, 40(4): 1250–1269.10.1137/S0363012900377791Suche in Google Scholar

[6] Grasselli M. A stability result for the HARA class with stochastic interest rates. Insurance: Mathematics and Economics, 2003, 33(3): 611–627.10.1016/j.insmatheco.2003.09.003Suche in Google Scholar

[7] Pang T. Stochastic Portfolio Optimization with Log Utility. International Journal of Theoretical and Applied Finance, 2006, 9(6): 869–887.10.1142/S0219024906003858Suche in Google Scholar

[8] Li J, Wu R. Optimal investment problem with stochastic interest rate and stochastic volatility: Maximizing a power utility. Applied Stochastic Models in Business and Industry, 2009, 25(3): 407–420.10.1002/asmb.759Suche in Google Scholar

[9] Guan G, Liang Z. Optimal management of DC pension plan in a stochastic interest rate and stochastic volatility framework. Insurance: Mathematics and Economics, 2014, 57(6): 58–66.10.1016/j.insmatheco.2014.05.004Suche in Google Scholar

[10] Guan G, Liang Z. Optimal reinsurance and investment strategies for insurer under interest rate and inflation risks. Insurance Math. Econom., 2014, 55(1): 105–115.10.1016/j.insmatheco.2014.01.007Suche in Google Scholar

[11] Ho T S Y, Lee S B. Term structure movements and pricing interest rate contingent claims. Journal of Finance, 1986, 41(5): 1011–1029.10.1111/j.1540-6261.1986.tb02528.xSuche in Google Scholar

[12] Fleming W H, Soner H M. Controlled Markov Processes and Viscosity Solutions. Springer, Berlin, 1993.Suche in Google Scholar

[13] Gu A, Guo X, Li Z, et al. Optimal control of excess-of-loss reinsurance and investment for insurers under a CEV model. Insurance: Mathematics and Economics, 2012, 51(3): 674–684.10.1016/j.insmatheco.2012.09.003Suche in Google Scholar

Received: 2015-9-25
Accepted: 2015-11-10
Published Online: 2016-6-25

© 2016 Walter de Gruyter GmbH, Berlin/Boston

Heruntergeladen am 20.9.2025 von https://www.degruyterbrill.com/document/doi/10.21078/JSSI-2016-244-14/html
Button zum nach oben scrollen