Next Article in Journal
Technical Analysis on the Bitcoin Market: Trading Opportunities or Investors’ Pitfall?
Next Article in Special Issue
Variance and Interest Rate Risk in Unit-Linked Insurance Policies
Previous Article in Journal
Erratum: Hui Ye, Anthony Bellotti. Modelling Recovery Rates for Non-Performing Loans. Risks 7 (2019): 19
Previous Article in Special Issue
Stochastic Mortality Modelling for Dependent Coupled Lives
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Multivariate Collective Risk Model: Dependent Claim Numbers and Panjer’s Recursion

1
msg life Austria Ges.m.b.H., Ausstellungsstraße 50, 1020 Vienna, Austria
2
Department of Financial and Actuarial Mathematics, Vienna University of Technology, Wiedner Hauptstraße 8–10/E105-1, 1040 Vienna, Austria
*
Author to whom correspondence should be addressed.
Risks 2020, 8(2), 43; https://doi.org/10.3390/risks8020043
Submission received: 21 March 2020 / Revised: 25 April 2020 / Accepted: 27 April 2020 / Published: 2 May 2020
(This article belongs to the Special Issue Interplay between Financial and Actuarial Mathematics)

Abstract

:
In this paper, we discuss a generalization of the collective risk model and of Panjer’s recursion. The model we consider consists of several business lines with dependent claim numbers. The distributions of the claim numbers are assumed to be Poisson mixture distributions. We let the claim causes have certain dependence structures and prove that Panjer’s recursion is also applicable by finding an appropriate equivalent representation of the claim numbers. These dependence structures are of a stochastic non-negative linear nature and may also produce negative correlations between the claim causes. The consideration of risk groups also includes dependence between claim sizes. Compounding the claim causes by common distributions also keeps Panjer’s recursion applicable.

1. Introduction

The aggregation of risks in a portfolio is an important component in both insurance mathematics and financial risk management such as credit risk, operational risk, and insurance risk. In the banking industry, risk management becomes more and more important with the introduction of the Basel II and Basel III accords. For this purpose, there are several portfolio risk models, the CreditRisk + model (introduced by Credit Suisse First Boston (1997)) being a particular example with a lot of practical features. It does not require many assumptions and enables the recursive and exact computation of the loss distribution, i.e., it is not necessary to apply Monte Carlo methods that introduce a stochastic error. The collected contributions in Gundlach and Lehrbass (2004) offer a wide overview of CreditRisk + and related results. However, a shortcoming of CreditRisk + is the assumed independence of the default causes. Deshpande and Iyer (2009) also consider a dependence model for default causes in the usual CreditRisk + model that comprises linear combinations of risk factors and is a generalization of Giese (2004). However, both models do not introduce negative correlation. Instead, we consider a more general linear dependence structure in an extension of the CreditRisk + model that allows us to model negative correlation between default numbers and we can also comprise a stochastic component in the linear dependence. Moreover, we consider risk groups that model joint defaults (for both claim numbers and claim sizes).
Insurance companies have to comply with Solvency II requirements, and the model we consider is capable of aggregating and quantifying risks and can thus be used to determine the minimal capital requirements, which is the first pillar of Solvency II. There is also an increasing demand to reflect dependencies between the risks. Our model can be also used to aggregate operational risk.
In the (extended) CreditRisk + model as in the collective risk model, it is necessary to calculate the distribution of a random sum. In risk management, the collective risk model is diversified into several risk clusters or lines of business and thus random sums where each claim number is driven by a claim cause. Thus, we consider a sum of several collective risk models. We should note that the term claim is used in an insurance context, whereas the term default is used in the credit risk context. We will use both notions interchangeably. Since we model claim numbers by Poisson mixture distributions, it is meaningful to speak of claim cause intensities as the claim causes are modelled by the mixing distributions (and accordingly default causes). Hence, we consider the random sum
g G h = 1 N g X g , h ,
where { X g , h } h N are independent sequences of independent and identically distributed random variables for each g G and ( N g ) g G are N 0 -valued random variables independent of { ( X g , h ) g G } h N .
For independent ( N g ) g G and independent { ( X g , h ) g G } h N that are independent of ( N g ) g G , Panjer’s recursion is the main mathematical tool to evaluate the distribution of the random sum in Equation (1). For the rest of this paper, we assume that N 0 are the natural numbers including zero, i.e., N 0 = { 0 , 1 , 2 , 3 , } . The distribution of an N 0 -valued random variable N, denoted by { q n } n N 0 , belongs to a Panjer ( a , b , k ) class with a , b R and k N 0 if q 0 = q 1 = = q k 1 = 0 and
q n = a + b n q n 1 for all n N with n k + 1 .
Let { X h } h N be a sequence of N 0 -valued random variables. If f n = P [ X 1 = n ] and a f 0 1 , then, according to (Hess et al. 2002, Corollary 4.3), the distribution p n = P [ X 1 + + X N = n ] is given by the recursion
p n = 1 1 a f 0 P [ X 1 + + X k = n ] q k + j = 1 n a + b j n f j p n j
for all n N with initial condition
p 0 = q 0 if f 0 = 0 , k 0 q k f 0 k otherwise ,
which is the probability-generating function of the distribution of N at f 0 . The recursive evaluation of such compound distributions was introduced to actuarial science by Panjer (1981) and extended by Willmot and Panjer (1987). The distributions belonging to a Panjer ( a , b , k ) class were identified by Sundt and Jewell (1981); Willmot (1988), and Hess et al. (2002). Gerhold et al. (2010) treat questions of numerical stability of Panjer’s recursion. Sundt and Vernic (2009) contributed extensively to this topic.
Note also that Panjer’s recursion can be generalized to multivariate claim sizes, cf. Sundt (1999), and the assumptions on the claim sizes can be weakened to exchangeability, cf. Hashorva (2011). Further details can be found in Rudolph (2014). For simplicity of the presentation, we concentrate on the dependence between the claim numbers and allow dependence between the individual losses only within risk groups, cf. Definition 2 below.
There are further approaches analyzing portfolio risk. The paper of Soltes and Danko (2017) proposes to create a portfolio with minimal risk. In Zoričák et al. (2019), we find an interesting approach to predict bankruptcy for imbalanced datasets by one-class classification methods. A general paper is proposed by Jajuga (2016) to introduce model risk in finance from an investor’s perspective, this is a complement to the modeling of risk in this paper. An analysis of real data with respect to dependencies and predictability of stock market returns is given in Fabiánová and Glova (2015).
The remainder of this paper is organized as follows: in Section 2, we introduce a non-negative linear dependence structure with a stochastic component of the claim cause intensities between risk factors and prove alternative representations of the claim numbers with the same joint distribution allowing for the application of Panjer’s recursion. It should be pointed out that this structure also allows for negative dependence, see Example 1. In addition, many models comprise dependence structures, but then the question arises how to calculate the distribution.
In Section 3, we consider a different dependence structure: we mix the default cause intensities with common mixture distributions by letting a parameter of the distribution be random. Then, we prove that an alternative representation with the same joint distribution exists that also allows for an application of Panjer’s recursion.
Finally, in Section 4, we conclude with some interesting examples that show that the choice of correlation and of the mixing distribution has an impact on the distribution of the total portfolio loss.
In the Appendix A, we recall some definitions and present a few useful results on compound and mixture distributions and the corresponding probability-generating functions. These results are important for the presentation of our results.

2. Construction of Dependent Claim Numbers by Linear Combinations

In this section, we present one of our major results. The starting point is the following. In an extension of the CreditRisk + model to be found in Gerhold et al. (2010), the authors consider a random claim number N = N 1 + + N m , where, for each i { 1 , , m } , the claim number N i has a Poisson mixture distribution and the mixing random variable is a gamma-distributed claim cause intensity Λ i . They only stipulate that the claim numbers are conditionally independent given Λ 1 , , Λ m . The conditional distribution of the random claim numbers given Λ 1 , , Λ m is assumed to satisfy
L N i Λ 1 , , Λ m = a . s . L N i Λ i = a . s . Poisson ( λ i Λ i ) , i { 1 , , m } ,
where λ i 0 . The claim cause intensities are assumed to be independent.
In contrast, we develop a generalization of this model by admitting dependence between risk group specific claim cause intensities. Dependent claim cause intensities provide the dependence structure. We construct dependence structures by linear combinations of several non-negative risk factors. These dependence structures are chosen stochastically allowing for risk groups. This means that risk clusters (or lines of business) within one risk group default together, i.e., the default of one risk cluster will immediately imply the default of all the other risks in this risk group. In this setting, we can also relax the assumption of gamma-distributed risk factors and consider τ -tempered α -stable distributions too.
In the course of our considerations, we need the following definition of a multivariate Poisson distribution which is motivated in Lindskog and McNeil (2003) and can be also found in (Sundt and Vernic 2009, Chp. 20.1) with a different notation.
Definition 1.
Let m N , G P ( { 1 , , m } ) a collection of subsets of { 1 , , m } with Ø G , and λ = ( λ g ) g G [ 0 , ) G Poisson parameters. For each g G , define the vector c g = ( c 1 , g , , c m , g ) { 0 , 1 } m by
c i , g = 1 g ( i ) .
Let ( N g ) g G be independent random variables with N g Poisson ( λ g ) for every g G . Then, the distribution of the N 0 m -valued random vector
N = ( N 1 , , N m ) = g G c g N g
is called the m-variate Poisson distribution, we use the notation MPoisson( G , λ , m ).
This notation allows us to consider risk groups. In the credit risk interpretation, the obligors in a non-empty group g { 1 , , m } default together with intensity λ g , but independently of the other groups. An additional form of dependence of the claim numbers can come from the linear and stochastic dependence of the intensities, see Assumption 1 below.
Remark 1.
If G = { { 1 } , , { m } } or, more generally, if λ g = 0 for all g G with | G | 2 , then N 1 , , N m in Equation (6) are independent.
We consider individual claim sizes to take values in a vector space V. Throughout the paper, the reader may assume that V equals R or a finite-dimensional Euclidean space and that V = V . Figure 1 below is a numerical illustration with V = R 2 . However, unless specified otherwise, all theoretical results hold for a separable real Banach space V with topological dual V .
Remark 2.
Multi-variate claim sizes have numerous applications, let us mention some examples:
(a) 
A two-component model is of interest for budgeting purposes: The first component represents the loss payments in the current budget period, the second one the actuarial reserves for later claim payments.
(b) 
In an insurance context, the components can represent different types of claim payments. For a portfolio of health insurance contracts, this can be costs of medical treatments and allowances for missing income of the insured. For a portfolio of personal liability or automobile collision insurances, these can be claims for bodily injuries and property damages.
(c) 
For liquidity considerations, the components can refer to the time period when a claim occurs or a counterparty defaults. Furthermore, the size of the loss given default can depend on the time of the default, in particular when a loan or a mortgage is amortized during its life span and not at maturity.
(d) 
In the context of stochastic claims reserving (see Wüthrich and Merz (2008) for a textbook presentation), the components can represent the development years. Here, the default probability (or intensity) refers to the claims originating from the initial insured period; the claims may be reported at a later year and payments may be spread out during the remaining years of the model.
Corresponding to Definition 1, we generalize the compound Poisson distribution (see Definition A1(b)) to the m-variate case.
Definition 2.
Suppose the random variable N has an m-variate Poisson distribution MPoisson G , λ , m with representation (6). Let V be a vector space, and X = ( X g ) g G , where each X g is a probability distribution on V m such that X g ( U 1 , g × × U m , g ) = 1 for all g G where
U i , g = { 0 } i f i g , V i f i g .
Let ( X g , h ) h N be independent i.i.d. sequences of V m -valued random variables with X g , 1 X g for all g G independent of ( N g ) g G . Then, the distribution of
S = ( S 1 , , S m ) = g G h = 1 N g X g , h
is called compound m-variate Poisson distribution CMPoi ( G , λ , m , X ) .
Remark 3.
If V = R and X = δ c g , i.e., the distribution concentrated in c g given by Equation (5), then CMPoi ( G , λ , m , X ) = MPoisson G , λ , m . The distribution X g describes the random losses vector of the obligors in g defaulting together.
As a direct consequence of the stochastic representation in Equation (7), the independence of ( h = 1 N g X g , h ) g G and the characteristic function for compound Poisson distributions, we obtain the following lemma.
Lemma 1.
The characteristic function of a random vector S with the compound m-variate Poisson distribution CMPoi ( G , λ , m , X ) is given by
φ S ( y ) = E e i y , S = g G exp λ g ( 1 φ X g ( y ) ) , y ( V ) m .
Remark 4.
For notational convenience, we write φ S for φ L S and analogously for other distributions in the following.
To state the additivity of compound m-variate Poisson distributions in Lemma 2 below, the following notation is useful:
Definition 3.
Let X 1 , , X n be probability distributions on V m . Then, their convex combination with relative weights α 1 , , α n 0 is a probability distribution on V m defined in terms of α : = α 1 + + α n by
Convex ( α i , X i ) i { 1 , , n } = 1 α i = 1 n α i X i i f α > 0 , δ 0 i f α = 0 .
Remark 5.
If α > 0 and if the joined support of the distributions X 1 , , X n is a finite or countably infinite set like N 0 d , then Convex ( α i , X i ) i { 1 , , n } can be computed in a numerically stable way by calculating, for each atom in the support, the convex combination of the corresponding probabilities of the atom.
Lemma 2.
Let m , n N and let S i CMPoi ( G i , λ i , m , X i ) for i { 1 , , n } be independent random vectors with compound m-variate Poisson distributions according to Definition 2, where λ i = ( λ g , i ) g G i and X i = ( X g , i ) g G i . Then, S 1 + + S n CMPoi ( G , λ , m , X ) with set of groups G = G 1 G n , intensities λ = ( λ g ) g G and group claim size distributions X = ( X g ) g G given by
λ g = i I g λ g , i a n d X g = Convex ( λ g , i , X g , i ) i I g ,
respectively, where I g : = { i { 1 , , n } g G i } for g G .
Proof. 
Using the independence of S 1 , , S n and Lemma 1, it follows with some rearrangement that
φ S ( y ) = i = 1 n φ S i ( y ) = g G exp i I g λ g , i + i I g λ g , i φ X g , i ( y ) , y ( V ) m .
Using Definition 3, Equations (9) and Remark A1, the result follows. □
Remark 6.
Note that Lemma 2 implies the infinite divisibility of the compound m-variate Poisson distribution.
We now provide the conditional characteristic function of the distribution of a compound m-variate Poisson random vector that is mixed with a random vector of dependent intensities, whose components naturally have to be non-negative. For simplicity, we first introduce an assumption on the dependence scenario.
Assumption 1.
Let J Ø be an arbitrary finite set, m N , and assume G P ( { 1 , , m } ) with G Ø and Ø G . Let A j = ( a g , l j ) g G , l { 0 , , n } [ 0 , ) G × { 0 , , n } for j J and let J be a random variable with values in J . Define the random matrix A J = j J 1 { J = j } A j . Let ( R 1 , , R n ) be non-negative random variables and let R 0 be a non-negative constant. Define Λ g = l = 0 n a g , l J R l for g G .
Remark 7.
The random variable J can be interpreted as selecting a dependence scenario for claim cause intensities ( Λ g ) g G . Accordingly, the random variables R 1 , , R n can be interpreted as non-idiosyncratic risk factors, the constant R 0 as idiosyncratic “risk factor” or idiosyncratic component.
Lemma 3.
Let Assumption 1 be satisfied, let λ g , j 0 for g G and j J , and
L S J , R 1 , , R n = a . s . CMPoi G , ( λ g , J Λ g ) g G , m , X .
Then, the characteristic function of the conditional distribution of the random vector S given ( J , R ) = ( J , R 1 , , R n ) satisfies
φ S | ( J , R ) ( y ) = a . s . g G exp λ g , J ( 1 φ X g ( y ) ) l = 0 n a g , l J R l , y ( V ) m .
Proof. 
As a consequence of Equation (10), the representation Λ g = l = 0 n a g , l J R l for g G in Assumption 1 and the characteristic function given in Lemma 1 the result follows. □
Remark 8.
To give a stochastic representation of the conditional distribution of S in Equation (10), consider random variables ( N g ) g G conditionally independent given J , R 1 , , R n such that
L ( N g | J , R 1 , , R n ) = a . s . L ( N g | J , Λ g ) = a . s . Poisson λ g , J Λ g ,
where λ g , j 0 for g G and j J . With sequences { X g , h } h N for g G as in Definition 2, independent of J , R 1 , , R n , define S by Equation (7).
Remark 9.
In Lemma 3, the non-negative parameters λ g , j with g G and j J might look like redundant notation. However, when taking this framework as an extension of the CreditRisk + model, these parameters are useful to control the expected value E N g for g G since in the CreditRisk + we have E Λ g = 1 that only models the structure of dependence. The dependence of λ g , j of j J is sustained to also have the possibility to model characteristics of scenarios.
Remark 10.
This model can be also considered over a time-horizon. Thus, it can be also used to model the impacts of one risk gradually spreading onto other risks. Therefore, we need to introduce sequences of ( X g , h , 1 , , X g , h , d ) h N independent i.i.d. sequences of V m -valued random variables with X g , 1 , r X g , r for all g G r with r { 1 , , d } and G r a risk group with respect to time. For each time step, we need to define corresponding risk groups to model a gradual spread. These risk groups need to be constructed according to potential distribution. A further result on such a multi-period model can be found in (Rudolph 2014, Lemma 6.21).

2.1. Alternative Representation

If we allow for dependent claim cause intensities and risk groups, the independence between the random sums S i = h = 1 N i X i , h for i { 1 , , m } is lost, hence a convolution is not applicable. While the independence is also lost in our setting, it is nevertheless possible to find an adequate representation of the claim numbers with the same distribution such that a variant of Panjer’s recursion can be applied, as will be seen later. These claim cause intensities may be also considered as a multivariate gamma distribution. There are several approaches to construct multivariate gamma distributions, and one of them includes a linear dependence between the marginal distributions, cf. (Johnson et al. 2000, Chp. 48.3.4).
Based on Lemma 3, we derive a theorem that provides us with a general structure of a stochastically chosen non-negative linear dependence scenario:
Theorem 1.
Let the assumptions of Lemma 3 be satisfied. On the other hand, for every j J , consider n + 1 independent sequences of i.i.d. random variables { B j , l , h } h N , independent of { J = j } , R 1 , , R n , such that, for all g G ,
P [ B j , l , 1 = g ] = p g , j , l
where p g , j , l [ 0 , 1 ] with g G p g , j , l = 1 satisfies the condition
p g , j , l k G λ k , j a k , l j = λ g , j a g , l j
for each j J , g G , and l { 0 , , n } . Let { Q j , l } l { 0 , , n } be random variables conditionally independent given { J = j } , R 1 , , R n which satisfy
L ( Q j , l | 1 { J = j } , R 1 , , R n ) = a . s . L ( Q j , l | R l ) = a . s . Poisson R l g G λ g , j a g , l j .
Assume further that ( Q j , 0 , , Q j , n ) and { B j , l , h } h N with l { 0 , , n } are independent. Let { X l , g , h } h N be independent sequences for l { 0 , , n } with identical distribution as X g for g G independent of all previous random variables. Finally, define the V m -valued random vector M by
M = j J 1 { J = j } l = 0 n h = 1 Q j , l X l , B j , l , h , h .
Then, ( M , R 1 , , R n ) and ( S , R 1 , , R n ) have the same distribution.
Proof. 
For the proof, we apply Remark A1. By Lemma 3, the characteristic function of the distribution of ( S , R ) = ( S , R 1 , , R n ) is given by Equation (11).
Now consider the characteristic function of the distribution of the random vector ( M , R ) = ( M 1 , , M m , R 1 , , R n ) . By partitioning J, we get for all y ( V ) m , z R n
φ ( M , R ) ( y , z ) = j J E 1 { J = j } e i z , R exp i l = 0 n h = 1 Q j , l y , X l , B j , l , h , h .
Using the fact that the sequences { X l , g , h } h N and { B j , l , h } h N are i.i.d. and for l { 0 , , n } independent of { J = j } , R 1 , , R n and ( Q j , 0 , , Q j , n ) for each j J and that { B j , l , h } h N has the distribution given in Equation (13) and thus applying the law of total probability, we observe that
φ ( M , R ) ( y , z ) = j J E 1 { J = j } e i z , R l = 0 n g G p g , j , l φ X g ( y ) Q j , l .
Using the conditional independence of the random variables Q j , 0 , , Q j , n given { J = j } , R 1 , , R n for j J and Equation (15),
φ ( M , R ) ( y , z ) = j J E 1 { J = j } e i z , R l = 0 n E g G p g , j , l φ X g ( y ) Q j , l | R l = j J E 1 { J = j } e i z , R l = 0 n exp R l g G λ g , j a g , l j 1 g G p g , j , l φ X g ( y ) .
Noting Equation (14) and the distributive law
R l g G λ g , j a g , l j k G p k , j , l φ X g ( y ) = k G R l g G λ g , j a g , l j p k , j , l φ X g ( y )
yields
φ ( M , R ) ( y , z ) = j J E 1 { J = j } e i z , R l = 0 n exp R l g G λ g , j a g , l j 1 φ X g ( y ) ,
and rearranging yields Equation (11), which completes the proof. □
Remark 11.
The random sums M i for i { 1 , , m } do not in general have a counting distribution in a Panjer ( a , b , k ) class because the distributions of the random variables Q j , l for j J and l { 0 , , n } are not generally in a Panjer ( a , b , k ) class.
Remark 12.
In the above theorem, the case g G λ g , j a g , l j = 0 for a j J and an l { 0 , , n } has not been considered separately. If such a sum is zero, then it follows that L ( Q j , l | R l ) = a . s . Poisson ( 0 ) = δ 0 . Hence, the corresponding term in the sum (16) can be omitted. Thus, the distribution of B j , l , 1 may be chosen arbitrarily.
Remark 13.
Of course, there are also other methods to construct dependence between claim cause intensities: Luo and Shevchenko (2010) propose in a bivariate set-up to model the dependence by a t-copula. They further use simulation and calibration methods. However, in such a case, recursion methods do not seem to be possible. In contrast, this is possible with our approach.
Remark 14.
Using this random and non-negative linear dependence and the concept of risk groups, we are able to compute the expectations and covariances for the claim cause intensities and hence claim numbers. Recall that Λ g = l = 0 n a g , l J R l for g G . Let J and ( R 1 , , R n ) be independent. We observe for g , k G
E Λ g J = l = 0 n a g , l J E R l
and, if the risk factors R 1 , , R n are in L 2 and uncorrelated,
cov Λ g , Λ k | J = l , p = 1 n a g , l J a k , p J cov R l , R p = l = 1 n a g , l J a k , l J Var R l .
Thus, we compute the covariance between Λ g and Λ k for g , k G
cov Λ g , Λ k = E cov Λ g , Λ k | J + cov E Λ g J , E Λ k J = l = 1 n E a g , l J a k , l J Var R l + l , p = 0 n cov a g , l J , a k , p J E R l E R p .
By the Poisson mixture distribution of the claim numbers ( N g ) g G given in Equation (12), the expectation of N g for g G is
E N g = E E [ N g | J , Λ g ] = E λ g , J Λ g = l = 0 n E λ g , J a g , l J E R l .
For Λ g in L 2 , the variance of N g for g G is given by
Var N g = E Var N g | J , Λ g + Var E [ N g | J , Λ g ] = E λ g , J Λ g + Var λ g , J Λ g ,
and the covariance between N g and N k for g , k G with g k is
cov N g , N k = E cov N g , N k | J , R 1 , , R n + cov E [ N g | J , Λ g ] , E [ N k | J , Λ g ] = cov λ g , J Λ g , λ k , J Λ k
due to the conditional independence of the ( N g ) g G given J , R 1 , , R n and the conditional distribution of ( N g ) g G .
It might not be obvious from this remark how this construction may also provide negative dependence. This is different from other dependence modelling, since it usually only comprises positive dependence. We give the following example (see also Section 4.1 below).
Example 1.
Let Assumption 1 be satisfied and assume G = { { 1 } , { 2 } } . Let R 1 and R 2 be two non-negative and non-degenerate risk factors, possibly dependent with unbounded support. Then, let N = 2 and J = { 0 , 1 } and λ g , j 0 for g = 1 , 2 and j = 0 , 1 . Consider a J -valued random variable J independent of ( R 1 , R 2 ) and
A j = 0 j 0 0 0 1 j f o r j = 0 , 1 .
Then, the claim cause intensities are given by Λ 1 = J R 1 and Λ 2 = ( 1 J ) R 2 . Thus,
cov Λ 1 , Λ 2 = E J R 1 ( 1 J ) R 2 E J R 1 E ( 1 J ) R 2 = E R 1 P [ J = 1 ] E R 2 P [ J = 0 ]
because J ( 1 J ) = 0 and J and ( R 1 , R 2 ) are independent. Therefore, with this antithetic choice, we obtain negative correlation between the claim cause intensities since the expected values are positive. Let { X l , g , h } h N with g , l = 1 , 2 be independent of Q j , l and { B j , l , h } h N with X l , g , h = ( X 1 , l , g , h , X 2 , l , g , h ) . Consequently, the equivalent representations are given by Equation (16) with B j , l , h = ( B 1 , j , l , h , B 2 , j , l , h )
M 1 = 1 { J = 1 } h = 1 Q 1 , 1 X 1 , 1 , B 1 , 1 , 1 , h , h a n d M 2 = 1 { J = 0 } h = 1 Q 0 , 2 X 2 , 2 , B 2 , 0 , 2 , h , h
because Poisson ( 0 ) and Ber 0 are δ 0 -degenerate distributions. Hence, the total loss is S ^ = M 1 + M 2 . If we assume R l to be gamma-distributed for l = 1 , 2 , then we can apply a simplified version of Calculation Method 1 given below.
If we choose Λ 1 = J R 1 / E J and Λ 2 = ( 1 J ) R 2 / E 1 J for gamma-distributed R 1 and R 2 , it is possible to obtain every value in [ 1 , 0 ) as correlation by an appropriate choice of the corresponding variances, see also (Schmock 2020, Example 6.38, Version March 25th 2020). In order to obtain perfect negative correlation, the risk factors need to have variance equal to zero, i.e., they have to be degenerate.

2.2. Evaluation of the Loss Distribution

In this section, we give an algorithm that shows how Panjer’s recursion is applicable in our model with the dependence structure given in this section. It might be suggested that the calculation of the distribution of the total loss S ^ = S 1 + + S m arising from S = ( S 1 , , S m ) given in Equation (7) with conditional representation specified in Lemma 3 requires at least n convolutions (one for each risk factor). However, our algorithm, based on the ideas of (Gerhold et al. 2010, sct. 5.5), circumvents these convolutions by an iterated application of Panjer’s recursion (cf. Equations (23) and (26) below) and a convex combination (cf. Equation (24) below). In applications of the algorithm to gamma-distributed risk factors, we use the well-known fact that a negative binomial distribution is a compound Poisson distribution where the severity distribution is a logarithmic distribution, cf. also Ammeter (1948). The convex combination uses the fact that the convolution of compound Poisson distributions is again a compound Poisson distribution.
We now adapt the algorithm in (Gerhold et al. 2010, sct. 5.5) to the evaluation of the total loss S ^ = S 1 + + S m in Equation (7).
Calculation Method 1.Consider the setting of Theorem 1 and assume in addition that R 1 , , R n are independent and have infinitely divisible distributions. Then, Q j , 0 , , Q j , n are independent for each j J . An application of Theorem 1 yields for the total loss S ^ given as the sum of the m components in Equation (7):
S ^ = S 1 + + S m = i = 1 m g G h = 1 N g X i , g , h = d j J 1 { J = j } l = 0 n h = 1 Q j , l i = 1 m X i , l , B j , l , h , h ,
where { ( X 1 , g , h , , X m , g , h ) } h N with g G as well as { ( X 1 , l , g , h , , X m , l , g , h ) } h N with g G and l { 0 , , n } are independent sequences of i.i.d. V m -valued random vectors with corresponding distribution X g . To calculate the distribution of the right-hand side of Equation (20) , we start with the inner sum and define the distributions of the group losses by
X g s = L i = 1 m X i , g , g G ,
where ( X 1 , g , , X m , g ) X g . If the components of ( X 1 , g , , X m , g ) are independent, then X g s can be computed by convolution. For certain classes of multivariate distributions (like the multivariate logarithmic or negative multinomial distribution), the distribution of the sum of its components is known in closed form. In any case, we assume that X g s is known or can be computed for every g G . Since the risk group g in the inner sum of the right-hand side of Equation (20) is selected with probability p g , j , l determined by Equation (14), we obtain the discrete mixture distribution
X j , l : = L i = 1 m X i , l , B j , l , 1 , 1 = Convex ( λ g , j a g , l j , X g s ) g G
for each j J and l { 0 , , n } , see Definition 3.
By (Sundt and Vernic 2009, Corollary 4.1) and the infinite divisibility of the distributions of R 1 , , R n , the random variable Q j , l with conditional distribution specified in Equation (15) has a compound Poisson distribution, cf. Definition A1(b), hence there exist a Poisson intensity ν j , l 0 and a distribution K j , l on N such that Q j , l CPoi ( ν j , l , K j , l ) for every j J and l { 1 , , n } ; see Section 2.3 below for more details and examples. Define the compound distribution
X j , l c = Compound ( K j , l , X j , l ) , j J , l { 1 , , n } ,
see Definition A1(a). If K j , l is a Panjer ( a , b , k ) distribution and X j , l is concentrated on N 0 d , then X j , k c can be evaluated with Panjer’s recursion, which is numerically stable in many cases, cf. (Gerhold et al. 2010, Theorem 4.5). For each j J , define the discrete mixture distribution
X j = Convex ( R 0 λ g , j a g , 0 j , X g s ) g G , ( ν j , l , X j , l c ) l { 1 , , n } ,
and the Poisson intensity
μ j = R 0 g G λ g , j a g , 0 j + l = 1 n ν j , l .
Then,
X j c : = L l = 0 n h = 1 Q j , l i = 1 m X i , l , B j , l , h , h = CPoi μ j , X j , j J .
The final step is to calculate the discrete mixture distribution
L ( S ^ ) = Convex ( P [ J = j ] , X j c ) j J .

2.3. Examples for Distributions of the Risk Factors

2.3.1. Gamma Distribution

Remark 15.
As a specific case in Calculation Method 1, fix a risk factor l { 1 , , n } and assume that R l has a gamma distribution, so R l Gamma ( α l , β l ) with shape parameter α l > 0 and rate parameter β l > 0 . Define, for j J ,
q j , l = g G q g , j , l w i t h q g , j , l = λ g , j a g , l j β l + k G λ k , j a k , l j .
By Lemma A1 (with m = 1 and T 1 in that framework) and comparison of Equations (15) and (A3), it follows that L ( Q j , l ) = NegBin ( α l , q j , l ) for each j J . According to Lemma A2 (with m = 1 and T 1 in that framework, see also Definitions A1–A3), we have the representation
NegBin ( α l , q j , l ) = CPoi α l ln ( 1 q j , l ) , Log ( q j , l ) .
Hence, ν j , l = α l ln ( 1 q j , l ) in Equations (24) and (25) and K j , l = Log ( q j , l ) in Equation (23), which for q j , l > 0 is the Panjer ( q j , l , q j , l , 1 ) -distribution and makes Panjer’s algorithm applicable for the calculation of X j , l c .

2.3.2. Generalized Tempered Stable Distribution

We consider another and more general case as usual in the CreditRisk + model that has already been considered in Gerhold et al. (2010). We allow one or several of the risk factors to have a τ -tempered α -stable distribution instead of a gamma distribution. This is a very flexible family of distributions and also allows for modelling of heavy tails. Stable distributions are denoted by S α ( σ , β , μ ) with α ( 0 , 2 ] , scale parameter σ > 0 , skewness parameter β [ 1 , 1 ] and shift parameter μ R , cf. e.g., (Samorodnitsky and Taqqu 1994, p. 9). If α ( 0 , 1 ) , β = 1 , and μ = 0 , then the support of S α ( σ , β , μ ) is on the non-negative real line (cf. (Samorodnitsky and Taqqu 1994, p. 15)). We can generalize such a distribution by a change of measure of Y S α ( σ , 1 , 0 ) as in (Gerhold et al. 2010, sct. 5.3) by introducing additionally a tempering parameter τ 0 and r N 0 and to obtain the distribution function
F α , σ , τ , r ( y ) = E Y r e τ Y 1 { Y y } E Y r e τ Y , y R .
In case r = 0 , this is a τ -tempered α -stable distribution (cf. (Rosiński 2007, Theorem 4.1)).
In this context, it is interesting to introduce the notation of the extended negative binomial distribution, cf. (Hess et al. 2002, p. 287):
Definition 4.
Let k N , 0 < p 1 , and α ( k , k + 1 ) . A random variable N has the extended negative binomial distribution ExtNegBin α , k , p if P [ N = n ] = 0 for n { 0 , 1 , , k 1 } and
P [ N = n ] = n + α 1 n p n ( 1 p ) α j = 0 k 1 j + α 1 j p j f o r n N w i t h n k ,
where the generalized binomial coefficient is given by
n + α 1 n = l = 1 n α + l 1 l = ( 1 ) n α n .
The probability-generating function of this distribution for the important case k = 1 is given by
G N ( z ) = 1 ( 1 p z ) α 1 ( 1 p ) α for z C with | z | 1 p ,
cf. e.g., (Gerhold et al. 2010, eq. (2.3) with an adjustment to our notation.
Remark 16.
As another case in Calculation Method 1, fix a risk factor l { 1 , , n } and assume that R l has a τ-tempered α-stable distribution, i.e., R l F α l , σ l , τ l , 0 . Then, the distribution of the random variables Q j , l for j J in Equation (15) can be written as a compound Poisson distribution and the application of (Gerhold et al. 2010, Lemma 5.10) provides a means of converting this distribution into a random sum with distributions in a Panjer ( a , b , k ) class, namely
Q j , l CPoi ( δ j , l , K j , l ) , j J
where
K j , l = ExtNegBin ( α l , 1 , q j , l ) i f q j , l > 0 , δ 0 i f q j , l = 0 ,
with parameters
δ j , l = γ α l , σ l τ l + g G λ g , j a g , l j α l τ l α l , γ α l , σ l = σ l α l cos ( α l π / 2 )
and
q j , l = g G λ g , j a g , l j τ l + g G λ g , j a g , l j i f τ l + g G λ g , j a g , l j > 0 , 0 o t h e r w i s e .
Then, Panjer’s recursion can be applied using (Gerhold et al. 2010, Theorem 4.5 and Algorithm 5.12) to this random sum. Altogether, Calculation Method 1 can be applied.
Remark 17.
In case the distribution of the risk factors is given in Equation (30) with r N , it is possible to apply (Gerhold et al. 2010, Algorithm 5.18).

3. Dependent Claim Numbers by Continuous Mixtures

In Section 2, we constructed dependence between the group-specific claim cause intensities ( Λ g ) g G by random non-negative linear combinations of the risk factors R 1 , , R n and by joint claims within risk groups. Now, we construct dependence between the risk factors by continuous mixture distributions. By choosing | J | = 1 and A the identity matrix, we can have this type of dependence directly. We turn one parameter of the distributions of the risk factors into a random variable, and, if it appears in several risk factors, we obtain dependence. As in Section 2, the independence between the components of the sum S given in Equation (7) is lost and there is also a need for an alternative to convolutions. Hence, we establish an alternative representation of the claim numbers with equal joint distribution if the claim cause intensities have either a gamma distribution or a τ -tempered α -stable distribution.

3.1. Continuous Compounding

In this subsection, we mix the risk factors R 1 , , R n with random variables T 1 , , T n that do not themselves have a mixture distribution. We consider a rather general result on mixture distributions to construct dependent claim numbers. As in the previous section, we let the group-specific claim cause intensities be random and non-negatively linearly dependent.
Theorem 2.
Let Assumption 1 be satisfied. Let T 1 , , T n be strictly positive random variables. Let α l , β l > 0 for all l { 1 , , n } and λ g , j 0 for all g G and j J . Then, we define two V m -valued random vectors S and M as follows:
(a) 
Let R 1 , , R n be random variables with conditional gamma distributions, conditionally independent given J , T 1 , , T n and let their shape parameters be randomized, i.e., for each l { 1 , , n }
L R l J , T 1 , , T n = a . s . L R l T l = a . s . Gamma ( α l T l , β l ) .
Then, let
L S J , R 1 , , R n , T 1 , , T n = a . s . CMPoi G , ( λ g , J Λ g ) g G , m , X .
(b) 
For each j J , consider random variables P j , 1 , , P j , n , which are conditionally independent given { J = j } , T 1 , , T n , such that for each l { 1 , , n }
L P j , l | 1 { J = j } , T 1 , , T n = a . s . L ( P j , l | T l ) = a . s . Poisson α l ln ( 1 q j , l ) T l
with q j , l given by Equation (28) and a random variable P j , 0 independent of { J = j } , T 1 , , T n , P j , 1 , , P j , n such that
L ( P j , 0 ) = Poisson R 0 g G λ g , j a g , 0 j .
Let { B j , l , h , k } h , k N with l { 0 , , n } be independent double-indexed sequences independent of { J = j } , P j , 0 , , P j , n , T 1 , , T n , consisting of i.i.d. random vectors such that for l { 0 , , n }
P [ B j , l , 1 , 1 = g ] = p g , j , l for g G ,
with p g , j , l as in Theorem 1. Let further { Y j , l , h } h N with l { 0 , , n } be independent sequences, independent of { J = j } , T 1 , , T n , P j , 0 , , P j , n , consisting of i.i.d. random variables such that
L ( Y j , l , 1 ) = Log ( q j , l ) for l { 1 , , n } , δ 1 for l = 0 ,
cf. Definition A3. Let { X l , g , h , k } h , k N be independent double-indexed sequences for l { 0 , , n } with identical distribution given by X g for g G independent of all previous random variables. Define the random vector M by
M = j J 1 { J = j } l = 0 n h = 1 P j , l k = 1 Y j , l , h X l , B j , l , h , k , h , k .
Then, ( M , T 1 , , T n ) and ( S , T 1 , , T n ) have the same distribution.
Remark 18.
We consider the same stochastic representation as in Remark 8 except for the conditional distribution of ( N g ) g G , which is given by
L ( N g | J , R 1 , , R n , T 1 , , T n ) = a . s . L ( N g | J , Λ g ) = a . s . Poisson ( λ g , J Λ g )
with ( Λ g ) g G as in Assumption A1.
Remark 19.
In setting (b) of Theorem 2, note that by Lemma A2 with m = 1 and T 1 in that framework
L h = 1 P j , l Y j , l , h | T l = a . s . NegBin α l T l , q j , l , j J , l { 1 , , n } ,
but the representation (38) is more convenient for Algorithm 2 below.
Proof. 
We apply Remark A1. Fix y ( V ) m , z R n . Let T = ( T 1 , , T n ) . First, we compute the characteristic function at ( y , z ) of the distribution of ( S , T ) . By an analogous argumentation as in the proof of Lemma 3, we obtain
φ ( S , T ) ( y , z ) = E e i z , T g G exp λ g , J ( 1 φ X g ( y ) ) l = 0 n R l a g , l J .
Conditioning on J , T 1 , , T n , using the conditional independence of R 1 , , R n given J , T 1 , , T n and using Equation (33) and the Laplace transform of the gamma distribution yields
φ ( S , T ) ( y , z ) = E [ e i z , T exp R 0 g G λ g , J a g , 0 J ( 1 φ X g ( y ) ) × l = 1 n β l β l + g G λ g , J a g , l J ( 1 φ X g ( y ) ) α l T l ] .
Now, consider the characteristic function at ( y , z ) of the distribution of the random vector ( M , T ) = ( M 1 , , M m , T 1 , , T n ) . Using the definition of M and partitioning J, we obtain
φ ( M , T ) ( y , z ) = j J E 1 { J = j } e i z , T exp i l = 0 n h = 1 P j , l k = 1 Y j , l , h y , X l , B j , l , h , k , h , k .
For each j J , using the fact that { B j , l , h , k } h , k N , { Y j , l , h } h N , and { X l , g , h , k } h , k N for l { 0 , , n } are i.i.d., independent, and independent of { J = j } , P j , 0 , , P j , n , T 1 , , T n and that { B j , l , h , k } h , k N and { X l , g , h , k } h , k N are independent of { Y j , l , h } h N gives
φ ( M , T ) ( y , z ) = j J E 1 { J = j } e i z , T l = 0 n E E exp i y , X l , B j , l , 1 , 1 , 1 , 1 Y j , l , 1 P j , l .
Conditioning on { J = j } , T 1 , , T n , using the conditional independence of P j , 1 , , P j , n given { J = j } , T 1 , , T n , the independence of P j , 0 of { J = j } , P j , 1 , , P j , n , T 1 , , T n for each j J , and using Equations (34) and (35) provides
φ ( M , T ) ( y , z ) = j J E [ 1 { J = j } e i z , T G P j , 0 ( y ) l = 1 n G P j , l ( y ) ]
with
G P j , 0 ( y ) = exp R 0 1 G Y j , 0 , 1 ( φ X B j , 0 , 1 , 1 ( y ) ) g G λ g , j a g , 0 j
and for l { 1 , , n }
G P j , l ( y ) = exp α l ln ( 1 q j , l ) T l 1 G Y j , l , 1 ( φ X B j , l , 1 , 1 ( y ) ) ,
where G Y j , l , 1 is the probability-generating function of the distribution of Y j , l , 1 . By Equation (36) and the law of total probability, we obtain for l { 0 , , n }
φ X B j , l , 1 , 1 ( y ) = g G p g , j , l φ X g ( y )
and by Equation (37) for l = 0
G P j , 0 ( y ) = exp R 0 1 k G p k , j , 0 φ X k ( y ) g G λ g , j a g , 0 j
and for l { 1 , , n }
G P j , l ( y ) = exp α l ln ( 1 q j , l ) T l 1 ln 1 q j , l g G p g , j , l φ X g ( y ) ln ( 1 q j , l ) .
In case q j , l = 0 for some j , l with l { 1 , , n } , it follows that ln ( 1 q j , l ) = 0 , and we have
exp α l ln ( 1 q j , l ) T l 1 G Y j , l , 1 ( φ X B j , l , 1 , 1 ( y ) ) = 1 .
If q j , 0 = 0 , this implies g G λ g , j a g , 0 j = 0 , and we have
exp R 0 1 G Y j , 0 , 1 ( φ X B j , 0 , 1 , 1 ( y ) ) g G λ g , j a g , 0 j = 1 .
If q j , l > 0 for l { 1 , , n } , then a simplification yields
exp α l ln ( 1 q j , l ) T l 1 ln 1 q j , l g G p g , j , l φ X g ( y ) ln ( 1 q j , l ) = exp α l T l ln ( 1 q j , l ) ln 1 g G p g , j , l φ X g ( y )
and if q j , 0 > 0 , we observe by Equation (36)
exp R 0 1 k G p k , j , 0 φ X k ( y ) g G λ g , j a g , 0 j = exp R 0 g G λ g , j a g , 0 j ( 1 φ X g ( y ) ) .
Altogether, we obtain
φ ( M , T ) ( y , z ) = j J E [ 1 { J = j } e i z , T exp R 0 g G λ g , j a g , 0 j ( 1 φ X g ( y ) ) × l = 1 n β l β l + g G λ g , j a g , l j ( 1 φ X g ( y ) ) α l T l ] ,
which completes the proof. □
Remark 20.
If T l for l { 1 , , n } has a beta distribution, then the distribution of P j , l for l { 1 , , n } and j J is a Poisson-beta distribution, also known as a general Waring distribution. (Hesselager 1996, Theorem 1, Example 3) provides a recursive algorithm for a compound distribution with such a counting distribution. Finally, an n-fold convolution becomes necessary. In this case, T 1 , , T n should be independent.
Now, we formulate a similar but different result to Theorem 2. Here, we consider τ -tempered α -stable distributions as distributions for the claim cause intensities.
Lemma 4.
Let Assumption 1 be satisfied. Let T 1 , , T n be non-negative random variables and T 0 a non-negative constant and σ l > 0 , τ l 0 and α l ( 0 , 1 ) for each l { 1 , , n } and λ g , j 0 for each g G and j J . Then, we define the two V m -valued random vectors S and M as follows:
(a) 
Let R 0 be a non-negative constant and let R 1 , , R n be random variables with conditional τ-tempered α-stable distributions, conditionally independent given J , T 1 , , T n and let their parameters be random, i.e., for each l { 1 , , n }
L R l J , T 1 , , T n = a . s . L R l T l = a . s . F α l , σ l , τ l T l , 0 .
Let R l = R l T l for l { 0 , , n } . Then, let
L S J , R 1 , , R n , T 1 , , T n = a . s . CMPoi G , ( λ g , J Λ g ) g G , m , X
(b) 
For each j J , consider random variables P j , 1 , , P j , n which are conditionally independent given { J = j } , T 1 , , T n such that for each l { 1 , , n }
L ( P j , l | 1 { J = j } , T 1 , , T n ) = a . s . L ( P j , l | T l ) = a . s . Poisson δ j , l T l α l ,
where δ j , l is given by Equation (31), and for each j J a random variable P j , 0 , which is independent of { J = j } , T 1 , , T n , P j , 1 , , P j , n such that
L P j , 0 = Poisson R 0 T 0 g G λ g , j a g , 0 j .
Let further { Y j , l , h } h N with l { 0 , , n } be independent sequences independent of { J = j } , T 1 , , T n , P j , 0 , , P j , n , consisting of i.i.d. random variables such that
L Y j , l , 1 = ExtNegBin α l , 1 , q j , l i f q j , l > 0 , δ 0 i f q j , l = 0 ,
for l { 1 , , n } where q j , l is given by Equation (32), and let Y j , 0 , 1 δ 1 . Let { B j , l , h , k } h , k N with l { 0 , , n } be independent double-indexed sequences independent of { ( Y j , 0 , h , , Y j , n , h ) } h N and { J = j } , T 1 , , T n , P j , 0 , , P j , n for j J , consisting of i.i.d. random variables such that
P [ B j , l , 1 , 1 = g ] = p g , j , l f o r g G ,
with p g , j , l [ 0 , 1 ] and g G p g , j , l = 1 satisfies p g , j , l k G λ k , j a k , l j = λ g , j a g , l j . Let { X l , g , h , k } h , k N be independent double-indexed sequences for l { 0 , , n } with identical distribution given by X g for g G independent of all previous random variables. Define the random vector M by
M = j J 1 { J = j } l = 0 n h = 1 P j , l k = 1 Y j , l , h X l , B j , l , h , k , h , k .
Then, ( M , T 1 , , T n ) and ( S , T 1 , , T n ) have the same distribution.
Remark 21.
We consider the same stochastic representation as in Remark 8 except for the conditional distribution of ( N g ) g G , which is given by
L ( N g | J , R 1 , , R n , T 1 , , T n ) = a . s . L ( N g | J , Λ g ) = a . s . Poisson ( λ g , J Λ g )
with ( Λ g ) g G as in Assumption A1.
Proof. 
We apply Remark A1. Fix y ( V ) m , z R n . Let T = ( T 1 , , T n ) . First, we compute the characteristic function at ( y , z ) of the distribution of ( S , T ) . By an analogous argumentation as in the proof of Lemma 3, we obtain
φ ( S , T ) ( y , z ) = E e i z , T g G exp λ g , J ( 1 φ X g ( y ) ) l = 0 n R l T l a g , l J .
Conditioning on J , T 1 , , T n , using the conditional independence of R 1 , , R n given J , T 1 , , T n and Equation (41) (and hence (Gerhold et al. 2010, eq. (5.25)) provides
φ ( S , T ) ( y , z ) = E [ e i z , T exp R 0 T 0 g G λ g , J a g , 0 J ( 1 φ X g ( y ) ) × l = 1 n exp γ α l , σ l T l α l τ l + g G λ g , J a g , l J ( 1 φ X g ( y ) ) α l ( T l τ l ) α l ] .
Now consider the characteristic function at ( y , z ) of the distribution of the random vector ( M , T ) = ( M 1 , , M m , T 1 , , T n ) . Using the definition of M and partitioning J yields
φ ( M , T ) ( y , z ) = j J E 1 { J = j } e i z , T exp i l = 0 n h = 1 P j , l k = 1 Y j , l , h y , X l , B j , l , h , k , h , k .
For each j J , using that { Y j , l , h } h N , { B j , l , h , k } h , k N , and { X l , g , h , k } h , k N are i.i.d. and independent of the random variables ( P j , 0 , , P j , n ) and l { 0 , , n } and that { Y j , l , h } h N are independent of { B j , l , h , k } h , k N and { X l , g , h , k } h , k N , we observe
φ ( M , T ) ( y , z ) = j J E 1 { J = j } e i z , T l = 0 n E E exp i y , X l , B j , l , 1 , 1 , 1 , 1 Y j , l , 1 P j , l .
By Equation (45) and due to the law of total probability, we obtain
φ ( M , T ) ( y , z ) = j J E 1 { J = j } e i z , T l = 0 n E g G p g , j , l φ X g ( y ) Y j , l , 1 P j , l .
Conditioning on { J = j } , T 1 , , T n , using the conditional independence of the random variables P j , 1 , , P j , n given { J = j } , T 1 , , T n and the independence of P j , 0 of { J = j } , T 1 , , T n for each j J and using Equations (42) and (43) yields
φ ( M , T ) ( y , z ) = j J E 1 { J = j } e i z , T G P j , 0 ( y ) l = 1 n G P j , l ( y )
with
G P j , 0 ( y ) = exp R 0 T 0 g G λ g , j a g , 0 j ( 1 φ X g ( y ) ) ,
where we have used already Equation (44) for l = 0 and simplified the term, and for l { 1 , , n }
G P j , l ( y ) = exp δ j , l T l α l 1 G Y j , l , 1 g G p g , j , l φ X g ( y ) ,
where G Y j , l , 1 denotes the probability-generating function of the distribution of Y j , l , 1 . If q j , l = 0 for some j , l , then, by Equation (44), G Y j , l , 1 ( z ) = 1 and hence
exp δ j , l T l α l 1 G Y j , l , 1 g G p g , j , l φ X g ( y ) = 1 .
If q j , l > 0 , by Equation (44) and q j , l g G p g , j , l φ X g ( y ) = g G λ g , j a g , l j φ X g ( y ) τ l + g G λ g , j a g , l j , we have
exp δ j , l T l α l 1 τ l τ l + g G λ g , j a g , l j α l 1 + 1 g G λ g , j a g , l j φ X g ( y ) τ l + g G λ g , j a g , l j α l 1 τ l τ l + g G λ g , j a g , l j α l = exp ( γ α l , σ l T l α l τ l + g G λ g , j a g , l j α l τ l α l × τ l α l + ( τ l + g G λ g , j a g , l j ( 1 φ X g ( y ) ) ) α l ( τ l + g G λ g , j a g , l j ) α l τ l α l ) ,
hence, since q j , l = 0 implies g G λ g , j a g , l j = 0 ,
φ ( M , T ) ( y , z ) = j J E [ 1 { J = j } e i z , T exp R 0 T 0 g G λ g , j a g , 0 j ( 1 φ X g ( y ) ) × l = 0 n exp γ α l , σ l T l α l τ l + g G λ g , j a g , l j ( 1 φ X g ( y ) ) α l τ l α l ] .
Unfortunately, due to the parameter α l ( 0 , 1 ) for l { 1 , , n } , it is not possible to choose, as before, the random variables T l with a gamma distribution or a τ -tempered α -stable distribution. By (Steutel and van Harn 2004, Example VI.12.8), T l α l is not even infinitely divisible if T l is gamma-distributed. It is only known that powers of τ -tempered α -stable distributed random variables T l α l (then the distribution of T l is by (Steutel and van Harn 2004, Proposition VI.5.7) and (Steutel and van Harn 2004, Proposition VI.5.26) a generalized gamma convolution) are infinitely divisible if α l > 1 , cf. (Steutel and van Harn 2004, Theorem VI.5.18) and (Steutel and van Harn 2004, Proposition VI.5.19(i). If α l ( 0 , 1 ) , then T l α l is only infinitely divisible if the corresponding characteristic function has no zeros (cf. (Lukacs 1970, Theorem 8.4.1)). Hence, we need to assume that T l α l has a gamma distribution or a τ -tempered α -stable distribution.

3.2. Evaluation of the Loss Distribution

In a constellation such as in Theorem 2, the evaluation of the distribution of the random sum S is an adaptation of the algorithm in (Gerhold et al. 2010, sct. 5.5) under certain assumptions and slighty different from the other dependence scenario.
Calculation Method 2.Consider the setting of Theorem 2. In addition, assume that the random variables T 1 , , T n are independent. Furthermore, similar to the approach of (Giese 2004, sct. 10.2), assume that T l Gamma ( σ l , ν l ) with shape parameter σ l > 0 and rate parameter ν l > 0 for every l { 1 , , n } . We apply Theorem 2 and obtain for the total loss S ^ = S 1 + + S m given as in Equation (7):
S ^ = i = 1 m g G h = 1 N g X i , g , h = d j J 1 { J = j } l = 0 n h = 1 P j , l k = 1 Y j , l , h i = 1 m X i , l , B j , l , h , k , h , k ,
where { ( X 1 , g , h , , X m , g , h ) } h N with g G and { ( X 1 , l , g , h , k , , X m , l , g , h , k ) } h , k N with g G and i { 1 , , m } are independent sequences of i.i.d. V m -valued random vectors with corresponding distribution X g . To calculate the distribution of the right-hand side of Equation (47), first define the group loss distribution X g s by Equation (21) and then the discrete mixture distributions X j , l by the right-hand side of Equation (22) for each j J and l { 0 , , n } .
In contrast to Calculation Method 1, the next step is to compute, for each j J , the compound logarithmic distribution
X j , l cl = L k = 1 Y j , l , 1 i = 1 m X i , l , B j , l , 1 , k , 1 , k = CLog ( q j , l , X j , l ) f o r l { 1 , , n } , X j , l f o r l = 0 .
If the support of X j , l is contained in N 0 d , then Panjer’s algorithm is applicable and numerically stable because Log ( q j , l ) = Panjer ( q j , l , q j , l , 1 ) for the case q j , l > 0 .
Let
s j , l = α l ln ( 1 q j , l ) ν l α l ln ( 1 q j , l ) , j J , l { 1 , , n } .
Since T 1 , , T n have gamma distributions, it follows from Lemma A1 (with m = 1 and T 1 in that framework) that
L ( P j , l ) = NegBin ( σ l , s j , l ) , j J , l { 1 , , n } .
According to Lemma A2 (with m = 1 and T 1 in that framework), cf. Equation (29),
P j , l CPoi ( ν j , l , K j , l ) w i t h ν j , l = σ l ln ( 1 s j , l ) a n d K j , l = Log ( s j , l ) .
Hence, we can proceed as in Calculation Method 1, where we have to replace Equation (23) by
X j , l c = CLog ( s j , l , X j , l cl ) , j J , l { 1 , , n } ,
and then use Equations (24)–(27).
Remark 22.
It should be pointed out that the random variables T 1 , , T n need not be independent. This is crucial for the construction of dependence. For instance, we could let T ˜ 1 , , T ˜ p be independent random variables and T ˜ 0 a non-negative constant. Define another random matrix B K = k K 1 { K = k } B k with K Ø an arbitrary finite set and B k [ 0 , ) n × ( p + 1 ) for k K and K a K -valued random variable. Then, let ( T 1 , , T n ) = B K ( T ˜ 0 , , T ˜ p ) . Then, an application of Theorem 1 should be inserted because P j , 0 , , P j , n are not independent.
Remark 23.
The algorithms do not only work this way if R 1 , , R n and T 1 , , T n , respectively, are gamma-distributed. They also work if R 1 , , R n or T 1 , , T n have a τ-tempered α-stable distribution. By an application of (Gerhold et al. 2010, Lemma 5.10), the severity distribution of the corresponding compound Poisson distribution of CPoi ( δ j , l , K j , l ) with parameters given in 16 is an extended negative binomial distribution. Thus, apply (Gerhold et al. 2010, Algorithm 5.12). A special case of this class of distributions is e.g., the inverse Gaussian distribution, which has also been described in (Sundt and Vernic 2009, p. 91). (Gerhold et al. 2010, Example 5.21) also present how to evaluate such a distribution.
Remark 24.
Using the convex combination for the evaluation of the portfolio distribution provides a claim size with a logarithmic distribution that might require high computational effort in the evaluation of Panjer’s recursion. This convex combination also requires a high computational effort in the evaluation of Panjer’s recursion with the Poisson distribution. However, this approach is preferable since the evaluation of Panjer’s recursion for a negative binomial distribution followed by several convolutions clearly demands more resources.

4. Numerical Illustration

In this section, we give examples of the impacts of different dependence structures in our model on the distribution of the portfolio loss. More precisely, this means we consider the impacts of different correlations. We will see that, given certain constraints originating from the extended CreditRisk + model, the resulting distributions differ significantly.

4.1. Dependence Modelling with Two Gamma-Distributed Risk Factors

We give an example of the application of Calculation Method 1 that implements the dependence structures given in Theorem 1. We consider four cases of dependence structures, and we use the following assumptions and parameters: in each case, we consider two independent and gamma-distributed risk factors, i.e., R i Gamma ( α i , β i ) with α i , β i > 0 , hence E R i = α i / β i and Var ( R i ) = α i / β i 2 for i = 1 , 2 . We consider two lines of business, meaning that G = { { 1 } , { 2 } } , and we simplify notation by writing g = 1 and g = 2 instead of g = { 1 } and g = { 2 } . A further constraint of our example is that the claim cause intensities satisfy E Λ i = E R i for i = 1 , 2 in each case of correlation. In each of these four cases, the distribution of the random sum S ^ from Equation (20) can be determined explicitly by computing the characteristic functions.
(a)
Positively correlated claim cause intensities: Take | J | = 1 , omit the index j, let λ 1 , λ 2 0 , take any a [ 0 , 1 ] , and define
A = 0 1 0 0 a ˜ 1 a with a ˜ = a E R 2 E R 1 = a α 2 β 1 α 1 β 2 .
By Equation (17),
Var Λ 1 = Var ( R 1 ) , Var Λ 2 = a ˜ 2 Var ( R 1 ) + ( 1 a ) 2 Var ( R 2 ) ,
and cov Λ 1 , Λ 2 = a ˜ Var R 1 . Hence, depending on a [ 0 , 1 ] , the correlation
corr ( Λ 1 , Λ 2 ) = a ˜ Var R 1 a ˜ 2 Var R 1 + ( 1 a ) 2 Var R 2
can take every value in [ 0 , 1 ] . Similar to Equation (11) using Equation (21), we obtain for every y V
φ S ^ ( y ) = E exp λ 1 R 1 ( 1 φ X 1 s ( y ) ) a ˜ λ 2 R 1 ( 1 φ X 2 s ( y ) ) × E exp ( 1 a ) λ 2 R 2 ( 1 φ X 2 s ( y ) ) = β 1 β 1 + ( λ 1 + a ˜ λ 2 ) ( 1 φ X 1 ( y ) ) α 1 β 2 β 2 + ( 1 a ) λ 2 ( 1 φ X 2 s ( y ) ) α 2
with claim size distribution X 1 : = Convex ( λ 1 , X 1 s ) , ( a ˜ λ 2 , X 2 s ) , see Definition 3. Using Definition A1(d) and Equation (A2), we see that L ( S ^ ) is the convolution of the two compound negative binomial distributions
CNegBin α 1 , λ 1 + a ˜ λ 2 β 1 + λ 1 + a ˜ λ 2 , X 1 and CNegBin α 2 , ( 1 a ) λ 2 β 2 + ( 1 a ) λ 2 , X 2 s .
(b)
Independent claim cause intensities: Take a = 0 in Case (a), hence X 1 = X 1 .
(c)
Comonotone claim cause intensities: Take a = 1 in Case (a). Then, the convolution is not necessary and L ( S ^ ) is a compound negative binomial distribution.
(d)
Negatively correlated claim cause intensities (cf. Example 1): We choose two dependence scenarios, i.e., J = { 1 , 2 } , and let J be a random variable independent of ( R 1 , R 2 ) such that p 1 : = P [ J = 1 ] ( 0 , 1 ) and p 2 = P [ J = 2 ] = 1 p 1 . Let λ i , j 0 for i , j = 1 , 2 . With b , c 0 define c i : = c E R i for i = 1 , 2 as well as
A 1 = 0 b 0 c 2 0 0 and A 2 = c 1 0 0 0 0 b .
Similar to Equation (11) using Equation (21), we observe that for every y V
φ S ^ ( y ) = p 1 E exp c 2 λ 2 , 1 R 0 ( 1 φ X 2 s ( y ) ) b λ 1 , 1 R 1 ( 1 φ X 1 s ( y ) ) + p 2 E exp c 1 λ 1 , 2 R 0 ( 1 φ X 1 s ( y ) ) b λ 2 , 2 R 2 ( 1 φ X 2 s ( y ) ) = p 1 exp c 2 λ 2 , 1 R 0 ( 1 φ X 2 s ( y ) ) β 1 β 1 + b λ 1 , 1 ( 1 φ X 1 s ( y ) ) α 1 + p 2 exp c 1 λ 1 , 2 R 0 ( 1 φ X 1 s ( y ) ) β 2 β 2 + b λ 2 , 2 ( 1 φ X 2 s ( y ) ) α 2 .
Hence, L ( S ^ ) is a mixture of
CPoi ( c 2 λ 2 , 1 R 0 , X 2 s ) CNegBin α 1 , b λ 1 , 1 β 1 + b λ 1 , 1 , X 1 s
and
CPoi ( c 1 λ 1 , 2 R 0 , X 1 s ) CNegBin α 2 , b λ 2 , 2 β 2 + b λ 2 , 2 , X 2 s ,
each of which is a convolution of a compound Poisson distribution with a compound negative binomial distributions. Note that, for i = 1 , 2
Λ i = c R 0 E R i + 1 { J = i } ( b R i c R 0 E R i ) , E Λ i J = c R 0 + 1 { J = i } ( b c R 0 ) E R i , E Λ i = c R 0 ( 1 p i ) + b p i E R i , cov ( Λ i , Λ j | J ) = b 2 cov R i , R j 1 { J = i } 1 { J = j } , cov ( Λ i , Λ j ) = E cov ( Λ i , Λ j | J ) + cov E Λ i J , E Λ j J ,
hence
Var ( Λ i ) = b 2 p i Var ( R i ) + p 1 p 2 ( b c R 0 ) 2 ( E R i ) 2 ,
cov ( Λ 1 , Λ 2 ) = p 1 p 2 ( b c R 0 ) 2 E R 1 E R 2 ,
where the last two equations can also be derived from Equation (17). To get c R 0 ( 1 p i ) + b p i = 1 for i = 1 , 2 such that Λ 1 and Λ 2 according to Equation (51) have the same expectations as in the previous cases, we can can take any b [ 0 , 2 ] and R 0 > 0 and define p 1 = p 2 = 1 / 2 and c = ( 2 b ) / R 0 , or we can take any R 0 > 0 and p 1 ( 0 , 1 ) and define b = 1 , p 2 = 1 p 1 and c = 1 / R 0 .
Remark 25.
In the case of independent, identically gamma-distributed risk factors R 1 and R 2 , meaning that α 1 = α 2 , β 1 = β 2 (we will drop these indices), we can not only have the same expectation but also the same variance of Λ 1 and Λ 2 in Cases (b), (c) and (d): Due to Equation (48), Var ( Λ i ) = α / β 2 in Cases (b) and (c). Due to Equation (52) with the choice b ( 0 , 2 ) , p 1 = p 2 = 1 / 2 and R 0 c = 2 b as discussed above,
Var ( Λ i ) = α 2 β 2 ( b 2 + 2 α ( b 1 ) 2 ) , i { 1 , 2 } ,
in Case (d). The quadratic equation b 2 + 2 α ( b 1 ) 2 = 2 has the solutions
b ± = 2 α ± 2 + 2 α 1 + 2 α .
If α ( 0 , 1 ) , then b < 0 and the solution cannot be used. If α 1 , then b [ 0 , 1 ) . However, b + ( 0 , 2 ) for all α > 0 , hence b + always works. Thus, given α , R 0 > 0 , one or two choices for the matrices A 1 and A 2 in Equation (50) are determined.
For the numerical example in the setting of Remark 25, see Figure 1 and Figure 2, we specify the remaining parameters as follows: we choose the risk factors such that E R 1 = E R 2 = 1 and Var ( R 1 ) = Var ( R 2 ) = 1 / 2 , e.g., α = β = 2 . In addition, we take R 0 = 1 . Thus, the two sets of parameters for the matrices A 1 and A 2 in Equation (50) are given by
b ± = 4 ± 6 5 and c ± = 2 b ± = 6 6 5 .
As a sanity check, Equation (48) for a ˜ = a { 0 , 1 } as well as Equation (52) for the above parameter sets give Var ( Λ i ) = 1 / 2 for i = 1 , 2 . Equation (49) for a { 0 , 1 } gives corr ( Λ 1 , Λ 2 ) = a , and by Equations (52) and (53) for the above parameters,
corr ( Λ 1 , Λ 2 ) = 14 4 6 25 0 . 1681 for b + 1 . 2899 , 0 . 9519 for b 0 . 3101 .
As Poisson parameters for Figure 1 and Figure 2, we choose λ i = λ i , j = 20 for i , j = 1 , 2 ; for the distribution of the (group) claim sizes, cf. Equation (21), we take X 1 s = δ ( 1 , 0 ) and X 2 s = δ ( 0 , 1 ) in Figure 1, and X 1 s = X 2 s = δ 1 in Figure 2.
An application of Theorem 1 gives us an equivalent representation such that we can apply Algorithm 1. For a better comparison, we put the probability mass functions into one graph, see Figure 2 with the first 120 values. We see interesting differences between the probability mass functions. Taking the probability mass function with independent claim cause intensities as an initial point of comparison, we observe that the probability mass function for negatively correlated claim cause intensities is a bit less, and, much less, respectively, light-tailed with a (much) taller maximum, whereas the probability mass function for comonotone, i.e., positively correlated claim cause intensities has a mass with a smaller maximum. We also observe that the probability mass function for the claim cause intensities with positive correlation has a heavier tail than the other probability mass functions.
While the probability mass functions for the Poisson case and the case with corr ( Λ 1 , Λ 2 ) 0.9519 look quite similar in Figure 2, the two-dimensional plot in Figure 1, which makes a difference between the losses of the two business lines, illustrates the tremendous difference between them.

4.2. Dependence Modelling with Two τ -Tempered α -stable-Distributed Risk Factors

For an illustration of Calculation Method 1 that implements the dependence structures given in Theorem 1, here with respect to Remark 16, we consider two risk factors meaning n = 2 . This implies that the considered risk factors have a τ -tempered α -stable distribution, i.e., R l F α l , σ l , τ l , 0 with l = 1 , , n . Following the notation of Remark 16, we need to compute a Panjer’s recursion for an extended negative binomial distribution, followed by a computation of Panjer’s recursion for a Poisson distribution and finally the distribution of a discrete mixture distribution.
We consider three different dependence structures: independent risk factors, positively correlated risk factors, and negatively correlated risk factors. We use the following assumptions and parameters: in each case, we consider two independent and τ -tempered α -stable-distributed risk factors, i.e., R l F α l , σ l , τ l , 0 with α l ( 0 , 1 ) , σ l > 0 and τ l > 0 . We consider two lines of business, meaning that G = { { 1 } , { 2 } } , and we simplify notation by writing g = 1 and g = 2 instead of g = { 1 } and g = { 2 } . The remaining parameters are chosen as follows: λ g , l = 20 , α l = 0 . 5 , σ l = 5 and τ l = 10 for g G and l = 1 , 2 . Furthermore, we let R 0 = 0 .
For the case of independent risk factors we take | J | = 1 and
A = 0 1 0 0 0 1 ,
and for the case of positive correlation, also with | J | = 1 , we choose
A = 0 1 0 0 1 0 .
In the case of negative correlation, we need two matrices:
A 1 = 0 1 0 0 0 0
and
A 2 = 0 0 0 0 0 1 .
The distribution of the required random variable J is P [ J = 1 ] = P [ J = 2 ] = 1 / 2 . For simplicity, we choose a degenerate distribution of the claim sizes X { 1 } s = X { 2 } s = δ 1 . The result of the computation is depicted in Figure 3.

4.3. Dependence Modeling with Three Gamma-Distributed Risk Factors

In this example, we extend the setting of the previous subsection and also apply Calculation Method 1. We want to study the impacts of the comonotone and two negative correlation structures between the claim cause intensities and compare with the case of independence. In each one of the four dependence structures, we consider three independent, gamma-distributed risk factors, i.e., R i Gamma ( α i , β i ) with α i , β i > 0 , hence E R i = α i / β i and Var ( R i ) = α i / β i 2 for i = 1 , 2 , 3 . We again consider m = 2 lines of business, but we allow for simultaneous claims in both of them, meaning that G = { { 1 } , { 2 } , { 1 , 2 } } . For comparison of the different dependence cases in our example, we impose the constraint that the claim cause intensities always satisfy
E [ Λ { i } ] = E R i for i = 1 , 2 and E [ Λ { 1 , 2 } ] = E R 3 .
(a)
Three positively correlated claim cause intensities: Take | J | = 1 , omit the index j, let λ { 1 } , λ { 2 } , λ { 1 , 2 } 0 , take any a , b , c [ 0 , 1 ] with b + c 1 , and define
A = 0 1 0 0 0 a ˜ 1 a 0 0 b ˜ c ˜ 1 b c with a ˜ = a E R 2 E R 1 , b ˜ = b E R 3 E R 1 and c ˜ = c E R 3 E R 2 ,
with the convention that ( Λ { 1 } , Λ { 2 } , Λ { 1 , 2 } ) = A ( R 0 , R 1 , R 2 , R 3 ) .
(b)
Independent claim cause intensities: Take a = b = c = 0 in Case (a).
(c)
Comonotone claim cause intensities: Take a = b = 1 and c = 0 in Case (a).
(d)
Positive correlation between the individual claim cause intensities Λ { 1 } and Λ { 2 } , and both of them negatively correlated to the joined claim cause intensity Λ { 1 , 2 } : We take two dependence scenarios, i.e., J = { 1 , 2 } , and let J be a random variable independent of ( R 1 , R 2 , R 3 ) such that p 1 : = P [ J = 1 ] ( 0 , 1 ) and p 2 = P [ J = 2 ] = 1 p 1 . Let λ g , j 0 for g G and j J . With b , c 0 , we define b ˜ = b E R 2 / E R 1 , c i = c E R i for i = 1 , 2 , 3 as well as
A 1 = 0 b 0 0 0 b ˜ 0 0 c 3 0 0 0 and A 2 = c 1 0 0 0 c 2 0 0 0 0 0 0 b .
Then, Λ { i } = c R 0 E R i + 1 { J = 1 } ( b R 1 / E R 1 c R 0 ) E R i for i = 1 , 2 and Λ { 1 , 2 } = c R 0 E R 3 + 1 { J = 2 } ( b R 3 c R 0 E R 3 ) . To satisfy the conditions in Equation (56), we need again that c R 0 ( 1 p i ) + b p i = 1 for i = 1 , 2 as in Case (d) of Section 4.1.
In the case of independent, identically gamma-distributed risk factors R 1 , R 2 , R 3 we can, in addition to the same expectations, also have the same variances of Λ { 1 } , Λ { 2 } and Λ { 1 , 2 } in Cases (b), (c) and (d) above; for details, see Remark 25.
For the numerical illustration in Figure 4 and the comparison with Figure 2, we take R 0 = 1 and independent R 1 , R 2 , R 3 Gamma ( α , β ) with α = β = 2 . Furthermore, we take λ g = λ g , j = 10 for g G and j = 1 , 2 . For g = { 1 , 2 } , we get claims in both business lines simultaneously, hence each business line gets 20 claims in expectation, which corresponds to the numerical example of the previous subsection. In Case (d), we take p 1 = p 2 = 1 / 2 and b ± , c ± as in (54), hence the correlation corr ( Λ { i } , Λ { 1 , 2 } ) for i = 1 , 2 is given by the right-hand side of Equation (55).

4.4. Dependence Modelling with Two Risk Factors and Mixture

For an illustration of Calculation Method 2, we choose the following contraints. Consider the conditions of Theorem 2. We assume two risk factors and a gamma-distributed random variable T with L T = a . s . Gamma σ , μ with σ > 0 and μ > 0 . Since such a single mixture distribution introduces a very strong measure of dependence, it is not useful to consider also negative correlation. For negative correlation, refer to (Rudolph 2014, Remark 4.11). For instance, an antithetic choice would suffice. Thus, it is sufficient to let J = { 1 } and we define G = { { 1 } , { 2 } } . The distribution of the conditionally independent risk factors is given by
L R g J , T = a . s . L R g T = a . s . Gamma α g T , β g , g = 1 , 2 .
We consider two dependence matrices: For a case, we call independent (since we do not impose further dependence)
A = 0 1 0 0 0 1
and another case that imposes further (very strong) dependence
A = 0 1 1 0 1 1 .
For simplicity, we choose a degenerate distribution of the claim sizes X { 1 } s = X { 2 } s = δ 1 . We let λ g = 20 , α g = 7 and β g = 2 for g = 1 , 2 . Furthermore, we choose σ = 0 . 3 and μ = 1 and R 0 = 0 . The result of the computation is shown in Figure 5.

Author Contributions

Supervision, U.S.; writing–original draft preparation, C.R.; writing–review and editing, C.R. and U.S. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Christian Doppler Research Association (CDG). The authors gratefully acknowledge the fruitful collaboration and support by UniCredit Bank Austria, Austrian Federal Financing Agency, COR & FJA (now msg life) and Österreichische Kontrollbank through CDG and the Christian Doppler Laboratory for Portfolio Risk Management (PRisMa Lab). The second author also gratefully acknowledges financial support from funds of the Österreichische Nationalbank (Anniversary Fund, project number 16549).

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

In this appendix, we provide for further reference some general representations of distributions when compounded or mixed that we needed in the preceding sections. Additionally, we introduce a dependence scenario between claim cause intensities and prove a corresponding formula for the probability-generating function of the distributions of Poisson mixture vectors.
Remark A1.
A probability-generating function uniquely determines a distribution, cf. (Kallenberg 2002, Theorem 5.3). The result is proven for characteristic functions but carries over to probability-generating functions. If V is a separable Banach space, see (Ledoux and Talagrand 1991, p. 39).
Definition A1
(Notation for compound distributions). Let N be an N 0 -valued random variable and { X h } h N an i.i.d. sequence of random variables, independent of N. Let K and F denote the distributions of N and X 1 , respectively.
(a) 
Let the distribution of the random sum S = h = 1 N X h be Compound ( K , F ) .
(b) 
If N has a Poisson distribution with expectation μ 0 , then S has a compound Poisson distribution and we write S CPoi μ , F .
(c) 
If N has a logarithmic distribution with parameter p [ 0 , 1 ) , cf. Definition A3, then S has a compound logarithmic distribution and we write S CLog p , F .
(d) 
If N has a negative binomial distribution NegBin ( α , p ) with parameters α > 0 and p [ 0 , 1 ) , cf. Definition A2, then S has a compound negative binomial distribution and we write S CNegBin ( α , p , F ) .
Results pointing out the relation between a compound Poisson distribution and a Poisson mixture distribution are already well known. It is interesting to consider compound distributions with negative binomial distributions instead. For this purpose, we give a definition of the negative multinomial distribution, cf. (Sibuya et al. 1964, sct. 2), which is a multivariate generalization of the negative binomial distribution.
Definition A2.
A random vector N = ( N 1 , , N m ) has an m-variate negative multinomial distribution with parameters α > 0 and p 1 , , p m 0 with p 0 : = 1 p 1 p m > 0 if
P [ N 1 = n 1 , , N m = n m ] = Γ α + i = 1 m n i Γ ( α ) i = 1 m n i ! p 0 α i = 1 m p i n i
for all n 1 , , n m N 0 . Denote this distribution on N 0 m by NegMult α ; p 1 , , p m . In the case m = 1 , we write NegBin ( α , p 1 ) . Note that NegBin ( α , 0 ) = δ 0 .
According to (Wishart 1949, eq. (2.1), the corresponding probability-generating function is
G N ( z 1 , , z m ) = p 0 1 i = 1 m p i z i α ,
defined for all ( z 1 , , z m ) C m with i = 1 m p i | z i | < 1 . For further reading on the negative multinomial distribution, we refer to Sibuya et al. (1964); Wishart (1949), or (Johnson et al. 1997, Chp. 36). This negative multinomial distribution can be also seen as a generalization in the sense of a generalization of the geometric distribution, see also Mai et al. (2013).
The following lemma generalizes the fact that a Poisson-gamma mixture is a negative binomial distribution and will be useful for further applications. It also generalizes (Sibuya et al. 1964, sct. 3.b) by allowing T to be random.
Lemma A1.
Let α , β > 0 and λ 1 , , λ m 0 . Let T be a strictly positive random variable and let Λ be a random variable such that L Λ T = a . s . Gamma ( α T , β ) . Let N = ( N 1 , , N m ) be a random vector with conditionally independent components given ( Λ , T ) such that
L N i Λ , T = a . s . Poisson ( λ i Λ ) , i { 1 , , m } .
Then, L N T = a . s . NegMult α T ; p 1 , , p m with p i = λ i / ( β + d = 1 m λ d ) for every i { 1 , , m } .
Proof. 
Apply Remark A1 and compute the probability-generating function of the conditional distribution of N given T. Conditioning on ( Λ , T ) and using the conditional independence of N 1 , , N m given ( Λ , T ) and the conditional distribution of N given ( Λ , T ) yields for z [ 0 , 1 ] m
G N | T ( z ) = a . s . E i = 1 m z i N i | T = a . s . E i = 1 m E z i N i | Λ , T | T = a . s . E i = 1 m e λ i Λ ( 1 z i ) | T .
Using the conditional distribution of Λ given T, we observe that
G N | T ( z ) = a . s . β + i = 1 m λ i ( 1 z i ) β α T = p 0 1 i = 1 m p i z i α T
with p 0 = β / ( β + d = 1 m λ d ) , see also Equation (A2).
A further result also exists on mixtures of the negative multinomial distribution. It is a generalization of a result that can be found in (Sibuya et al. 1964, sct. 3.d). Let us first introduce a generalization of the logarithmic distribution to multivariate dimensions taken from (Sibuya et al. 1964, sct. 3.d).
Definition A3.
Let p 1 , , p m 0 with 0 < p 0 : = 1 p 1 p m < 1 . A random vector ( L 1 , , L m ) has a multivariate logarithmic distribution if
P [ L 1 = l 1 , , L m = l m ] = ( l 1 ) ! ln p 0 i = 1 m p i l i l i ! , ( l 1 , , l m ) N 0 m { 0 } ,
where l = i = 1 m l i . We denote this distribution by MultLog p 1 , , p m . In the case m = 1 , we just write Log p and define Log 0 = δ 1 .
According to (Sibuya et al. 1964, eq. (3.2), its probability-generating function is given by
G L ( z ) = ln ( 1 i = 1 m p i z i ) ln p 0 , z [ 0 , 1 ] m .
Then, the aforementioned generalization is the following. It is a generalization of (Schröter 1995, p. 95/96) and a special case of it with m = 1 and T 1 can be also found in Ammeter (1948).
Lemma A2.
Let λ > 0 , p 0 ( 0 , 1 ) and p 1 , , p m 0 with p 0 + + p m = 1 . Let T be a strictly positive random variable and let N be a random variable such that
L N T = a . s . Poisson ( λ T ln p 0 ) .
Let { L h } h N be a sequence independent of ( N , T ) , consisting of i.i.d. m-dimensional random vectors with L 1 MultLog p 1 , , p m . Define S = h = 1 N L h . Then,
L S T = a . s . NegMult λ T ; p 1 , , p m .
Proof. 
Using Remark A1, we prove the claim by considering the probability-generating function G S | T of the conditional distribution of S given T for z ( 0 , 1 ] m with L h = ( L 1 , h , , L m , h )
G S | T ( z ) = a . s . E i = 1 m z i h = 1 N L i , h | T = a . s . E exp h = 1 N i = 1 m L i , h ln z i | T .
The sequence { L h } h N is i.i.d. and independent of ( N , T ) . Hence, using also Equation (A4), we observe
G S | T ( z ) = a . s . E E exp i = 1 m L i , 1 ln z i N | T = a . s . E ln ( 1 i = 1 m p i z i ) ln p 0 N | T .
Using Equation (A5) and a simplification provide
G S | T ( z ) = a . s . exp ( λ T ( ln p 0 ln 1 i = 1 m p i z i = a . s . 1 i = 1 m p i z i p 0 λ T ,
which completes the proof. □

References

  1. Ammeter, Hans. 1948. A generalization of the collective theory of risk in regard to fluctuating basic-probabilities. Scandinavian Actuarial Journal 31: 171–98. [Google Scholar] [CrossRef]
  2. Credit Suisse First Boston. 1997. CreditRisk+: A Credit Risk Management Framework. Available online: http://www.csfb.com/institutional/research/assets/creditrisk.pdf (accessed on 15 April 2014).
  3. Deshpande, Amogh, and Srikanth K. Iyer. 2009. The CreditRisk+ model with general sector correlations. Central European Journal of Operations Research 17: 219–28. [Google Scholar] [CrossRef]
  4. Fabiánová, Lucia, and Jozef Glova. 2015. Dependency and predictability of stock market returns: An empirical study of German capital market. Journal of Applied Economic Sciences 10: 1020–22. [Google Scholar]
  5. Gerhold, Stefan, Uwe Schmock, and Richard Warnung. 2010. A generalization of Panjer’s recursion and numerically stable risk aggregation. Finance and Stochastics 14: 81–128. [Google Scholar] [CrossRef]
  6. Giese, Götz. 2004. Dependent risk factors. In CreditRisk+ in the Banking Industry. Berlin: Springer, pp. 153–66. [Google Scholar]
  7. Gundlach, Matthias, and Frank Lehrbass. 2004. CreditRisk+ in the Banking Industry. Berlin: Springer. [Google Scholar]
  8. Hashorva, Enkelejd. 2011. A convolution identity for exchangeable risks. Albanian Journal of Mathematics 5: 41–43. [Google Scholar]
  9. Hess, Klaus Th., Anett Liewald, and Klaus D. Schmidt. 2002. An extension of Panjer’s recursion. ASTIN Bulletin: The Journal of the IAA 32: 283–97. [Google Scholar] [CrossRef] [Green Version]
  10. Hesselager, Ole. 1996. A recursive procedure for calculation of some mixed compound Poisson distributions. Scandinavian Actuarial Journal 1996: 54–63. [Google Scholar] [CrossRef]
  11. Jajuga, Krzysztof. 2016. Model risk in finance—Investor perspective. Transformations in Business and Economics 15: 289–304. [Google Scholar]
  12. Johnson, Norman Lloyd, Samuel Kotz, and Narayanaswamy Balakrishnan. 1997. Discrete Multivariate Distributions. New York: John Wiley & Sons Inc. [Google Scholar]
  13. Johnson, Norman Lloyd, Samuel Kotz, and Narayanaswamy Balakrishnan. 2000. Continuous Multivariate Distributions. In Models and Applications, 2nd ed. New York: Wiley-Interscience, vol. 1. [Google Scholar]
  14. Kallenberg, Olav. 2002. Foundations of Modern Probability, 2nd ed. New York: Springer. [Google Scholar]
  15. Ledoux, Michel, and Michel Talagrand. 1991. Probability in Banach Spaces. Berlin: Springer. [Google Scholar]
  16. Lindskog, Filip, and Alexander J. McNeil. 2003. Common Poisson shock models: Applications to insurance and credit risk modelling. ASTIN Bulletin: The Journal of the IAA 33: 209–38. [Google Scholar] [CrossRef] [Green Version]
  17. Lukacs, Eugene. 1970. Characteristic Functions, 2nd ed. New York: Hafner Publishing Co. [Google Scholar]
  18. Luo, Xiaolin, and Pavel V. Shevchenko. 2010. The t-copula with multiple parameters of degrees of freedom: Bivariate characteristics and application to risk management. Quantitative Finance 10: 1039–54. [Google Scholar] [CrossRef] [Green Version]
  19. Mai, Jan-Frederik, Matthias Scherer, and Natalia Shenkman. 2013. Multivariate geometric distributions, (logarithmically) monotone sequences, and infinitely divisible laws. Journal of Multivariate Analysis 115: 457–80. [Google Scholar] [CrossRef]
  20. Panjer, Harry H. 1981. Recursive evaluation of a family of compound distributions. ASTIN Bulletin: The Journal of the IAA 12: 22–26. [Google Scholar] [CrossRef] [Green Version]
  21. Rosiński, Jan. 2007. Tempering stable processes. Stochastic Processes and Their Applications 117: 677–707. [Google Scholar] [CrossRef] [Green Version]
  22. Rudolph, Cordelia. 2014. A Generalization of Panjer’s Recursion for Dependent Claim Numbers and an Approximation of Poisson Mixture Models. Ph.D. Thesis, Vienna University of Technology, Vienna, Austria. [Google Scholar]
  23. Samorodnitsky, Gennady, and M. S. Taqqu. 1994. Stable Non-Gaussian Random Processes. Stochastic Models with Infinite Variance. New York: Chapman & Hall. [Google Scholar]
  24. Schmock, Uwe. 2020. Modelling Dependent Credit Risks with Extensions of CreditRisk+ and Application to Operational Risk. Lecture Notes. Vienna: Vienna University of Technology. [Google Scholar]
  25. Schröter, Klaus Jürgen. 1995. Verfahren zur Approximation der Gesamtschadenverteilung. Karlsruhe: Verlag Versicherungswirtschaft. [Google Scholar]
  26. Sibuya, Masaaki, Isao Yoshimura, and Ryoichi Shimizu. 1964. Negative multinomial distribution. Annals of the Institute of Statistical Mathematics 16: 409–26. [Google Scholar] [CrossRef]
  27. Soltes, Vincent, and Jakub Danko. 2017. Proposal of creation of a portfolio with minimal risk. Investment Management and Financial Innovations 14: 107–15. [Google Scholar] [CrossRef] [Green Version]
  28. Steutel, Fred W., and Klaas van Harn. 2004. Infinite Divisibility of Probability Distributions on the Real Line. New York: Marcel Dekker Inc., vol. 259. [Google Scholar]
  29. Sundt, Bjørn. 1999. On multivariate Panjer recursions. ASTIN Bulletin: The Journal of the IAA 29: 29–45. [Google Scholar] [CrossRef] [Green Version]
  30. Sundt, Bjørn, and William S. Jewell. 1981. Further results on recursive evaluation of compound distributions. ASTIN Bulletin: The Journal of the IAA 12: 27–39. [Google Scholar] [CrossRef] [Green Version]
  31. Sundt, Bjørn, and Raluca Vernic. 2009. Recursions for Convolutions and Compound Distributions with Insurance Applications. Berlin: Springer. [Google Scholar]
  32. Willmot, Gordon E. 1988. Sundt and Jewell’s familiy of discrete distributions. ASTIN Bulletin: The Journal of the IAA 18: 17–29. [Google Scholar] [CrossRef] [Green Version]
  33. Willmot, Gordon E., and Harry H. Panjer. 1987. Difference equation approaches in evaluation of compound distributions. Insurance: Mathematics and Economics 6: 43–56. [Google Scholar] [CrossRef]
  34. Wishart, John. 1949. Cumulants of multivariate multinomial distributions. Biometrika 36: 47–58. [Google Scholar] [CrossRef]
  35. Wüthrich, Mario V., and Michael Merz. 2008. Stochastic Claims Reserving Methods in Insurance. Chichester: John Wiley & Sons Inc. [Google Scholar]
  36. Zoričák, Martin, Peter Gnip, Peter Drotár, and Vladimír Gazda. 2019. Bankruptcy prediction for small- and medium-sized companies using severely imbalanced datasets. Economic Modelling 84: 165–76. [Google Scholar] [CrossRef]
Figure 1. Highly smoothed probability mass functions of the N 0 2 -valued portfolio loss S ^ corresponding to Figure 2 for the Poisson case (central truncated graph), the comonotone case (with a hump along the diagonal) and the case of corr ( Λ 1 , Λ 2 ) 0.9519 (two symmetric humps near the corners). We take V = R 2 . The (group) claim sizes equal ( 1 , 0 ) and ( 0 , 1 ) in these cases, i.e., X 1 s = δ ( 1 , 0 ) and X 2 s = δ ( 0 , 1 ) .
Figure 1. Highly smoothed probability mass functions of the N 0 2 -valued portfolio loss S ^ corresponding to Figure 2 for the Poisson case (central truncated graph), the comonotone case (with a hump along the diagonal) and the case of corr ( Λ 1 , Λ 2 ) 0.9519 (two symmetric humps near the corners). We take V = R 2 . The (group) claim sizes equal ( 1 , 0 ) and ( 0 , 1 ) in these cases, i.e., X 1 s = δ ( 1 , 0 ) and X 2 s = δ ( 0 , 1 ) .
Risks 08 00043 g001
Figure 2. Probability mass functions of the N 0 -valued portfolio loss S ^ for different dependence structures for the two claim cause intensities Λ 1 and Λ 2 satisfying E Λ 1 = E Λ 2 = 1 and Var Λ 1 = Var Λ 2 = 1 / 2 as discussed in Section 4.1. We take V = R . The claim size equals 1 in all cases, i.e., X 1 s = X 2 s = δ 1 . For comparison, also the case Λ 1 = Λ 2 1 , yielding a Poisson distribution with expectation 40, is shown.
Figure 2. Probability mass functions of the N 0 -valued portfolio loss S ^ for different dependence structures for the two claim cause intensities Λ 1 and Λ 2 satisfying E Λ 1 = E Λ 2 = 1 and Var Λ 1 = Var Λ 2 = 1 / 2 as discussed in Section 4.1. We take V = R . The claim size equals 1 in all cases, i.e., X 1 s = X 2 s = δ 1 . For comparison, also the case Λ 1 = Λ 2 1 , yielding a Poisson distribution with expectation 40, is shown.
Risks 08 00043 g002
Figure 3. Probability mass functions of the N 0 -valued portfolio loss S ^ for different dependence structures for the two claim cause intensities Λ 1 and Λ 2 . We take V = R . The claim size equals 1 in all cases, i.e., X 1 s = X 2 s = δ 1 .
Figure 3. Probability mass functions of the N 0 -valued portfolio loss S ^ for different dependence structures for the two claim cause intensities Λ 1 and Λ 2 . We take V = R . The claim size equals 1 in all cases, i.e., X 1 s = X 2 s = δ 1 .
Risks 08 00043 g003
Figure 4. Probability mass functions of the portfolio loss S ^ for four different dependence structures for the three claim cause intensities Λ { 1 } , Λ { 2 } and Λ { 1 , 2 } satisfying E Λ g = 1 and Var Λ g = 1 / 2 for all g G = { { 1 } , { 2 } , { 1 , 2 } } as discussed in Section 4.3. The group claim size, cf. Equation (21), equals 1 for g = { 1 } , { 2 } and 2 for g = { 1 , 2 } , i.e., X { 1 } s = X { 2 } s = δ 1 and X { 1 , 2 } s = δ 2 . The given correlation refers to corr ( Λ { i } , Λ { 1 , 2 } ) for i = 1 , 2 , see the right-hand side of Equation (55). For comparison with Section 4.1, the corresponding probability mass functions from Figure 2 are indicated as solid lines.
Figure 4. Probability mass functions of the portfolio loss S ^ for four different dependence structures for the three claim cause intensities Λ { 1 } , Λ { 2 } and Λ { 1 , 2 } satisfying E Λ g = 1 and Var Λ g = 1 / 2 for all g G = { { 1 } , { 2 } , { 1 , 2 } } as discussed in Section 4.3. The group claim size, cf. Equation (21), equals 1 for g = { 1 } , { 2 } and 2 for g = { 1 , 2 } , i.e., X { 1 } s = X { 2 } s = δ 1 and X { 1 , 2 } s = δ 2 . The given correlation refers to corr ( Λ { i } , Λ { 1 , 2 } ) for i = 1 , 2 , see the right-hand side of Equation (55). For comparison with Section 4.1, the corresponding probability mass functions from Figure 2 are indicated as solid lines.
Risks 08 00043 g004
Figure 5. Probability mass functions of the N 0 -valued portfolio loss S ^ for different dependence structures for the two claim cause intensities Λ 1 and Λ 2 satisfying E Λ 1 = E Λ 2 = 1.05 . We take V = R . The claim size equals 1 in all cases, i.e., X 1 s = X 2 s = δ 1 .
Figure 5. Probability mass functions of the N 0 -valued portfolio loss S ^ for different dependence structures for the two claim cause intensities Λ 1 and Λ 2 satisfying E Λ 1 = E Λ 2 = 1.05 . We take V = R . The claim size equals 1 in all cases, i.e., X 1 s = X 2 s = δ 1 .
Risks 08 00043 g005

Share and Cite

MDPI and ACS Style

Rudolph, C.; Schmock, U. Multivariate Collective Risk Model: Dependent Claim Numbers and Panjer’s Recursion. Risks 2020, 8, 43. https://doi.org/10.3390/risks8020043

AMA Style

Rudolph C, Schmock U. Multivariate Collective Risk Model: Dependent Claim Numbers and Panjer’s Recursion. Risks. 2020; 8(2):43. https://doi.org/10.3390/risks8020043

Chicago/Turabian Style

Rudolph, Cordelia, and Uwe Schmock. 2020. "Multivariate Collective Risk Model: Dependent Claim Numbers and Panjer’s Recursion" Risks 8, no. 2: 43. https://doi.org/10.3390/risks8020043

APA Style

Rudolph, C., & Schmock, U. (2020). Multivariate Collective Risk Model: Dependent Claim Numbers and Panjer’s Recursion. Risks, 8(2), 43. https://doi.org/10.3390/risks8020043

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