Next Article in Journal
A Moiré Removal Method Based on Peak Filtering and Image Enhancement
Next Article in Special Issue
Dominations in Intutionistic Fuzzy Directed Graphs with Applications towards Influential Graphs
Previous Article in Journal
Loss-Averse Supply Chain Coordination with a Revenue-Sharing Contract
Previous Article in Special Issue
Hyperconnectedness and Resolvability of Soft Ideal Topological Spaces
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Calculating Insurance Claim Reserves with an Intuitionistic Fuzzy Chain-Ladder Method

by
Jorge De Andrés-Sánchez
Social and Business Research Laboratory, Universitat Rovira i Virgili, Campus de Bellissens, 43204 Reus, Spain
Mathematics 2024, 12(6), 845; https://doi.org/10.3390/math12060845
Submission received: 21 February 2024 / Revised: 10 March 2024 / Accepted: 11 March 2024 / Published: 13 March 2024
(This article belongs to the Special Issue Advances and Applications of Soft Computing)

Abstract

:
Estimating loss reserves is a crucial activity for non-life insurance companies. It involves adjusting the expected evolution of claims over different periods of active policies and their fluctuations. The chain-ladder (CL) technique is recognized as one of the most effective methods for calculating claim reserves in this context. It has become a benchmark within the insurance sector for predicting loss reserves and has been adapted to estimate variability margins. This variability has been addressed through both stochastic and possibilistic analyses. This study adopts the latter approach, proposing the use of the CL framework combined with intuitionistic fuzzy numbers (IFNs). While modeling with fuzzy numbers (FNs) introduces only epistemic uncertainty, employing IFNs allows for the representation of bipolar data regarding the feasible and infeasible values of loss reserves. In short, this paper presents an extension of the chain-ladder technique that estimates the parameters governing claim development through intuitionistic fuzzy regression, such as symmetric triangular IFNs. Additionally, it compares the results obtained with this method with those derived from the stochastic chain ladder by England and Verrall.

1. Introduction

The estimation of loss reserves is a fundamental process in the management of insurance companies. It consists of setting a prudent value on claims not yet made on active policies, which will ultimately impact the financial statements and the required capital to continue with the current insurance portfolio [1]. Thus, a prudent estimation of these provisions, which ultimately requires the application of so-called actuarial judgement, needs to use a value of maximum reliability but, at the same time, estimate the possible variability around that expected value [2]. The final estimated value for reserves, although it should tend to overestimate them and cover possible unfavorable deviations from their expected value, should not be excessive [3].
Within claim-reserving methods, the actuarial literature often distinguishes between deterministic and stochastic methods. While the former provides a point value of reserves that can be considered the “expected” or maximum confidence value, stochastic methods allow the variability around that reasonable value to be measured [2]. To this commonly accepted typology, we can add fuzzy methods [4].
Among the various applications that fuzzy set theory (FST) has had in insurance mathematics, we can outline the modeling of uncertain and vague parameters with possibility distributions [5,6]. In these applications, fuzzy modeling allows the quantification of epistemic uncertainty, that is, a measure of the reliability with which a certain variable A takes a specific value x [7]. In the context of determining loss reserves, this vagueness may be induced, first, by the imprecision of some of the data available to the insurance company [8]. An additional source to consider is the scarcity of the sample used for reserve calculation; since it is not advisable to use data too far from the present, it can bias estimates due to factors such as changes in judicial practices and public awareness of liability issues [9].
The literature on the variability of loss reserves starts from a scheme used in practice that allows loss reserves to be obtained as a point value. A very common scheme for calculating the value of reliable mathematical reserves is the chain-ladder (CL) method, or variants of this method, such as the London CL or the Bornhuetter–Ferguson methods [10]. The chain-ladder method has been the subject of adaptations that allow modeling of the variability of reserves stochastically [11,12] but also with possibility distributions [13,14,15].
As shown in Table 1, among the most commonly used schemes for quantifying loss reserves, in addition to the CL method, we can highlight the geometric separation method [16] and methods that model incremental claims in a two-way manner, such as in [17]. The methodology for adjusting the parameters governing the evolution of claims over time can be performed heuristically [4,14,15,18] or with fuzzy regression methods that apply both the principle of minimum fuzziness [19,20] and the fuzzy least-squares approach [21,22].
All of the methods reviewed in Table 1 model the uncertainty of parameters with type-one fuzzy numbers, i.e., simple fuzzy numbers (FNs), that is, through possibility distributions that allow introducing epistemic uncertainty about the real value of the parameters [27]. However, FNs do not allow the introduction of negative information about these parameters that the evaluator might have, i.e., about what the parameters “are not”. This paper extends fuzzy loss reserving to the use of bipolar information, i.e., imprecise estimations about the values that the parameters of interest can take and about those they cannot take. Bipolarity does not introduce additional uncertainty but provides new information [28].
Our paper uses the chain-ladder scheme to capture the dynamics of claiming processes and the concept of intuitionistic fuzzy numbers (IFNs) that model uncertain quantities [29] within Atanassov’s theory of intuitionistic fuzzy sets [30,31]. Thus, this work expands the practical applications of IFNs, which are relatively scarce in finance and insurance. Among such applications, we can highlight the following:
  • Capital budgeting [32,33,34,35];
  • Option pricing [36,37,38];
  • Productivity measurements [39,40,41];
  • Actuarial field: while Uzhga-Rebrov and Grabusts [42] use intuitionistic fuzzy values to address environmental risk analysis, Andrés-Sánchez [43] does so to price the life contingencies of people with impaired life expectancies.
This paper falls within the fourth domain, specifically in the field of claim reserving. In this regard, our aim is threefold. First, we demonstrate that the estimation of stochastic loss reserves can be interpreted as estimates made through possibility distributions with type-one fuzzy numbers. Subsequently, we introduce intuitionistic fuzzy regression in claim-reserving calculations. Although fuzzy regression has been applied in several areas of actuarial analysis, such as mortality adjustment [44,45], the use of intuitionistic fuzzy regression in actuarial science is nonexistent. We do so by employing the intuitionistic fuzzy regression method [46], which extends the possibilistic regression models [47,48,49]. Similarly, we compare the results of the proposed method to those obtained with the stochastic chain-ladder (SCL) method [12].

2. Intuitionistic Fuzzy Numbers

2.1. Fuzzy Numbers and Intuitionistic Fuzzy Numbers

Definition 1.
A fuzzy set (FS) in a referential set  X A ¨ , is defined as follows [50]:
A ¨ = x , μ A x ,   x X ,
where  μ A : X 0 , 1  is the membership function of  A ¨ .
Definition 2.
The fuzzy set  A ¨  can be represented through level sets or  α -cuts,  A α [50]:
A α = x | μ A x α ,   0 < α 1 .
Definition 3.
A fuzzy number (FN),  A ¨ , is a fuzzy subset of the real line [51] such that
i. 
is normal, i.e., x | μ A x = 1 ;
ii. 
is convex, i.e., x 1 , x 2 R , 0   λ 1 ,   μ A λ x 1 , 1 λ x 2 min μ A x 1 , μ A x 2 .
Remark 1.
As a consequence, the α-cuts of  A ¨  and  A α  are confidence intervals:
A α = x | μ A x α ,   0 < α 1 = A _ α ,   A ¯ α ,
where  A _ α  is an increasing function of  α  and  A ¯ α is a decreasing function.
Remark 2.
The membership function of  A ¨ ,   μ A x is also called the possibility distribution function.
Fuzzy set theory commonly relies on fuzzy numbers (FNs) to represent imprecise quantities [51]. Specifically, triangular fuzzy numbers are very common in practical applications because the grading of the membership level is linear. This approach is reasonable because it applies the principle of parsimony when dealing with vague information [52].
Definition 4.
A symmetric triangular fuzzy number (STFN) is a particular case in which a triangular fuzzy number (TFN) can be represented by the couple  A ¨ = A , r A ,  r A 0 . Then, the membership function is
μ A x = 1 x A r A x A < r A 1 x = A 0 o t h e r w i s e ,
with the following being its α-cut representation:
A α = A _ α ,   A ¯ α = A r A 1 α , A + r A 1 α ,   0 α 1 .
Within TFNs, shapes are of special interest when the available information about the reference variable is scarce and can be summarized in a center and plausible deviations from it [53,54]. Symmetric triangular fuzzy numbers (STFNs) allow for a good balance between comprehensiveness in capturing the available information and the use of the parsimony principle [54]. In insurance modeling, the usefulness of STFNs has been shown in several papers [13,14,55].
Definition 5.
Let us take a continuous random variable  A  and a family of confidence intervals  A α , such that  P x A α 1 α  and P(·) is a probability measure. Therefore, an equivalent fuzzy quantity  A ¨  has the following α-cut,  A α  [56]:
A α = A _ α ,   A ¯ α = x | F A x = α 2 ,   x | F A x = 1 α 2 ,
where F(·) is the distribution function.
Remark 3.
Consequently, the possibility distribution function of  A ¨  equivalent to  A  is  μ A x = s u p α | x ϵ A α .
It should be emphasized that the interpretation of probabilistic confidence intervals as α-level sets of possibility distributions has been widely argued in the literature [53,54,56,57,58,59,60,61]. Buckley [58] justified the transformation of a set of probabilistic confidence intervals into fuzzy numbers with the fact that in subsequent calculations, more information is used than simple point estimates or confidence intervals.
Definition 6.
The intuitionistic fuzzy set (IFS)  A ~  defined in a referential set  X  is
A ~ = x , μ A x , v A x ,   x X ,
where  μ A : X 0 , 1  measures the membership of  x  in  A ~  and  v A : X 0 , 1  is nonmembership. These functions must accomplish  0 μ A x + v A x 1 .
Remark 4.
The degree of hesitancy,  h A x , of  A ~  is  h A x = 1 μ A x v A x .
Remark 5.
An IFS generalizes the concept of an FS such that if  h A x = 0   x ,  A ~  is a conventional FS  A ¨ .
Definition 7.
An IFN can be expressed using  α , β -levels or  α , β -cuts, as  A α , β :
A α , β = x | μ A x α ,   v A x β , 0 α + β 1 , α , β 0 , 1 .
Remark 6.
A α , β  can be decoupled into two level sets [62], such as  A α = x | μ A x α   a n d   A β * = x | v A x β ,  in such a way that
A α , β = A α = x | μ A x α   ,   A β * = x | v A x β , 0 α + β 1 , α , β 0 , 1 .
Definition 8.
An intuitionistic fuzzy number (IFN) is an IFS defined on real numbers, such that
i. 
is normal, i.e.,  x | μ A x = 1 v A x = h A x = 0 ;
ii. 
μ A x  is convex,  x 1 , x 2 R , 0 λ 1 , μ A λ x 1 , 1 λ x 2 min μ A x 1 , μ A x 2 ;
iii. 
v A x  is concave,  x 1 , x 2 R , 0 λ 1 , v A λ x 1 , 1 λ x 2 max v A x 1 , v A x 2 .
Remark 7.
The  α , β -cuts of  A ~  and  A α , β  can be decoupled as follows:  A α = x | μ A x α = A _ α ,   A ¯ α   a n d   A β * = v A x β = A * _ β , A * ¯ β ,   0 α + β 1 ,   α , β 0 , 1 .
Remark 8.
Thus, from Remark 7, an  α , β -level of  A α , β  can be represented as
A α , β = A α = A _ α , A ¯ α , A * β = A * _ β , A * ¯ β , 0 α + β 1 , α , β ( 0 , 1 ) ,
where  A _ α  and  A * ¯ β  increase with their arguments,  α  and  β , respectively. Similarly,  A ¯ α  and  A * _ β  are decreasing with respect to these arguments.
Remark 9.
In an IFN,  μ A x  can be interpreted as the lower possibility distribution function of the quantity of interest A, and  μ A * x = 1 v A x  is the upper distribution function of that quantity.
The functions μ A * x and μ A x can be interpreted as bipolar possibility distribution measurements, in such a way that μ A * x accounts for the potential possibility and μ A x quantifies the real possibility of A being x [28].
Definition 9.
A symmetric triangular intuitionistic fuzzy number (STIFN) is a particular case of a triangular intuitionistic fuzzy number (TIFN) that can be denoted as  A ~ = A , r A , r A *   , with membership and nonmembership functions:
μ A x = 1 x A r A x A < r A 1 x = A 0 o t h e r w i s e ,
and
υ A x = x A r A * x A < r A 0 x = A 1 o t h e r w i s e ,
where  r A r A * . Figure 1 depicts the shape of an STIFN and the relationship between the embedded functions  μ A x  (the actual possibility distribution function),  v A x ,  μ A * x  (the potential possibility distribution function), and  h A x .
Remark 10.
The level sets  A α , β    of a TIFN can be decoupled into
A α = A _ α ,   A ¯ α = A r A 1 α , A + r A 1 α ,
A * β = A * _ β , A * ¯ β = A r A * β , A + r A * β .
Thus, STIFNs are an extension of STFNs such that if r A = r A * , we deal with conventional TFNs [63]. Thus, the use of symmetrical TFNs based on the principle of parsimony to justify their use can be extended to the use of STIFNs.

2.2. Intuitionistic Fuzzy Number Arithmetic

The fuzzy loss-reserving methods in Table 1 calculate provisions based on the assumption that the parameters governing the claiming process are determined by fuzzy numbers. Performing arithmetic operations with FNs requires the application of Zadeh’s extension principle, which can be implemented through α-cuts [64].
Zadeh’s extension principle and its compatibility with α -cuts arithmetic can be extended to the evaluation of functions defined in real numbers when the parameters are IFNs instead of FNs [65]. This paper considers the case of continuous and differentiable functions y = f x 1 , x 2 , , x n , such that the values of the input variables are given as the means of IFNs A ~ ( i ) ,   i = 1 , 2 , , n . This generates an IFN B ~ , B ~ = f A ~ ( 1 ) , A ~ ( 2 ) , , A ~ ( n ) . Thus, the membership and nonmembership functions of B ~ are as follows:
μ B y = max y = f x 1 , x 2 , , x n min μ A ( 1 ) x 1 , μ A ( 2 ) x 2 , , μ A ( n ) x n ,
v B y = min y = f x 1 , x 2 , , x n max v A ( 1 ) x 1 , v A ( 2 ) x 2 , , v A ( n ) x n .
Therefore, if A ~ ( i ) ,   i = 1 , 2 , , n are FNs, it is only necessary to obtain μ B y using the usual max/min principle. However, we can fit B ~ thoughout B α , β from A ( i ) α , β by calculating B α , β = f A ( 1 ) α , β , A ( 2 ) α , β , , A ( n ) α , β . Thus, given that f is continuous, the α , β -cuts of B ~ are defined as B α , β = B α = B _ α , B ¯ α , B * β = B * _ β , B * ¯ β , 0 α + β 1 , α , β 0 , 1 , where
B _ α = inf y | y = f x 1 , , x n , x i A i α , B ¯ α = sup y | y = f x 1 , , x n , x i A i α ,
B * _ β = inf y | y = f x 1 , , x n , x i A ( i ) * β , B * ¯ β = sup y | y = f x 1 , , x n , x i A ( i ) * β .
Following [66], when f monotonically increases with respect to x i ,   i = 1 , 2 , m and monotonically decreases in x i ,   i = m + 1 , m + 2 , , n ,   m n , B α = B _ α , B ¯ α is as follows:
B _ α = f A 1 _ α , A 2 _ α , , A m _ α , A m + 1 ¯ α , , A n ¯ α   and B ¯ α = f A 1 ¯ α , A 2 ¯ α , , A m ¯ α , A m + 1 _ α , , A n _ α .
By analogy, the β-cut representation of B * β = B * _ β , B * ¯ β is
B * _ β = f A ( 1 ) * _ β , A ( 2 ) * _ β , , A m * _ β , A m + 1 * ¯ β , , A ( n ) * ¯ β   and B * ¯ β = f A ( 1 ) * ¯ β , A ( 2 ) * ¯ β , , A ( m ) * β ¯ , A ( m + 1 ) * _ β , , A ( n ) * _ β .
The linear combination of STIFNs is also an STIFN. Therefore, from the STIFNs A ~ ( i ) = A ( i ) , r A ( i ) , r A ( i ) *   , B ~ = B ,   r B ,   r B *   , where [46]
B = i = 1 n λ i A ( i ) ,     r B = i = 1 n λ i · r A ( i ) ,     r B * = i = 1 n λ i · r A ( i ) * .
The evaluation of nonlinear functions using STIFNs does not produce a new STIFN. Despite this limitation, we feel that maintaining a linear shape is relevant. Following the argument in [67] justifying the use of approximating linear fuzzy numbers, complicated forms of IFNs may cause drawbacks in processing imprecise information modeled by these fuzzy structures, and the interpretation of the results becomes more difficult. The same argument, based on the parsimony principle, can be used to maintain the symmetrical structure of the input data.
Thus, we evaluate the approximation to B ~ = f A ~ ( 1 ) , A ~ ( 2 ) , , A ~ ( n ) with an STIFN B T ~ = B , r B , r B * when the inputs are A ~ ( i ) = A ( i ) , r A ( i ) , r A ( i ) * , i = 1 , 2 , , n . To do this, we rely on the results of [51] to approximate an LR fuzzy number to the result of a nonlinear function of LR fuzzy numbers, which is based on the linear approximation of α-cuts with a Taylor expansion. In the field of FNs, this methodology produces an STFN that approximates the functions of STFNs, as shown in several actuarial applications [13,14,15,55].
The extremes of the α -cuts B α = B _ α , B ¯ α in (19) are approximated by means of a Taylor expansion to the first grade from α = 1 to any α ∈ [0, 1]. To do this, we use the gradient f A = f x 1 A , f x 2 A , , f x n A , such that A = A 1 , A 2 , , A n . Therefore, B _ α can be developed as follows:
B _ α B T _ α = f A + i = 1 m f x i A · r A i α 1 i = m + 1 n f x i A · r A i α 1 = f A i = 1 n f x i A · r A i 1 α .
Analogously, we develop B ¯ α as follows:
B ¯ α B T ¯ α = f A i = 1 m f x i A · r A i α 1 + i = m + 1 n f x i A · r A i α 1 = f A + i = 1 n f x i A · r A i 1 α .
Similarly, B * β = B * _ β , B * ¯ β , whose exact values are given in (20), can also be determined via Taylor expansion to the first grade from β = 0 to β ( 0 , 1 ] . Therefore, for B * _ β , we state
B * _ β B T * _ β = f A i = 1 m f x i A · r A ( i ) * β + i = m + 1 n f x i A · r A ( i ) * β = f A i = 1 n f x i A · r A i * β .
In the same manner, we expand B * ¯ β as follows:
B * ¯ β B T * ¯ β = f A + i = 1 m f x i A · r A ( i ) * β i = m + 1 n f x i A · r A ( i ) * β = f A + i = 1 n f x i A · r A ( i ) * β .
Consequently, from (22), (23), (24), and (25), we find that
B = f A ,     r B = i = 1 n f x i A · r A i ,     r B * = i = 1 n f x i A · r A ( i ) * .
Analogous to [68], we evaluate the relative error measurement in the bounds of B α , β , whose exact value can be calculated with (19)–(20), by those of its symmetrical triangular approximation, B α , β T , which must be stated by applying (26) in (11) and (12). Thus, the deviations in B α = B _ α , B ¯ α are
ε _ α = B _ α B T _ α B _ α , ε ¯ α = B ¯ α B T ¯ α B ¯ α ,
and for B * β = B * _ β , B * ¯ β ,
ε * _ β = B * _ β B T * _ β B * _ β , ε * ¯ β = B * ¯ β B T * ¯ β B * ¯ β .

2.3. Intuitionistic Linear Regression with the Minimum Fuzziness Principle and Symmetric Coefficients

Within the fuzzy field, there are a large number of regression methodologies that can be divided into two main groups. In the first type, we can group those based on the minimum fuzziness principle (MFP), and the second includes those based on the minimization of the distance between observations and predictions, such as models that can be labeled fuzzy least-squares models [69].
This dichotomy between the minimum fuzziness principle and distance minimization is also observed in intuitionistic fuzzy regression models. For example, [46] extended the minimum fuzziness principle to an intuitionistic regression setting, and [70] developed a least-squares intuitionistic fuzzy regression methodology. Our paper uses the first approach. Therefore, our model is based on fuzzy regression models with symmetric parameters.
Let us suppose that the equation to be fitted has a dependent factor dependent on m real-value explanatory variables x i   , i = 0, 1, 2, …, m, where x 0 = 1 and x i ϵ R , i = 1, 2, …, n. The outcome is then a linear function of intuitionistic coefficients A ~ ( i ) = A ( i ) , r A ( i ) , r A ( i ) *   ,   i = 0 , 1 , , m and, thus, an STIFN Y ~ = Y , r Y , r Y *   . This is obtained from (21), as follows:
Y = j = 0 m A ( i ) x i ,   r Y = i = 0 n r A ( i ) x i ,   r Y * = i = 0 n r A ( i ) * x i .
Moreover, both the observations of the input variables and the output variable are crisp, which is a common hypothesis in intuitionistic fuzzy regression models. Thus, for the jth observation, the outcome is the crisp number y j , generated by the crisp income ( 1 , x 1 j ,   x 2 j , , x i j , , x i m ) . Therefore, y j is a possible value of a TIFN Y j ~ = Y j , r Y j , r Y j *   , whose membership function μ Y j y j and nonmembership function ν Y j y j in (11)–(12) are determined from (29)
Y j = i = 0 m A ( i ) x i j ,   r Y j = i = 0 m r A ( i ) x i j ,   r Y j * = i = 0 m r A ( i ) * x i j .
The objective is to fit for A ~ ( i ) ,   i = 0 , 1 , 2 , , m , an STIFN estimate a ~ ( i ) = a ( i ) , r a ( i ) , r a ( i ) *   ,   i = 0 , 1 , 2 , , m that simultaneously maximizes the membership of the observations in the fitted system and minimizes the uncertainty of that system. Therefore, to find a ~ ( i ) , the following multiple-objective programming problem must be implemented:
minimize A ( i ) , r A ( i ) r A ( i ) * , i = 0 , 1 , , n α ,   β ,   z 1 = j = 1 n   r Y j   ,   z 2 = j = 1 n   r Y j *   ,
which is subject to
μ Y j y j   α ,   ν Y j y j β ,   j = 1 , 2 , , n , r A i , r A i * 0 ,   i = 0 , 1 , ,   m .
0 α + β 1 , α , β 0 , 1 .
To solve (31), we implement the following steps:
Step 1: We state a minimum reachable value α = g and β = 1 − g h. Like in possibilistic regression models, g ∈ [0, 1) scales the total fuzziness of the estimated system. If g = 0, the uncertainty of the system is minimal; conversely, the inclusiveness of the observations may be low. On the other hand, a higher g causes all observations to be included with greater intensity, and the predictions of the fitted system are less specific [49].
The value of h ∈ [0, 1 − g ) reflects the level of hesitancy in the system. For h = 0, the actual and potential possibility of a particular value are identical; therefore, we have a conventional possibilistic regression. At this step, we decouple (31) as follows:
minimize A ( i ) , r A ( i ) , i = 0 , 1 , , n z 1 = j = 1 n   r Y j     ,
subject to
μ Y j y j   g ,   j = 1 ,   2 , ,   n ,   r A ( i ) 0 ,   i = 0, 1 , , m ,
and
minimize A ( i ) , r A ( i ) * , i = 0 , 1 , , n z 2 = j = 1 n   r Y j *   ,
subject to
ν Y j y j 1 g h ,   = 1 ,   2 , ,   n ,   r A ( i ) * 0 ,   i = 0 , 1 , , m .
Step 2: We initially state for (32)–(33) that g = h = 0 . This implies the minimum fuzziness level and no hesitancy. Thus, we adjust a possibilistic regression model and r A ( i ) = r A ( i ) * . This leads us to obtain the estimates of A ( i ) and r A ( i )   , which we denote as a ( i ) ( 0 ) and r a ( i ) ( 0 ) , respectively, where i = 0, 1, …, m. Thus, we must solve the following:
minimize A ( i ) , r A ( i ) , i = 0 , 1 , , n z 1 = z 2 = i = 0 m r A ( i ) j = 1 n x i j   ,
subject to
i = 0 m A i x i j i = 0 m r A i x i j y j i = 0 m A i x i j + i = 0 m r A i x i j ,   j = 1 , 2 , , n
A ( i ) , r A ( i ) 0 ,   i = 0 , 1 , , m .
Step 3: To fit the centers in (34), a ( i ) ( 0 ) for A i and r a ( i ) ( 0 ) for r A i i = 0, 1, …, m, the literature proposes two alternatives:
  • Alternative 1. The values of a ( i ) ( 0 ) and r a ( i ) ( 0 ) are those that are solved in a unique step (34). In this case, the centers a ( i ) ( 0 ) are those that are obtained in a quantile regression at the median, identical to [71].
  • Alternative 2. The value a ( i ) ( 0 ) in the first step is obtained by using ordinary least squares [72]. However, there is no reason why any other method, such as the maximum likelihood or weighted least-squares methods, cannot be used. In the second step, r a ( i ) ( 0 ) is obtained by solving (34) and taking into account that this linear programming problem is as follows, after independently stating a ( i ) ( 0 ) :
minimize r A ( i ) , i = 0 , 1 , , n z 1 = z 2 = i = 0 m r A ( i ) j = 1 n x i j   ,
subject to
i = 0 n r A i x i j y j i = 0 n a ( i ) ( 0 ) x i j i = 0 n r A i x i j ,   j = 1 , 2 , , n
r A ( i ) 0 ,   i = 0 , 1 , , m .
Step 4: We establish the optimal value of g based on this criterion. This value optimizes what these authors refer to as the credibility of the system [73]. To achieve this, we define the estimation of Y j ~ , obtained from the parameters adjusted in Step 1 and Step 3 as y ~ j ( 0 ) = y j ( 0 ) , r y j ( 0 ) , r y j ( 0 ) , i.e., y ~ j ( 0 ) is an STFN where y j ( 0 ) = i = 0 m a ( i ) x i j , and r y j ( 0 ) = i = 0 m r a ( i ) ( 0 ) x i j . So,
g = 1 2 1 γ ( 0 ) δ ( 0 ) γ ( 0 ) < δ ( 0 ) 0 o t h e r w i s e ,
where γ ( 0 ) = j = 1 n μ y j ( 0 ) y j r y j ( 0 )   , δ ( 0 ) = j = 1 n 1 μ y j ( 0 ) y j r y j ( 0 ) and then we state that
r a ( i ) = r a ( i ) ( 0 ) 1 g .
Step 5: We subsequently proceed to obtain the estimates of r a ( i ) * . To achieve this, the decision maker must determine the degree of hesitancy in the system, where h ϵ [ 0 , 1 g ) . In the case where h = 0 , there is no hesitancy; if h 1 g , the level of hesitancy tends to be at its maximum. Thus,
r a ( i ) * = r a ( i ) ( 0 ) 1 g h .

3. An Intuitionistic Chain Ladder for Claim Reserving

3.1. Claim Reserving with the Chain-Ladder Method and Stochastic Variability and a Probability–Possibility Transformation

The historical data illustrating the evolution of claims are typically presented in a run-off triangle format, similar to Table 2 [10]. In this table, Ci,j represents the accumulated claim cost of insurance contracts originating in the ith development period (i = 0, 1, …, n) during the jth claiming period (j = 0, 1, …, n). Therefore, the accumulated claims Ci,j, i = 1, 2, …, n; j = ni + 1, ni + 2, …, n are unknown and must be fitted.
An alternative way to present historical data consists of the run-off triangle of incremental claims, in a way similar to Table 3. Table 3 can be obtained from Table 2 by taking into account that Si,j = Ci,jCi,j−1, i = 0, 1, 2, …, n − 1, j = 1, 2, …, ni, and Si,0 = Ci,0. Therefore, the incremental claims Si,j, i = 1, 2, …, n; j = ni + 1, ni + 2, …, n are unknown and must be fitted.
The triangle of accumulated claims (Table 2) is the input of several common methods to fit claim reserves, such as the chain-ladder method or the Bornhuetter–Ferguson method. The key concept of the CL method is the so-called link ratio between development year j and j + 1, F j , which allows us to obtain the cumulative claims of the (j + 1)th development period from those of the jth period:
C i , j + 1 = F j · C i , j F j = C i , j + 1 C i , j ,
where the available observations of F j are as follows:
f i , j = C i , j + 1 C i , j ,   i = 0 , 1 , , n j 1
To obtain an average value of F j , f j , we consider the widely used CL, which provides an unbiased estimator of F j [11]. Thus, the average development factor for the jth year is
f j = i = 0 n j 1 C i , j + 1 i = 0 n j 1 C i , j .
The terminal value of accumulated claims for the ith origin year Ci,n, i = 1, 2, …, n is approximated by c i . n , as follows:
c i . n = C i , n i j = n i n 1 f j ,
and c i . n is an increasing function of development factors, since the partial derivative c i . n f j is
c i . n f j = C i , n i k = n i k j n 1 f k .
Thus, the reserves are linked with the origin year i = 1, 2, …, n, R O i :
R O i = c i . n C i , n i = C i , n i j = n i n 1 f j 1 .
So, the overall provisions, R, are
R = i = 1 n R O i .
The classical chain-ladder method is deterministic. However, this methodology is flexible enough to generate stochastic estimates of fluctuations by applying the SCL method, which is implemented in the following six steps:
  • Obtain the estimates of the observations C i , j , c i , j , i = 0, 1, …, n; j < ni by using f j backwards from C i , n i .
  • Calculate an estimate of observed incremental claims (Table 2) by stating s i , j = c i , j c i , j 1 , in the case of s i , 0 = c i , 0 .
  • Calculate the descaled Pearson residuals due to fitting the real incremental claims in Table 2, S i , j , with s i , j :
    m i , j = S i , j s i , j s i , j , i = 0 , 1 , , n ; j n i .
  • Resample m i , j , i = 0, 1, …, n; j n i . Therefore, we find m i , j b ,   i = 0, 1, …, n; j n i .
  • Calculate the incremental claims s i , j b = s i , j + s i , j m i , j b , i = 0, 1, …, n; j n i . This implies adjusting a new Table 3.
  • From Table 3, in the above step, we can resample the accumulated claims and construct Table 2. This new table allows us to obtain the development factors (40) and reserves (44) and (45). These six steps can be implemented B times in such a way that predictions of claiming reserves can be obtained as confidence intervals.
Note that Steps 1–6 allow B simulations of loss reserves to be obtained for every origin year R O i ( b i ) ,   b i = 1, 2, …, B and the whole reserve R ( b )   a n d   b = 1 , 2 , , B . Without losing generality, let us suppose that in all the cases, R O i ( b i ) R O i ( b i + 1 ) . Then, the reserve R O i is contained with a probability 1 α in the interval R O i α , such that
R O i α = R O i _ α , R O i ¯ α = R O i r o u n d B · α 2 , R O i r o u n d B · 1 α 2 ,
which can be interpreted as the α-cuts of a possibilistic estimate of the reserves of the ith year R O ¨ i .
Therefore, we can estimate a confidence interval for the overall reserves in two ways. A conservative estimate, R α , is R α = i = 1 n R O i α , such that R α can be considered the α-cuts of the possibility distribution R ¨ :
R α = R _ α , R α ¯ = i = 1 n R O i _ α , i = 1 n R O i ¯ α
A more specific approximation of the overall reserves implies inducing a confidence interval with a probability level 1 α ,   R α from R ( b ) R b + 1 ,  b = 1, 2, …, B by calculating the following:
R α = R _ α , R ¯ α = R r o u n d B · α 2 , R r o u n d B · 1 α 2 .
Therefore, from the probabilistic confidence interval R α , we can induce a possibility distribution function R ¨ by considering Definition 5.

3.2. An Intuitionistic Fuzzy Chain-Ladder Method

3.2.1. Fitting Symmetrical Intuitionistic Triangular Fuzzy Development Factors

Let us express relation (39), in which, from a known accumulated claim amount in the jth development period, we must obtain the accumulated quantity in the (j + 1)th development period, which is uncertain because the development factor is an STIFN. So,
C ~ i , j + 1 = F ~ ( j ) · C i , j
where C ~ i , j + 1 = C i , j + 1 , r C i , j + 1 , r C i , j + 1 * and F ~ ( j ) = F ( j ) , r F ( j ) , r F ( j ) * . Therefore, from (21), C i , j + 1 = F ( j ) · C i , j , r C i , j + 1 = C i , j · r F ( j ) and r C i , j + 1 * = C i , j · r F ( j ) * .
To fit F ~ ( j ) by means of f ~ ( j ) = f ( j ) , r f ( j ) , r f ( j ) * , we consider the data in Table 2. The couples (y,x) are defined as C i , j + 1 , C i , j , i = 0 , 1 , , n j 1 . Therefore, f ~ ( j ) = f ( j ) , r f ( j ) , r f ( j ) * is first fitted, where f ( j ) ( 0 ) = f ( j ) and r f ( j ) ( 0 ) = r f ( j ) = r f ( j ) * are the optimum values of the arguments in the programming problem (32)–(33) for g = h = 0 . Therefore, we must solve the version of (34) for relation (49):
minimize F ( j ) , r F ( j ) z 1 = z 2 = r F ( j ) i = 0 n j 1 C i , j   ,
which is subject to
F j C i , j C i , j r F j C i , j + 1 F j C i , j + C i , j r F j ,   i = 0 , 1 , , n j 1 .
r F j 0 .
By dividing the inclusion constraints in (50) by C i , j , i = 0 , 1 , , n j 1 , the independent terms turn into (40), f i , j = C i , j + 1 C i , j . Likewise, the cost function of (50) has only one argument. Therefore, the linear pro-gramming problem becomes
minimize F ( j ) , r F ( j ) z 1 = z 2 = r F ( j ) ,  
which is subject to
F j r F j f i , j F j + r F j ,   i = 0 , 1 , , n j 1 .
r F j 0 .
To solve (51), we can follow Alternatives 1 and 2 in Section 2.3. By using Alternative 1, the solution of that linear programming problem allows us to obtain f ( j ) ( 0 ) as the result of the quantile regression at the median and, simultaneously, r f ( j ) ( 0 ) . Models (51) can be implemented by Alternative 2 in Section 2.3 by prefixing f ( j ) ( 0 ) with the CL formula (41). In this case, the linear programming problem (51) becomes
minimize r F ( j ) z 1 = z 2 = r F ( j )   ,
which is subject to
r F j f i , j f ( j ) ( 0 ) r F j ,   i = 0 , 1 , , n j 1 .
r F j 0 ,
and so,
r f ( j ) ( 0 ) = max i = 0 , 1 , , n i 1 f i , j f ( j ) ( 0 ) .
Then, the empirical estimates of F ~ ( j ) = F ( j ) , r F ( j ) , r F ( j ) * , f ~ ( j ) = f ( j ) , r f ( j ) , r f ( j ) * , j = 0, 1, …, n − 1, are obtained as follows:
i.
Considering that f ( j ) = f ( j ) ( 0 ) ,
ii.
r f ( j ) is obtained by calculating g in Step 4 of Section 2.3 with (36) and (37).
iii.
Finally, r f ( j ) * is adjusted in Step 5 of Section 2.3 by subjectively stating the degree of system hesitancy, h , and using (38).
The structure of the data in Table 2 leads us to obtain the last development factor f ~ ( n 1 ) with only the pair C 0 , n , C 0 , n 1 . Therefore, it is easy to verify that f ( j ) = f ( j ) ( 0 ) = C 0 , n C 0 , n 1 , but this approach also leads to the unrealistic conclusion that it is a certain parameter, i.e., r f ( j ) = r f ( j ) * = 0 . Mack [11], in his stochastic free-distribution modeling of reserves over the CL model, addresses this issue based on the intuition that the absolute uncertainty of the development factors tends to decrease over time, as does the expected value of these factors. Thus, the standard deviation of the development factor F n 1 is estimated as the minimum of the standard deviation of F n 3 and F n 2 and the ratio between the variance of F n 3 and the standard deviation of F n 2 . Taking this idea into consideration and considering that k times the standard deviation of random quantities can be interpreted as the radius of an equivalent STFN [54,56], f ~ ( n 1 ) = f ( n 1 ) , r f ( n 1 ) , r f ( n 1 ) *   , where
f ( n 1 ) = C 0 , n C 0 , n 1 ,  
and
r f ( n 1 ) = min r f ( n 3 ) 2 r f ( n 2 ) , r f ( n 3 ) , r f ( n 2 ) ,   r f ( n 1 ) * = min r f ( n 3 ) * 2 r f ( n 2 ) * , r f ( n 3 ) * , r f ( n 2 ) * .

3.2.2. Fitting Reserves with Symmetric Triangular Intuitionistic Fuzzy Development Factors

To state the reserves, we must estimate the terminal value of the claims in every origin year i = 1, 2, …, n, c ~ i , n , which can be expressed through its α , β -cuts as follows:
c i , n α , β = c i , n α = c i , n _ α , c i , n ¯ α , c i , n * β = c i , n * _ β , c i , n * ¯ β , 0 α + β 1 , α , β ( 0 , 1 ) .
Specifically, c i , n α , β is obtained from C i , n i and the α , β -cuts of f ~ ( j ) , f ( j ) α , β ,   j = ni, ni + 1, …, n − 1 by adapting (42) to an intuitionistic setting:
c i , n α , β = C i , n i j = n i n 1 f ( j ) α , β ,
and thus, c i , n α = c i , n _ α , c i , n ¯ α is obtained considering that (42) is an increasing function of development factors:
c i , n α = C i , n i j = n i n 1 f j _ α , C i , n i j = n i n 1 f j ¯ α = C i , n i j = n i n 1 f ( j ) r f ( j ) 1 α   , C i , n i j = n i n 1 f ( j ) + r f ( j ) 1 α   .
Similarly, c i , n * β = c i , n * _ β , c i , n * ¯ β is calculated as follows:
c i , n * β = C i , n i j = n i n 1 f ( j ) * _ β , C i , n i j = n i n 1 f ( j ) * ¯ β = C i , n i j = n i n 1 f ( j ) r f ( j ) * β , C i , n i j = n i n 1 f ( j ) + r f ( j ) * β .
Note that c ~ i , n is not an STIFN. However, by using derivatives (43) and (22)–(25), we can approximate c ~ i , n c ~ i , n T = c i , n , , r c i , n , r c i , n * , where the center is
c i , n = C i , n i j = n i n 1 f ( j ) ,
and the radii are
r c i , n = C i , n i j = n i n 1 k = n i k j n 1 f k r f ( j ) , r c i , n * = C i , n i j = n i n 1 k = n i k j n 1 f k r f ( j ) * .
Therefore, by using (44), we can obtain the intuitionistic reserves for the ith origin year R O ~ i through R O i α , β = R O i α = R O i _ α , R O i ¯ α , R O i * β = R O i * _ β , R O i * ¯ β , 0 α + β 1 , α , β ( 0 , 1 ) by calculating:
R O i α , β = c i , n α , β C i , n i .
Then,
R O i α = C i , n i j = n i n 1 f j r f j 1 α   1 , C i , n i j = n i n 1 f j + r f j 1 α   1 , R O i * β = C i , n i j = n i n 1 f ( j ) r f ( j ) * 1 α   1 , C i , n i j = n i n 1 f j + r f ( j ) * 1 α   1 .
The intuitionistic fuzzy estimate of reserves of the ith year is not an STIFN. However, an STIFN approximate R O ~ i R O ~ i T = R O i , r R O i , r R O i * is obtained by the following:
R O ~ i T = c i , n T C i , n i = c i , n C i , n i , r c i , n , r c i , n * ,
and so, considering (60) and (61),
R O i = c ~ i , n T C i , n i = C i , n i j = n i n 1 f j 1 r R O i = C i , n i j = n i n 1 k = n i k j n 1 f k r f ( j ) , r R O i * = C i , n i j = n i n 1 k = n i k j n 1 f k r f ( j ) * .
Similarly, an intuitionistic fuzzy estimate of the overall reserve R ~ is obtained with (45) through R α , β = R α = R _ α , R ¯ α , R β * = R * _ β , R * ¯ β , 0 α + β 1 , α , β ( 0 , 1 ) . By implementing R α , β = i = 1 n R O i α , β . Then,
R α = i = 1 n R O i _ α , i = 1 n R O i ¯ α , R β * = i = 1 n R O i * _ β , i = 1 n R O i * ¯ β .
Therefore, an STIFN approximate to R ~ R ~ T = R , r R , r R * is obtained simply as follows:
R ~ T = i = 1 n R O ~ i T = i = 1 n R O i , i = 1 n r R O i , i = 1 n r R O i * .

4. Empirical Application

4.1. Estimating Loss Reserves with Deterministic and Stochastic Chain-Ladder Method

Below, we present an empirical application based on the run-off triangle of accumulated claims shown in Table 4. These data were utilized in [74,75]. Table 4 also illustrates the development factors found using (41). Thus, we observe that a crisp development factor f 0 = 1.899 is estimated, indicating that the accumulated claims from development years zero to one increase on average by 89.90% for all origin years. Similarly, we can interpret the estimates of the development factors f 1 , f 2 , f 3 , and f 4 .
Table 5 presents the individual reserves obtained for each of the origin years i = 1, 2, …, 5 and the total reserves with a deterministic CL. Thus, we can observe that as the origin year increases, the reserve to be allocated increases, as claims from more development years are pending. It can be noted that in both Table 4 and Table 5, we obtain the expected values of the link ratios and reserves, but we do not have any estimation of their variability. This analysis is carried out in Table 6 and Table 7, where reserves are estimated using the SCL method, and the obtained possibilistic confidence intervals are interpreted as possibility distributions, using Definition 5 of Section 2.
Table 6 displays a table of incremental claims analogous to Table 3 that is deduced from Table 4. Table 4 also shows the theoretical table of incremental claims that are deduced from the development factors of the chain-ladder method. The difference between the observed and theoretical tables of incremental claims through descaled Pearson residuals allows the implementation of the SCL method, described in Section 3.1, to fit the variability of reserves by origin year and total reserves.
Table 7 presents the results obtained with B = 5000 bootstrapping resamples. The confidence intervals were calculated with Equations (46)–(48). The upper endpoints of the confidence intervals obtained for probability levels α = 0, 0.01, 0.05, and 0.1 are the 100%, 99.5%, 97.5%, and 95% estimated percentiles for the reserves, respectively. These quantiles are commonly used to estimate extreme claim scenarios.
Within the total reserves, we distinguished two confidence intervals: R α (47) and R α (48). The former arises from adding the confidence intervals associated with the reserves of each origin year. Thus, as shown in Table 7, for a confidence level of 100% ( α = 0 ) , we obtained the overall reserves R 0 = [63.68, 97.77] + [526.48, 628.62] + [1507.29, 1677.46] + [2720.66, 2997.42] + [4641.72, 5105.40] = [9459.82, 10,506.67]. The most prudent reserve value would be 10506.67 since it arises from the sum of the estimates of the value that accumulates 100% probability of the reserves from each origin year. In contrast, R α arises from the application of (45) in each of the B = 5000 simulations, making it a narrower confidence interval. Table 7 shows that R 0 = [9533.03, 10,481.08], so the prudent value for the reserve based on this confidence interval is 10,481.08.
Table 7 also shows that the reserves of each origin year, and the total reserves can be estimated through a possibility distribution by gathering and fitting successive confidence intervals (46)–(48) from α = ε (ε 0) to α = 1 . By applying Definition 5 and Remark 3, we can adjust the reserves of the fifth period to a possibility distribution R O ¨ 5 whose core is 4826.23 and support [4641.72, 5105.40]. Likewise, Table 6 also shows that we can obtain a possible estimate of overall reserves R ¨ = i = 1 5 R O ¨ i whose center is 9899.31, supporting [9459.82, 10,506.67].

4.2. Estimating Loss Reserves with a Symmetric Triangular Intuitionistic Fuzzy Chain Ladder

Next, we put into work the methodology developed in Section 3.2, which allows us to estimate the claim reserves with STIFNs, with the data of the run-off triangle in Table 4. We also compare the results obtained with those of the bootstrap estimates using the SCL method in Table 6, which we reinterpret as α-cuts of possibility distributions. Therefore, to obtain the estimates of F ~ ( 0 ) = F ( 0 ) , r F ( 0 ) , r F ( 0 ) *   and f ~ ( 0 ) = f ( 0 ) , r f ( 0 ) , r f ( 0 ) *   , we solve the linear programming problem (51), whose constraints are built up with the link ratios of each origin year i = 1, 2, …, 5, as shown in Table 8:
minimize F ( 0 ) , r F ( 0 ) z 1 = z 2 = r F ( 0 ) ,  
which is subject to
F 0 r F 0 1.853 F 0 + r F 0 ,
F 0 r F 0 1.889 F 0 + r F 0 ,
F 0 r F 0 1.923 F 0 + r F 0 ,
F 0 r F 0 1.928 F 0 + r F 0 ,
F 0 r F 0 1.890 F 0 + r F 0 ,
r F 0 0 .
Table 9 shows the STIFNs adjusted to development factors for j = 0, 1, 2, 3, 4. Thus, if the estimate F 0 is not prefixed with (41), we obtain f ( 0 ) ( 0 ) = 1.891 and r f ( 0 ) ( 0 ) = 0.038 , and (36) and (37) allow us to obtain an optimum uncertainty level for membership functions of development factor g = 0.14 . Thus, from (38), r f ( 0 ) = 0.044 . The degree of system hesitancy, h , must be estimated subjectively by the decision maker. This may be linked, for example, to the perceived reliability of the data or the predictability of the insurance environment. The calculations in this numerical application are performed with h = 0.1 , so we obtain r f ( 0 ) * = 0.049 .
In Section 3.2, we also state that f ( 0 ) ( 0 ) can be predefined with the deterministic CL shown in Table 4. Therefore, f ( 0 ) ( 0 ) = 1.899   , and by using (51), r f ( j ) ( 0 ) = 0.046 . Equations (36) and (37) allow us to obtain an optimum uncertainty degree g = 0 . Therefore, from (38) r f ( 0 ) = 0.046 and by using the hesitancy level h = 0.1 , r f ( 0 ) * = 0.051 .
Note that the spreads of f ~ ( 4 ) , r f ( 4 ) and r f ( 4 ) * cannot be obtained from the sample in Table 4 since only one individual link ratio exists. To fit these spreads, we use (55) and then set the following:
r f ( 4 ) = min 0.0035 2 0.0067 ; 0.0035 ; 0.0067 = 0.0035
r f ( 4 ) * = min 0.0042 2 0.0081 ; 0.0042 ; 0.0081 = 0.0042 .
Table 10 shows the estimates of the overall loss reserves of the intuitionistic claim reserves calculated with the two alternatives proposed in Section 3.2, R ~ . Thus, first, we compute the “exact” α , β -cuts of both methods. This involves using (58) and (59) to determine the terminal accumulated claims; (63) to find R O ~ i ; and i = 1, 2, …, 5 and (66) to determine the total value. Table 10 also shows the STIFN approximation of the total reserve, R ~ T , which is obtained using the sequential use of (60), (61), (65), and (67). Table 10 also shows the errors calculated with (27) and (28). Their values suggest that the symmetric triangular approximation is almost perfect. Note that the maximum error lies in the β -cuts of the nonmembership function at β = 1 and does not exceed 0.15% in any case.
The results of Table 7, which come from bootstrapping resamples, can be interpreted as α-level sets of possibility distributions. Therefore, they can be compared with fuzzy intuitionistic estimates, which are constructed through two possibility distributions. In other words, the probabilistic intervals obtained with bootstrapping and the α , β -cuts can be interpreted by the actuary in a similar manner. Thus, according to Table 7, the value of reserves that includes 100% of their possible values could be given as 10,506.67 if we sum the 100th percentile of the reserves associated with all origin years, and 10,481.08 if we consider the 100th percentile of the bootstrap simulations of overall reserves. These results are similar and comparable to those obtained with the membership function of the overall reserves obtained in Table 10. We can observe in the α -cuts of the reserves, R α T , that if they are calculated with Alternative 1, their prudent estimate can range between 10,391.15 (at the 0.25-cut) and 10,565.45 (at the 0-cut). The conclusions we can draw from the fit obtained with Alternative 2 are similar, as the upper end of the 0-cut is 10,563.24 and that of the 0.25-cut is 10,397.25.
The β -cuts of the nonmembership functions complement the information provided by the α-cuts of the membership functions, introducing the existence of bipolarity. Thus, in Table 10, Alternative 1 for estimating the development factors offers an upper bound at the 0-cut for the loss reserves of R T ¯ 0 = 10,565.45 and an upper limit of the 1-cut of the nonmembership function, R T * ¯ 1 = 10,688.71. Thus, the use of IFNs in reserve estimation allows us to first obtain an estimation of the most extreme possible scenario (10,565.45), whose adjustment does not use subjective information at any time but rather uses only run-off triangle data. That is, the meaning of the estimation is analogous to that obtained with stochastic simulation or what we would obtain with the use of possibilistic regression. However, the use of IFNs also allows us to obtain an estimation of the scenario that we could classify as potentially more extreme through the higher value of the 1-cut of the nonmembership function. The quantification of this scenario requires the participation of the decision maker, who must indicate a perceived degree of hesitancy, which in this numerical application was h = 0.1 .
Using other types of modeling for the underlying link ratios in the run-off triangle, such as LR or adaptive functions, allows us to obtain the same <0,1>-cut and <1,0>-cut as our method if the IFN estimates have the same centers and radii. However, the rest of the α , β -cuts, which can be assimilated to structured simulations of the variables involved in the analysis, would change, varying their amplitude. In the case of using adaptive functions, the linear functions used in this paper can be considered as a baseline, with an order of m = 1. From this baseline, m < 1 implies a dilation of the α , β -levels, and thus, they will incorporate more uncertainty. In contrast, m > 1 indicates a contraction of the results compared to those obtained with STIFNs. Thus, the α , β -cuts will have a smaller width.
Table 11 shows the estimated reserves associated with the five origin years and the overall reserves through the STIFNs. The use of this type of IFN can be very useful for applying the actuarial judgement required to set a definite crisp value for the loss reserves. In the case of total reserves, if we take Alternative 2 from Section 3.2 as a reference for decision making, the most reliable value is 9899.31, which coincides with (40). Possible deviations of up to 663.93 are estimated, and deviations exceeding 765.15 are considered not possible. Regarding deviations between 663.93 and 765.15, there is hesitancy about their feasibility. When considering only the data in the run-off triangle, the conclusion must be that they are not possible. On the other hand, they are possible based on the degree of hesitancy perceived by the decision maker.

5. Conclusions and Further Research

The determination of insurance loss reserves must be prudent, necessitating the quantification of their expected value and potential deviations from that value. To ascertain the most plausible value, a statistical method such as the chain-ladder (CL) method is utilized to estimate the expected claim evolution. Subsequently, it is necessary to estimate possible deviations from these values with greater reliability. The contributions of this work include providing tools for estimating and interpreting such values using fuzzy set theory and intuitionistic fuzzy set theory.
The first contribution of this work is to show that the information obtained through stochastic models such as bootstrapping and the use of conventional fuzzy numbers are similar. In fact, we can reinterpret the value and variability of reserves obtained with the stochastic CL (SCL) methodology with possibility distributions. Therefore, both instruments capture epistemic uncertainty.
The second and main contribution of our work is the generalization of developments in claim reserving with fuzzy numbers to the use of intuitionistic fuzzy numbers (IFNs). This tool allows the introduction of bipolar information about possible reserve variability into the estimation, i.e., both “positive” information about feasible parameter values and negative information about those that cannot be taken in any case.
This work assumes that the parameters governing the evolution of claims are symmetrical and triangular IFNs (STIFNs). Special attention is given to the approximation of each IFN to be of the same nature as the results that arise from its functional handling. Linear shapes often provide effective resolution in practical applications of fuzzy set theory. Moreover, symmetry often allows for a good balance between parsimony and comprehensiveness in capturing available information and facilitates interpretability of the results by end-users who may not necessarily have knowledge of fuzzy logic. The value of loss reserves when the development factors are estimated using the STIFN technique can be easily approximated through the most likely scenario, obtained with conventional chain-ladder methodology, and by evaluating the deviations from this value with the gradient function of the terminal value of claims from each origin year in the spreads of the membership and nonmembership functions of the link ratios.
The results provided by the proposed method can be very useful in actuarial practice since they can be interpreted very intuitively by the person responsible for establishing reserves, as there is no need for knowledge of fuzzy set theory. While the center of an STIFN quantifies, in a very synthetic way, the most reliable value of reserves, the two spreads provide an approximation of the maximum deviations from this value, the maximum achievable deviation, and the first not-achievable deviation. On the other hand, representing reserves through α , β -cuts allows for the structuring of simulations on their appropriate value in multiple scenarios, which can be of great help to decision makers.
Certainly, the limitation of using STIFNs is that they do not account for asymmetry in the link ratios, and similarly, they do not allow for the introduction of more refined calibration of possibility distributions, such as adaptive membership functions. This latter issue implies that introducing nuances, such as concentration and dilation, is not possible. However, the proposed scheme can be adapted to accommodate more sophisticated forms of membership and nonmembership functions.
Our extension of intuitionistic regression can be applied in other financial and actuarial contexts where possibilistic regression has already been used, such as, for example, estimating the implied moments of options [76,77]. A natural extension of this work would involve introducing intuitionistic uncertainty into the analysis of non-life insurance claims, expanding the results obtained with fuzzy numbers to calculate discounted reserves [78], the discounted values of non-life insurance liabilities [79], or the terminal values of an insurance company [80,81].

Funding

This research benefited from the Research Project of the Spanish Science and Technology Ministry “Sostenibilidad, digitalizacion e innovacion: nuevos retos en el derecho del seguro” (PID2020-117169GB-I00).

Institutional Review Board Statement

Not applicable.

Data Availability Statement

The data used in the paper have been referenced.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Hindley, D. Introduction. In Claims Reserving in General Insurance. International Series on Actuarial Science; Hindley, D., Ed.; Cambridge University Press: Cambridge, UK, 2017; pp. 1–15. [Google Scholar]
  2. Hindley, D. Stochastic Reserving Methods. In Claims Reserving in General Insurance. International Series on Actuarial Science; Hindley, D., Ed.; Cambridge University Press: Cambridge, UK, 2017; pp. 146–319. [Google Scholar]
  3. Andrés Sánchez, J. Calculating insurance claim reserves with fuzzy regression. Fuzzy Sets Syst. 2006, 157, 3091–3108. [Google Scholar] [CrossRef]
  4. Bastos, I.S.; Vana, L.B.; Novo, C.C. Estimating IBNR claim reserves using Gaussian Fuzzy Numbers. Context. Rev. Cont. Econ. Gest. 2023, 21, e83343. [Google Scholar] [CrossRef]
  5. Shapiro, A.F. Fuzzy logic in insurance. Insur. Math. Econ. 2004, 35, 399–424. [Google Scholar] [CrossRef]
  6. Derrig, R.A.; Ostaszewski, K.M. Fuzzy Set Theory. In Encyclopedia of Actuarial Science; Teugels, J.F., Sundt, B., Asmussen, S., Eds.; John Willey and Sons Ltd.: Chichester, NH, USA, 2006. [Google Scholar]
  7. Dubois, D.; Prade, H. The three semantics of fuzzy sets. Fuzzy Sets Syst. 1997, 90, 141–150. [Google Scholar] [CrossRef]
  8. Hindley, D. Data. In Claims Reserving in General Insurance. International Series on Actuarial Science; Hindley, D., Ed.; Cambridge University Press: Cambridge, UK, 2017; pp. 16–39. [Google Scholar]
  9. Straub, E. Nonlife Insurance Mathematics; Springer: Berlin/Heidelberg, Germany, 1997; pp. 102–115. [Google Scholar]
  10. Schmidt, K.D.; Zocher, M. Additive Method. In Handbook on Loss Reserving. EAA Series; Radtke, M., Schmidt, K.D., Schnaus, A., Eds.; Springer: Cologne, Germany, 2016. [Google Scholar]
  11. Mack, T. Distribution-free Calculation of the Standard Error of Chain Ladder Reserve Estimates. Astin Bullet. 1993, 23, 213–225. [Google Scholar] [CrossRef]
  12. England, P.D.; Verrall, R. Analytic and bootstrap estimates of prediction errors in claims reserving. Insur. Math. Econ. 1999, 25, 281–293. [Google Scholar] [CrossRef]
  13. Andrés-Sanchez, J.; Gómez, A.T. Applications of fuzzy regression in actuarial analysis. J. Risk Insur. 2003, 70, 665–699. [Google Scholar] [CrossRef]
  14. Heberle, J.; Thomas, A. Combining chain-ladder claims reserving with fuzzy numbers. Insur. Math. Econ. 2014, 55, 96–104. [Google Scholar] [CrossRef]
  15. Heberle, J.; Thomas, A. The fuzzy Bornhuetter–Ferguson method: An approach with fuzzy numbers. Ann. Actuar. Sci. 2016, 10, 303–321. [Google Scholar] [CrossRef]
  16. Taylor, G.C. Separation of inflation and other effects from the distribution of nonlife insurance claim delays. Astin Bullet. 1977, 10, 219–230. [Google Scholar] [CrossRef]
  17. Kremer, E. IBNR-claims and the two-way model of ANOVA. Scand. Actuar. J. 1982, 1, 47–55. [Google Scholar] [CrossRef]
  18. Yan, C.; Liu, T.; Dong, Q.; Liu, W. Payments Per Claim Method Based on Fuzzy Numbers. In Proceedings of the 14th International Conference on Natural Computation, Fuzzy Systems and Knowledge Discovery (ICNC-FSKD), Huangshan, China, 28–30 July 2018; pp. 643–648. [Google Scholar]
  19. Kim, J.H.; Kim, J. Fuzzy regression towards a general insurance application. J. Appl. Math. Inform. 2014, 32, 343–357. [Google Scholar] [CrossRef]
  20. Woundjiagué, A.; Mbele Bidima, M.L.D.; Waweru Mwangi, R. An Estimation of a Hybrid Log-Poisson Regression Using a Quadratic Optimization Program for Optimal Loss Reserving in Insurance. Adv. Fuzzy Syst. 2019, 1393946. [Google Scholar] [CrossRef]
  21. Apaydin, A.; Baser, F. Hybrid fuzzy least-squares regression analysis in claims reserving with geometric separation method. Insur. Math. Econ. 2010, 47, 113–122. [Google Scholar] [CrossRef]
  22. Woundjiagué, A.; Mbele Bidima, M.L.D.; Waweru Mwangi, R. A fuzzy least-squares estimation of a hybrid log-Poisson regression and its goodness of fit for optimal loss reserves in insurance. Int. J. Fuzzy Syst. 2019, 21, 930–944. [Google Scholar]
  23. Andrés-Sánchez, J. Claim reserving with fuzzy regression and Taylor’s geometric separation method. Insur. Math. Econ. 2007, 40, 145–163. [Google Scholar] [CrossRef]
  24. Yan, C.; Liu, Q.; Liu, J.; Liu, W.; Li, M.; Qi, M. Payments per claim model of outstanding claims reserve based on fuzzy linear regression. Int. J. Fuzzy Syst. 2019, 21, 1950–1960. [Google Scholar] [CrossRef]
  25. Andrés-Sánchez, J. Claim reserving with fuzzy regression and the two ways of ANOVA. Appl. Soft Comput. 2012, 12, 2435–2441. [Google Scholar] [CrossRef]
  26. Baser, F.; Apaydin, A. Calculating insurance claim reserves with hybrid fuzzy least squares regression analysis. Gazi Univ. J. Sci. 2010, 23, 163–170. [Google Scholar]
  27. Dubois, D.; Prade, H. An overview of the asymmetric bipolar representation of positive and negative information in possibility theory. Fuzzy Sets Syst. 2009, 160, 1355–1366. [Google Scholar] [CrossRef]
  28. Dubois, D.; Prade, H. Gradualness, uncertainty and bipolarity: Making sense of fuzzy sets. Fuzzy Sets Syst. 2012, 192, 3–24. [Google Scholar] [CrossRef]
  29. Mitchell, H.B. Ranking-intuitionistic fuzzy numbers. Int. J. Uncertain. Fuzziness Knowl.-Based Syst 2004, 12, 377–386. [Google Scholar] [CrossRef]
  30. Atanassov, K.T. Intuitionistic fuzzy sets. Fuzzy Sets Syst. 1986, 20, 87–96. [Google Scholar] [CrossRef]
  31. Atanassov, K.T. More on intuitionistic fuzzy sets. Fuzzy Sets Syst. 1989, 33, 37–45. [Google Scholar] [CrossRef]
  32. Kumar, G.; Bajaj, R.K. Implementation of intuitionistic fuzzy approach in maximizing net present value. Int. J. Math. Comput. Sci. 2014, 8, 1069–1073. [Google Scholar]
  33. Kahraman, C.; Çevik Onar, S.; Öztayşi, B. Engineering economic analyses using intuitionistic and hesitant fuzzy sets. J. Intell. Fuzzy Syst. 2015, 29, 1151–1168. [Google Scholar] [CrossRef]
  34. Boltürk, E.; Kahraman, C. Interval-valued and circular intuitionistic fuzzy present worth analyses. Informatica 2022, 33, 693–711. [Google Scholar] [CrossRef]
  35. Haktanır, E.; Kahraman, C. Intuitionistic fuzzy risk adjusted discount rate and certainty equivalent methods for risky projects. Int. J. Prod. Econ. 2023, 257, 108757. [Google Scholar] [CrossRef]
  36. Wu, L.; Liu, J.F.; Wang, J.T.; Zhuang, Y.M. Pricing for a basket of LCDS under fuzzy environments. SpringerPlus 2016, 5, 1–12. [Google Scholar] [CrossRef]
  37. Ersen, H.Y.; Tas, O.; Kahraman, C. Intuitionistic fuzzy real-options theory and its application to solar energy investment projects. Eng. Econ. 2018, 29, 140–150. [Google Scholar] [CrossRef]
  38. Ersen, H.Y.; Tas, O.; Ugurlu, U. Solar energy investment valuation with intuitionistic fuzzy trinomial lattice real option model. IEEE Trans. Eng. Manag. 2023, 70, 2584–2593. [Google Scholar] [CrossRef]
  39. Puri, J.; Yadav, S.P. Intuitionistic fuzzy data envelopment analysis: An application to the banking sector in India. Expert Syst. Appl. 2015, 42, 4982–4998. [Google Scholar] [CrossRef]
  40. Arya, A.; Yadav, S.P. Development of intuitionistic fuzzy data envelopment analysis models and intuitionistic fuzzy input–output targets. Soft Comp. 2019, 23, 8975–8993. [Google Scholar] [CrossRef]
  41. Davoudabadi, R.; Mousavi, S.M.; Mohagheghi, V. A new decision model based on DEA and simulation to evaluate renewable energy projects under interval-valued intuitionistic fuzzy uncertainty. Renew. Energ. 2021, 164, 1588–1601. [Google Scholar] [CrossRef]
  42. Uzhga-Rebrov, O.; Grabusts, P. Methodology for Environmental Risk Analysis Based on Intuitionistic Fuzzy Values. Risks 2023, 11, 88. [Google Scholar] [CrossRef]
  43. Andrés-Sánchez, J. Pricing Life Contingencies Linked to Impaired Life Expectancies Using Intuitionistic Fuzzy Parameters. Risks 2024, 12, 29. [Google Scholar] [CrossRef]
  44. Koissi, M.C.; Shapiro, A.F. Fuzzy formulation of the Lee–Carter model for mortality forecasting. Insur. Math. Econ. 2006, 39, 287–309. [Google Scholar] [CrossRef]
  45. Szymański, A.; Rossa, A. The modified fuzzy mortality model based on the algebra of ordered fuzzy numbers. Biom. J. 2021, 63, 671–689. [Google Scholar] [CrossRef]
  46. Parvathi, R.; Malathi, C.; Akram, M.; Atanassov, K. Intuitionistic fuzzy linear regression analysis. Fuzzy Optim Decis Making 2013, 12, 215–229. [Google Scholar] [CrossRef]
  47. Tanaka, H.; Ishibuchi, H. A possibilistic regression analysis based on linear programming. In Fuzzy Regression Analysis; Kacprzuk, J., Fedrizzi, M., Eds.; Physica-Verlag: Heildelberg, Germany, 1992; pp. 47–60. [Google Scholar]
  48. Lee, H.; Tanaka, H. Upper and lower approximation models in interval regression using regression quantile techniques. Eur. J. Oper. Res. 1999, 116, 653–666. [Google Scholar] [CrossRef]
  49. Savic, D.; Predrycz, W. Fuzzy linear models: Construction and evaluation. In Fuzzy Regression Analysis; Kacprzuk, J., Fedrizzi, M., Eds.; Physica-Verlag: Heildelberg, Germany, 1992; pp. 91–100. [Google Scholar]
  50. Zadeh, L.A. Fuzzy Sets. Inf. Control 1965, 8, 338–353. [Google Scholar] [CrossRef]
  51. Dubois, D.; Prade, H. Fuzzy numbers: An overview. Read. Fuzzy Sets Intell. Syst. 1993, 112–148. [Google Scholar]
  52. Kreinovich, V.; Kosheleva, O.; Shahbazova, S.N. Why Triangular and Trapezoid Membership Functions: A Simple Explanation. In Recent Developments in Fuzzy Logic and Fuzzy Sets. Studies in Fuzziness and Soft Computing; Shahbazova, S., Sugeno, M., Kacprzyk, J., Eds.; Springer: Cologne, Germany, 2020; Volume 391, pp. 25–51. [Google Scholar]
  53. Mauris, G.; Lasserre, V.; Foulloy, L. A fuzzy approach for the expression of uncertainty in measurement. Measurement 2001, 29, 165–177. [Google Scholar] [CrossRef]
  54. Dubois, D.; Folloy, L.; Mauris, G.; Prade, H. Probability–possibility transformations, triangular fuzzy sets, and probabilistic inequalities. Reliab. Comput. 2004, 10, 273–297. [Google Scholar] [CrossRef]
  55. Andrés-Sánchez, J.; Gonzalez-Vila, L.G.V. The valuation of life contingencies: A symmetrical triangular fuzzy approximation. Insur. Math. Econ. 2017, 72, 83–94. [Google Scholar] [CrossRef]
  56. Mauris, G. Possibility distributions: A unified representation of usual direct-probability-based parameter estimation methods. Int. J. Approx. Reason. 2011, 52, 1232–1242. [Google Scholar] [CrossRef]
  57. Couso, I.; Montes, S.; Gil, P. The necessity of the strong α-cuts of a fuzzy set. Int. J. Uncertain. Fuzziness Knowl. Based Syst. 2001, 9, 249–262. [Google Scholar] [CrossRef]
  58. Buckley, J.J. Fuzzy statistics: Regression and prediction. Soft Comput. 2005, 9, 769–775. [Google Scholar] [CrossRef]
  59. Sfiris, D.S.; Papadopoulos, B.K. Non-asymptotic fuzzy estimators based on confidence intervals. Inf. Sci. 2014, 279, 446–459. [Google Scholar] [CrossRef]
  60. Dubois, D.; Prade, H. Practical methods for constructing possibility distributions. Inte. J. Intell. Syst. 2016, 31, 215–239. [Google Scholar] [CrossRef]
  61. Adjenughwure, K.; Papadopoulos, B. Fuzzy-statistical prediction intervals from crisp regression models. Evolv. Syst. 2020, 11, 201–213. [Google Scholar] [CrossRef]
  62. Yuan, X.H.; Li, H.X.; Zhang, C. The theory of intuitionistic fuzzy sets based on the intuitionistic fuzzy special sets. Inf. Sci. 2014, 277, 284–298. [Google Scholar] [CrossRef]
  63. Kumar, P.S.; Hussain, R.J. A method for solving unbalanced intuitionistic fuzzy transportation problems. Notes Intuition. Fuzzy Sets 2015, 21, 54–65. [Google Scholar]
  64. Mukherjee, A.K.; Gazi, K.H.; Salahshour, S.; Ghosh, A.; Mondal, S.P. A brief analysis and interpretation on arithmetic operations of fuzzy numbers. Res. Contr. Optim. 2023, 13, 100312. [Google Scholar] [CrossRef]
  65. Bayeg, S.; Mert, R. On intuitionistic fuzzy version of Zadeh’s extension principle. Notes Intuition. Fuzzy Sets 2021, 27, 9–17. [Google Scholar] [CrossRef]
  66. Buckley, J.J.; Qu, Y. On using α-cuts to evaluate fuzzy equations. Fuzzy Sets Syst. 1990, 38, 309–312. [Google Scholar] [CrossRef]
  67. Grzegorzewski, P.; Pasternak-Winiarska, K. Natural trapezoidal approximations of fuzzy numbers. Fuzzy Sets Syst. 2014, 250, 90–109. [Google Scholar] [CrossRef]
  68. Andrés-Sánchez, J.; González-Vila Puchades, L. Life settlement pricing with fuzzy parameters. Appl. Soft Comput. 2023, 148, 110924. [Google Scholar] [CrossRef]
  69. Chukhrova, N.; Johannssen, A. Fuzzy regression analysis: Systematic review and bibliography. Appl. Soft Comput. 2019, 84, 105708. [Google Scholar] [CrossRef]
  70. Arefi, M.; Taheri, S.M. Least-Squares Regression Based on Atanassov’s Intuitionistic Fuzzy Inputs–Outputs and Atanassov’s Intuitionistic Fuzzy Parameters. IEEE Trans. Fuzzy. Syst. 2014, 23, 1142–1154. [Google Scholar] [CrossRef]
  71. Koenker, R.; Bassett Jr, G. Regression quantiles. Econometrica 1978, 46, 33–50. [Google Scholar] [CrossRef]
  72. Ishibuchi, H.; Nii, M. Fuzzy regression using asymmetric fuzzy coefficients and fuzzified neural networks. Fuzzy Sets Syst. 2001, 119, 273–290. [Google Scholar] [CrossRef]
  73. Chen, F.; Chen, Y.; Zhou, J.; Liu, Y. Optimizing h value for fuzzy linear regression with asymmetric triangular fuzzy coeffi-cients. Eng. Appl. Artif. Intell. 2016, 47, 16–24. [Google Scholar] [CrossRef]
  74. The Faculty and Institute of Actuaries. Claims Reserving Manual, 2nd ed.; The Faculty and Institute of Actuaries: London, UK, 1997. [Google Scholar]
  75. Schmidt, K.D.; Zocher, M. The Bornhuetter-Ferguson. Variance 2008, 2, 85–110. Available online: https://www.casact.org/sites/default/files/2021-07/Bornhuetter-Ferguson-Schmidt-Zocher.pdf (accessed on 9 March 2024).
  76. Muzzioli, S.; Gambarelli, L.; De Baets, B. Option implied moments obtained through fuzzy regression. Fuzzy Optim. Decis. Making 2020, 19, 211–238. [Google Scholar] [CrossRef]
  77. Muzzioli, S.; Ruggieri, A.; De Baets, B. A comparison of fuzzy regression methods for the estimation of the implied volatility smile function. Fuzzy Sets Syst. 2015, 266, 131–143. [Google Scholar] [CrossRef]
  78. Andrés-Sánchez, J. Fuzzy claim reserving in nonlife insurance. Comput. Sci. Inf. Syst. 2014, 11, 825–838. [Google Scholar]
  79. Cummins, D.J.; Derrig, R.A. Fuzzy financial pricing of property-liability insurance. N. Am. Actuar. J. 1997, 1, 21–40. [Google Scholar] [CrossRef]
  80. Mircea, I.; Covrig, M. A discrete time insurance model with reinvested surplus and a fuzzy number interest rate. Procedia Econ. Financ. 2015, 32, 1005–1011. [Google Scholar] [CrossRef]
  81. Ungureanu, D.; Vernic, R. On a fuzzy cash flow model with insurance applications. Decis. Econ. Financ. 2015, 38, 39–54. [Google Scholar] [CrossRef]
Figure 1. Triangular intuitionistic fuzzy numbers.
Figure 1. Triangular intuitionistic fuzzy numbers.
Mathematics 12 00845 g001
Table 1. A revision of contributions to claim reserve modeling with fuzzy mathematics.
Table 1. A revision of contributions to claim reserve modeling with fuzzy mathematics.
Method to Fit Fuzzy ParametersNote ExtensionsTaylor’s Separation MethodTwo-Way Methods
Heuristically[4,14,15,18]------
FR-MFP[3,13][23,24][19,20,25]
FR-FLS[26][21][22]
Note: FR-MFP stands for fuzzy regression with the minimum fuzziness principle, and FR-FLS stands for fuzzy least squares.
Table 2. Run-off triangle of accumulated claims.
Table 2. Run-off triangle of accumulated claims.
Development/Payment Period
i|j01j = nin − 1n
Occurrence/Origin Period0C0,0C0,1C0,jC0,n−1C0,n
1C1,0C1,1C1,jC1,n−1
iCi,0Ci,1Ci,n−i
n − 1Cn−1,0Cn−1,1
nCn,0
Table 3. Run-off triangle of incremental claims.
Table 3. Run-off triangle of incremental claims.
Development/Payment Period
01j = nin − 1n
Occurrence/Origin Period0S0,0S0,1S0,jS0,n−1S0,n
1S1,0S1,1S1,jS1,n−1
iSi,0Si,1Si,ni
n − 1Sn−1,0Sn−1,1
nSn,0
Table 4. Run-off triangle of accumulated claims used in this paper.
Table 4. Run-off triangle of accumulated claims used in this paper.
i|j012345
0100118552423298833353403
111132103277434223844
21265243332333977
3149028733883
417253261
51889
f j 1.8991.3291.2321.1201.020
Source: Faculty and Institute of Actuaries [74].
Table 5. Deterministic loss reserves obtained by using the chain-ladder method.
Table 5. Deterministic loss reserves obtained by using the chain-ladder method.
R O 1 R O 2 R O 3 R O 4 R O 5 R
78.38567.931584.672842.104826.239899.31
Table 6. Run-off triangles of observed incremental claims and theoretical incremental claims with chain-ladder development factors (41).
Table 6. Run-off triangles of observed incremental claims and theoretical incremental claims with chain-ladder development factors (41).
Observed Incremental Claims Theoretical Incremental Claims
i|j012345i|j012345
01001854568565347680957.27861.02598.44561.04357.2468
11113990671648422 11103.37992.43689.78646.66411.76
212651168800744 21278.491149.95799.26749.30
3149013831010 31538.061383.41961.53
417251536 41716.811544.19
51889 51889
Source: Own elaboration from the Faculty and Institute of Actuaries [74].
Table 7. Estimates of reserves with bootstrapping confidence intervals and chain-ladder development factors.
Table 7. Estimates of reserves with bootstrapping confidence intervals and chain-ladder development factors.
α R O 1 α R O 2 α R O 3 α R O 4 α
1[78.38, 78.38][567.93, 567.93][1584.67, 1584.67][2842.10, 2842.10]
0.25[76.94, 79.83][564.14, 571.63][1577.02, 1592.89][2830.92, 2852.36]
0.5[75.08, 81.62][560.41, 576.65][1569.64, 1602.17][2820.37, 2866.04]
0.75[71.78, 87.40][554.77, 583.26][1558.79, 1613.60][2804.07, 2883.77]
0.1[69.09, 91.11][548.46, 589.15][1547.12, 1625.22][2787.26, 2901.24]
0.05[68.06, 92.87][544.98, 593.35][1539.65, 1633.24][2775.83, 2915.09]
0.01[66.74, 94.64][537.02, 602.96][1522.32, 1651.93][2751.60, 2942.85]
0[63.68, 97.77][526.48, 628.62][1507.29, 1677.46][2720.66, 2997.42]
α R O 5 α R α R α
1[4826.23, 4826.23][9899.31, 9899.31][9899.31, 9899.31]
0.25[4808.51, 4838.96][9857.53, 9935.67][9866.55, 9930.80]
0.5[4792.94, 4859.91][9818.44, 9986.38][9834.38, 9972.34]
0.75[4768.92, 4887.95][9758.33, 10,055.97][9786.68, 10,028.49]
0.1[4739.83, 4918.17][9691.76, 10,124.88][9733.84, 10,078.13]
0.05[4721.56, 4937.74][9650.07, 10,172.29][9702.70, 10,107.32]
0.01[4682.60, 4994.47][9560.29, 10,286.87][9644.19, 10,193.22]
0[4641.72, 5105.40][9459.82, 10,506.67][9533.03, 10,481.08]
Note: (a) 1 − α represents the confidence level of the probabilistic confidence interval, which can be interpreted as the α-cut of the equivalent fuzzy number; (b) R α is the overall reserve calculated by summing the confidence intervals i = 1 5 R O i α and R α , which are the confidence intervals of the reserves, by applying bootstrapping.
Table 8. Run-off triangle of individual link ratios, f i , j , i = 0, 1, …, 4; j = 0, 1, …, ni − 1.
Table 8. Run-off triangle of individual link ratios, f i , j , i = 0, 1, …, 4; j = 0, 1, …, ni − 1.
i|j012345
01.8531.3061.2331.1161.020
11.8891.3191.2341.123
21.9231.3291.230
31.9281.352
41.890
5
Table 9. Symmetric triangular intuitionistic fuzzy number estimation of development factors with h = 0.1.
Table 9. Symmetric triangular intuitionistic fuzzy number estimation of development factors with h = 0.1.
Parameters of Intuitionistic Fuzzy Regression (Alternative 1)
j01234
f ( j ) 1.8911.3291.2321.1201.020
r f ( j ) ( 0 ) 0.0380.0230.0020.004---
g 0.1400.1790.4570.500---
r f ( j ) 0.0440.0280.0030.0070.003
r f ( j ) * 0.0490.0310.0040.0090.004
Parameters of Intuitionistic Fuzzy Regression (Alternative 2)
j01234
f ( j ) 1.89951.32911.23211.12001.0204
r f ( j ) ( 0 ) 0.04630.02290.00200.0038---
g 0.00000.00000.41790.4274---
r f ( j ) 0.04630.02290.00350.00670.0035
r f ( j ) * 0.05150.02550.00420.00810.0042
Table 10. α , β -cuts of the intuitionistic fuzzy estimates of overall reserves with the two methodologies proposed in this paper.
Table 10. α , β -cuts of the intuitionistic fuzzy estimates of overall reserves with the two methodologies proposed in this paper.
Alternative 1Alternative 2
α β R _ α R ¯ α R * _ β R * ¯ β R _ α R ¯ α R * _ β R * ¯ β
109868.259868.259868.259868.259899.319899.319899.319899.31
0.750.259694.5210,043.129663.9210,074.169733.8610,065.829708.7210,091.30
0.50.59521.9410,219.169461.1710,281.669569.4610,233.409519.5210,284.69
0.250.759350.5010,396.369260.0010,490.779406.1010,402.049331.7010,479.50
019180.1910,574.729060.3810,701.509243.8010,571.769145.2610,675.72
R ~ T = ( 9868.25 ,   697.21 ,   820.46 ) R ~ T = ( 9899.31 ,   663.93 ,   765.15 )
α β R T _ α R T ¯ α R T * _ β R T * ¯ β R T _ α R T ¯ α R T * _ β R T * ¯ β
109868.259868.259868.259868.259899.319899.319899.319899.31
0.750.259693.9410,042.559663.1310,073.369733.3310,065.299708.0210,090.60
0.50.59519.6410,216.859458.0110,278.489567.3510,231.279516.7310,281.88
0.250.759345.3410,391.159252.9010,483.599401.3610,397.259325.4510,473.17
019171.0410,565.459047.7810,688.719235.3810,563.249134.1610,664.46
α β ε _ α ε ¯ α ε * _ β ε * ¯ β ε _ α ε ¯ α ε * _ β ε * ¯ β
100.00%0.00%0.00%0.00%0.00%0.00%0.00%0.00%
0.750.250.01%0.01%0.01%0.01%0.01%0.01%0.01%0.01%
0.50.50.02%0.02%0.03%0.03%0.02%0.02%0.03%0.03%
0.250.750.06%0.05%0.08%0.07%0.05%0.05%0.07%0.06%
010.10%0.09%0.14%0.12%0.09%0.08%0.12%0.11%
Note: The errors ε * ¯ β , ε ¯ α , and ε * _ β ,   ε * ¯ β are obtained with (27) and (28).
Table 11. Symmetrical triangular intuitionistic fuzzy estimates of reserves.
Table 11. Symmetrical triangular intuitionistic fuzzy estimates of reserves.
Alternative 1Alternative 2
R O ~ 1 T (78.38, 12.28, 15.05)(78.38, 13.34, 16.11)
R O ~ 2 T (567.93, 43.40, 53.90)(567.93, 42.66, 51.62)
R O ~ 3 T (1584.67, 66.38, 82.21)(1584.67, 66.72, 80.69)
R O ~ 3 T (2842.10, 200.87, 236.11)(2842.10, 179.74, 207.03)
R O ~ 4 T (4795.16, 374.28, 433.19)(4826.23, 361.47, 409.70)
R ~ (9868.25, 697.21, 820.46)(9899.31, 663.93, 765.15)
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Andrés-Sánchez, J.D. Calculating Insurance Claim Reserves with an Intuitionistic Fuzzy Chain-Ladder Method. Mathematics 2024, 12, 845. https://doi.org/10.3390/math12060845

AMA Style

Andrés-Sánchez JD. Calculating Insurance Claim Reserves with an Intuitionistic Fuzzy Chain-Ladder Method. Mathematics. 2024; 12(6):845. https://doi.org/10.3390/math12060845

Chicago/Turabian Style

Andrés-Sánchez, Jorge De. 2024. "Calculating Insurance Claim Reserves with an Intuitionistic Fuzzy Chain-Ladder Method" Mathematics 12, no. 6: 845. https://doi.org/10.3390/math12060845

APA Style

Andrés-Sánchez, J. D. (2024). Calculating Insurance Claim Reserves with an Intuitionistic Fuzzy Chain-Ladder Method. Mathematics, 12(6), 845. https://doi.org/10.3390/math12060845

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