1. Introduction
1.1. Model descriptions
This paper studies a bi-dimensional continuous-time compound risk model, where an insurance company concurrently engages in two lines of business. Both lines are exposed to severe accidents such as car accidents, fire disasters, or catastrophes such as earthquakes, floods, and hurricanes. It is reasonable to assume that one of these risks may lead to numerous claims for both two lines of business at the same time or for only one line of business (if some catastrophe does not affect any line of business, it is then not considered as a risk). To characterize this scenario, we assume that the arrival times of successive risks are 
$0 \lt \sigma_1\leq\sigma_2\leq\cdots$, which constitute a counting process
\begin{equation*}
N(t)=\sup\{n\geq 1:\sigma_n\leq t\},\quad t\geq 0
\end{equation*} with a mean function 
$\lambda(t)=EN(t)=\sum_{i=1}^{\infty}P\left(\sigma_i\leq t\right)$ being finite. We adopt non-negative random vectors 
$\left(\tau_i^{(1)}, \tau_i^{(2)}\right)$, 
$i\geq 1$, to represent the numbers of claims resulting from each line of business caused by catastrophic risks. Then, the discounted surplus process of the insurance firm at time t is as follows:
\begin{equation}
\left (
\begin{matrix}
U_1(t) \\
U_2(t)
\end{matrix}
\right )
=\left (
\begin{matrix}
x_1 \\
x_2
\end{matrix}
\right )+\left (
\begin{matrix}
\int_{0}^{t}e^{-rs}dC_1(s) \\
\int_{0}^{t}e^{-rs}dC_2(s)
\end{matrix}
\right )-\sum_{i=1}^{N(t)}\left (
\begin{matrix}
\sum_{j=1}^{\tau_i^{(1)}}X_{ij}^{(1)}e^{-r\sigma_{i}} \\
\sum_{j=1}^{\tau_i^{(2)}}X_{ij}^{(2)}e^{-r\sigma_{i}}
\end{matrix}
\right ),\;t\geq 0,
\end{equation} where, for 
$l=1,2$, we denote by xl the initial reserve, 
$r\geq 0$ the constant force of interest, 
$X_{ij}^{(l)}$ the jth claim triggered by the ith risk from the lth line of business, 
$\{C_l(t),t\geq 0\}$ the premium accumulation process of the lth line of business whose paths are non-negative, non-decreasing, right continuous, and 
$C_l(t) \lt \infty$ almost surely for any t > 0.
In the above bi-dimensional continuous-time compound risk model (1.1), we suppose that
•
$\left( \tau_{i}^{(1)},\tau_i^{(2)}\right),i\geq 1$, are non-negative, independent and identically distributed random vectors with generic random pair $(\tau_1,\tau_2)$ satisfying $P(\tau_1=0) \gt 0$, $P(\tau_2=0) \gt 0$, and $ P(\tau_1=0, \tau_2=0)=0$;• for
$l=1,2$, $\left\{X_{ij}^{(l)},i\geq 1,j\geq 1\right\}$ is a non-negative identically distributed array with mutually independent rows and $\left\{X_{ij}^{(l)},j\geq 1\right\}$, $i\geq 2$ are independent copies of $\left\{X_{1j}^{(l)},j\geq 1\right\}$;•
$\{(C_1(t),$ $C_2(t));t\geq 0\}$ and $\{N(t);t\geq 0\}$ are mutually independent, and they are independent of $\left\{X_{ij}^{(1)},i\geq 1,j\geq 1\right\}$, $\left\{X_{ij}^{(2)},\right.$ $\left.i\geq 1,j\geq 1\right\}$, and $\left\{\left( \tau_{i}^{(1)},\tau_i^{(2)}\right),i\geq 1\right\}$, where the later three sequences are arbitrarily dependent.
1.2. Brief review
The asymptotic analysis for ruin probabilities in compound risk models is an enduring topic that can be observed from numerous references.
Tang et al. [Reference Tang, Su, Jiang and Zhang24] gave the concept of compound risk models first. According to them, the compound risk models are more realistic than the classical risk models and natural extensions of the classical ones. They established the precise large deviations for this kind of model. Since then, many researchers have further explored this topic and advanced the research. Yang and Wang [Reference Yang and Wang27] investigated the precise large deviations for dependent dominatedly varying-tailed random variables, and the results were utilized to derive asymptotic bounds for finite-time ruin probabilities in compound renewal risk models. Yang et al. [Reference Yang, Wang, Leipus and Šiaulys29] and Zong [Reference Zong31] obtained asymptotic formulae for finite-time ruin probabilities of nonstandard compound renewal risk models with negatively dependent claims, respectively. The former discussed two cases, the first one considered the distributions of the numbers of claims to be dominatedly varying-tailed and the second one considered the distributions of the numbers of claims to be in the maximum domain of attraction of the Gumbel distribution, and the claim sizes are light-tailed; the later studied the case where both the distributions of the numbers of claims and the claim sizes are heavy-tailed but without a dominating relationship among their tails. Aleškevičiene et al. [Reference Aleškevičiene, Leipus and Šiaulys1] considered the random sums of independent and identically distributed random variables and applied the main result to study the asymptotic behavior of the finite-time ruin probability. Leipus and Šiaulys [Reference Leipus and Šiaulys19] studied the modified compound discrete-time risk model and obtained the asymptotics of the finite-horizon ruin probability in such a model for a subclass of heavy-tailed claim sizes and numbers of claims. For some recent research for ruin probabilities in compound risk models, the reader is referred to Yang et al. [Reference Yang, Wang and Chen28], who developed a second-order asymptotic formula that provides a more accurate and efficient solution compared to traditional methods for calculating ruin probabilities over an infinite time and to Liu and Gao [Reference Liu and Gao21], who analyzed a nonstandard compound renewal risk model incorporating stochastic return on investments and derived asymptotic formulae for finite-time ruin probabilities.
As the insurance industry continues to develop, insurance operations are becoming increasingly diversified. Traditional unidimensional models are sometimes becoming pale in interpreting realistic actuarial issues. Chan et al. [Reference Chan, Yang and Zhang4] promoted it by introducing bi-dimensional ones. For bi-dimensional risk models, various kinds of ruin probabilities are defined, and some relevant results are derived. In this paper, we mainly explore the sum-ruin probability, which is defined in the following way.
\begin{equation}
\psi_{\rm sum}(x_1,x_2;t)=P\left(T_{\rm sum}\leq t|(U_1(0),U_2(0))=(x_1,x_2)\right),
\end{equation}where T sum is the corresponding time of ruin defined by
\begin{equation*}
T_{\rm sum}=\inf\{s\geq 0|U_1(s)+U_2(s) \lt 0\}.
\end{equation*} T sum, first defined in Chan et al. [Reference Chan, Yang and Zhang4], is a rather appealing choice to measure the time when an insurance firm ruins. This measure is generally considered a vital actuarial quantity in the measurement of the solvency of an insurance company and has natural interpretations. For instance, 
$\{T_{\rm sum}\leq t\}$ indicates that the sum of 
$U_1(s)$ and 
$U_2(s)$ is negative for one or more instances within the time interval 
$[0,t]$. For the papers investigating 
$\psi_{\rm sum}(x_1,x_2;t)$, we refer the reader to Gao and Yang [Reference Gao and Yang16] for pairwise strongly quasi-asymptotically independent (QAI) claims, Cheng and Yu [Reference Cheng and Yu10] for different claim-number processes and strongly subexponential claims, Cheng [Reference Cheng9] for uniform asymptotics for 
$\psi_{\rm sum}(x_1,x_2;t)$ with two arbitrarily dependent claim-number processes, Chen et al. [Reference Chen, Wang, Cheng and Yan8] for the claims that follow the dependence structure appearing in Ko and Tang [Reference Ko and Tang18], for Chen et al. [Reference Chen, Li and Cheng7] which partly extended the results of Chen et al. [Reference Chen, Wang, Cheng and Yan8] by introducing the stochastic returns.
1.3. Motivation
Inspired by the papers mentioned above, this paper aims to study the asymptotic formula for 
$\psi_{\rm sum}(x_1,x_2;t)$ in compound risk model (1.1) mainly based on the following considerations:
• We find that almost all papers concerning compound risk models are unidimensional and little work focuses on bi-dimensional ones. However, it is meaningful to investigate bi-dimensional risk models because of the increasing complexity of insurance products.
•
$\psi_{\rm sum}(x_1,x_2;t)$ is rarely studied in bi-dimensional risk models for mathematical complexity. This pushes us to investigate $\psi_{\rm sum}(x_1,x_2;t)$ in a bi-dimensional risk model for it is an important tool to measure the solvency of an insurance company.• In a major catastrophe, the number of insurance claims may be beyond expectations, and a larger number of claims implies more casualties and also indicates larger claim amounts. This highlights the importance of assuming arbitrary dependency between the claims sequence and the numbers of claims sequence.
Consequently, in this paper, we are going to investigate the asymptotic formula for the sum-ruin probability in a bi-dimensional compound risk model, where the tails of claims are dominated by the tails of the numbers of claims and no extra dependence structures are equipped between the numbers of claims and the claims, which has not been considered as far as we know.
The remainder of this paper is organized as follows. Section 2 shows some preliminaries and presents our main results. Section 3 conducts some numerical studies, and the results are presented by some figures. Section 4 gives some necessary lemmas and the proofs of the main results.
2. Preliminaries and main results
Throughout the paper, we denote by 1A the indicator function of a set A and by 
$x^{+}=x1_{(x\geq 0)}$, 
$x^{-}=-x1_{(x \lt 0)}$ the positive and negative parts of x, respectively. All limit relations hold as 
$x\to\infty$ or 
$(x_1,x_2)\to(\infty,\infty)$ unless otherwise noted. For two univariate or bivariate functions f and g, we write 
$f=o(g)$ if 
$\lim f/g=0$; 
$f=O(g)$ if 
$\limsup f/g \lt \infty$; 
$f\asymp g$ if both 
$f=O(g)$ and 
$g=O(f)$; 
$f\lesssim g$ if 
$\limsup f/g\leq 1$ and 
$f\sim g$ if 
$\lim f/g=1$. For two real numbers a and b, we write 
$a\vee b=\max\{a,b\}$ and 
$a\wedge b=\min\{a,b\}$.
2.1. Heavy-tailed distributions
Heavy-tailed distributions are a superior choice to model the claims or the numbers of claims. Recall a random variable ξ or its distribution function Fξ is said to be heavy-tailed if 
$\int_{-\infty}^{\infty}e^{ux}F_{\xi}(dx)=\infty$ for any u > 0. We are particularly interested in the class of distribution functions with consistent variation, denoted by 
$\mathscr{C}$. Recall that a random variable ξ or its distribution function Fξ supported on 
$(-\infty,\infty)$ is consistently varying-tailed if
\begin{equation*}
\lim_{y\uparrow 1}\limsup\frac{\overline{F_{\xi}}(xy)}{\overline{F_{\xi}}(x)}=1\quad \mbox{or}\quad \lim_{y\downarrow 1}\liminf\frac{\overline{F_{\xi}}(xy)}{\overline{F_{\xi}}(x)}=1.
\end{equation*} We need the following subclasses of the heavy-tailed distribution class to continue our research. A random variable ξ or its distribution function Fξ is called long-tailed, denoted by 
$\xi\in\mathscr{L}$ or 
$F_{\xi}\in\mathscr{L}$, if for any fixed y > 0,
\begin{equation*}
\overline{F_{\xi}}(x+y)\sim\overline{F_{\xi}}(x);
\end{equation*} is called dominatedly varying-tailed, denoted by 
$\xi\in\mathscr{D}$ or 
$F_{\xi}\in\mathscr{D}$ if
\begin{equation*}
\overline{F_{\xi}}(xy)=O\left(\overline{F_{\xi}}(x)\right)
\end{equation*} holds for some 
$0 \lt y \lt 1$; is called regularly varying-tailed, denoted by 
$\xi\in\mathscr{R}_{-\alpha}$ or 
$F_{\xi}\in\mathscr{R}_{-\alpha}$ for some α > 0 if
\begin{equation*}
\overline{F_{\xi}}(xy)\sim y^{-\alpha}\overline{F_{\xi}}(x)
\end{equation*}for any y > 0. In conclusion, we have the following relations,
\begin{equation*}
\mathscr{R}_{-\alpha}\subset\mathscr{C}\subset\mathscr{L}\cap\mathscr{D}\subset\mathscr{L}.
\end{equation*}For a more comprehensive discussion about the heavy-tailed distribution class or its subclasses and their applications to insurance and finance, the reader is referred to Bingham et al. [Reference Bingham, Goldie, Teugels and Teugels3], Cline and Samorodnitsky [Reference Cline and Samorodnitsky13], and Embrechts et al. [Reference Embrechts, Klüppelberg and Mikosch15], among others.
At the end of this subsection, we recall an elementary quantity. We set
\begin{equation*}
\overline{F_{\xi}}_*(y)=\liminf\frac{\overline{F_{\xi}}(xy)}{\overline{F_{\xi}}(x)}
\end{equation*}and define
\begin{equation*}
M^{+}_{F_{\xi}}=-\lim_{y\to\infty}\frac{\log\overline{F_{\xi}}_*(y)}{\log y}
\end{equation*}as the upper Matuszewska index of distribution function Fξ. The reader is referred to Bingham et al. [Reference Bingham, Goldie, Teugels and Teugels3] for a formal discussion.
2.2. Dependence structures and main results
To describe the dependency of random variables, various dependence structures have been introduced over the past few decades. To learn more about popular ones, the reader is referred but not limited to the following sources: Ko and Tang [Reference Ko and Tang18], Liu [Reference Liu20], Geluk and Tang [Reference Geluk and Tang17], Chen and Yuen [Reference Chen and Yuen5], Asimit and Badescu [Reference Asimit and Badescu2], Chen and Yuen [Reference Chen and Yuen6], and Cheng and Cheng [Reference Cheng and Cheng12]. Among these papers, Chen and Yuen [Reference Chen and Yuen5] proposed a dependence structure named quasi-asymptotic independence.
Definition 2.1. For two non-negative random variables ξ 1 and ξ 2 with distributions F 1 and F 2, respectively, they are said to be QAI if the following relation
\begin{equation}
\lim\frac{P\left(\xi_1 \gt x, \xi_2 \gt x\right)}{\overline{F_{1}}(x)+\overline{F_{2}}(x)}=0
\end{equation}holds.
Note that (2.1) can be expressed in terms of copula,
\begin{equation*}
\lim_{(u,v)\to(0^{+},0^{+})}\frac{\hat{C}(u,v)}{u+v}=0,\quad with\quad u=\overline{F_1}(x)\quad and \quad v=\overline{F_2}(x),
\end{equation*} where 
$\hat{C}(\cdot,\cdot)$ is the survival copula, which connects with the corresponding copula 
${C}(\cdot,\cdot)$ through the equality 
$\hat{C}(u,v)=u+v-1+C(1-u,1-v)$. For a formal and comprehensive discussion about copula, the reader is referred to Nelson [Reference Nelsen22]. We point out that Definition 2.1 provides us with a wide range of random variables. Many bivariate copulas satisfy (2.1) such as the Frank copula, the Johnson-Kotz iterated Farlie-Gumbel-Morgenstern copula, the Ali-Mikhail-Haq copula, and so on.
In this paper, we assume that the number of claims τ 1 and τ 2 are QAI. Concurrently, we apply extended negative dependence, which first appeared in Liu [Reference Liu20], to characterize the dependency among the claims.
Definition 2.2. For random variables 
$\{\xi_n,n\geq 1\}$, if there exists a finite constant M > 0 such that for each 
$n\geq 1$ and all 
$x_i\in(-\infty, \infty)$, 
$1\leq i\leq n$,
\begin{equation}
P\left(\bigcap_{i=1}^n\left\{\xi_i \gt x_i\right\}\right)\leq M\prod_{i=1}^nP\left(\xi_i \gt x_i\right),
\end{equation} then the random variables 
$\{\xi_n,n\geq 1\}$ are said to be extended upper negatively dependent (EUND); analogously, if there exists a finite constant M > 0 such that for each 
$n\geq 1$ and all 
$x_i\in(-\infty, \infty)$, 
$1\leq i\leq n$,
\begin{equation}
P\left(\bigcap_{i=1}^n\left\{\xi_i\leq x_i\right\}\right)\leq M\prod_{i=1}^nP\left(\xi_i\leq x_i\right),
\end{equation} then the random variables 
$\{\xi_n,n\geq 1\}$ are said to be extended lower negatively dependent (ELND). Further, the random variables 
$\{\xi_n,n\geq 1\}$ are called extended negatively dependent (END) if they are both EUND and ELND.
The constant M is called the dominating constant. When dealing with END random variables, we always assume that the dominating constants provided in (2.2) and (2.3) are identical. The extended negative dependence structure indeed contains a wide range of random variables. If M = 1 in (2.2) and (2.3), then the random variables 
$\{\xi_n,n\geq 1\}$ are negatively upper dependent and negatively lower dependent, respectively. And they are called positively dependent (PD) if the inequalities (2.2) and (2.3) both hold in the reverse direction when M = 1. An ND sequence obviously must be an END sequence. Nevertheless, it is possible for some PD sequences to find a corresponding positive constant M such that both (2.2) and (2.3) hold, {for example, the Ali-Mikhail-Haq copula with parameter 
$\theta\in(0,1]$, and Example 4.1 in [Reference Liu20]. The END structure is notably more extensive than the ND structure as it can not only represent a negative dependence structure but also, to a certain extent, a positive one.
The END structure and some other dependence structures appearing in the papers mentioned at the beginning of this subsection have been extensively studied. These structures hold significance for research due to their breaking of the classical independent hypothesis and thus bridging the gap between theoretical research and practical applications. Moreover, these concepts have demonstrated their value in various other fields, including precise large deviations, probability inequalities, asymptotic analysis on random sums, and so on.
We would like to share our main results in this part.
Theorem 2.1.
Consider the bi-dimensional compound risk model (1.1). For 
$l=1,2$, let 
$\left\{X_{1j}^{(l)},j\geq 1\right\}$ be an END sequence and Fl, Gl be the distribution functions of 
$X_{11}^{(l)}$ and τl, respectively. Suppose that 
$G_l\in\mathscr{C}$, 
$\mu_l:=EX_{11}^{(l)} \lt \infty$ and
\begin{equation}
x\overline{F_l}(x)=o\left(\overline{G_{l}}(x)\right).
\end{equation} Let τ 1, τ 2 be QAI. If T > 0 satisfies 
$\lambda(T) \gt 0$ and
\begin{equation}
E\left[N^{\beta+1}(T)\right] \lt \infty \quad {\textit{for some}}\quad \beta \gt M^{+}_{H},
\end{equation} where H is any distribution whose tail satisfies 
$\overline{H}(x)\sim\overline{G_1}\left(\frac{x}{\mu_1}\right)+\overline{G_2}\left(\frac{x}{\mu_2}\right)$. Then, we have
\begin{equation}
\psi_{\rm sum}(x_1,x_2;T)\sim\sum_{l=1}^2\int_0^T\overline{G_l}\left(\frac{x_1+x_2}{\mu_l}e^{rs}\right)\lambda(ds).
\end{equation}We will give a corollary in the following to provide convenience to numerical studies, which simplifies (2.6).
Corollary 2.1. Let the conditions of Theorem 2.1 be valid, if G 1, 
$G_2\in\mathscr{R}_{-\alpha}$ for some α > 0, and 
$\{N(t),t\geq 0\}$ is a Poisson process with rate λ > 0, then we have
\begin{equation}
\psi_{\rm sum}(x_1,x_2;T)\sim\frac{\lambda\left(1-e^{-\alpha rT}\right)}{\alpha r}\sum_{l=1}^2\mu_l^{\alpha}\overline{G_l}(x_1+x_2).
\end{equation}3. Numerical examples
This section will give some numerical examples to examine the accuracy of (2.7). Concretely speaking, we will use the crude Monte Carlo method to estimate the ruin probability defined in (1.2) and compare the simulated values with the asymptotic values on the RHS of (2.7). In what follows, for 
$l=1,2$, we present some settings.
• sample size
$m=100,000$;•
$x_1+x_2$ varies from $20,000$ to $40,000$ with step 100;•
$C_l(t)=t$, for any $t\geq 0$;• the rate λ = 1, T = 10, and r = 0.03;
• for each
$n=1,2,\cdots$, $P(\tau_1=n)=P(\tau_2=n)=\frac{6}{\pi^2n^2}$ and τ 1 and τ 2 are connected by Frank copula as follows.
\begin{equation*}
C(u_1,u_2)=-\log\left(1+\frac{\left(e^{-u_1}-1\right)\left(e^{- u_2}-1\right)}{e^{-1}-1}\right),
\end{equation*}then, in this case,
$G_1,G_2\in\mathscr{R}_{-1}$.
We will then give the details about how we estimate the empirical values of (1.2) in several steps. For any 
$1\leq k\leq m$,
• step 1: independently generate a sequence of inter-arrival times
$\left\{\vartheta_{k1},\vartheta_{k2},\cdots, \right\}$ and determine
\begin{equation*}
\tilde{N}_{k}(T)=\sup\left\{n\geq 1:\sum_{i=1}^n\vartheta_{ki}\leq T\right\};
\end{equation*}• step 2: generate
$\tilde{N}_{k}(T)$ pairs corresponding numbers of claims $\left(\tau_{kj}^{(1)},\tau_{kj}^{(2)}\right)$, $1\leq j\leq \tilde{N}_{k}(T)$;• step 3: for every pair of
$\left(\tau_{kj}^{(1)},\tau_{kj}^{(2)}\right)$, $1\leq j\leq \tilde{N}_{k}(T)$, generate two sequences of claims $\left\{X^{(1)}_{kij},1\leq i\leq \tau_{kj}^{(1)}\right\}$ and $\left\{X^{(2)}_{kij},1\leq i\leq \tau_{kj}^{(2)}\right\}$;• step 4: compute the quantity
$U_1\left(\sigma_{kj}\right)+U_2\left(\sigma_{kj}\right)$ in the following way,
(3.1)\begin{equation}
U_1\left(\sigma_{kj}\right)+U_2\left(\sigma_{kj}\right)=x_1+x_2+2\left(\frac1{r}-\frac1{r}e^{-r\left(\sigma_{kj}\right)}\right)-\sum_{l=1}^{j}\left(\sum_{i=1}^{\tau_{kl}^{(1)}}X^{(1)}_{kij}+\sum_{i=1}^{\tau_{kl}^{(2)}}X^{(2)}_{kij}\right),
\end{equation}where
$\sigma_{kj}=\sum_{i=1}^{j}\vartheta_{ki}$ and j = 1, 2, ⋯, $\tilde{N}_{k}(T).$
After m iterations of such loops, we may give empirical values of (1.2) as follows:
\begin{equation*}
\tilde{\psi}_{\rm sum}(x_1,x_2;T)=\frac1{m}\sum_{k=1}^{m}1_{\left(\bigwedge_{j=1}^{\tilde{N}_{k}(T)}\left(U_1\left(\sigma_{kj}\right)+U_2\left(\sigma_{kj}\right) \lt 0 \right)\right)}.
\end{equation*}We will subsequently introduce our numerical examples.
3.1. Independent claims
In the following examples, the claims from the same line of business are independent of each other.
3.1.1. Numerical example 1
Let 
$F_1(x)=F_2(x)=1-e^{-\sqrt{x}}$, x > 0, be Weibull distributions, which means that the insurance firm has to face heavy-tailed claims. When the initial reserve varies, we will repeat steps 1–4 for m times and compute the proportion of ruining samples according to (3.1). Figure 1 presents the result of this example.
Figure 1. Asymptotic and empirical values of 
$\psi_{\rm sum}(x_1,x_2;T)$ with independent Weibull distributed claims.
3.1.2. Numerical example 2
Let 
$F_1(x)=F_2(x)=1-e^{-x}$, x > 0, be exponential distributions, which means that the insurance firm faces light-tailed claims in the sense that the claims have finite exponential moments. In a manner analogous to Numerical example 1, the result is presented in Figure 2.
Figure 2. Asymptotic and empirical values of 
$\psi_{\rm sum}(x_1,x_2;T)$ with independent exponentially distributed claims.
3.2. Dependent claims
In the following examples, the claims from each line of business are jointly distributed by Clayton copula of the form
\begin{equation*}
C(u_1,u_2,\cdots,u_n)=\left(u_1^{-1}+u_2^{-1}+\cdots+u_n^{-1}-n+1\right)^{-1},
\end{equation*} where 
$n=\tau_{kj}^{(1)}$ or 
$\tau_{kj}^{(2)}$ were indicated in the beginning of this section. According to Example 4.2 of [Reference Liu20] and Remark 3.1 of [Reference Ko and Tang18], these claims are EDN random variables.
3.2.1. Numerical example 3
Let 
$F_1(x)=F_2(x)=1-e^{-\sqrt{x}}$, x > 0, be Weibull distributions. Similarly, the result is shown in Figure 3.
Figure 3. Asymptotic and empirical values of 
$\psi_{\rm sum}(x_1,x_2;T)$ with dependent Weibull distributed claims.
3.2.2. Numerical example 4
Let 
$F_1(x)=F_2(x)=1-e^{-x}$, x > 0, be exponential distributions, which means that the insurance firm faces light-tailed claims. Figure 4 shows the result of this example.
Figure 4. Asymptotic and empirical values of 
$\psi_{\rm sum}(x_1,x_2;T)$ with dependent exponential distributed claims.
Through Figures 1–4, we find that the obtained asymptotic results perform well when the initial wealth is large enough. However, we find that the empirical values fluctuate around the asymptotic values, which is not a surprise since the larger the initial wealth is the smaller the ruin probability 
$\psi_{\rm sum}(x_1,x_2; T)$ becomes and the more fluctuation the empirical value exhibits.
4. Proofs of the main results
4.1. Some lemmas
In preparation for proving Theorem 2.1, we need some necessary lemmas that not only aid in establishing the main results but also have their own merits.
The first lemma obtains a result about tail probabilities of random sums, where the distributions of random numbers of summands dominate the distributions of summands. Aleškevičiene textitet al. [Reference Aleškevičiene, Leipus and Šiaulys1], Robert and Segers [Reference Robert and Segers23], Zhang et al. [Reference Zhang, Cheng and Wang30], and Cheng [Reference Cheng11] have explored this topic in different ways. In contrast to them, the following lemma discusses the random sums of two arbitrarily dependent END sequences, which is powerful in deriving (2.6).
Lemma 4.1. Let 
$\{\xi,\xi_i, i\geq 1\}$, 
$\{\eta,\eta_i, i\geq 1\}$ be two sequences of non-negative identically distributed END random variables with positive means, and ω 1, ω 2 be two non-negative integer-valued QAI random variables. If 
$F_{\omega_1}$, 
$F_{\omega_2}\in\mathscr{C}$,
\begin{equation}
x\overline{F_\xi}(x)=o\left(\overline{F_{\omega_1}}(x)\right), \quad x\overline{F_\eta}(x)=o\left(\overline{F_{\omega_2}}(x)\right),
\end{equation}then it holds that
\begin{equation}
P\left(\sum_{i=1}^{\omega_1}\xi_i+\sum_{j=1}^{\omega_2}\eta_j \gt x\right)\sim \overline{F_{\omega_1}}\left(\frac{x}{E\xi}\right)+\overline{F_{\omega_2}}\left(\frac{x}{E\eta}\right).
\end{equation}Proof. We denote by H 1 and H 2 the distributions of 
$\sum_{i=1}^{\omega_1}\xi_i$ and 
$\sum_{j=1}^{\omega_2}\eta_j$, respectively. According to Theorem 3.1 of [Reference Cheng11], we have
\begin{equation*}
\overline{H_1}(x)\sim \overline{F_{\omega_1}}\left(\frac{x}{E\xi}\right)\quad {\rm{and}}\quad \overline{H_2}(x)\sim\overline{F_{\omega_2}}\left(\frac{x}{E\eta}\right).
\end{equation*} This implies that both H 1 and H 2 belong to 
$\mathscr{C}$. Therefore, because of Theorem 3.1 of [Reference Chen and Yuen5], if we want to prove (4.2), it is sufficient to prove that 
$\sum_{i=1}^{\omega_1}\xi_i$ and 
$\sum_{j=1}^{\omega_2}\eta_j$ are QAI. For every 
$0 \lt \delta \lt 1$, we do the following decompositions.
\begin{equation}
\begin{aligned}
P\left(\sum_{i=1}^{\omega_1}\xi_i \gt x,\sum_{j=1}^{\omega_2}\eta_j \gt x\right)&\leq P\left(\sum_{i=1}^{\omega_1}\xi_i \gt x,\sum_{j=1}^{\omega_2}\eta_j \gt x,\left\{\omega_1\leq \frac{(1-\delta)x}{E\xi}\right\}\bigcup \left\{\omega_2\leq \frac{(1-\delta)x}{E\eta}\right\}\right)\\
&\quad +P\left(\omega_1 \gt \frac{(1-\delta)x}{E\xi}, \omega_2 \gt \frac{(1-\delta)x}{E\eta}\right)\\
&:=I_1+I_2.
\end{aligned}
\end{equation} Since 
$F_{\omega_1}$ and 
$F_{\omega_1}$ belong to 
$\mathscr{C}\subset\mathscr{D}$ and ω 1, ω 2 are QAI, it is obvious that
\begin{equation}
\begin{aligned}
I_2&\leq P\left(\omega_1 \gt \frac{(1-\delta)x}{E\xi\vee E\eta}, \omega_2 \gt \frac{(1-\delta)x}{E\xi\vee E\eta}\right)\\
&=o\left(\overline{F_{\omega_1}}\left(\frac{(1-\delta)x}{E\xi\vee E\eta}\right)+\overline{F_{\omega_2}}\left(\frac{(1-\delta)x}{E\xi\vee E\eta}\right)\right)\\
&=o\left(\overline{H_1}(x)+\overline{H_2}(x)\right).
\end{aligned}
\end{equation} Before we estimate I 1, we need to make some preparations. According to (4.1) and Lemma 4.4 of [Reference Embrechts, Klüppelberg and Mikosch14], there exists some non-decreasing slowly varying function L(x) with 
$L(x)\to\infty$ such that 
$\overline{F_{\xi}}(x)=o\left(\frac{\overline{F_{\omega_1}}(x)}{xL(x)}\right)$ and 
$\overline{F_{\eta}}(x)=o\left(\frac{\overline{F_{\omega_2}}(x)}{xL(x)}\right)$ (a function L(x) is said to be slowly varying if 
$L(xy)\sim L(x)$ for any y > 0). We define two distributions A 1 and A 2 by
\begin{equation*}
A_l(x)=\left(1-\frac{\overline{F_{\omega_l}}(x)}{xL(x)}\wedge 1\right)1_{(x \gt 1)},\quad l=1,2.
\end{equation*} It is easy to check that 
$A_l\in\mathscr{D}$ since 
$F_{\omega_l}\in\mathscr{C}$, 
$l=1,2$. Note that 
$\overline{F_{\xi}}(x)=o\left(\overline{A_1}(x)\right)$ and 
$\overline{F_{\eta}}(x)=o\left(\overline{A_2}(x)\right)$, then there exists some positive constant C (irrespective of x) such that for large x,
\begin{equation}
\begin{aligned}
I_1&\leq P\left(\sum_{i=1}^{[(1-\delta) x/E\xi]}\xi_i \gt x\right)+P\left(\sum_{j=1}^{[(1-\delta)x/E\eta]}\eta_j \gt x\right)\\
&\leq P\left(\sum_{i=1}^{[(1-\delta) x/E\xi]}(\xi_i-E\xi) \gt \frac{\delta}{2}x\right)+P\left(\sum_{j=1}^{[(1-\delta)x/E\eta]}(\eta_j-E\eta) \gt \frac{\delta}{2} x\right)\\
&\leq C\frac{(1-\delta) x}{E\xi}\overline{A_{1}}\left(\frac{\delta}{2}x\right)+C\frac{(1-\delta)x}{E\eta}\overline{A_{2}}\left(\frac{\delta}{2}x\right)\\
&=C\frac{(1-\delta)}{E\xi}\frac{\overline{F_{\omega_1}}\left(\frac{\delta}{2} x\right)}{\frac{\delta}{2}L\left(\frac{\delta}{2}x\right)}+C\frac{(1-\delta)}{E\eta}\frac{\overline{F_{\omega_2}}\left(\frac{\delta}{2}x\right)}{\frac{\delta}{2}L\left(\frac{\delta}{2}x\right)}\\
&=o(1)\left(\overline{H_1}(x)+\overline{H_2}(x)\right),
\end{aligned}
\end{equation} where the third step holds because of Lemma 5.1 of [Reference Cheng11] and the last step holds due to 
$L(x)\to\infty$ and 
${F_{\omega_l}}\in\mathscr{C}\subset\mathscr{D}$, 
$l=1,2$. Plugging (4.4) and (4.5) into (4.3) yields the conclusion. Thus, (4.2) is proved.
The following lemma helps prove Theorem 2.1, which can be proved by exploiting Lemma 4.1.
Lemma 4.2. Let the conditions of Theorem 2.1 be valid, then for any T satisfying 
$\lambda(T) \gt 0$ and (2.5), it holds that
\begin{equation}
P\left(\sum_{i=1}^{N(T)}\sum_{j=1}^{\tau_i^{(1)}}X_{ij}^{(1)}e^{-r\sigma_i}+\sum_{i=1}^{N(T)}\sum_{j=1}^{\tau_i^{(2)}}X_{ij}^{(2)}e^{-r\sigma_i} \gt x\right)\sim\sum_{l=1}^2\int_0^T\overline{G_l}\left(\frac{x}{\mu_l}e^{rs}\right)\lambda(ds).
\end{equation}Proof. For simplicity, we write 
$Y_i^{(1)}=\sum_{j=1}^{\tau_i^{(1)}}X_{ij}^{(1)}$, 
$Y_i^{(2)}=\sum_{j=1}^{\tau_i^{(2)}}X_{ij}^{(2)}$, 
$i\geq 1$. It is easy to see that 
$\left\{Y_i^{(1)}, i\geq 1\right\}$ and 
$\left\{Y_i^{(2)}, i\geq 1\right\}$ are two sequences of independent and identically distributed random variables. Because of Lemma 4.1, it is not hard to check that for 
$i\geq 1$, the survival function
\begin{equation}
P\left(Y_i^{(1)}+Y_i^{(2)} \gt x\right)\sim \overline{G_1}\left(\frac{x}{\mu_1}\right)+\overline{G_2}\left(\frac{x}{\mu_2}\right),
\end{equation} which suggests that the distribution of 
$Y_i^{(1)}+Y_i^{(2)}$ belongs to the class 
$\mathscr{C}$. Relation (4.6) amounts to the conjunction of
\begin{equation}
P\left(\sum_{i=1}^{N(T)}Y_i^{(1)}e^{-r\sigma_i}+\sum_{i=1}^{N(T)}Y_i^{(2)}e^{-r\sigma_i} \gt x\right)\lesssim\sum_{l=1}^2\int_0^T\overline{G_l}\left(\frac{x}{\mu_l}e^{rs}\right)\lambda(ds)
\end{equation}and
\begin{equation}
P\left(\sum_{i=1}^{N(T)}Y_i^{(1)}e^{-r\sigma_i}+\sum_{i=1}^{N(T)}Y_i^{(2)}e^{-r\sigma_i} \gt x\right) \gt rsim\sum_{l=1}^2\int_0^T\overline{G_l}\left(\frac{x}{\mu_l}e^{rs}\right)\lambda(ds).
\end{equation} On the one hand, for any fixed integer 
$n_0 \gt 0$, we split the probability into two parts.
\begin{equation}
\begin{aligned}
&\quad P\left(\sum_{i=1}^{N(T)}Y_i^{(1)}e^{-r\sigma_i}+\sum_{i=1}^{N(T)}Y_i^{(2)}e^{-r\sigma_i} \gt x\right)\\
&=\left(\sum_{n=1}^{n_0}+\sum_{n=n_0+1}^{\infty}\right)P\left(\sum_{i=1}^{n}Y_i^{(1)}e^{-r\sigma_i}+\sum_{i=1}^{n}Y_i^{(2)}e^{-r\sigma_i} \gt x,N(T)=n\right)\\
&:=I_1+I_2.
\end{aligned}
\end{equation} We write 
$\Omega_n=\left\{\left(y_1,y_2,\cdots,y_{n+1}\right):0 \lt y_1 \lt y_2 \lt \cdots \lt y_n\leq T \lt y_{n+1}\right\}$, then we use the law of total probability to obtain
\begin{equation}
\begin{aligned}
I_1&=\sum_{n=1}^{n_0}\int_{\Omega_n}P\left(\sum_{i=1}^{n}\left(Y_i^{(1)}+Y_i^{(2)}\right)e^{-ry_i} \gt x\right)P\left(\bigcap_{i=1}^{n+1}\left\{\sigma_i\in dy_i\right\}\right)\\
&\sim\sum_{n=1}^{n_0}\int_{\Omega_n}\sum_{i=1}^n\left(\overline{G_1}\left(\frac{xe^{ry_i}}{\mu_1}\right)+\overline{G_2}\left(\frac{xe^{ry_i}}{\mu_2}\right)\right)P\left(\bigcap_{i=1}^{n+1}\left\{\sigma_i\in dy_i\right\}\right)\\
&\leq \sum_{i=1}^{\infty}\sum_{n=i}^{\infty}\sum_{l=1}^2P\left(\tau_le^{-r\sigma_i} \gt \frac{x}{\mu_l}, N(T)=n\right)\\
&=\sum_{l=1}^2\int_0^T\overline{G_l}\left(\frac{x}{\mu_l}e^{rs}\right)\lambda(ds),
\end{aligned}
\end{equation}where (4.7) and Proposition 5.1 of [Reference Tang and Tsitsiashvili25] are used in the second step. As for the estimate of I 2, by utilizing Lemma 2.4 of [Reference Wang and Chen26] and (4.7), we know for the β presented in (2.5), there exists some constant C > 0 such that
\begin{equation}
\begin{aligned}
I_2&\leq \sum_{n=n_0+1}^{\infty}P\left(\sum_{i=1}^{n}\left(Y_i^{(1)}+Y_i^{(2)}\right) \gt x\right)P(N(T)=n)\\
&\lesssim C\sum_{n=n_0+1}^{\infty}n^{\beta+1}\left(\overline{G_1}\left(\frac{x}{\mu_1}\right)+\overline{G_2}\left(\frac{x}{\mu_2}\right)\right)P(N(T)=n)\\
&\leq C\left(\overline{G_1}\left(\frac{x}{\mu_1}\right)+\overline{G_2}\left(\frac{x}{\mu_2}\right)\right)E\left[N^{\beta+1}(T)1_{(N(T) \gt n_0)}\right]\\
&=o\left(\sum_{l=1}^2\int_0^T\overline{G_l}\left(\frac{x}{\mu_l}e^{rs}\right)\lambda(ds)\right),\quad \mathrm{as}\quad n_0\to\infty,
\end{aligned}
\end{equation} where 
$G_1,G_2\in\mathscr{C}$ and (2.5) are used in the last step. Plugging (4.11) and (4.12) into (4.10) yields (4.8). On the other hand, in a manner analogous to deriving (4.11), it holds that
\begin{equation*}
\begin{aligned}
&\quad P\left(\sum_{i=1}^{N(T)}Y_i^{(1)}e^{-r\sigma_i}+\sum_{i=1}^{N(T)}Y_i^{(2)}e^{-r\sigma_i} \gt x\right)\\
& \gt rsim\left(\sum_{n=1}^{\infty}-\sum_{n=n_0+1}^{\infty}\right)\int_{\Omega_n}\sum_{i=1}^n\left(\overline{G_1}\left(\frac{xe^{ry_i}}{\mu_1}\right)+\overline{G_2}\left(\frac{xe^{ry_i}}{\mu_2}\right)\right)P\left(\bigcap_{i=1}^{n+1}\left\{\sigma_i\in dy_i\right\}\right)\\
&:=J_1-J_2,
\end{aligned}
\end{equation*} where 
$J_1= \sum_{l=1}^2\int_0^T\overline{G_l}\left(\frac{x}{\mu_l}e^{rs}\right)\lambda(ds)$ is direct. While J 2 can be bounded from above as follows:
\begin{equation*}
\begin{aligned}
J_2&\leq \left(\overline{G_1}\left(\frac{x}{\mu_1}\right)+\overline{G_2}\left(\frac{x}{\mu_2}\right)\right)E\left[\left(N(T)\right)1_{(N(T) \gt n_0)}\right]\\
&=o\left(\sum_{l=1}^2\int_0^T\overline{G_l}\left(\frac{x}{\mu_l}e^{rs}\right)\lambda(ds)\right),\quad \mathrm{as}\quad n_0\to\infty,
\end{aligned}
\end{equation*}where the last step is obtained similarly to (4.12). Then, (4.9) is proved and so is (4.6). This completes the proof of Lemma 4.2.
4.2. Proofs of the main results
Proof of Theorem 2.1
The asymptotic upper bound of 
$\psi_{\rm sum}(x_1,x_2; T)$ is immediately obtained by Lemma 4.2, now we aim to prove the asymptotic lower bound of 
$\psi_{\rm sum}(x_1,x_2; T)$, i.e.,
\begin{equation}
\psi_{\rm sum}(x_1,x_2;T) \gt rsim \sum_{l=1}^2\int_0^T\overline{G_l}\left(\frac{x_1+x_2}{\mu_l}e^{rs}\right)\lambda(ds).
\end{equation}Let B be a positive constant. By utilizing Lemma 4.2, it is not hard to verify that
\begin{equation*}
\begin{aligned}
\psi_{\rm sum}(x_1,x_2;T)&\geq P\left(\sum_{l=1}^2\left(\sum_{i=1}^{N(T)}\sum_{j=1}^{\tau_i^{(l)}}X_{ij}^{(l)}e^{-r\sigma_i}-\int_0^Te^{-rs}C_l(ds)\right) \gt x_1+x_2\right)\\
& \gt rsim\sum_{l=1}^2\int_0^T\overline{G_l}\left(\frac{x_1+x_2+2B}{\mu_l}e^{rs}\right)\lambda(ds)P\left(\int_0^Te^{-rs}C_1(ds)\leq B, \int_0^Te^{-rs}C_2(ds)\leq B\right)\\
&\sim \sum_{l=1}^2\int_0^T\overline{G_l}\left(\frac{x_1+x_2}{\mu_l}e^{rs}\right)\lambda(ds)P\left(\int_0^Te^{-rs}C_1(ds)\leq B, \int_0^Te^{-rs}C_2(ds)\leq B\right)\\
&\sim \sum_{l=1}^2\int_0^T\overline{G_l}\left(\frac{x_1+x_2}{\mu_l}e^{rs}\right)\lambda(ds),\quad \mathrm{as} \quad B\to\infty,
\end{aligned}
\end{equation*} where the third step holds due to 
$G_l\in\mathscr{C}\subset\mathscr{L}\cap\mathscr{D}$, 
$l=1,2$. This confirms (4.13) and thus (2.6) is proved.
Proof of Corollary 2.1
It is not hard to verify that (2.7) immediately follows from Theorem 2.1 and Theorem 1.5 of [Reference Bingham, Goldie, Teugels and Teugels3]. To save space, we omit the details.
Acknowledgments
The authors wish to thank the editor and the referees for their careful reading and exceedingly valuable comments on the earlier version of this paper, which have helped us significantly improve the quality of the manuscript.
Open access
