1. Introduction
Consider an insurance company that invests an asset in a financial market. Such a situation carries two kinds of risk, as summarized by [Reference Norberg24] and [Reference Tang and Tsitsiashvili32]. One kind, called insurance risk in the literature, stems from insurance claims, and the other, called financial risk, stems from financial investments. These two types of risk should be dealt with carefully in conducting solvency assessments of insurance companies. In this paper, we use multidimensional risk models to accommodate the two types of risks in the discrete-time and continuous-time cases.
Major events, such as the COVID-19 crisis and natural catastrophes, significantly affect almost all economic sectors. This applies in particular to financial markets and the insurance industry, implying that it makes sense to consider interdependencies between financial and insurance risks. In this paper, we propose to describe the dependence structure using Assumption 1 within a framework of general multivariate regular variation (MRV). Under certain conditions, multiple dependence structures satisfy Assumption 1, such as those of Proposition 2 and the relation (14) in this paper. Moreover, Assumption 1 includes scale mixtures, which have a variety of interpretations in different applications. For instance, they can reflect both individual and common factors of a risk class, and they can describe loss variables. However, the implications of Assumption 1 are not confined to this, and one can find available information after carefully studying Assumption 1.
Since the pioneering work of [Reference Norberg24] and [Reference Tang and Tsitsiashvili32], various studies have been performed on discrete-time insurance risk models with financial and insurance risks. Some recent works include [Reference Chen6], [Reference Chen and Yuan9], [Reference Li and Tang21], [Reference Tang and Yang29], and [Reference Yang and Konstantinides33]. All of these works restrict attention to an insurance company with one business line. However, special attention is also given to those insurance companies with multiple business lines. In Section 4 of this paper, we introduce the multidimensional discrete-time risk model (16) with a dependence structure between financial and insurance risks satisfying Assumption 1. We then obtain a precise asymptotic expansion for the ruin probability that holds uniformly for the whole time horizon. Furthermore, we present a two-dimensional numerical example. The contribution of this numerical example is twofold. First, it provides a specific dependence structure that satisfies Assumption 1. Second, it shows the accuracy of the asymptotic expansions for the ruin probability.
For the continuous-time risk model with financial and insurance risks, there is abundant literature on the asymptotic estimation of ruin probabilities; see [Reference Asmussen and Albrecher3], [Reference Cheng, Konstantinides and Wang10], [Reference Heyde and Wang12], [15], [Reference Li19], [Reference Paulsen26], and [Reference Tang, Wang and Yuen28]. However, in almost all of these papers, it is assumed that financial and insurance risks are independent. There is no doubt that, in a similar or the same macroeconomic environment, it is unrealistic to assume that financial and insurance risks over the same period are independent. In Section 5 of this paper, under Assumption 1, we construct a multidimensional continuous-time risk model with a dependence structure between insurance and financial risks. Three ruin probabilities in the infinite-time horizon are then investigated.
The rest of the paper consists of four parts. Section 2 describes regular variation and MRV and presents some properties of MRV. Section 3 introduces Assumption 1. Under Assumption 1, in Section 4 we study the d-dimensional discrete-time risk model (16) (for $d\geq1$ ) and construct a two-dimensional numerical example satisfying Assumption 1. In Section 5 we consider a d-dimensional continuous-time risk model under Assumption 1.
2. Preliminaries
If there exist two positive functions $g(\cdot,\cdot)$ and $f(\cdot,\cdot)$ satisfying
then we say that $f(x,t)\lesssim g(x,t)$ holds if $l_2\leq1$ ; $f(x,t)\gtrsim g(x,t)$ holds if $l_1\geq1$ ; $f(x,t)\sim g(x,t)$ holds if $l_1=l_2=1$ ; and $f(x,t)\asymp g(x,t)$ holds if $0<l_1 \leq l_2< \infty$ . All limit relations, unless otherwise stated, hold as $x\rightarrow\infty$ . Throughout this paper, C denotes a generic positive constant without relation to x, which may vary with the context. Moreover, for any $x,y\in\mathbb{R}$ , we write $x\wedge y=\min\left\{x,y\right\}$ , $x\vee y=\max\left\{x,y\right\}$ , and $x^{+}=x\vee 0$ .
In this paper, random and real vectors, supposed to be d-dimensional, are expressed by bold English letters. For two real vectors $\boldsymbol{a}=(a_1, \ldots,a_d)$ and $\boldsymbol{b}=(b_1,\ldots,b_d)$ , we assume that operations between vectors such as $\boldsymbol{a}>\boldsymbol{b}$ , $\boldsymbol{a}\pm\boldsymbol{b}$ , $\boldsymbol{a}\boldsymbol{b}$ , and $\boldsymbol{a}^{-1}$ should be interpreted componentwise, and we let $[\boldsymbol{a},\boldsymbol{b}]=[a_1,b_1]\times\cdots\times[a_d,b_d]$ , $[\boldsymbol{a},\boldsymbol{\infty})=[a_1,\infty)\times\cdots\times[a_d,\infty)$ . Additionally, for $y\in\mathbb{R}$ , write $y\boldsymbol{a}=\left(ya_1,\ldots,ya_d\right)$ as usual. We also write $\mathbf{0}=(0,\ldots,0)$ , $\mathbf{1}=(1,\ldots,1)$ , and $\mathbf{e}_{k}=(0,\ldots,1,\ldots,0)$ , where the value 1 occurs in the kth spot.
Definition 1. A distribution F on $[0,\infty)$ satisfying $\bar{F}(x)>0$ for all $x\geq0$ is said to be of regular variation, and we write $F\in\mathcal{R}_{-\alpha}$ , if
for some $0<\alpha<\infty$ .
For some distribution function with $F\in\mathcal{R}_{-\alpha}$ , $0<\alpha<\infty$ , we get by [Reference Bingham, Goldie and Teugels5, Proposition 2.2.3] that there exist positive constants $C_F>1$ and $D_F$ such that, for all $x, xy\geq D_F$ and any $0<\varepsilon<\alpha$ ,
From the relation (2), it follows that for $p>\alpha$ ,
Definition 2. An $\mathbb{R}_+^d$ -valued random vector $\boldsymbol{Z}$ is said to follow a multivariate regularly varying distribution if there exist a Radon measure $\nu$ not identically 0 and a positive normalizing function $a(\!\cdot\!)$ monotonically increasing to $\infty$ , such that when $t\rightarrow\infty$ ,
where $\stackrel{v}{\longrightarrow}$ denotes vague convergence.
From the above definition, one can show that the Radon measure $\nu$ is homogeneous—namely, that there exists some $\alpha>0$ , denoting the MRV index, such that the equality
holds for every Borel set $B\subset[0,\infty]^d\backslash \{\textbf{0}\}$ . For details on this, see [Reference Resnick27, p. 178]. From [Reference Konstantinides17, p. 459], we obtain $a(t)\in\mathcal{R}_{1/\alpha}$ , and hence there exists some distribution $F\in\mathcal{R}_{-\alpha}$ such that
where $a^{\leftarrow}(t)$ is the generalized inverse function of a(t). Consequently, the relation (4) can be expressed as follows:
Hence, we write $\boldsymbol{Z}\in \textrm{MRV}(\alpha,F,\nu)$ . For more details on MRV, the reader is referred to [Reference Resnick27, Chapter 6] or [Reference Konstantinides17, Chapter 13].
Next we present a lemma that is important in the following proofs.
Lemma 1. Let random vector $\boldsymbol{Z}\in \boldsymbol{MRV}(\alpha,F,\nu)$ for some $\alpha>0$ , let $\boldsymbol{\xi}$ be a positive random vector with arbitrarily dependent components satisfying $\mathbb{E}\boldsymbol{\xi}^{\beta}<\boldsymbol{\infty}$ for some $\beta>\alpha$ , and let $\left\{\Delta_t,t\in\mathcal{T}\right\}$ be a set of random events such that $\lim_{t\rightarrow t_0}\mathbb{P}(\Delta_t)=0$ for some $t_0$ in the closure of the index set $\mathcal{T}$ . Furthermore, assume that $\left\{\boldsymbol{\xi},\left\{\Delta_t,t\in\mathcal{T}\right\}\right\}$ and $\boldsymbol{Z}$ are independent. Then
Proof. Since the random vector $\boldsymbol{Z}\in \textrm{MRV}(\alpha,F,\nu)$ , the marginal tail of $\boldsymbol{Z}$ is regularly varying from [Reference Konstantinides17, p. 458]. Therefore, we have $Z_k\in\mathcal{R}_{-\alpha}$ , $1\leq k\leq d$ . Then
where we used [Reference Tang and Yuan30, Lemma 6.2]. This completes the proof.
3. Main assumption
This paper proposes the following assumption to describe certain dependence structures between the random vectors $\mathbf{U}$ and $\mathbf{V}$ .
Assumption 1. Let $\left\{\boldsymbol{U}_i\boldsymbol{V}_i =(U_{1i}V_{1i}, \,\ldots, \, U_{di}V_{di} ), i\in\mathbb{N} \right\}$ , $d\geq1$ , be a sequence of independent and identically distributed (i.i.d.) nonnegative random vectors with generic vector $\boldsymbol{U}\boldsymbol{V} =(U_1V_1, \,\ldots, \, U_dV_d)\in\boldsymbol{MRV}(\alpha,F,\nu)$ such that $\nu\left((\textbf{1},\boldsymbol{\infty})\right)>0$ .
Remark 1. From Assumption 1, we can derive that
and
The relation (7) indicates that $U_1V_1,\ldots,U_dV_d$ reveal large joint movements and accordingly are asymptotically dependent. The relation (8) indicates that the tails of the products $U_1V_1,\ldots,U_dV_d$ are regularly varying and that $U_iV_i$ and $U_jV_j$ have comparable tails. Furthermore, note that by (5) and (7), for any $(b_1, \ldots, b_d)\in[0,\infty]^d\backslash \{\textbf{0}\}$ ,
Remark 2. Set $\mathbf{U}=U\mathbf{1}$ ; then the product $\mathbf{UV}$ can be rewritten as
which is called the scale mixture. This structure can be applied in various areas, including investments or insurance portfolios, and it has various interpretations in different applications. For example, [Reference Li and Sun18] developed a method to estimate the tail dependence of heavy-tailed scale mixtures of multivariate distributions. The paper [Reference Zhu and Li34] studied the asymptotic relations between the value-at-Risk and multivariate tail conditional expectation for heavy-tailed scale mixtures of multivariate distributions. Additionally, the class (9) of loss distributions is a subset of all multivariate regularly varying distributions considered in [Reference Joe and Li14], but it includes various loss distributions. In addition, the scale mixture can reflect both individual factors (via the $V_i$ ) and the common factor (via U) of a risk class. For instance, U is a common systemic risk factor associated with macroeconomic conditions, including regulatory bodies and economic conditions, while the quantities $V_k$ , $1\leq k\leq d$ , are individual risks explained as individual business risks and assets.
Remark 3. Since the structure of $\mathbf{UV}\in\textrm{MRV}(\alpha, F, \nu)$ in Assumption 1 makes no dependence assumption between random vectors $\mathbf{U}$ and $\mathbf{V}$ , it allows for a wide range of dependence structures between $\mathbf{U}$ and $\mathbf{V}$ , and many dependence structures satisfy Assumption 1 under certain conditions (such as those of Proposition 2 and the relation (14)). This further enhances the practical and theoretical interest of Assumption 1.
For Assumption 1, the following question arises naturally: what are the conditions under which the distribution of $\mathbf{UV}$ is MRV? This question has received an increasing amount of attention. The paper [Reference Basrak, Davis and Mikosch4] demonstrated the MRV of the product $\boldsymbol{UV}$ , when $\mathbf{U}$ is MRV and independent of $\mathbf{V}$ , as in Proposition 1 (see also [Reference Basrak, Davis and Mikosch4, Appendix], [Reference Fougeres and Mercadier11, Theorem 1], and [Reference Konstantinides and Li16, Lemma 3.1]).
Proposition 1. Let the random vector $\boldsymbol{U}\in {\boldsymbol{MRV}}(\alpha,F,\nu)$ for some $\alpha>0$ , and let $\boldsymbol{V}$ be a nonnegative random vector with dependent components satisfying $\mathbb{E}\boldsymbol{V}^p<\boldsymbol{\infty}$ for some $p>\alpha$ . Assume that $\boldsymbol{V}$ and $\boldsymbol{U}$ are independent. Then the following relation holds for every Borel set $K\subset[0,\infty]^d\backslash \{\textbf{0}\}$ :
where
The paper [Reference Chen, Wang and Zhang7] studied the tail asymptotics of the nonnegative random loss $\sum_{i=1}^d R_iS$ , where the stand-alone risk vector $\mathbf{R}=(R_1, \, \ldots, \, R_d)$ follows a multivariate regularly varying distributions with index $\alpha>0$ and is independent of S, representing the background risk factor. Furthermore, [Reference Chen, Wang and Zhang7] also assumed that $\mathbb{E}S^{\alpha+\delta}<\infty$ for some $\delta>0$ . Essentially, the conditions in [Reference Chen, Wang and Zhang7] imply that the random vector $(R_1S, \, \ldots, \, R_dS)$ is MRV by Proposition 1.
The hypothesis of independence between $\mathbf{U}$ and $\mathbf{V}$ in Proposition 1 may be too strong in certain settings. If this condition can be weakened, Assumption 1 will become broader and more meaningful. The paper [Reference Fougeres and Mercadier11] further improved Proposition 1 to the following result, Proposition 2, which seems meaningful in the context of actuarial risk theory.
Proposition 2. Let us assume that
on the Borel sets of $([0,\infty]^d\backslash \{\textbf{0}\}) \times [0,\infty]^d$ , where L represents a probability measure on $[0,\infty]^d$ and $\nu$ denotes a Radon measure on $[0,\infty]^d\backslash \{\textbf{0}\}$ not concentrated at $\infty$ . Suppose that, for some $\delta>0$ and any $i=1, \, 2, \, \ldots, \, d$ ,
and also that
where $ \mathbb{I}_{[E]}$ is the indicator of the Borel set E and $|| \cdot ||$ denotes some norm on $\mathbb{R}^d$ . Then the random vector $\mathbf{UV}$ follows a multivariate regularly varying distribution under the relations (10), (11), and (12). More precisely,
as $t\rightarrow\infty$ , where $\nu_L$ is the measure defined on $[0,\infty]^d\backslash \{\textbf{0}\}$ by
Assume that $(U_1, \ldots, U_d, V)$ , $d\geq1$ , is a positive random vector. Li [Reference Li19] introduced the following dependence structure among the components of $(U_1, \ldots, U_d, V)$ :
-
• There are some d-variate function $f\;:\; \, [0,\infty]^d\backslash \{\textbf{0}\} \mapsto (0, \infty)$ , some univariate function $h\;:\; \, (0,\infty) \mapsto (0,\infty)$ satisfying
(13) \begin{align} 0 < \inf_{0 < t < \infty} h(t) \leq \sup_{0 < t < \infty} h(t) <\infty, \end{align}and some distribution F supported on $(0,\infty)$ with an infinite upper endpoint, such that for every $\mathbf{b} \in [0,\infty]^d\backslash \{\textbf{0}\}$ the relation(14) \begin{align} \mathbb{P}(U_1>b_1u,\ldots,U_d>b_du | V = v) \sim f(\mathbf{b}) h(v) \bar{F}(u) \end{align}holds uniformly for $v\in(0,\infty)$ as $u\rightarrow\infty$ .
According to the relation (2.10) in [Reference Li19], taking $\mathbf{b}=\mathbf{e}_{k}$ , $1\leq k\leq d$ , in (14) yields that the relation
holds uniformly for $v\in(0,\infty)$ as $u\rightarrow\infty$ , where $\bar{F}_k(u)$ is the distribution function of $U_k$ , and $\mu=\mathbb{E}[h(V)]\in(0,\infty)$ . Therefore, the dependence structure specified in (14) implies that the marginal vector $(U_k,V)$ follows a parallel two-dimensional dependence structure as shown in (15). The structure (15) is defined in [Reference Asimit and Jones2] and further developed in [Reference Asimit and Badescu1] and [Reference Li, Tang and Wu22]; it includes a wide range of commonly used bivariate copulas.
Now define a new positive random variable $V^*$ with distribution
For simplicity, we introduce a new random variable $\xi$ with distribution function $F\in\mathcal{R}_{-\alpha}$ for $0<\alpha<\infty$ , and let $\xi$ and $V^*$ be independent of all other sources of randomness. If
holds for $\varepsilon>0$ , then by (13),
and for $\frac{\alpha}{\alpha+\varepsilon}<p<1$ and large $u>0$ ,
and
hold. Consequently, we have that by (14),
where $\mathbb{P}(\xi V^*>u)\sim \bar{F}(u)\mathbb{E}[(V^*)^{\alpha}]$ is due to the relation (4.4) in [Reference Tang, Wang and Yuen28]. Therefore, $(U_1V,\ldots,U_dV)$ satisfies Assumption 1.
4. The study of a d-dimensional discrete-time risk model under Assumption 1
This section considers an insurer who runs multiple lines of business and makes risky assets along individual lines. For every $i \in \mathbb{N}=\{1, \, 2, \, \ldots\}$ and integer $d\geq1$ , the real-valued random variable $X_{ki}$ , $1\leq k \leq d$ , denotes the net insurance loss (described by the aggregate claim amount minus the aggregate premium income) of the kth business of the insurer over the period i, and the positive random variable $\theta_{ki}$ denotes the discount factor, related to the return on the investment, of the kth business of the insurer over the same period. Let $\{(X_{1},\ldots,X_{d}), (X_{1i},\ldots,X_{di}), i\in \mathbb{N}\}$ and $\{(\theta_{1},\ldots,\theta_{d}), (\theta_{1i},\ldots,\theta_{di}), i\in \mathbb{N}\}$ be sequences of i.i.d. random vectors. The stochastic discounted value of total net insurance losses for the insurance company up to the time n can be described as
where multiplication over the empty set is understood to be 1, and $Y_{ki}=\prod_{j=1}^{i} \theta_{kj}$ , $1\leq k\leq d$ . The product $\boldsymbol{\rho}x = (\rho_1x, \, \ldots, \rho_dx)$ denotes the vector of initial reserves assigned to different businesses, with positive $\rho_1, \, \ldots, \, \rho_d$ such that $\sum_{k=1}^d \rho_k =1$ . For $n\in\mathbb{N}$ , the ruin probability can be defined as
which denotes the probability of the insurance company’s wealth process going below zero by time n. When $d=1$ , [Reference Fougeres and Mercadier11] considered the ruin probability of (16) if the random vector
is MRV, and [Reference Tang and Yuan31] studied the ruin probability of (16) if the random vector $(X_{1i}, \theta_{1i})$ , $i\in\mathbb{N}$ , follows a bivariate regular variation structure. When $d=2$ , according to [Reference Chen and Yang8], we can obtain that the random vector
still follows a bivariate-regular-variation structure if the pairs $(X_{1i}\theta_{1i}, X_{2i}\theta_{2i})$ , $i\in\mathbb{N}$ , follow some bivariate-regular-variation structure. In this section, we assume that $(X_{1i}^+\theta_{1i}, \ldots, X_{di}^+\theta_{di})$ , $i\in\mathbb{N}$ and $d\geq1$ , is a sequence of i.i.d. random vectors with generic vector $(X_{1}^+\theta_{1}, \ldots, X_{d}^+\theta_{d})\in \textrm{MRV}(\alpha, F,\nu)$ such that $\nu((\textbf{1},\boldsymbol{\infty}))>0$ ; then we study the asymptotic formula for the ruin probability for $n\in\mathbb{N}$ .
Theorem 1. Consider the d-dimensional discrete-time risk model (16). For each $i\in\mathbb{N}$ , assume that $\boldsymbol{X}_{i}^+ \boldsymbol{\theta}_{i}$ satisfies Assumption 1. If there is a constant $\beta>\alpha$ such that $\mathbb{E}\boldsymbol{\theta}^{\beta}<\textbf{1}$ , then it holds uniformly for $n\in\mathbb{N}$ that
4.1. Numerical results
This subsection presents a two-dimensional numerical example to examine the accuracy of Theorem 1. All computations are conducted in R. For simplicity, we assume that $\theta_{1}=\theta_{2}=\theta$ , and let $\theta$ be exponential with parameter $\lambda=2$ .
We first construct a random pair $(X_1, X_2)$ . Let the distribution function F of the random variable Z follow a Pareto distribution, with scale parameter $\kappa=4$ and shape parameter $\alpha=2$ ( $\alpha>0$ ); that is, $F(x) = 1-\left( \frac{\kappa}{\kappa+x} \right)^{\alpha} \in \mathcal{R}_{-\alpha}$ , $x>0$ . Suppose that the dependence structure between Z and $\theta$ is established by a Farlie–Gumbel–Morgenstern (FGM) copula,
Then, by [Reference Asimit and Badescu1, Example 2.2] or [Reference Li, Tang and Wu22, Example 3.2], the random pair $(Z, \theta)$ satisfies
with $h(t)=1+\gamma(1-2e^{-\lambda t})$ . Let $(X_1, X_2)=(\zeta_1Z, \zeta_2Z)$ , where $\zeta_k$ , $k=1,2$ , are uniform on [0,1] and independent of $(Z, \theta)$ , and $\mathbb{E}\zeta_k^q<\infty$ for some $q>\alpha$ .
We next verify that $(X_1\theta, X_2\theta)$ satisfies Assumption 1. Using [Reference Li19, Example 3.1], we obtain
where $b_1,\, b_2>0$ are constants, and
Moreover, for $\beta=3>\alpha$ satisfying $\mathbb{E}\theta^{\beta}=\frac{6}{\lambda^3}<1$ and $\alpha/\beta<p<1$ , we obtain that
where $\mu=\mathbb{E}\theta^{\alpha}+\gamma \mathbb{E}\theta^{\alpha}-2\gamma \mathbb{E}\left(\theta^{\alpha}e^{-\lambda\theta}\right)$ , and in the second and fourth steps we use the relation $\mathbb{P}(\theta>x^p)\leq x^{-p\beta} \mathbb{E}\theta^\beta=o(1)\bar{F}(x)$ from (3). This implies that $\mathbb{P}(X_1\theta>\cdot, X_2\theta>\cdot)$ possesses a bivariate regularly varying tail. Set $\gamma=0.5$ , so $\mu={11}/{16}$ and $\nu((1,\infty]\times(1,\infty])=\mu V(1,1)={11}/{96}>0$ ; then $(X_1\theta, X_2\theta)$ satisfies Assumption 1.
Finally, we estimate and compare the numerical results of the asymptotic formula (17) and the simulation of $\Psi(x,n)$ . Set $n=5$ ; then the asymptotic formula (17) becomes
Denote by $\Psi_2(x,n)$ the Monte Carlo simulation of $\Psi(x,n)$ . From [Reference Nelsen25, Exercise 3.23] we can generate an FGM random pair $(Z, \theta)$ , and then we get $(X_1, X_2, \theta)=(Z\zeta_1, Z\zeta_2, \theta)$ by generating two independent uniform (0,1) variates, $\zeta_1$ and $\zeta_2$ . We simulate $N=10^7$ samples of $\left((X_{11},X_{21}, \theta_{1}), \ldots,(X_{1n},X_{2n}, \theta_{n})\right)$ , and for each $k=1,\ldots,N$ , we denote by $\left( (X_{11}^{(k)},X_{21}^{(k)}, \theta_{1}^{(k)}), \ldots,(X_{1n}^{(k)},X_{2n}^{(k)}, \theta_{n}^{(k)}) \right)$ the independent copy of $\left((X_{11},X_{21}, \theta_{1}), \ldots,(X_{1n},X_{2n}, \theta_{n})\right)$ . Hence, $\Psi_2(x,n)$ is estimated by
From Table 1, we observe that the ratio $\Psi_1(x,n)/\Psi_2(x,n)$ approaches 1 as x becomes large. Fixing x, we notice that the ruin probabilities decay when $\rho_1$ decreases from $0.5$ to $0.1$ , meaning that a larger value of $(1-\rho_1)$ leads to a smaller ruin probability of the corresponding business line from the insurance company.
4.2. Some lemmas before the proof of Theorem 1
Clearly, one can derive the following relation for all $x>0$ :
Lemma 2. If the conditions of Theorem 1 hold, then for every fixed $\boldsymbol{b}\in[0,\infty]^d\backslash \{\textbf{0}\}$ and n we have
Proof. For any $1\leq i\leq n$ , we have by the conditions of Theorem 1 that $\boldsymbol{X}_{i}^+ \boldsymbol{\theta}_i\in {\textrm{MRV}}(\alpha,F,\nu)$ and
Since $\boldsymbol{X}_{i}^+ \boldsymbol{\theta}_i$ and $ \boldsymbol{Y}_{i-1}=\prod_{j=1}^{i-1}\boldsymbol{\theta}_{j} $ are independent, it follows from Proposition 1 that
Hence, it suffices to show that for every fixed $\boldsymbol{b}\in[0,\infty]^d\backslash \{\textbf{0}\}$ and n, the following asymptotic formula holds:
Now we aim to show the upper bound of $\mathbb{P}\left(\sum_{i=1}^n \boldsymbol{X}_{i} \boldsymbol{Y}_{i} > \boldsymbol{b}x \right)$ . For an arbitrary $\varepsilon>0$ , choose small $v_1>0$ such that $(1-v_1)^{-\alpha}\leq1+\varepsilon$ and $(1+v_1)^{-\alpha}\geq 1-\varepsilon$ hold. Then
For $K_1(x,n)$ , we have by the relation (21) that
For $K_2(x,n)$ we obtain, for $1\leq k\leq d$ and some $b_k>0$ ,
where we applied (20), Lemma 1, and (8) in the last step, using the independence between
when $i>j$ and the independence between $X_{kj}^+\theta_{kj}$ and $(Y_{k,j-1} , X_{li}^+\theta_{li}Y_{l,i-1})$ when $i<j$ . Hence, for large x, we get
Finally, we construct the lower bound of $\mathbb{P}\left(\sum_{i=1}^n \boldsymbol{X}_{i} \boldsymbol{Y}_{i} > \boldsymbol{b}x \right)$ . We have
so we find
For $K'_{\!\!1}(x,n)$ , we have by the relations (19) and (21) that
For $K'_{\!\!2}(x,n)$ , we obtain
where the last equality follows from Lemma 1, (20), and (8) by the independence between $X_{km}^+\theta_{km}$ and $(Y_{k,m-1} , |X_{ki}|Y_{ki})$ . For
we find that, by Chebyshev’s inequality,
where the last step is due to (3), Lemma 1, (20), and (8), using the independence between the product $X_{km}^+\theta_{km}$ and $\left( Y_{k,m-1} , |X_{ki}| \prod_{j=1,j\neq m}^{i}\theta_{kj} \right)$ . For $K'_{\!\!3}(x,n)$ , by Lemma 1 we have
Therefore, for large x, we obtain
This completes the proof.
Lemma 3. If the conditions of Theorem 1 hold, then for every fixed $\boldsymbol{b}\in[0,\infty]^d\backslash \{\textbf{0}\}$ we get
Proof. For any fixed $0<p<\beta$ and $0<\varepsilon<\beta-\alpha$ , we choose some $0<q<1$ satisfying
Then there exists some $n_1>0$ such that the following relation holds:
Since $X_{1i}^+\theta_{1i}$ is of regular variation from (8), applying [Reference Li20, Lemma 1] for large x we obtain
This completes the proof.
4.3. Proof of Theorem 1
Proof. First we show that for any $n_2$ sufficiently large, the relation (17) holds uniformly for $1\leq n\leq n_2$ . By Lemma 2, the following relations hold uniformly for $1\leq n\leq n_2$ :
and
Next, we study the uniformity of the relation (17) when $n>n_2$ . Clearly,
On the one hand, it holds uniformly for $n>n_2$ that, by Lemmas 2 and 3 and the relation (24),
On the other hand, it holds uniformly for $n>n_2$ that, for $v_3$ satisfying $(1-v_3)^{-\alpha}\leq 1+\varepsilon$ ,
For $I_1(x,n_2)$ , we have by Lemma 2 that
For $I_2(x,n)$ , we obtain by Lemmas 2 and 3 and the relation (24) that
This completes the proof.
5. The study of d-dimensional continuous-time risk model under Assumption 1
In this section we consider an insurance company with d lines of business, for $d\geq1$ . Let x be the initial reserve and let $\boldsymbol{\rho}$ be the allocation vector. For each $1\leq k\leq d$ , the price process of the investment portfolio in the kth business is expressed by a geometric Lévy process $\{ e^{L_k(t)} , \, t\geq0 \}$ ; namely, $\{ L_k(t), \, t\geq0 \}$ is a Lévy process that has stationary and independent increments, is stochastically continuous, and starts from 0. Then the insurer’s discounted risk process, $\mathbf{R}(t)=(R_1(t), \, \ldots, \, R_d(t))$ , is given by
with $t\geq0$ , where $A_{ki}$ , $1\leq k \leq d$ and $i\in\mathbb{N}$ , denotes the ith claim amount from the kth business, and $(c_1, \, \ldots, \, c_d)$ denotes the vector of constant premium rates. The successive claim arrival times are denoted by $0<\tau_1<\tau_2<\cdots$ , and the claim arrival process $\{ N(t); \, t\geq0 \}$ is a renewal process with the following finite renewal function:
The vectors $ (A_{1i}, \, \ldots, \, A_{di})$ , $i\in\mathbb{N}$ , form a sequence of i.i.d. nonnegative random vectors with a generic random vector $(A_1,\, \ldots, \, A_d)$ . The inter-arrival times are denoted by $\chi_1=\tau_1$ and $\chi_i=\tau_i-\tau_{i-1}$ for $i=2, \, 3, \, \ldots$ . The sequence $\{\chi_i, \, i\geq1\}$ is i.i.d. with generic random variable $\chi$ . We define the Laplace exponent for the Lévy process { $L_k(t), \, t\geq0$ }, $1\leq k\leq d$ , by the formula
If $\phi_k(s)$ is finite, then $\mathbb{E}e^{-sL_k(t)}=e^{t\phi_k(s)}<\infty$ for $t\geq0$ .
In the multivariate risk model case, ruin may appear in various situations. Thus, there are several versions of the probabilities of ruin based on various ruin sets. For the model (25) we adopt the following three ruin times:
which is the first instant when all line reserves become negative simultaneously,
which is the first instant when at least one of the businesses is below zero, and
which is the first instant at which all lines have at some point run into deficit, though not necessarily simultaneously. Here, $\inf \emptyset$ is understood as $\infty$ by convention. In this section we study the following three types of ruin probabilities:
where ‘ $\#$ ’ denotes either ‘max’, ‘min’, or ‘ult’. The ruin probability represents a significant indicator of the functional quality of the insurance company and has been studied extensively by [Reference Asmussen and Albrecher3], [Reference Ji13], [Reference Konstantinides and Li16], and [Reference Li19], and others.
Theorem 2. Consider the d-dimensional continuous-time risk model (25). Assume that the vector $\mathbf{A}_i e^{ \mathbf{L}(\tau_{i-1})-\mathbf{L}(\tau_i) } =\left( A_{1i}e^{L_1(\tau_{i-1})-L_1(\tau_i)}, \, \ldots, \, A_{di}e^{L_d(\tau_{i-1})-L_d(\tau_i)}\right)$ , $i\in\mathbb{N}$ , satisfies Assumption 1. If $\mathbb{E}L_k(1)<\infty$ for $1\leq k\leq d$ and there is a constant $\beta>\alpha$ such that the Laplace exponent of the Lévy process $\{ L_k(t), \ t\geq0 \}$ satisfies $\phi_k(\beta)<0$ , then the following asymptotic formulas hold:
Remark 4. Suppose that $L_k(t)$ , $1\leq k\leq d$ , follow the jump-diffusion process
where $r_k\in \mathbb{R}$ stands for the log return rate, $\sigma_k>0$ denotes the volatility, and N(t) represents the aforementioned claim arrival process. Here, $\{ W_k(t), \, t\geq0 \}$ is a Wiener process and $B_{ki}$ describes the jump sizes. Assume that for $1\leq k \leq d$ all the random sources $\{(A_{ki}, B_{ki}), i\geq1\}$ , $\{N(t), t\geq0\}$ , and $\{W_k(t), t\geq0\}$ are mutually independent and $(A_{1i}e^{-B_{1i}}, \, \ldots, \, A_{di}e^{-B_{di}})$ satisfies Assumption 1. We observe that $(A_{1i}e^{-B_{1i}}, \, \ldots, \, A_{di}e^{-B_{di}})$ and $(e^{-r_1\chi_i-\sigma_1W_1(\chi_i)},\,\ldots,\,e^{-r_d\chi_i-\sigma_dW_d(\chi_i)})$ are independent and
for $\beta>\alpha$ . Furthermore,
for any $\mathbf{b}\in[0,\infty]^d\setminus \{\mathbf{0}\}$ , which implies that the random vector $\boldsymbol{A}_i e^{\mathbf{L}(\tau_{i-1}) - \mathbf{L}(\tau_{i})}$ is MRV by Proposition 1.
5.1. Some lemmas before the proof of Theorem 2
Lemma 4. Let the conditions of Theorem 2 hold. Then for any $\mathbf{b}\in[0,\infty]^d\setminus \{\mathbf{0}\}$ and n we get
Proof. From the conditions of Theorem 2, we obtain by Hölder’s inequality that for any fixed $i\in\mathbb{N}$ ,
Hence by Proposition 1 we find
Indeed, in Lemma 2, for each $i=1, \ldots, n$ , take $\boldsymbol{X}_{i}=\boldsymbol{A}_{i}$ ,
then, applying Lemma 2 together with (35), we obtain
The proof is complete.
Lemma 5. Let the conditions of Theorem 2 hold. Then for any $\mathbf{b}\in[0,\infty]^d\setminus \{\mathbf{0}\}$ we get
Proof. Take $\boldsymbol{X}_{i}=\boldsymbol{A}_{i}$ and
for each $i=1, \ldots, n$ . Using Theorem 1 with $n=\infty$ and (35), we obtain
This completes the proof.
The following lemma is significant in the literature on the uniform estimates for the ruin probability of insurance risk models with an exponential Lévy process investment return (see for example [Reference Li19, Lemma 4.8] or [Reference Tang, Wang and Yuen28, Lemma 4.6]).
Lemma 6. Let $Z_k$ , $1\leq k\leq d$ , be an exponential functional of the Lévy process $\{ L_k(t), \, t \geq 0 \}$ defined as
If $\mathbb{E}L_k(1)<\infty$ , then for every $\beta>0$ satisfying $\phi_k(\beta)<0$ , we have $\mathbb{E}Z_k^{\beta}<\infty$ .
Proof. From [Reference Tang, Wang and Yuen28], the condition $\phi_k(\beta)<0$ implies that $\mathbb{E}L_k(1)>0$ . Since $0<\mathbb{E}L_k(1)<\infty$ , we can apply [Reference Maulik and Zwart23, Lemma 2.1] to complete the proof.
5.2. Proof of Theorem 2
Proof. From the conditions in Theorem 2 we may conclude that $\phi_k(\alpha)<0$ for $\alpha<\beta$ and $1\leq k\leq d$ . Then
Since
there exists some large $n_3$ such that, for any $\varepsilon>0$ ,
We first focus on the relation (31) for $\Psi_{\textrm{max}}(x, \, \infty)$ . By (25), (27), and (30), we have
Then, by Lemma 5, it holds that
Let
By Lemma 1, we have
For $I_1(x)$ , choose $v_4$ satisfying
By Lemma 4, $F\in\mathcal{R}_{-\alpha}$ , and (37), we have
From Markov’s inequality, Lemma 6, and the relation (3), we obtain
Then we find
Combining the relations (38), (39), and (41) yields the relation (31) for $\Psi_{\textrm{max}}(x, \, \infty)$ .
Now we establish the relation (32) for $\Psi_{\textrm{min}}(x,\infty)$ . We obtain by the inclusion–exclusion principle and Lemma 5 that
By the inclusion–exclusion principle and (40), we get
Combining (42) and (43) gives the relation (32).
To establish the relation (31) for $\Psi_{\textrm{ult}}(x,\infty)$ , we note that by Lemma 5,
Analogously to the reasoning for the lower bound of $\Psi_{\textrm{max}}(x, \, \infty)$ , it follows that
Combining the relations (44) and (45) gives the relation (31) for $\Psi_{\textrm{ult}}(x, \, \infty)$ . This completes the proof.
Acknowledgements
We are sincerely grateful to the anonymous reviewers, whose insightful comments helped us improve the quality of this paper. We also express our great appreciation to the editor for their kind work.
Funding information
This work was supported by the National Natural Science Foundation of China (Project No. 71271042).
Competing interests
There were no competing interests to declare which arose during the preparation or publication process of this article.