The paper deals with bonus–malus systems with different claim types and varying deductibles. The premium relativities are softened for the policyholders who are in the malus zone and these policyholders are subject to per claim deductibles depending on their levels in the bonus–malus scale and the types of the reported claims. We introduce such bonus–malus systems and study their basic properties. In particular, we investigate when it is possible to introduce varying deductibles, what restrictions we have and how we can do this. Moreover, we deal with the special case where varying deductibles are applied to the claims reported by policyholders occupying the highest level in the bonus–malus scale and consider two allocation principles for the deductibles. Finally, numerical illustrations are presented.
One of the main tasks of an actuary is to design a tariff structure that fairly distributes the total risk of potential losses among policyholders. To this end, he often has to grade all policyholders into risk classes such that all policyholders belonging to the same class pay the same premium. Rating systems penalizing policyholders responsible for one or more accidents by premium surcharges (or maluses), and rewarding claim-free policyholders by giving them discounts (or bonuses) are now in force in many developed countries. Such systems, which are often called bonus–malus systems, aim to assess individual risks better.
The amount of premium is adjusted each year on the basis of the individual claims experience. In practice, a bonus–malus scale consists of a finite number of levels and each of the levels has its own relative premium. After each year, the policyholder moves up or down according to the transition rules and the number of claims reported during the current year. Thus, bonus–malus systems also encourage policyholders to be careful. Note that the premium relativities are traditionally computed using a quadratic loss function. This method is proposed by Norberg in his pioneering work [14] on segmented tariffs. Alternatively, Denuit and Dhaene [2] use an exponential loss function to compute the relativities.
In most of the commercial bonus–malus systems used by insurance companies, knowing the current level and the number of claims during the current period suffices to determine the next level in the scale. Therefore, the future level depends only on the present and not on the past. The numbers of claims in different years are usually assumed to be independent. So the trajectory of each policyholder in the bonus–malus scale can be considered as a Markov chain. For details and more information concerning bonus–malus systems, we refer the reader to [3, 10, 19]. In particular, [3] presents a comprehensive treatment of the various experience rating systems and their relationships with risk classification.
As pointed out by many authors, traditional bonus–malus systems suffer from two considerable drawbacks:
The claim amounts are not taken into account. So a posteriori corrections depend only on the number of claims. In this case, policyholders who had accidents with small or large claims are penalized unfairly in the same way. In particular, this breeds bonus hunger when policyholders cover small claims themselves in order to avoid future premium increases.
At any time the policyholders may leave the insurance company without any further financial penalties. Thus, bonus–malus systems create the possibility of malus evasion, i.e. the situation when the policyholders leave the insurance company to avoid premium increase because of reported claims.
An alternative approach, which, at least theoretically, eliminates the second drawback, is proposed by Holtan [9]. He suggests the use of very high deductibles that may be borrowed by the policyholders in the insurance company. The deductibles are assumed to be constant for all policyholders, i.e. independent of the level they occupy in the bonus–malus system at the time of claim occurrence. Although technically acceptable, this approach obviously causes considerable practical problems. Practical consequences of Holtan’s proposal are investigated by Lemaire and Zi [11]. Particularly, it is shown that the introduction of high deductibles increases the variability of payments and the efficiency of the rating systems for most policyholders.
Bonus–malus systems with varying deductibles are introduced by Pitrebois, Denuit and Walhin [17]. Specifically, the a posteriori premium correction induced by the bonus–malus system is replaced by a deductible (in whole or in part). To each level of the bonus–malus system in the malus zone is assigned an amount of deductible, which is applied to the claims filed during the coverage period either annually or claim by claim. Relative premiums at high levels of bonus–malus systems are often very large and the systems can be softened by introducing deductibles. The insurance company compensates the reduced penalties in the malus zone with the deductibles paid by policyholders who report claims being in the malus zone. This can be commercially attractive since the policyholders are penalized only if they report claims in the future.
As pointed out in [17], combining bonus–malus systems with varying deductibles presents a number of advantages. Firstly, the policyholders will do all their best to prevent or at least decrease the losses. Secondly, even if the policyholder leaves the company after a claim, he has to pay for the deductible. Thirdly, relative premiums and amounts of deductibles may be tuned in an optimal way in order to attract the policyholders. The numerical illustrations show that the mixed case (reduced relative premiums combined with per claim deductibles) gives the best results. The amounts of deductibles are moderate in this case. Nevertheless, this approach does not eliminate the first drawback mentioned above.
To eliminate the first drawback, a few other approaches have been proposed. Bonus–malus systems involving different claim types are designed in [18]. Each claim type induces a specific penalty for the policyholder. In particular, claim amounts can be taken into account in this way. For some early results in this direction, see also [10, 15, 16].
The next approach is based on a generalization of results obtained by Dionne and Vanasse [4, 5], who propose a bonus–malus system that integrates a priori and a posteriori information on an individual basis. Specifically, the system is derived as a function of the years that the policyholder is in the portfolio, the number of claims and his individual characteristics. Frangos and Vrontos [6] expand the frame developed in [4, 5] and propose a generalized bonus–malus system that takes into consideration simultaneously the policyholder’s characteristics, the number of his claims and the exact amount of each claim. Particularly, it is assumed that we have all information about the claim frequency history and the claim amount history for each policyholder for the time period he is in the portfolio. In this generalized bonus–malus system, the premium is a function of the years that the policyholder is in the portfolio, his number of claims, the amount of each claim and the significant a priori rating variables. Therefore, the future premium depends on the past (all policyholder’s history) and is actually calculated individually for each policyholder. This approach is extended and developed in [12, 13, 20].
Another approach that takes into account claim amounts is proposed by Bonsdorff [1]. The author considers a general framework for a bonus–malus system based on the number of claims during the previous year and the total amount of claims during the previous year. The set of the bonus levels is some interval and the transitions between the levels are determined by these characteristics.
Gómez-Déniz, Hernández-Bastida and Fernández-Sánchez [8] obtain expressions that can be used to compute bonus–malus premiums based on the distribution of the total claim amount but not on the claims which produce the amounts.
Another modification of traditional bonus–malus systems, which take into account only the number of claims, is considered by Gómez-Déniz in the recent paper [7]. The author presents a statistical model, which distinguishes between two different claim types, incorporating a bivariate distribution based on the assumption of dependence.
The present paper deals with the case where both penalty types are used for the policyholders who are in the malus zone. Specifically, the premium relativities, which are computed using a quadratic loss function, are softened and the policyholders who are in the malus zone are subject to per claim deductibles. The mixed bonus–malus system combining premium relativities and deductibles is expected to be the most relevant in practice (see [3, 17]).
We try to eliminate both drawbacks mentioned above. To take into account claim amounts, we consider different claim types and use the multi-event bonus–malus systems introduced in [18]. To eliminate the second drawback, we introduce varying deductibles for the policyholders who are in the malus zone. The deductibles depend on the level of the policyholder in the bonus–malus scale and the types of the reported claims. Such bonus–malus systems present a number of advantages and seem to be very attractive for policyholders. Namely, policyholders reporting small and large claims are not penalized in the same way. This helps to avoid or at least decrease bonus hunger. Moreover, all advantages mentioned in [17] are also in force.
The rest of the paper is organized as follows. In Section 2, we describe the multi-event bonus–malus systems introduced in [18]. In Section 3, we introduce varying deductibles in such bonus–malus systems and study basic properties of such systems. In particular, we investigate when it is possible to introduce varying deductibles in the bonus–malus systems with different claim types, what restrictions we have and how we can do this. Section 4 deals with the special case where varying deductibles are applied to the claims reported by policyholders occupying the highest level in the bonus–malus scale. We consider two allocation principles for the deductibles, which seems to be natural and fair for policyholders. In Section 5, we consider an example of such a bonus–malus system, deal with exponentially distributed claim sizes and present numerical illustrations. Section 6 completes the paper.
Bonus–malus systems with different claim types
To take into consideration different claim types, we use the multi-event bonus–malus systems introduced in [18] (see also [3]).
Let us pick a policyholder at random from the portfolio and denote by N the number of claims reported by the policyholder during the year. In what follows, we assume that there is no a priori risk classification (or we work inside a specified rating cell). Denote by λ>00$]]> the a priori annual expected claim frequency. Let Θ be the (unknown) accident proneness of this policyholder, i.e. Θ represents the residual effect of unobserved characteristics. The risk profile of the portfolio is described by the distribution function FΘ of Θ. It is usually assumed that P[Θ≥0]=1 and E[Θ]=1. Thus, the actual (unknown) annual claim frequency of this policyholder is λΘ.
We assume that the number of claims N is mixed-Poisson distributed. To be more precise, the conditional probability mass function of N is given by
P[N=j|Θ=θ]=(λθ)jj!e−λθ,j≥0.
Hence, the unconditional probability mass function of N is given by
P[N=j]=∫0+∞P[N=j|Θ=θ]dFΘ(θ),j≥0.
We introduce m+1 different claim types reported by the policyholder. Each claim type induces a specific penalty for the policyholder, which will be described below. The type of the given claim is determined by the claim amount C. We assume that all claim amounts are independent and identically distributed, and claim amounts and claim frequencies are mutually independent. The claims are classified according to a multinomial scheme. Let 0<c1∗<c2∗<⋯<cm∗<∞. We choose these numbers so that all claims of size less than or equal to c1∗ are considered as claims of type 0, all claims of size from the interval (c1∗,c2∗] are considered as claims of type 1 and so on; finally, all claims of size greater than cm∗ are considered as claims of type m. Let
q0=P[C≤c1∗],q1=P[c1∗<C≤c2∗],q2=P[c2∗<C≤c3∗],…,qm=P[C>cm∗].{c_{m}^{\ast }}\big].\end{array}\]]]>
Thus, each time a claim is reported, it is classified in one of the m+1 possible categories with probabilities q0, q1, …, qm. It is clear that ∑i=0mqi=1. Moreover, it is natural to choose the ci∗ so that all qi are strictly positive.
Denote by Ni the number of claims of type i, 0≤i≤m. Therefore, for a given Θ, the random variables N0, N1, …, Nm are mutually independent and the corresponding conditional probability mass function is given by
P[Ni=j|Θ=θ]=(λθqi)jj!e−λθqi,j≥0,0≤i≤m.
The bonus–malus scale is assumed to have s+1 levels numbered from 0 to s. A higher level number indicates a higher premium. In particular, the policyholders who are at level 0 enjoy the maximal bonus. A specified level is assigned to a new policyholder. Each claim-free year is rewarded by a bonus point, i.e. the policyholder goes one level down. Each claim type entails a specific penalty expressed as a fixed number of levels per claim. It is natural to assume that larger claims entail more severe penalties.
We suppose that knowing the present level and the number of claims of each type filed during the present year suffices to determine the level to which the policyholder is transferred. So the bonus–malus system may be represented by a Markov chain.
Let pl0l(λθ;q) be the probability of moving from level l0 to level l for a policyholder with annual mean claim frequency λθ and vector of probabilities q=(q0,q1,…,qm)T, where 0≤l0≤s, 0≤l≤s and qi is the probability that the claim is of type i. Denote by P(λθ;q) the one-step transition matrix, i.e.
P(λθ;q)=p00(λθ;q)p01(λθ;q)…p0s(λθ;q)p10(λθ;q)p11(λθ;q)…p1s(λθ;q)⋮⋮⋱⋮ps0(λθ;q)ps1(λθ;q)…pss(λθ;q).
Taking the nth power of P(λθ;q) yields the n-step transition matrix with elements pl0l(n)(λθ;q). Here pl0l(n)(λθ;q) is the probability of moving from level l0 to level l in n transitions.
The transition matrix P(λθ;q) is assumed to be regular, i.e. there is some integer n0≥1 such that all entries of P((λθ;q))n0 are strictly positive. Consequently, the Markov chain describing the trajectory of a policyholder with expected claim frequency λθ and vector of probabilities q is ergodic and thus possesses a stationary distribution
π(λθ;q)=(π0(λθ;q),π1(λθ;q),…,πs(λθ;q))T.
Here πl(λθ;q) is the stationary probability for a policyholder with annual mean claim frequency λθ to be at level l, i.e.
πl(λθ;q)=limn→∞pl0l(n)(λθ;q).
The stationary probabilities πl(λθ;q) can be obtained directly (see [3, 19]). Indeed, since the matrix P(λθ;q) is regular, the matrix I−P(λθ;q)+E is invertible and π(λθ;q) is given by
πT(λθ;q)=eT(I−P(λθ;q)+E)−1,
where e is a column vector of 1s, E is an (s+1)×(s+1) matrix consisting of s+1 column vectors e and I is an (s+1)×(s+1) identity matrix.
Next, let L be the level occupied in the scale by a randomly selected policyholder and πl be the proportion of policyholders at level l once the steady state has been reached. Then
πl=P[L=l]=∫0+∞πl(λθ;q)dFΘ(θ),0≤l≤s.
It is easily seen that ∑l=0sπl=1.
We denote by rl the premium relativities associated with level l. The meaning is that a policyholder occupying level l pays the premium equal to λrlE[C]. Note that here and below we consider only net premiums. In what follows, we assume that E[C]<∞ and the distribution function FC of C is continuous.
To compute the premium relativities rl, we use a quadratic loss function as proposed in [14] (see also [3, 10]). To this end, we minimize the expected squared difference between the “true” relative premium Θ and the relative premium rL applicable to this policyholder after the stationary state has been reached, i.e. we minimize E[(Θ−rL)2]. The solution to this problem is given by
rl=∫0+∞θπl(λθ;q)dFΘ(θ)∫0+∞πl(λθ;q)dFΘ(θ)
(see, e.g., [3, pp. 185–186] for the details).
Varying deductibles in the bonus–malus systems with different claim types
The policyholder occupying level l in the bonus–malus systems described in Section 2 should pay λrlE[C]. We suppose that the relative premiums of the policyholders who are in the bonus zone, i.e. with rl≤1, are unchanged. So they pay premiums equal to λrlE[C] and are subject to no further penalties. Let s0=min{l:rl>1}1\}$]]>. The relative premiums of the policyholders who are in the malus zone, i.e. with rl>11$]]>, or equivalently occupying level l such that s0≤l≤s, are softened in the following way. Instead of paying the premium equal to λrlE[C], the policyholder who is at level l pays a reduced premium equal to (1−αl)λrlE[C] for some specified αl≥0 depending on level l. We suppose that αl satisfy the following assumption.
1≤(1−αs0)rs0≤(1−αs0+1)rs0+1≤⋯≤(1−αs)rs.
Note that Assumption 1 implies that the policyholder occupying a higher level in the bonus–malus scale pays a higher reduced premium, which is not less than the basic premium λE[C]. Moreover, from Assumption 1 we have 0≤αl≤1−1/rl for all l such that s0≤l≤s. The case αl=0 corresponds to the situation when the policyholder occupying level l pays premium equal to λrlE[C] and is subject to no further penalties. If αl=1−1/rl, then the policyholder pays only the basic premium λE[C] and has to pay something for claims in the future. To compensate the reduced premium, the policyholder is subject to per claim deductible, i.e. applied to each reported claim separately, equal to dl,i depending on level l occupied in the malus zone and claim type i. We impose the following natural restrictions on dl,i.
For all l such that s0≤l≤s, we have
0≤dl,0≤c1∗,0≤dl,1≤c1∗,0≤dl,2≤c2∗,…,0≤dl,m≤cm∗.
For every fixed l such that s0≤l≤s, we have
dl,0≤dl,1≤⋯≤dl,m.
For every fixed i such that 0≤i≤m, we have
ds0,i≤ds0+1,i≤⋯≤ds,i.
Assertion (i) of Assumption 2 means that if a claim reported is of type i, where 0≤i≤m, then the deductible applied to it is strictly less than the claim amount, i.e. the insurance company covers at least some part of the losses. Next, assertion (ii) implies that larger claims are subject to higher deductibles. Finally, assertion (iii) means that the higher level in the bonus–malus scale implies a higher deductible for each specified claim amount.
For policyholders occupying level l, the deductibles dl,0, dl,1, …, dl,m are found using the indifference principle (see [3, 17]): for this group of policyholders, the part αl of the penalties induced by the bonus–malus system is on average equal to the total amount of deductibles paid by these policyholders. Thus, the indifference principle for the policyholder occupying level l can be written in the following way:
λrlE[C]=(1−αl)λrlE[C]+λrl(E[C|C≤dl,0]P[C≤dl,0]+dl,0P[dl,0<C≤c1∗]+dl,1P[c1∗<C≤c2∗]+…+dl,mP[C>cm∗]),s0≤l≤s.{c_{m}^{\ast }}\big]\big),\hspace{1em}s_{0}\le l\le s.\end{array}\]]]>
On the left-hand side of (5), we have the premium paid by the policyholder occupying level l in the bonus–malus system described in Section 2. On the right-hand side of (5), we have the expected amount paid by this policyholder in the bonus–malus system with varying deductibles. This amount consists of the reduced premium and the expected amount of penalties induced by deductibles.
Equation (5) can be rewritten as
αlE[C]=E[C|C≤dl,0]P[C≤dl,0]+dl,0P[dl,0<C≤c1∗]+dl,1P[c1∗<C≤c2∗]+⋯+dl,mP[C>cm∗],s0≤l≤s,{c_{m}^{\ast }}\big],\hspace{1em}s_{0}\le l\le s,\end{array}\]]]>
which is equivalent to
αlE[C]=E[C|C≤dl,0]P[C≤dl,0]+dl,0(q0−P[C≤dl,0])+dl,1q1+⋯+dl,mqm,s0≤l≤s.
Thus, in order to introduce varying deductibles in the bonus–malus system, we have to choose αl, dl,0, dl,1, …, dl,m satisfying (6) and Assumptions 1 and 2 for all l such that s0≤l≤s. In what follows, we call any such combination of αl, dl,0, dl,1, …, dl,m, where s0≤l≤s, by a solution to (6).
The right-hand side of (6) is a non-decreasing function of each of the variablesdl,0,dl,1, …,dl,m.
The assertion of Lemma 1 is evident for variables dl,1, …, dl,m. We now show it for dl,0. Let 0≤dl,0′≤dl,0″≤c1∗. Then we get
E[C|C≤dl,0″]P[C≤dl,0″]+dl,0″P[dl,0″<C≤c1∗]−E[C|C≤dl,0′]P[C≤dl,0′]−dl,0′P[dl,0′<C≤c1∗]≥E[C|C≤dl,0′]P[C≤dl,0′]+dl,0′P[dl,0′<C≤dl,0″]P[C≤dl,0″]P[C≤dl,0″]+dl,0″P[dl,0″<C≤c1∗]−E[C|C≤dl,0′]P[C≤dl,0′]−dl,0′P[dl,0′<C≤c1∗]=dl,0′P[dl,0′<C≤dl,0″]+dl,0″P[dl,0″<C≤c1∗]−dl,0′P[dl,0′<C≤c1∗]=(dl,0″−dl,0′)P[dl,0″<C≤c1∗]≥0,
which proves the lemma. □
In addition, taking into account the continuity of FC gives that the right-hand side of (6) is continuous in dl,0, dl,1, …, dl,m.
For any fixedαlsatisfying Assumption1, we havedl,m≤min{αlE[C]/qm,cm∗},s0≤l≤s.
From Lemma 1, we conclude that for any fixed αl, the maximum possible value of dl,m is when dl,0=dl,1=⋯=dl,m−1=0. Therefore, dl,m≤αlE[C]/qm. Taking into account assertion (i) of Assumption 2 yields (7). □
Ifl1andl2are such thats0≤l1<l2≤s, thenαl1≤αl2.
By assertion (iii) of Assumption 2, we have dl1,0≤dl2,0, dl1,1≤dl2,1, …, dl1,m≤dl2,m. Consequently, from Lemma 1, it follows that
E[C|C≤dl1,0]P[C≤dl1,0]+dl1,0(q0−P[C≤dl1,0])+dl1,1q1+⋯+dl1,mqm≤E[C|C≤dl2,0]P[C≤dl2,0]+dl2,0(q0−P[C≤dl2,0])+dl2,1q1+⋯+dl2,mqm.
By (6), we have αl1≤αl2. □
For existence of a solution to (6), it is necessary thatαl≤min{1−1rl,f(c1∗,c2∗,…,cm∗)E[C]}for alls0≤l≤s,wheref(c1∗,c2∗,…,cm∗)=E[C|C≤c1∗]q0+c1∗q1+⋯+cm∗qm.
By Lemma 1, the maximum value of the right-hand side of (6) is attained when dl,0=c1∗, dl,1=c1∗, dl,2=c2∗, … , dl,m=cm∗ and equals to E[C|C≤c1∗]q0+c1∗q1+⋯+cm∗qm. Therefore, we get
αl≤E[C|C≤c1∗]q0+c1∗q1+⋯+cm∗qmE[C],s0≤l≤s.
Since
E[C|C≤c1∗]q0+c1∗q1+⋯+cm∗qm<E[C|C≤c1∗]q0+E[C|c1∗<C≤c2∗]q1+⋯+E[C|C>cm∗]qm=E[C],{c_{m}^{\ast }}\big]\hspace{0.1667em}q_{m}\\{} & \displaystyle \hspace{1em}=\mathbb{E}[C],\end{array}\]]]>
the right-hand side of (9) is less than 1. Thus, combining (9) and the inequality αl≤1−1/rl, which follows immediately from Assumption 1, gives (8). □
Let condition (8) hold. Then for any fixed l such thats0≤l≤s, we can always choosedl,0,dl,1, … ,dl,msatisfying assertions (i) and (ii) of Assumption2such that (6) is true.
The assertion of the lemma is evident.
Note only that if
min{1−1rl,f(c1∗,c2∗,…,cm∗)E[C]}=f(c1∗,c2∗,…,cm∗)E[C],
i.e.
αl=f(c1∗,c2∗,…,cm∗)E[C],
then the combination is unique: dl,0=c1∗, dl,1=c1∗, dl,2=c2∗, …, dl,m=cm∗. Otherwise, if
αl<f(c1∗,c2∗,…,cm∗)E[C],
then there are infinitely many such combinations.
Next, note that by Lemma 2, we can choose dl,0, dl,1, …, dl,m for any fixed l, but assertion (iii) of Assumption 2 may not hold.
The next theorem shows that we can always replace the bonus–malus system described in Section 2 by the bonus–malus system with varying deductibles.
There is always a solution to (6), i.e. a combination ofαl,dl,0,dl,1, … ,dl,m, wheres0≤l≤s, such that Assumptions1and2hold.
Consider now two cases.
1) If
min{1−1rs0,f(c1∗,c2∗,…,cm∗)E[C]}=f(c1∗,c2∗,…,cm∗)E[C],
let
αs0=f(c1∗,c2∗,…,cm∗)E[C].
By Lemma 2 and Remark 1 applied to l=s0, we get ds0,0=c1∗, ds0,1=c1∗, ds0,2=c2∗, … , ds0,m=cm∗. To satisfy Assumptions 1 and 2, for all l such that s0+1≤l≤s, we set αl=αs0, dl,0=ds0,0, dl,1=ds0,1, … , dl,m=ds0,m, which proves the theorem in the first case.
2) If
min{1−1rs0,f(c1∗,c2∗,…,cm∗)E[C]}=1−1rs0,
let
αs0=1−1rs0.
By Lemma 2 and Remark 1 applied to l=s0, we can always choose required ds0,0, ds0,1, …, ds0,m and there are infinitely many such combinations. We take any of them. Finally, to satisfy Assumptions 1 and 2, for all l such that s0+1≤l≤s, we also set αl=αs0, dl,0=ds0,0, dl,1=ds0,1, … , dl,m=ds0,m, which proves the theorem in the second case. □
The considerations given above show that the replacement of the bonus–malus system described in Section 2 (without deductibles) by the bonus–malus system with varying deductibles is not unique. So the insurance company can consider different replacements and choose one that seems to be more attractive from the policyholders’ point of view.
We now consider two special cases. Another one is considered in Section 4.
We now suppose that deductibles are applied only to claims of type m, i.e. dl,i=0 for all s0≤l≤s and 0≤i≤m−1. Therefore, (6) can be rewritten as
αlE[C]=dl,mqm.
Taking into account Theorem 1, we take any αs0 such that
0<αs0≤min{1−1rs0,cm∗qmE[C]}.
By (10), ds0,m=αs0E[C]/qm. To satisfy Assumptions 1 and 2, we set αl=αs0 and dl,m=ds0,m for all s0+1≤l≤s. Thus, we get a simple replacement to the bonus–malus system.
Let αl=1−1/rl for all s0≤l≤s. By Theorem 1, if
f(c1∗,c2∗,…,cm∗)E[C]<1−1rs0,
then (6) has no suitable solution. So the bonus–malus system cannot be replaced in this way.
The case where deductibles are applied only to the claims reported by policyholders occupying the highest level in the bonus–malus scale
We now consider the special case where dl,i=0 for all s0≤l≤s−1 and 0≤i≤m. This case is of great interest when the premium relativity for policyholders who are at level s is high enough and varying deductibles are applied to soften the bonus–malus system for such policyholders. To meet the conditions of Assumption 1, we require that (1−αs)rs≥rs−1, which gives αs≤1−rs−1/rs. So we can choose any positive αs such that
αs≤min{1−rs−1rs,f(c1∗,c2∗,…,cm∗)E[C]}
and then find ds,0, ds,1, …, ds,m from equation (6) applied to l=s, i.e. from
αsE[C]=E[C|C≤ds,0]P[C≤ds,0]+ds,0(q0−P[C≤ds,0])+ds,1q1+⋯+ds,mqm.
If inequality (11) is strict, then there are infinitely many solutions to (12). We now consider two allocation principles for the deductibles.
The first principle is when the deductibles are proportional to the average claims of each type. This principle seems to be natural and fair for policyholders, but that is not always possible, which is easily seen from the next theorem.
Letαsbe such that inequality (11) is strict and setx0=min{c1∗E[C|c1∗<C≤c2∗],c2∗E[C|c2∗<C≤c3∗],…,cm∗E[C|C>cm∗]}.{c_{m}^{\ast }}]}\hspace{-0.1667em}\bigg\}.\]]]>IfαsE[C]>E[C|C≤x0E[C|C≤c1∗]]P[C≤x0E[C|C≤c1∗]]+x0(E[C|C≤c1∗](q0−P[C≤x0E[C|C≤c1∗]])+E[C|c1∗<C≤c2∗]q1+⋯+E[C|C>cm∗]qm),\mathbb{E}\big[C\hspace{0.1667em}|\hspace{0.1667em}C\le x_{0}\hspace{0.1667em}\mathbb{E}\big[C\hspace{0.1667em}|\hspace{0.1667em}C\le {c_{1}^{\ast }}\big]\big]\hspace{0.2222em}\mathbb{P}\big[C\le x_{0}\hspace{0.1667em}\mathbb{E}\big[C\hspace{0.1667em}|\hspace{0.1667em}C\le {c_{1}^{\ast }}\big]\big]\\{} & \displaystyle \hspace{1em}+x_{0}\big(\mathbb{E}\big[C\hspace{0.1667em}|\hspace{0.1667em}C\le {c_{1}^{\ast }}\big]\hspace{0.1667em}\big(q_{0}-\mathbb{P}\big[C\le x_{0}\hspace{0.1667em}\mathbb{E}\big[C\hspace{0.1667em}|\hspace{0.1667em}C\le {c_{1}^{\ast }}\big]\big]\big)\\{} & \displaystyle \hspace{1em}+\mathbb{E}\big[C\hspace{0.1667em}|\hspace{0.1667em}{c_{1}^{\ast }}{c_{m}^{\ast }}\big]\hspace{0.1667em}q_{m}\big),\end{array}\]]]>then we cannot allocate deductibles proportionally to the average claims.
Otherwise, ifαsE[C]≤E[C|C≤x0E[C|C≤c1∗]]P[C≤x0E[C|C≤c1∗]]+x0(E[C|C≤c1∗](q0−P[C≤x0E[C|C≤c1∗]])+E[C|c1∗<C≤c2∗]q1+⋯+E[C|C>cm∗]qm),{c_{m}^{\ast }}\big]\hspace{0.1667em}q_{m}\big),\end{array}\]]]>then that is possible and (12) has a unique solution of such kind expressed byds,0=xE[C|C≤c1∗],ds,1=xE[C|c1∗<C≤c2∗],…,ds,m=xE[C|C>cm∗],{c_{m}^{\ast }}\big],\end{array}\]]]>where x is a unique positive solution to the equationαsE[C]=E[C|C≤xE[C|C≤c1∗]]P[C≤xE[C|C≤c1∗]]+x(E[C|C≤c1∗](q0−P[C≤xE[C|C≤c1∗]])+E[C|c1∗<C≤c2∗]q1+⋯+E[C|C>cm∗]qm).{c_{m}^{\ast }}\big]\hspace{0.1667em}q_{m}\big).\end{array}\]]]>
Let x be a proportionality coefficient, i.e. we have (16). To meet the assertion (i) of Assumption 2, we require that
xE[C|C≤c1∗]≤c1∗,xE[C|c1∗<C≤c2∗]≤c1∗,xE[C|c2∗<C≤c3∗]≤c2∗,…,xE[C|C>cm∗]≤cm∗,{c_{m}^{\ast }}\big]\le {c_{m}^{\ast }},\end{array}\]]]>
which is equivalent to
x≤min{c1∗E[C|c1∗<C≤c2∗],c2∗E[C|c2∗<C≤c3∗],…,cm∗E[C|C>cm∗]}.{c_{m}^{\ast }}]}\bigg\}.\]]]>
So we suppose that x∈[0,x0], where x0 is given by (13), and substitute (16) into (12), which yields (17).
The left-hand side of (17) is a positive constant. By Lemma 1, the right-hand side of (17) is an increasing function of x on [0,x0]. By the continuity of FC, it is continuous in x on [0,x0]. Moreover, this function is equal to 0 as x=0. Therefore, if (15) holds, then (17) has a unique solution x∈[0,x0]. Otherwise, if (14) is true, then there is no solution x∈[0,x0] to (17), which completes the proof. □
Equation (17) is not solvable analytically in the general case. To find x, we should use numerical methods.
The second principle implies that large claims are penalized by means of deductibles strictly and small claims are not penalized at all (if that is possible) or at least not so strictly. Anyway (12) must hold. Firstly, we check if it is possible to penalize only claims of type m. To this end, we set ds,0=ds,1=⋯=ds,m−1=0, ds,m=cm∗ and substitute this into (12). If
αsE[C]≤cm∗qm,
that is possible and the desired allocation is given by
ds,0=ds,1=⋯=ds,m−1=0andds,m=αsE[C]/qm.
Otherwise, if
αsE[C]>cm∗qm,{c_{m}^{\ast }}q_{m},\]]]>
we also have to penalize at least claims of type m−1. We set ds,0=ds,1=⋯=ds,m−2=0, ds,m−1=cm−1∗, ds,m=cm∗ and substitute this into (12). If
αsE[C]≤cm−1∗qm−1+cm∗qm,
then the desired allocation is given by
ds,0=ds,1=⋯=ds,m−2=0,ds,m−1=αsE[C]−cm∗qmqm−1andds,m=cm∗.
Otherwise, we also have to penalize at least claims of type m−2. So we set ds,0=ds,1=⋯=ds,m−3=0, ds,m−2=cm−2∗, ds,m−1=cm−1∗, ds,m=cm∗, substitute this into (12) and continue in this way until we get the desired allocation. Note that such an allocation is always possible provided that inequality (11) holds.
Numerical illustrations
We now consider the scale with 4 levels (numbered from 0 to 3) and 4 claim types (numbered from 0 to 3 with probabilities q0, q1, q2 and q3). If no claims are reported during the current year, the policyholder moves one level down. Each claim of type 0 is penalized by one level, each claim of type 1 is penalized by 2 levels, each claim of type 2 or 3 is penalized by 3 levels. Note that penalties for claims of types 2 and 3 are different due to varying deductibles. For instance, if 2 claims of type 0 are reported during the year, then the policyholder moves 2 levels up; if 1 claim of type 0 and 1 claim of type 1 are reported, then the policyholder moves 3 levels up, i.e. goes at the highest level anyway.
The elements of the one-step transition matrix for a policyholder with annual mean claim frequency λθ and vector of probabilities q=(q0,q1,q2,q3)T are calculated using formula (1). Thus, the one-step transition matrix is given by
P(λθ;q)=e−λθλθq0e−λθλθq1e−λθ+(λθq0)22e−λθ1−e−λθ−λθ(q0+q1)e−λθ−(λθq0)22e−λθe−λθ0λθq0e−λθ1−e−λθ−λθq0e−λθ0e−λθ01−e−λθ00e−λθ1−e−λθ.
Next, the stationary probabilities πl(λθ;q), 0≤l≤3, are calculated using formula (2). A standard computation shows that
π0(λθ;q)=e−3λθΔ,π1(λθ;q)=e−2λθ−e−3λθΔ,π2(λθ;q)=e−λθ−e−2λθ−λθq0e−3λθΔ,π3(λθ;q)=1−e−λθ−2λθq0e−2λθ+λθ(q0−q1)e−3λθ−(λθq0)22e−3λθΔ,
where
Δ=1−2λθq0e−2λθ−λθq1e−3λθ−(λθq0)22e−3λθ.
It is easily seen that
e−3λθ≥0,e−2λθ−e−3λθ≥0ande−λθ−e−2λθ−λθq0e−3λθ≥0
for all λ>00$]]>, θ≥0 and 0<q0<1. Moreover, it is evident that ∑l=03πl(λθ;q)=1. So to see that πl(λθ;q), 0≤l≤3, are indeed probabilities, we must show that
1−e−λθ−2λθq0e−2λθ+λθ(q0−q1)e−3λθ−(λθq0)22e−3λθ≥0
for all λ>00$]]>, θ≥0, q0>00$]]> and q1>00$]]> such that q0+q1<1.
It is easy to check that the minimum value of the left-hand side of (18) with respect to q0≥0 and q1≥0 such that q0+q1≤1 is attained as q0=1, q1=0 and equal to
1−e−λθ−2λθe−2λθ+λθe−3λθ−(λθ)22e−3λθ.
Introduce the function
h(y)=1−e−y−2ye−2y+ye−3y−y22e−3y,y≥0.
Taking the derivative yields
h′(y)=e−3y(ey−1)2+2y(2e−2y−e−3y)+3y22e−3y≥0,y≥0.
Therefore, h(y) is non-decreasing and its minimum value is attained as y=0 and equals to 0. Consequently, the minimum value of (19) is also 0 as λθ=0, which gives (18).
In this section, we deal with exponentially distributed claim sizes with mean μ>00$]]>, i.e. the distribution function FC of C is equal to FC(y)=1−e−y/μ, y≥0. Thus, we have
q0=1−e−c1∗/μ,q1=e−c1∗/μ−e−c2∗/μ,q2=e−c2∗/μ−e−c3∗/μ,q3=1−e−c3∗/μ;E[C|C≤c1∗]=1q0∫0c1∗yμe−y/μdy=μ(1−e−c1∗/μ)−c1∗e−c1∗/μ1−e−c1∗/μ=μ−c1∗e−c1∗/μ1−e−c1∗/μ,E[C|c1∗<C≤c2∗]=1q1∫c1∗c2∗yμe−y/μdy=(c1∗+μ)e−c1∗/μ−(c2∗+μ)e−c2∗/μe−c1∗/μ−e−c2∗/μ=μ+c1∗e−c1∗/μ−c2∗e−c2∗/μe−c1∗/μ−e−c2∗/μ,E[C|c2∗<C≤c3∗]=1q2∫c2∗c3∗yμe−y/μdy=(c2∗+μ)e−c2∗/μ−(c3∗+μ)e−c3∗/μe−c2∗/μ−e−c3∗/μ=μ+c2∗e−c2∗/μ−c3∗e−c3∗/μe−c2∗/μ−e−c3∗/μ,E[C|C>c3∗]=1q3∫c3∗+∞yμe−y/μdy=(c3∗+μ)e−c3∗/μe−c3∗/μ=μ+c3∗;{c_{3}^{\ast }}\big]=\frac{1}{q_{3}}{\int _{{c_{3}^{\ast }}}^{+\infty }}\frac{y}{\mu }\hspace{0.1667em}{e}^{-y/\mu }\hspace{0.1667em}\mathrm{d}y=\frac{({c_{3}^{\ast }}+\mu ){e}^{-{c_{3}^{\ast }}/\mu }}{{e}^{-{c_{3}^{\ast }}/\mu }}=\mu +{c_{3}^{\ast }};\]]]>f(c1∗,c2∗,c3∗)=μ−μe−c1∗/μ−c1∗e−c1∗/μ+c1∗(e−c1∗/μ−e−c2∗/μ)+c2∗(e−c2∗/μ−e−c3∗/μ)+c3∗e−c3∗/μ=μ−μe−c1∗/μ+(c2∗−c1∗)e−c2∗/μ+(c3∗−c2∗)e−c3∗/μ.
In addition, we suppose that FΘ(θ)=1−e−θ, θ≥0.
Let λ=0.1, μ=2, c1∗=1, c2∗=2, c3∗=4. Then we have q0≈0.3935, q1≈0.2387, q2≈0.2325, q3≈0.1353. The corresponding values of πl and rl, 0≤l≤3, are calculated using formulas (3) and (4), respectively, and are given in Tables 1–9.
Bonus–malus system with varying deductibles for Example 3 where deductibles are applied only to claims reported by policyholders occupying the highest level, α3=0.05 and the first principle is used
l
πl
rl
λrlE[C]
αl
(1−αl)λrlE[C]
dl,0
dl,1
dl,2
dl,3
3
0.0508
2.1844
0.4369
0.05
0.4150
0.0230
0.0730
0.1420
0.3004
2
0.0591
1.8899
0.3780
0
0.3780
0
0
0
0
1
0.0716
1.6543
0.3309
0
0.3309
0
0
0
0
0
0.8185
0.8050
0.1610
0
0.1610
0
0
0
0
Bonus–malus system with varying deductibles for Example 3 where deductibles are applied only to claims reported by policyholders occupying the highest level, α3=0.13 and the first principle is used
l
πl
rl
λrlE[C]
αl
(1−αl)λrlE[C]
dl,0
dl,1
dl,2
dl,3
3
0.0508
2.1844
0.4369
0.13
0.3801
0.0598
0.1903
0.3699
0.7827
2
0.0591
1.8899
0.3780
0
0.3780
0
0
0
0
1
0.0716
1.6543
0.3309
0
0.3309
0
0
0
0
0
0.8185
0.8050
0.1610
0
0.1610
0
0
0
0
We first consider the case when deductibles are applied only to the claims reported by policyholders occupying the highest level in the bonus–malus system. By (11), we have
α3≤min{1−r2r3,f(c1∗,c2∗,c3∗)μ}≈{1−1.88992.1844,1.4254892}≈min{0.1348,0.7127}=0.1348.
By (13), we get x0=0.6667. Taking α3=0.05 and α3=0.13 shows that (15) is true in both cases. Hence, we can allocate deductibles proportionally to the average claims. By (17), the proportionality coefficient x=0.050066 if α3=0.05 and x=0.130443 if α3=0.13. The corresponding values of deductibles are calculated using (16) and given in Tables 1 and 2.
If we apply the second principle considered in Section 4, we get α3μ≤c3∗q3 for both values of α3. The desired allocation is presented in Tables 3 and 4.
Bonus–malus system with varying deductibles for Example 3 where deductibles are applied only to claims reported by policyholders occupying the highest level, α3=0.05 and the second principle is used
l
πl
rl
λrlE[C]
αl
(1−αl)λrlE[C]
dl,0
dl,1
dl,2
dl,3
3
0.0508
2.1844
0.4369
0.05
0.4150
0
0
0
0.7389
2
0.0591
1.8899
0.3780
0
0.3780
0
0
0
0
1
0.0716
1.6543
0.3309
0
0.3309
0
0
0
0
0
0.8185
0.8050
0.1610
0
0.1610
0
0
0
0
Bonus–malus system with varying deductibles for Example 3 where deductibles are applied only to claims reported by policyholders occupying the highest level, α3=0.13 and the second principle is used
l
πl
rl
λrlE[C]
αl
(1−αl)λrlE[C]
dl,0
dl,1
dl,2
dl,3
3
0.0508
2.1844
0.4369
0.05
0.3801
0
0
0
1.9212
2
0.0591
1.8899
0.3780
0
0.3780
0
0
0
0
1
0.0716
1.6543
0.3309
0
0.3309
0
0
0
0
0
0.8185
0.8050
0.1610
0
0.1610
0
0
0
0
Bonus–malus system with varying deductibles for Example 3 where deductibles are applied only to claims of type 3
l
πl
rl
λrlE[C]
αl
(1−αl)λrlE[C]
dl,0
dl,1
dl,2
dl,3
3
0.0508
2.1844
0.4369
0.24
0.3320
0
0
0
3.5467
2
0.0591
1.8899
0.3780
0.13
0.3288
0
0
0
1.9212
1
0.0716
1.6543
0.3309
0.06
0.3110
0
0
0
0.8867
0
0.8185
0.8050
0.1610
0
0.1610
0
0
0
0
Bonus–malus system with varying deductibles for Example 3 where deductibles are applied only to claims of type 3. Larger values of αl
l
πl
rl
λrlE[C]
αl
(1−αl)λrlE[C]
dl,0
dl,1
dl,2
dl,3
3
0.0508
2.1844
0.4369
0.26
0.3233
0
0
0
3.8423
2
0.0591
1.8899
0.3780
0.25
0.2835
0
0
0
3.6945
1
0.0716
1.6543
0.3309
0.24
0.2514
0
0
0
3.5467
0
0.8185
0.8050
0.1610
0
0.1610
0
0
0
0
Tables 5 and 6 give examples of bonus–malus system in the case where deductibles are applied only to claims of type 3. In Table 6, we take larger values of αl, so we get higher values of the deductibles.
If we also apply deductibles to claims of type 2 for the same values of αl, the values of deductibles dl,3 are not so high (see Table 7). Moreover, in this case, we can take larger values of αl such that Assumptions 1 and 2 hold (see Table 8).
Table 9 presents an example of bonus–malus system in the case where deductibles are applied to claims of types 1, 2 and 3. On the one hand, policyholders who are in the malus zone pay not such high premiums. On the other hand, the corresponding values of deductibles are moderate. This bonus–malus system seems to be more attractive.
Bonus–malus system with varying deductibles for Example 3 where deductibles are applied to claims of types 2 and 3
l
πl
rl
λrlE[C]
αl
(1−αl)λrlE[C]
dl,0
dl,1
dl,2
dl,3
3
0.0508
2.1844
0.4369
0.26
0.3233
0
0
1.1
1.9522
2
0.0591
1.8899
0.3780
0.25
0.2835
0
0
1.1
1.8044
1
0.0716
1.6543
0.3309
0.24
0.2514
0
0
1.1
1.6566
0
0.8185
0.8050
0.1610
0
0.1610
0
0
0
0
Bonus–malus system with varying deductibles for Example 3 where deductibles are applied to claims of types 2 and 3. Larger values of αl
l
πl
rl
λrlE[C]
αl
(1−αl)λrlE[C]
dl,0
dl,1
dl,2
dl,3
3
0.0508
2.1844
0.4369
0.45
0.2403
0
0
1.7
3.7291
2
0.0591
1.8899
0.3780
0.40
0.2268
0
0
1.6
3.1620
1
0.0716
1.6543
0.3309
0.35
0.2151
0
0
1.5
2.5949
0
0.8185
0.8050
0.1610
0
0.1610
0
0
0
0
Bonus–malus system with varying deductibles for Example 3 where deductibles are applied to claims of types 1, 2 and 3
l
πl
rl
λrlE[C]
αl
(1−αl)λrlE[C]
dl,0
dl,1
dl,2
dl,3
3
0.0508
2.1844
0.4369
0.45
0.2403
0
0.7
1.5
2.8383
2
0.0591
1.8899
0.3780
0.40
0.2268
0
0.5
1.4
2.6239
1
0.0716
1.6543
0.3309
0.35
0.2151
0
0.3
1.3
2.4096
0
0.8185
0.8050
0.1610
0
0.1610
0
0
0
0
Let λ=0.1, μ=2, c1∗=0.3, c2∗=1.2, c3∗=2.8. Then we have q0≈0.1393, q1≈0.3119, q2≈0.3022, q3≈0.2466. The corresponding values of πl and rl as well as examples of bonus–malus systems with varying deductibles are given in Tables 10–12.
Bonus–malus system with varying deductibles for Example 4 where deductibles are applied to claims of types 2 and 3
l
πl
rl
λrlE[C]
αl
(1−αl)λrlE[C]
dl,0
dl,1
dl,2
dl,3
3
0.0653
2.0731
0.4146
0.20
0.3317
0
0
0.30
1.2544
2
0.0717
1.7925
0.3585
0.15
0.3047
0
0
0.25
0.9102
1
0.0679
1.6263
0.3253
0.10
0.2927
0
0
0.20
0.5659
0
0.7951
0.7869
0.1574
0
0.1574
0
0
0
0
Bonus–malus system with varying deductibles for Example 4 where deductibles are applied to claims of types 1, 2 and 3
l
πl
rl
λrlE[C]
αl
(1−αl)λrlE[C]
dl,0
dl,1
dl,2
dl,3
3
0.0653
2.0731
0.4146
0.24
0.3151
0
0.10
0.60
1.0847
2
0.0717
1.7925
0.3585
0.22
0.2796
0
0.10
0.55
0.9838
1
0.0679
1.6263
0.3253
0.20
0.2602
0
0.05
0.50
0.9461
0
0.7951
0.7869
0.1574
0
0.1574
0
0
0
0
Bonus–malus system with varying deductibles for Example 4 where deductibles are applied to claims of types 1, 2 and 3. Larger values of αl
l
πl
rl
λrlE[C]
αl
(1−αl)λrlE[C]
dl,0
dl,1
dl,2
dl,3
3
0.0653
2.0731
0.4146
0.45
0.2280
0
0.20
0.9
2.2937
2
0.0717
1.7925
0.3585
0.40
0.2151
0
0.15
0.8
2.0740
1
0.0679
1.6263
0.3253
0.35
0.2114
0
0.10
0.7
1.8543
0
0.7951
0.7869
0.1574
0
0.1574
0
0
0
0
Conclusion
The bonus–malus systems with different claim types and varying deductibles eliminate both drawbacks of the traditional bonus–malus systems mentioned in Section 1 and present a number of advantages, namely:
Policyholders reporting small and large claims are not penalized in the same way. This helps to avoid or at least decrease bonus hunger.
Policyholders will do all their best to prevent or at least decrease the losses.
Even if a policyholder leaves the company after a claim, he has to pay for the deductible.
Relative premiums and amounts of deductibles may be tuned in an optimal way in order to attract policyholders.
Section 5 gives a few examples of such bonus–malus systems. The numerical illustrations show that use of both penalty types (premium relativities and varying deductibles) in this way seems indeed attractive and fair for policyholders. On the one hand, policyholders who are in the malus zone pay not such high premiums. On the other hand, the corresponding values of deductibles are moderate. In addition, it is fair that only policyholders reporting (large) claims are subject to further penalties, i.e. deductibles.
Acknowledgments
The author is deeply grateful to the anonymous referees for careful reading and valuable comments and suggestions, which helped to improve the earlier version of the paper.
ReferencesBonsdorff, H.: On asymptotic properties of bonus–malus systems based on the number and on the size of the claims. 2005, 309–320 (2005). MR2164049. doi:10.1080/03461230510009826Denuit, M., Dhaene, J.: Bonus–malus scales using exponential loss functions. 25, 13–27 (2001)Denuit, M., Maréchal, X., Pitrebois, S., Walhin, J.-F.: . John Wiley & Sons, Chichester (2007). MR2384837. doi:10.1002/9780470517420Dionne, G., Vanasse, C.: A generalization of actuarial automobile insurance rating models: the negative binomial distribution with a regression component. 19, 199–212 (1989)Dionne, G., Vanasse, C.: Automobile insurance ratemaking in the presence of asymmetrical information. 7, 149–165 (1992)Frangos, N.E., Vrontos, S.D.: Design of optimal bonus–malus systems with a frequency and a severity component on an individual basis in automobile insurance. 31, 1–22 (2001). MR1945629. doi:10.2143/AST.31.1.991Gómez-Déniz, E.: Bivariate credibility bonus–malus premiums distinguishing between two types of claims. 70, 117–124 (2016). MR3543037. doi:10.1016/j.insmatheco.2016.06.009Gómez-Déniz, E., Hernández-Bastida, A., Fernández-Sánchez, M.P.: Computing credibility bonus–malus premiums using the total claim amount distribution. 43, 1047–1061 (2014). MR3331161Holtan, J.: Bonus made easy. 24, 61–74 (1994)Lemaire, J.: . Kluwer Academic Publisher, Boston (1995)Lemaire, J., Zi, H.: High deductibles instead of bonus–malus. Can it work?24, 75–88 (1994)Mahmoudvand, R., Hassani, H.: Generalized bonus–malus systems with a frequency and a severity component on an individual basis in automobile insurance. 39, 307–315 (2009). MR2749888. doi:10.2143/AST.39.1.2038066Mert, M., Saykan, Y.: On a bonus–malus system where the claim frequency distribution is geometric and the claim severity distribution is Pareto. 34, 75–81 (2005). MR2212712Norberg, R.: A credibility theory for automobile bonus system. 1976, 92–107 (1976). MR0428666. doi:10.1080/03461238.1976.10405605Pinquet, J.: Allowance for cost of claims in bonus–malus systems. 27, 33–57 (1997)Pinquet, J.: Designing optimal bonus–malus systems from different types of claims. 28, 205–220 (1998)Pitrebois, S., Denuit, M., Walhin, J.-F.: Bonus–malus systems with varying deductibles. 35, 261–274 (2005). MR2143217. doi:10.2143/AST.35.1.583175Pitrebois, S., Denuit, M., Walhin, J.-F.: Multi-event bonus–malus scales. 73, 517–528 (2006)Rolski, T., Schmidli, H., Schmidt, V., Teugels, J.: . John Wiley & Sons, Chichester (1999). MR1680267. doi:10.1002/9780470317044Tzougas, G., Vrontos, S., Frangos, N.: Optimal bonus–malus systems using finite mixture models. 44, 417–444 (2014). MR3389581. doi:10.1017/asb.2013.31