Next Article in Journal
Monte Carlo Simulation of the Moments of a Copula-Dependent Risk Process with Weibull Interwaiting Time
Next Article in Special Issue
Transformations of Telegraph Processes and Their Financial Applications
Previous Article in Journal
The Empirical Analysis of the Core Competencies of the Company’s Resource Management Risk. Preliminary Study
Previous Article in Special Issue
Multivariate General Compound Point Processes in Limit Order Books
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Merton Investment Problems in Finance and Insurance for the Hawkes-Based Models

by
Anatoliy Swishchuk
Department of Mathematics & Statistics, University of Calgary, Calgary, AB T2N 1N4, Canada
Risks 2021, 9(6), 108; https://doi.org/10.3390/risks9060108
Submission received: 6 May 2021 / Revised: 24 May 2021 / Accepted: 26 May 2021 / Published: 3 June 2021
(This article belongs to the Special Issue Stochastic Modelling in Financial Mathematics)

Abstract

:
We show how to solve Merton optimal investment stochastic control problem for Hawkes-based models in finance and insurance (Propositions 1 and 2), i.e., for a wealth portfolio X ( t ) consisting of a bond and a stock price described by general compound Hawkes process (GCHP), and for a capital R ( t ) (risk process) of an insurance company with the amount of claims described by the risk model based on GCHP. The main approach in both cases is to use functional central limit theorem for the GCHP to approximate it with a diffusion process. Then we construct and solve Hamilton–Jacobi–Bellman (HJB) equation for the expected utility function. The novelty of the results consists of the new Hawkes-based models and in the new optimal investment results in finance and insurance for those models.

1. Introduction

Merton optimal investment and consumption stochastic problem is one of the most studied classical problem in finance (Merton 1969, 1971, 1990; Bjork 2009; Karatzas and Shreve 1998). In this paper, we will show how to solve the Merton optimal investment stochastic control problem for Hawkes-based models in finance and insurance, i.e., for a wealth portfolio X ( t ) consisting of a bond and a stock price described by general compound Hawkes process (GCHP) (Swishchuk 2020; Swishchuk and Huffman 2020; Swishchuk 2017b), and for a capital R ( t ) of an insurance company with the amount of claims described by risk model based on GCHP (Swishchuk 2018; Swishchuk et al. 2020).
Namely, we will show how to solve the following two portfolio investment problems.
(1) Merton portfolio optimization problem in finance (Merton 1969, 1971) aims to find the optimal investment strategy for the investor with those two objects of investment, namely risk-less asset (e.g., investment grade government bonds)), with a fixed rate of interest r , and a number of risky assets (e.g., stocks) whose price GCHP. In this way, in our case, we suppose that B t and S t follows the following dynamics, respectively:
B t = B 0 exp { r t } S t = S 0 exp { G ( t ) } ,
where G ( t ) : = i = 1 N ( t ) a ( X i ) is the GCHP, X i is a discrete-time Markov chain (MC) with finite or infinite states, r > 0 is the interest rate, and a ( x ) is a continuous and finite function on X . We note, that the justification of using HP in finance may be found in (Da Fonseca and Zaatour 2013), and using GCHP that based on HP N ( t ) and a ( X i ) may be found in (Swishchuk and He 2019; Swishchuk and Huffman 2020). The model for a stock price S ( t ) that based on GCHP is a new and original in this paper. The investor starts with an initial amount of money, say X 0 = x , and wishes to decide how much money to invest in risky and risk-less assets to maximize the final wealth X t at the maturity T ;
(2) Merton portfolio optimization problem in insurance aims to find an optimal investment for the capital R ( t ) of an insurance company at time t ( R ( t ) is actually the risk model based on general compound Hawkes process (GCHP) (Swishchuk 2018; Swishchuk et al. 2020), when an investor decides to invest some capital A ( t ) in risky assets (e.g., stocks) and the rest, ( R ( t ) A ( t ) ) in risk-free assets (e.g., bonds or bank account). We note, that the risk model R ( t ) , based on GCHP, has the following representation:
R ( t ) = u + c t i = 1 N ( t ) a ( X i ) .
Here R ( 0 ) = u is an insurance company’s initial capital, c > 0 is the premium rate, N ( t ) is the Hawkes process, X i is a discrete-time finite or infinite state Markov chain with state space X = { 1 , 2 , 3 , . . . , N } or X = { 1 , 2 , . . . , N , . . . } , , respectively, and a ( x ) is a continuous and finite function on X . We note, that the justification for the Hawkes-based risk model in the form of the above equation may be found in (Swishchuk et al. 2020).
The investor starts with an initial capital, say R ( 0 ) = u , and wishes to decide how much money to invest in risky and risk-less assets to maximize the capital R ( t ) .
We solve both problems using diffusion approximation for GCHP G ( t ) : = i = 1 N ( t ) a ( X i ) (see Section 4 and Section 5). We note, that it is not a simplification of the initial models: the resulting models contain all the parameters of the initial models, including the parameters of Hawkes process. Furthermore, significance and insights of the results are discussed in Remark 6 and Remark 7 (Some Insights into the Results). In this case both problems can be solved explicitly. However, we cannot say this if we would like to dal with initial models and going without diffusion approximation (see Section 6. Discussion). We believe that these two problems for those two different models in finance and insurance are considered in the literature for the first time, because none of author’s 9 papers in the References contain similar or even close results.
The novelties of the paper are the following ones: (1) we consider a new model for the stock price S ( t ) in the form S t = S 0 exp { G ( t ) } , where G ( t ) : = i = 1 N ( t ) a ( X i ) is the GCHP; we call it Hawkes-based model for the stock price (or geometric general compound Hawkes process, similar to geometric Brownian motion); (2) solution of Merton investment problem for this Hawkes-based model; (3) solution of Merton investment problem for the Hawkes-based risk model.
The structure of the paper is the following one. Literature review is presented in Section 2. Section 3 is devoted to the definitions and properties of Hawkes process and general compound Hawkes processes, and LLN (Law of Large Numbers) and FCLT (Functional Central Limit Theorem) for them. Section 4 deals with Merton investment problem in finance for the stock price described by GCHP, and Section 5 deals with Merton investment problem in insurance for the risk model based on GCHP. Section 6 contains some discussions and describes the future work, and Section 7 concludes the paper.

2. Literature Review

The Merton optimal investment and consumption stochastic problem in finance was first considered in the seminal papers of Merton, (Merton 1969, 1971). General description of the problem and coverage of most today’s problems and methods may be found in (Fleming and Rishel 1975; Fleming and Soner 1993; Bjork 2009; Merton 1990; Karatzas and Shreve 1998).
The first papers on stochastic optimal control in insurance appeared relatively recently, e.g., we would like to mention the papers written by (Martin-Löf 1994; Brockett and Xia 1995; Browne 1995), to name a few. Since then many papers and books were written on this topic including (Asmussen and Taksar 1997; Højgaard and Taksar 1999; Schmidli 2008). Financial control methods applied in insurance, such as, e.g., in the management and the control of the specific risk insurance companies, are described in (Hipp 2004). The paper (Dufresne and Gerber 1991) considered a model in risk theory based on the compound Poisson process perturbed by diffusion. The paper (Grandits 2003) studied the Cramer-Lundberg approximation in the optimal investment case. Asymptotic ruin probabilities and optimal investment were investigated in (Gaier et al. 2003). Optimal risk distribution control model with application to insurance was studied in (Taksar 2000). Applications of stochastic processes in insurance and finance may be found in (Rolski et al. 1998). Many mathematical methods and aspects in risk theory may be found in (Bühlmann 1970; Gerber 1979; Grandell 1990).
We note that an alternative approach to consumption-portfolio optimization problem based on martingale methods were developed by (Karatzas et al. 1986; Pliska 1986; Cox and Huang 1989). Applications of martingale methods to the basic optimization problems can be found in (Cox and Huang 1989; Karatzas 1997; Korn and Korn 2000). Investment problem is also closely associated with risk management problems, such as, e.g., insurance/reinsurance and risk prevention. Probably Arrow’s 1963 paper (Arrow 1963) was the first one that drown attention to risk management with insurance. How insurance can be used as a risk prevention tool was shown by (Ehrlich and Becker 1972). Some early contributions to insurance/reinsurance problems may be found in Louberge (Louberge 1998; Dionne 2001).
Hawkes process was first introduced in (Hawkes 1971a, 1971b). Good introduction into Hawkes processes and their properties may be found in (Laub et al. 2015). GCHP and regime-switching GCHP were first introduced in (Swishchuk 2017b) and studied in details using real data in (Swishchuk et al. 2019; Swishchuk and He 2019; Swishchuk 2020; Swishchuk and Huffman 2020). Risk model based on GCHP was first introduced in (Swishchuk 2017a) and described in details in (Swishchuk 2018). Applications of the risk model based on GCHP to empirical data and optimal investment problem were considered in (Swishchuk et al. 2020).

3. General Compound Hawkes process

This section contains main definitions and results on one-dimensional Hawkes and general compound Hawkes processes which we will use in our paper. For the completeness, we present them in three subsections below.

3.1. Hawkes Process

Definition 1.
(One-dimensionalHawkes Process) (Hawkes 1971a, 1971b). The one-dimensional Hawkes process is a point process N ( t ) which is characterized by its intensity λ ( t ) with respect to its natural filtration:
λ ( t ) = λ + 0 t μ ( t s ) d N ( s ) ,
where λ > 0 , and the response function or self-exciting function μ ( t ) is a positive function and satisfies 0 + μ ( s ) d s : = μ ^ < 1 .
If ( t 1 , t 2 , . . . , t k ) denotes the observed sequence of past arrival times of the point process up to time t , the Hawkes conditional intensity is
λ ( t ) = λ + t k < t μ ( t t k ) .
The function μ ( t ) is sometimes also called the excitation function, and the constant λ is called the background intensit .
To avoid the trivial case, we suppose that μ ( t ) 0 , which is, a homogeneous Poisson process. Therefore, the Hawkes process is a non-Markovian extension of the Poisson process.
We note, that the Hawkes process is a self-exciting simple point process first introduced by A. Hawkes in 1971 (Hawkes 1971a, 1971b). Thus, the future evolution of a self-exciting point process is influenced by the timing of past events.
Except for some very special cases (e.g., exponential self-exiting function μ ( t ) ), the Hawkes process is non-Markovian . In this way, the Hawkes process has a long memory and depends on the entire past history .
Among many applications of the Hawkes process, we mention applications in finance, insurance, neuroscience, seismology, genome analysis, to name a few.
The above equation for λ ( t ) has the following interpretation: the events occur according to an intensity with a background intensity λ which increases by μ ( 0 ) at each new event then decays back to the background intensity value according to the function μ ( t ) .
Therefore, choosing μ ( 0 ) > 0 leads to a jolt in the intensity at each new event. This feature is often called a self-exciting feature: an arrival causes the conditional intensity function λ ( t ) in (1) and (2) to increase then the process is said to be self-exciting.
The following LLN and CLT for HP may be found in (Bacry et al. 2013). The convergences are considered in weak sense for the Skorokhod topology.
LLN for HP (Bacry et al. 2013). Let 0 < μ ^ : = 0 + μ ( s ) d s < 1 . Then
N ( t ) t t + λ 1 μ ^ .
Remark 1.
By LLN N ( t ) λ 1 μ ^ t for large t .
FCLT for HP (Bacry et al. 2013). Under LLN and 0 + s μ ( s ) d s < + conditions
P N ( t ) λ t / ( 1 μ ^ ) λ t / ( 1 μ ^ ) 3 < y t + Φ ( y ) ,
where Φ ( · ) is the c.d.f. of the standard normal distribution.
Remark 2.
By FCLT N ( t ) λ 1 μ ^ t + λ / ( 1 μ ^ ) 3 ) W ( t ) for large t , where W ( t ) is a standard Wiener process (see Bacry et al. 2013).
Remarks 1 and 2 above give the ideas about the averaged and diffusion approximated HP on a large time interval.

3.2. General Compound Hawkes Process

Definition 2.
(General Compound Hawkes Process). General compound Hawkes Process is defined as (Swishchuk 2020; Swishchuk and Huffman 2020; Swishchuk 2017b)
S ( t ) = S ( 0 ) + i = 1 N ( t ) a ( X i ) .
Here, X i is a discrete-time finite or infinite state Markov chain with state space X = { 1 , 2 , . . . , N } , or X = { 1 , 2 , . . . , N , . . . } , , respectively, a ( x ) is a continuous and bounded function on X , and N ( t ) is a Hawkes process with intensity λ ( t ) > 0 , independent of X i .
This general model is rich enough to:
• incorporate non-exponential distribution of inter-arrival times of orders in HFT or claims in insurance (hidden in N ( t ) )
• incorporate the dependence of orders or claims (via MC X i )
• incorporate clustering of of orders in HFT or claims (properties of N ( t ) )
• incorporate order or claim price changes different from one single number (in a ( X i ) ).
This model is also very general to include:
-in finance:
• compound Poisson process: S t = S 0 + k = 1 N ( t ) X k , where N ( t ) is a Poisson process and a ( X k ) = X k are i.i.d.r.v.
• compound Hawkes process (Swishchuk et al. 2019): S t = S 0 + k = 1 N ( t ) X k , where N ( t ) is a Hawkes process and a ( X k ) = X k are i.i.d.r.v.
• compound Markov renewal process: S t = S 0 + k = 1 N ( t ) a ( X k ) , where N ( t ) is a renewal process and X k is a Markov chain;
-in insurance:
• classical Cramer-Lundberg model: a ( X i ) = X i , X i are i.i.d.r.v., and μ ( t ) = 0 (then N ( t ) is a poisson process);
• Sparre-Andersen model: a ( X i ) = X i , X i are i.i.d.r.v., μ ( t ) = 0 , and N ( t ) is a renewal process;
• Markov-modulated model: a ( X i ) = X i , X i are i.i.d.r.v., λ ( t ) = λ ( X ( t ) ) , where X ( t ) is a MC; we call this model regime-switching risk model based on GCHP (Swishchuk 2020, 2017b).

3.3. LLN and FCLT for GCHP

Lemma 1. (LLN for GCHP)
(Swishchuk 2020; Swishchuk and Huffman 2020; Swishchuk 2017b). Let μ ^ : = 0 + μ ( s ) d s < 1 , and Markov chain X i is ergodic with stationary probabilities π i * . Then the GCHP S n t satisfies the following weak convergence in the Skorokhod topology:
S ( n t ) n n + a * · λ 1 μ ^ t ,
or
S ( t ) t t + a * · λ 1 μ ^ .
Here: a * is defined as a * : = i X π i * a ( i ) , where π i are ergodic probabilities for Markov chain X i .
Theorem 1. (FCLT (or Jump-Diffusion Limit) for GCHP)
(Swishchuk 2020; Swishchuk and Huffman 2020; Swishchuk 2017b). Let X k be an ergodic Markov chain and with ergodic probabilities ( π 1 * , π 2 * , . . . , π n * ) . Let also S t be LGCHP, and 0 < μ ^ : = 0 + μ ( s ) d s < 1 a n d 0 + μ ( s ) s d s < + .
Then
S ( n t ) N ( n t ) · a * n n + σ λ / ( 1 μ ^ ) W ( t ) ,
in weak sense for the Skorokhod topology, where W ( t ) is a standard Wiener process, σ is defined as:
( σ ) 2 : = i X π i * v ( i )
v ( i ) = b ( i ) 2 + j X ( g ( j ) g ( i ) ) 2 P ( i , j ) 2 b ( i ) j X ( g ( j ) g ( i ) ) P ( i , j ) , b = ( b ( 1 ) , b ( 2 ) , . . . , b ( n ) ) , b ( i ) : = a ( i ) a * , g : = ( P + Π * I ) 1 b ,
P is a transition probability matrix for X k , , i.e., P ( i , j ) = P ( X k + 1 = j | X k = i ) , Π * denotes the matrix of stationary distributions of P , and g ( j ) is the jth entry of g . The expressions above are valid for both finite and infinite state Markov chain, that is why we used the notation j X , where X is the states space of a MC X i .
Remark 3.
The formulas for a * and σ look much simpler in the case of two-state Markov chain X i = { δ , + δ } :
a * : = δ ( 2 π * 1 ) a n d ( σ * ) 2 : = 4 δ 2 1 p + π * ( p p ) ( p + p 2 ) 2 π * ( 1 π * ) ,
( p , p ) are transition probabilities of Markov chain X k , and π 1 * = π * , π 2 * = 1 π * .
From FCLT for HP, Section 3, and from Theorem 1 above follow the following FCLT for GCHP (pure jump diffusion limit).
Theorem 2. (FCLT (or Pure Diffusion Limit) for GCHP
(Swishchuk 2020; Swishchuk and Huffman 2020; Swishchuk 2017b; Swishchuk et al. 2020). Let X k be an ergodic Markov chain and with ergodic probabilities ( π 1 * , π 2 * , . . . , π n * ) . Let also S t be LGCHP, and 0 < μ ^ : = 0 + μ ( s ) d s < 1 a n d 0 + μ ( s ) s d s < + .
Then
S ( t ) a * λ 1 μ ^ t t t + σ ¯ N ( 0 , 1 ) ,
in weak sense for the Skorokhod topology, where N ( 0 , 1 ) is the standard normal c.d.f., and σ ¯ is defined as:
σ ¯ = ( σ * ) 2 + a * λ ( 1 μ ^ ) 3 2 ,
where σ * = σ λ / ( 1 μ ^ ) , σ and a * are defined in Theorem 1 and Lemma above, respectively.
Remark 4.
From Theorem 2 it follows that S ( t ) can be approximated by the pure diffusion process:
S ( t ) S ( 0 ) + a * λ 1 μ ^ t + σ ¯ W ( t ) ,
where W ( t ) is a standard Wiener process. This Remark 4 gives the idea about the pure diffusion approximation of GCHP on a large time interval.
Remark 5.
We note, that the rate of convergence in the Theorem 2 is C ( T ) / t , 0 t T , where C ( T ) > 0 is a constant (Swishchuk et al. 2020). Thus, the error of approximation for S ( t ) in Remark 4 is small for large t .

4. Merton Investment Problem in Finance for the Hawkes-Based Model

Let us consider Merton portfolio optimization problem. We suppose that B t and S t follows the following dynamics, respectively:
B t = B 0 exp { r t } S t = S 0 exp { G ( t ) } ,
where G ( t ) : = i = 1 N ( t ) a ( X i ) is the GCHP, X i a discrete-time finite or infinite state Markov chain with state space X = { 1 , 2 , 3 , . . . , N } or X = { 1 , 2 , . . . , N , . . . } , , respectively, r > 0 is the interest rate.
We note, that the justification of using HP in finance may be found in (Da Fonseca and Zaatour 2013), and using GCHP that based on HP N ( t ) and a ( X i ) may be found in (Swishchuk and He 2019; Swishchuk and Huffman 2020). The model for a stock price S ( t ) in (1) that based on GCHP is a new and original in this paper.
The investor starts with an initial amount of money, say X 0 = x , and wishes to decide how much money to invest in risky and risk-less assets to maximize the expected utility of the terminal wealth X t at the maturity T , , i.e., X T .
We denote by n ( t ) : = ( n B ( t ) , n S ( t ) ) an investor portfolio, where n B ( t ) and n S ( t ) are the amounts in cash invested in the bonds and the risky assets, respectively. The value X ( t ) at time t of such portfolio is
X ( t ) = n B ( t ) + n S ( t ) .
We suppose that our portfolio is admissible, i.e., X ( t ) 0 , a.s., 0 t T , and self-financing, i.e.,
d X ( t ) = n B ( t ) d B ( t ) B ( t ) + n S ( t ) d S ( t ) S ( t ) .
Suppose that G ( t ) : = i = 1 N ( t ) a ( X i ) follows FCLT when t + , (Swishchuk 2018; Swishchuk et al. 2020), thus G ( t ) can be approximated as (see Section 3, Remark 4)
G ( t ) a * λ 1 μ ^ t + σ ¯ W ( t ) ,
where μ ^ : = 0 + μ ( s ) d s < 1 , a * is a average of a ( x ) over stationary distribution of MC X i , λ is a background intensity, σ ¯ > 0 is defined in Section 3, Theorem 2. For exponential decaying intensity μ ^ = α / β .
Thus, S ( t ) in (1) can be presented in the following way using (2):
S ( t ) = S ( 0 ) e a * λ 1 μ ^ t + σ ¯ W ( t ) .
Using I t o ^ formula we can get from (3):
d S ( t ) = S ( t ) [ ( a * λ 1 μ ^ + σ ¯ 2 2 ) d t + σ ¯ d W ( t ) ] .
Then the change of the wealth process X t can be rewritten in the following way, taking into account (1)–(4):
d X t = n B ( t ) d B ( t ) B ( t ) + n S ( t ) d S ( t ) S ( t ) = n B ( t ) r d t + n S ( t ) [ ( a * λ 1 μ ^ + σ ¯ 2 2 ) d t + σ ¯ d W ( t ) ] = r X ( t ) d t + n S ( t ) ( a * λ 1 μ ^ + σ ¯ 2 2 r ) d t + n S ( t ) σ ¯ d W ( t ) .
Let π ( t ) : = n S ( t ) / X ( t ) be the portion of wealth invested in the assets/stocks at time t .
Then, from (1)–(5), we have the following expression for d X ( t ) :
d X t = r X ( t ) d t + n S ( t ) ( a * λ 1 μ ^ + σ ¯ 2 2 r ) d t + n S ( t ) σ ¯ d W ( t ) . = d X ( t ) = X ( t ) [ ( r + π t ( a * λ 1 μ ^ + σ ¯ 2 2 r ) ) d t + π t σ ¯ d W ( t ) ) ) ]
Finally, after replacing X ( t ) with X π ( t ) , to stress the dependence of X ( t ) on π t , from (6) we have the following equation for d X π ( t ) :
d X π ( t ) = X π ( t ) [ ( r + π ( t ) ( a * λ 1 μ ^ + σ ¯ 2 2 r ) ) d t + π ( t ) σ ¯ d W ( t ) ) ) ] .
Our main goal is to solve the following optimization problem:
max π E [ U ( X T π ) | X 0 = x ] ,
meaning to maximize the wealth/value function or performance criterion
E [ U ( X T π ) | X 0 = x ] , where U ( x ) is a utility function.
To find optimal π , we follow the standard procedure in this case (see Bjork 2009; Karatzas and Shreve 1998). For the utility function we take the logarithmic one, U ( x ) = log ( x ) . Therefore, we have to maximize max π E [ log ( X T π ) | X 0 = x ] , . Solving the Equation (7) and maximizing non-martingale term in the exponent for the solution, we can find the optimal investment solution π * ( t ) :
π * ( t ) = a * λ 1 μ ^ + σ ¯ 2 2 r σ ¯ 2 ,
where
σ ¯ = ( σ * ) 2 + a * λ ( 1 μ ^ ) 3 2 ,
and σ * and a * are defined in Section 3, Theorem 2.
Thus, we have arrived to the following proposition:
Proposition 1.
Let the conditions of Theorem 2, Section 3, are satisfied. Then the optimal investment solution for the Merton portfolio optimization problem is presented by π * ( t ) in (8) with σ ¯ in (9).
Remark 6.
(Some Insights into the Results). As we can see from the expression for π * ( t ) , the optimal investment solution depends on all parameters of the Hawkes-based model, namely, Hawkes process’s parameters λ and μ ^ , Markov chain X i and function a ( x ) through a * . For example, if a * increases then π * ( t ) increase, and if σ ¯ increases then π * ( t ) decreases, which follows from (8). The latter is obvious: in a very volatile market we should avoid a risk associated with investing in stocks. Furthermore, also obvious that if r increases then π * ( t ) decreases: it’s better to invest in bonds than in stocks.

5. Merton Investment Problem in Insurance for the Hawkes-Based Risk Model

Let us consider R ( t ) as the risk model based on GCHP, namely,
R ( t ) = u + c t i = 1 N ( t ) a ( X i ) .
Here, X i (claim sizes) is a discrete-time finite or infinite state Markov chain with state space X = { 1 , 2 , 3 , . . . , N } or X = { 1 , 2 , . . . , N , . . . } , , respectively, a ( x ) is a bounded and continuous function on X = { 1 , 2 , . . . , N } - space state for X i , and N ( t ) is a Hawkes process with intensity λ ( t ) > 0 , independent of X i , and satisfying:
λ ( t ) = λ + 0 t μ ( t s ) d N ( t ) .
Here, μ ( t ) is self-exiting function.
We note, that the justification for the Hawkes-based risk model in the form of the above Equation (10) may be found in (Swishchuk et al. 2020).
As long as we will consider optimization for a first insurer, thus we will focus on problems with infinite planning horizon.
Let A ( t ) be an amount invested in a risky asset, and suppose that the price S ( t ) of the risky asset follows GBM, i.e.,
d S ( t ) = S ( t ) ( a d t + b d W ( t ) ) ,
where a is a real constant, b > 0 .
Further, the leftover, R ( t ) A ( t ) > 0 , is invested in a bank account (or bonds) with interest rate r > 0 , thus
d ( R ( t ) A ( t ) ) = r ( R ( t ) A ( t ) ) d t .
By β ( t ) : = A ( t ) / S ( t ) we define the number of assets held at time t . Thus, the position of the insurer has the following evolution:
d R ( t ) = r R ( t ) d t + β ( t ) d S ( t ) r β ( t ) S ( t ) d t + c d t d ( G ( t ) ) .
Therefore, the dynamics for R ( t ) is (taking into account all above equations for S ( t ) , d ( R ( t ) A ( t ) ) and dR(t)):
d R ( t ) = r R ( t ) d t + A ( t ) ( a d t + b d W ( t ) ) r A ( t ) d t + c d t d ( G ( t ) )
Here, G ( t ) = i = 1 N ( t ) a ( X i ) .
Let π ( t ) : = A ( t ) / R ( t ) be the fraction of the total wealth R ( t ) invested in the risky assets.
Then, we can rewrite the equation for d R ( t ) in the following way (we use notation R π ( t ) to stress dependence of R ( t ) from π ):
d R π ( t ) = r R π ( t ) d t + A ( t ) ( a d t + b d W ( t ) ) r A ( t ) d t + c d t d ( G ( t ) ) = r R π ( t ) d t + A ( t ) ( ( a r ) d t + b d W ( t ) ) + c d t d ( G ( t ) ) = R π ( t ) [ ( r + π ( t ) ( a r ) ) d t + π ( t ) b d W ( t ) ] + c d t d ( G ( t ) ) ,
where W ( t ) is a standard Brownian motion.
As for the control at time t we will take the function π ( t ) , , i.e., the fraction of the total wealth R π ( t ) which should be invested in risky assets.
We will show how to find the optimal strategy π ( t ) , which maximizes our expected utility function, E [ U ( R π ( t ) ) | R π ( 0 ) = r ] , where U ( r ) is a special utility function.
We suppose that G ( t ) = i = 1 N ( t ) a ( X i ) follows FCLT when t + , (Swishchuk 2018; Swishchuk et al. 2020), thus G ( t ) can be approximated as (see Section 3, Remark 4)
G ( t ) a * λ 1 μ ^ t + σ ¯ W 1 ( t ) ,
where μ ^ : = 0 + μ ( s ) d s < 1 , a * is a average of a ( x ) over stationary distribution of MC X i , λ is a background intensity, σ ¯ is defined in Section 3. For exponential decaying intensity μ ^ = α / β .
Here, W 1 ( t ) is a Wiener process independent of W ( t ) (the case for correlated W 1 ( t ) with W ( t ) , i.e., such that [ W ( t ) , W 1 ( t ) ] = ρ d t , can be considered as well with some modifications).
We suppose that
c > a * λ 1 μ ^
(safety loading condition (SLC)).
After substituting (12) into (11) for dR(t) we get:
d R π ( t ) = [ R π ( t ) ( r + ( a r ) π ( t ) ) + ( c a * λ 1 μ ^ ) ] d t + R 2 ( t ) b 2 π 2 ( t ) + σ ¯ 2 d W 2 ( t ) ,
where W 2 ( t ) is a standard Wiener process independent of W ( t ) and W 1 ( t ) .
Generator for R π ( t ) in (13) is (here, R π ( 0 ) = x )
A π = [ x ( r + ( a r ) π ( t ) ) + ( c a * λ 1 μ ^ ) ] x + [ ( x 2 b 2 π 2 ( t ) + σ ¯ 2 ) / 2 ] x 2
Thus, we have to maximize E r [ U ( R π ( t ) ) ] , where U ( r ) is a utility function.
The HJB equation has the following form:
r v ( t , r ) + sup π [ A π v ( t , r ) ] = 0 ,
where v ( t , r ) = sup π E r [ U ( R π ( t ) ) ] .
We take the exponential utility function U ( r ) = e p r , p > 0 .
Solving the HJB equation we get the optimal control π ( t ) :
π ( t ) = ( a r ) x p b 2
where x = R π ( 0 ) > 0 , and π ( t ) depends on p . After finding
p = θ + θ 2 + σ ¯ 2 ( a r ) 2 / b 2 σ ¯ 2 ,
where θ : = x r + ( c a * λ 1 μ ^ ) , we can finally find that
π ( t ) = σ ¯ 2 ( a r ) x b ( θ b + θ 2 b 2 + σ ¯ 2 ( a r ) 2 ) ,
where
σ ¯ = ( σ * ) 2 + a * λ ( 1 μ ^ ) 3 2 ,
σ * = σ λ / ( 1 μ ^ ) , σ and a * are defined in Theorem 1 and Lemma, Section 3, respectively.
Thus, we arrived to the following proposition:
Proposition 2.
Let the conditions of Theorem 2, Section 3, are satisfied. Then the optimal investment solution for the Merton investment problem in insurance is presented by π ( t ) in (14) with σ ¯ in (15).
As we can see, the optimal control π ( t ) = π does not depend on t , thus is a constant, and contains all initial parameters of the risk model based on GCHP.
Remark 7.
(Some Insights into the Results). As we can see from the expression for π ( t ) , the optimal control depends not only from interest rate r , but also from all parameters of the Hawkes-based model, namely, Hawkes process’s parameters λ and μ ^ , Markov chain X i and function a ( x ) through a * , and the asset’s parameters a and b . For example, from (14) we can see that π ( t ) decreases if σ ¯ increases, where σ ¯ is, as we could call it, a ’volatility of all volatilities’, and thus it is not a good idea to invest in stocks in highly volatile market; also, if, for example, self-exiting function μ ( t ) = β exp ( α t ) is exponential, then we can see from (14) that π ( t ) depends on parameters α and β in the following way (here: μ ^ = α / β ): if α increases then π ( t ) is also increases, and when β increases then π ( t ) decreases. Furthermore, π ( t ) has a significant dependence on parameter θ , as we could call it, ’rate of increase of company’s capital’, θ : = x r + ( c a * λ 1 μ ^ ) , ( θ > 0 , which follows from SLC): if θ increases, then π ( t ) decreases (see (14)). The latter is understandable: if the capital of a company increases due to interest r and premium c , then there is no need to take a risk investing in stocks, and as a result, π ( t ) should be decreased.
Corollary 1.
Merton Investment Problem for Poisson-based Risk Model in Insurance.
The optimal control π P ( t ) for Poisson Risk Model (in this case, N ( t ) is a Poisson process, a ( x ) = x , X i are i.i.d.r.vs) is (which follows from (14)):
π P ( t ) = λ E [ X i ] 2 ( a r ) x b ( b ( x r + ( c λ E X i ) ) + b 2 ( x r + ( c λ E X i ) ) + λ E [ X i ] 2 ( a r ) 2 ) ) .
Here: σ ¯ = λ E [ X i ] 2 .

6. Discussion

Even though there is the advantage of the solution of Merton problems in finance and insurance that associated with the approximation of the initial complicated model by much simpler diffusion model, there are some limitations related to proposed portfolio design and to improving the proposed approximation method. Thus, the future work will be devoted to numerical example based on real data, simulations and considering the case without diffusion approximation for R ( t ) , by creating HJB equation for initial risk model R ( t ) = R ( 0 ) + c t i = 1 N ( t ) a ( X i ) and by solving this HJB equation for initial risk model base on GCHP. Probably we cannot avoid here numerical methods, because the HJB equation in this case cannot be solved exactly with a close form solution.

7. Conclusions

We described in this paper how to solve Merton optimal investment stochastic control problem for Hawkes-based models in finance and insurance (Propositions 1 and 2), i.e., for a wealth portfolio X ( t ) consisting of a bond and a stock price described by general compound Hawkes process (GCHP), and for a capital R ( t ) of an insurance company with the amount of claims described by the risk model based on GCHP. The novelty of the results consists of the new Hawkes-based models and in the new optimal investment results in finance and insurance for those models. Specifically: (1) we considered a new model for the stock price S ( t ) in the form of S t = S 0 exp { G ( t ) } , where G ( t ) : = i = 1 N ( t ) a ( X i ) is the GCHP; we call it Hawkes-based model for the stock price (or geometric general compound Hawkes process, similar to the geometric Brownian motion); (2) solution of Merton investment problem for this Hawkes-based model; (3) solution of Merton investment problem for the Hawkes-based risk model. We also gave some insights into the obtained results (see Remarks 6 and 7).

Funding

This research received no external funding.

Acknowledgments

The author thanks to NSERC for continuing support, and also to three anonymous referees for very valuable and constructive comments, remarks and suggestions. All remaining errors are mine.

Conflicts of Interest

The authors declare no conflict of interest.

Abbreviations

The following abbreviations are used in this manuscript:
LLNLaw of Large Numbers
FCLTFunctional Central Limit Theorem
GCHPGeneral Compound Hawkes Process
HJBHamilton–Jaocobi–Bellman Equation

References

  1. Arrow, Kenneth J. 1963. Uncertainty and the welfare economics of medical care. American Economic Review 53: 941–73. [Google Scholar]
  2. Asmussen, Søren, and Michael Taksar. 1997. Controlled diffusion models for optimal dividend payout. Insurance: Mathematics and Economics 20: 1–15. [Google Scholar]
  3. Bacry, Emmanuel, Sylvain Delattre, Marc Hoffmann, and Jean-François Muzy. 2013. Some limit theorems for Hawkes processes and application to financial statistics. Stochastic Processes and Their Applications 123: 2475–99. [Google Scholar] [CrossRef]
  4. Brockett, Patrick, and Xia Xiaohua. 1995. Operations research in insurance: A review. Transactions of Society of Actuaries XLVII: 7–80. [Google Scholar]
  5. Browne, Sid. 1995. Optimal investment policies for a firm with a random risk process: Exponential utility and minimizing the probability of ruin. Mathematics of Operations Research 20: 937–58. [Google Scholar] [CrossRef] [Green Version]
  6. Bühlmann, Hans. 1970. Mathematical Methods in Risk Theory. New York: Springer. [Google Scholar]
  7. Bjork, Tomas. 2009. Arbitrage Theory in Continuous Time, 3rd ed. Oxford: Oxford University Press. [Google Scholar]
  8. Cox, John C., and Chi-fu Huang. 1989. Optimal consumption and portfolio poly- cies when asset prices follow a diffusion process. Journal of Economic Theory 49: 33–83. [Google Scholar] [CrossRef] [Green Version]
  9. Da Fonseca, José, and Riadh Zaatour. 2013. Hawkes Process: Fast Calibration, Application to Trade Clustering, and Diffusive Limit. Journal of Futures Markets 34: 548–79. [Google Scholar] [CrossRef]
  10. Dionne, Georges, ed. 2001. Handbook of Insurance. Dordrecht: Kluwer Academic Publishers. [Google Scholar]
  11. Dufresne, Francois, and Hans U. Gerber. 1991. Risk theory for the compound Poisson process that is perturbed by diffusion. Insurance: Mathematics and Economics 10: 51–59. [Google Scholar] [CrossRef]
  12. Ehrlich, Isaac, and Gary S. Becker. 1972. Market insurance, self-insurance and self protection. Journal of Political Economy 80: 623–48. [Google Scholar] [CrossRef]
  13. Fleming, Wendell H., and Raymond W. Rishel. 1975. Deterministic and Stochastic Optimal Control. New York: Springer. [Google Scholar]
  14. Fleming, Wendell H., and Halil Mete Soner. 1993. Controlled Markov Processes and Viscosity Solutions. New York: Springer. [Google Scholar]
  15. Gaier, Johanna, Peter Grandits, and Walter Schachermayer. 2003. Asymptotic ruin probabilities and optimal investment. Annals of Applied Probability 13: 1054–76. [Google Scholar] [CrossRef]
  16. Gerber, Hans. 1979. An Introduction to Mathematical Risk Theory. Philadelphia: S.S. Huebner Foundation Monographs, University of Pennsylvania. [Google Scholar]
  17. Grandits, Peter. 2003. An analogue of the Cramer-Lundberg approximation in the optimal investment case. Applied Mathematics and Optimization 50: 1–20. [Google Scholar]
  18. Grandell, Jan. 1990. Aspects of Risk Theory. New York: Springer. [Google Scholar]
  19. Hawkes, Alan G. 1971a. Spectra of some self-exciting and mutually-exciting point processes. Biometrika 58: 83–90. [Google Scholar] [CrossRef]
  20. Hawkes, Alan G. 1971b. Point spectra of some mutually-exciting point processes. Journal of the Royal Statistical Society: Series B (Methodological) 33: 438–43. [Google Scholar] [CrossRef]
  21. Hipp, Christian 2004. Stochasric optimal control in insurance. In Stochastic Methods in Finance. Frittelli, M., and W Runggaldier, eds. Lecture Notes in Math. Berlin/Heidelberg: Springer, pp. 127–64.
  22. Højgaard, Bjarne, and Michael Taksar. 1999. Controlling risk exposure and div- idends pay-out schemes: Insurance company example. Mathematical Finance 2: 153–82. [Google Scholar] [CrossRef]
  23. Karatzas, Ioannis, John P. Lehoczky, Suresh P. Sethi, and Steven E. Shreve. 1986. Explicit solution of a general consumption/investment problem. Mathematics of Operations Research 11: 261–94. [Google Scholar] [CrossRef]
  24. Karatzas, Ioannis. 1998. Lectures of Mathematics of Finance. Providence: American Mathematics Society, p. 8. [Google Scholar]
  25. Karatzas, Ioannis, and Steven Shreve. 1998. Probability Theory and Mathematical Modelling. In Methods of Mathematical Finance. Berlin/Heidelberg: Springer, vol. 39. [Google Scholar]
  26. Korn, Ralf, and Elke Korn. 2000. Graduate Studies in Mathematics. In Option Pricing and Portfolio Optimization. Providence: American Mathematical Society, vol. 31. [Google Scholar]
  27. Laub, Patrick J., Thomas Taimre, and Philip K. Pollett. 2015. Hawkes Processes. arXiv arXiv:1507.02822v1. [Google Scholar]
  28. Louberge, Henri. 1998. Risk and insurance economics 25 years after. The Geneva Papers on Risk and Insurance-Issues and Practice 23: 1973–98. [Google Scholar] [CrossRef] [Green Version]
  29. Martin-Löf, Anders. 1994. Lectures on the use of control theory in insurance. Scandinavian Actuarial Journal, 1–25. [Google Scholar] [CrossRef]
  30. Merton, Robert. 1969. Lifetime portfolio selection under uncertainty: The continuous-time case. The Review of Economics and Statistics 51: 247–57. [Google Scholar] [CrossRef] [Green Version]
  31. Merton, Robert. 1971. Optimum consumption and portfolio rules in a continuous-time model. Journal of Economic Theory 3: 373–413. [Google Scholar] [CrossRef] [Green Version]
  32. Merton, Robert. 1990. Continuous-Time Finance. Cambridge: Basil Blackwell Inc. [Google Scholar]
  33. Pliska, Stanley. 1986. Stochastic calculus model of continuous trading: Optimal portfolios. Mathematics of Operation Research 11: 371–82. [Google Scholar] [CrossRef] [Green Version]
  34. Rolski, Tomasz, Hanspeter Schmidli, Volker Schmidt, and Jozef L. Teugels. 1998. Stochastic Processes for Insurance and Finance. Wiley Series in Probability and Statistics; Hoboken: Wiley. [Google Scholar]
  35. Schmidli, Hanspeter. 2008. Stochastic Control in Insurance. Berlin/Heidelberg: Springer. [Google Scholar]
  36. Swishchuk, Anatoliy. 2017a. Risk processes based on general compound Hawkes processes. Paper presented at IME 2017 Congress, Vienna, Austria, July 3–5; Available online: https://fam.tuwien.ac.at/events/ime2017/IME2017_abstracts.pdf (accessed on 26 May 2021).
  37. Swishchuk, Anatoliy. 2017b. General compound Hawkes processes in LOB. arXiv arXiv:1706.07459v2. [Google Scholar]
  38. Swishchuk, Anatoliy. 2018. Risk Model Based on Compound Hawkes Process. Wilmott 2018. [Google Scholar] [CrossRef]
  39. Swishchuk, Anatoliy, Remillard Bruno, Elliott Robert, and Chavez-Casillas Jonathan. 2019. Compound Hawkes processes in limit order books. In Financial Mathematics, Volatility and Covariance Modelling, 1st ed. Edited by Chevallier Julien, Goutte Stephane, Guerreiro David, Saglio Sophie and Sanhaji Bilel. Routledge: Taylor and Francis Group, vol. 2. [Google Scholar]
  40. Swishchuk, Anatoliy, and Qiyue He. 2019. Quantitative and Comparative Analyses of Limit Order Books with General Compound Hawkes Processes. Risks 7: 110. [Google Scholar]
  41. Swishchuk, Anatoliy. 2020. Modelling of limit order books by general compound Hawkes processes with implementations. Methodology and Computing in Applied Probability 23: 299–428. [Google Scholar]
  42. Swishchuk, Anatoliy, and Aiden Huffman. 2020. General compound Hawkes processes in limit order books. Risks 8: 28. [Google Scholar] [CrossRef] [Green Version]
  43. Swishchuk, Anatoliy, Rudi Zagst, and Gabriela Zeller. 2020. Hawkes processes in insurance: Risk model, application to empirical data and optimal investment. Insurance: Mathematics and Economics. [Google Scholar] [CrossRef]
  44. Taksar, Michael. 2000. Optimal Risk/Dividend Distribution Control Models: Applications to Insurance. Mathematical Methods of Operations Research 1: 1–42. [Google Scholar] [CrossRef]
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Swishchuk, A. Merton Investment Problems in Finance and Insurance for the Hawkes-Based Models. Risks 2021, 9, 108. https://doi.org/10.3390/risks9060108

AMA Style

Swishchuk A. Merton Investment Problems in Finance and Insurance for the Hawkes-Based Models. Risks. 2021; 9(6):108. https://doi.org/10.3390/risks9060108

Chicago/Turabian Style

Swishchuk, Anatoliy. 2021. "Merton Investment Problems in Finance and Insurance for the Hawkes-Based Models" Risks 9, no. 6: 108. https://doi.org/10.3390/risks9060108

APA Style

Swishchuk, A. (2021). Merton Investment Problems in Finance and Insurance for the Hawkes-Based Models. Risks, 9(6), 108. https://doi.org/10.3390/risks9060108

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