Next Article in Journal
An Efficient Framework for Adequacy Evaluation through Extraction of Rare Load Curtailment Events in Composite Power Systems
Next Article in Special Issue
Lorenz Surfaces Based on the Sarmanov–Lee Distribution with Applications to Multidimensional Inequality in Well-Being
Previous Article in Journal
Imaging Noise Suppression: Fourth-Order Partial Differential Equations and Travelling Wave Solutions
Previous Article in Special Issue
Frequency and Severity Dependence in the Collective Risk Model: An Approach Based on Sarmanov Distribution
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Sarmanov Distribution with Beta Marginals: An Application to Motor Insurance Pricing

by
Catalina Bolancé
,
Montserrat Guillen
*,† and
Albert Pitarque
Department Econometrics, Riskcenter-IREA, Universitat de Barcelona, E08034 Barcelona, Spain
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2020, 8(11), 2020; https://doi.org/10.3390/math8112020
Submission received: 26 October 2020 / Accepted: 9 November 2020 / Published: 13 November 2020
(This article belongs to the Special Issue Multivariate Sarmanov Distributions and Applications)

Abstract

:
Background: The Beta distribution is useful for fitting variables that measure a probability or a relative frequency. Methods: We propose a Sarmanov distribution with Beta marginals specified as generalised linear models. We analyse its theoretical properties and its dependence limits. Results: We use a real motor insurance sample of drivers and analyse the percentage of kilometres driven above the posted speed limit and the percentage of kilometres driven at night, together with some additional covariates. We fit a Beta model for the marginals of the bivariate Sarmanov distribution. Conclusions: We find negative dependence in the high quantiles indicating that excess speed and night-time driving are not uniformly correlated.

1. Introduction

We analyse a bivariate model based on the Sarmanov distribution with marginal Beta distributions. These marginals are specified based on a generalized linear model (Beta-GLM) or Beta regression as defined by Ferrari and Cribari-Neto [1]. The objective is to fit data defined in the ( 0 , 1 ) interval.
Many authors have analysed bivariate Beta distributions (see, for example, [2,3,4,5]). However, these distributions pose several difficult challenges: their generalization to higher dimensions and their specification as a generalized linear model are not straightforward. The Sarmanov distribution provides a way to address these challenges.
Originally, the Sarmanov distribution in its bivariate form was introduced by Sarmanov [6], its multivariate version was suggested by Lee [7] and was generalized by Bairamov et al. [8]. Its use to model the bivariate behaviour of random variables with a marginal B e t a ( α , β ) distribution was proposed by Gupta and Wong [3]. These authors defined the five parameter bivariate Beta distribution from what is known as Morgenstern’s distribution [9] with marginal Beta, which is a particular case of the Sarmanov distribution.
The bivariate Sarmanov distribution is characterized by its flexibility in the marginal distributions and, furthermore, given that its functional form establishes that the marginals are clearly separated from the dependency model, the specification in terms of a bivariate generalized linear model turns out to be natural. Generalizing from two dimensions to higher dimensions is simple—(see [10] for an example of a trivariate Sarmanov distribution specified as a generalized linear model with Negative Binomial marginals).
In this work, we show an application of the bivariate Sarmanov distribution with Beta marginals generalised linear model to predict two of the most relevant telematics variables in motor insurance [11]. Telematics variables are obtained from GPS/inertial devices installed in vehicles and they provide an abundant source of information to motor insurers. In our case study, a bivariate model is specified, for the proportion of kilometres driven above the posted speed limit and the proportion of kilometres driven at night. These two variables seem to be related, but researchers have not yet been able to find a good way to understand their association. The explanatory variables are the characteristics of the insured policyholder and the vehicle. The database used in our application has already been analysed in various works published in statistical, transport and risk analysis journals (see [11,12,13,14,15,16,17]). In all previous studies, the two telematics variables that we analyse here were used as predictors of the accident rate, and they were assumed to be uncorrelated.
In Section 2, the new bivariate Sarmanov model is specified and the particular case with marginal Beta-GLM with a domain in the ( 0 , 1 ) interval is analysed; the estimation method is also discussed. The results of our case study are shown in Section 3. Finally, Section 4 contains the conclusions.

2. The Models

Let ( Y 1 , Y 2 ) be a bivariate random vector that, for convenience, is defined in ( 0 , 1 ) 2 . Its distribution depends on a set of k quantitative or binary covariates, whose values are represented by the vector x j = x 1 j , , x k j , j = 1 , 2 , where x 1 j = 1 is a constant term. The relationship between Y j and the covariates is given by the linear predictor x j β j , where β j = β 1 j , , β k j , j = 1 , 2 , are vectors of parameters to be estimated. To simplify the notation, the covariates are assumed to be the same for j = 1 and j = 2 , and so the vector of explanatory variables is denoted as x. The bivariate probability density function (pdf) associated with the Sarmanov distribution is:
f Y 1 , Y 2 ( y 1 , y 2 | x β 1 , x β 2 ) = f 1 ( y 1 | x β 1 ) f 2 ( y 2 | x β 2 ) × 1 + ω ϕ 1 ( y 1 | x β 1 ) ϕ 2 ( y 2 | x β 2 ) , y 1 , y 2 ( 0 , 1 )
where ω is the dependence parameter and ϕ j , j = 1 , 2 , are bounded kernel functions. For the function defined in (1) to be a pdf, the following conditions must hold:
0 1 ϕ j ( y j | x β j ) f j ( y j | x β j ) d y j = 0 , j = 1 , 2
and
1 + ω ϕ 1 ( y 1 | x β 1 ) ϕ 2 ( y 2 | x β 2 ) 0 , ( y 1 , y 2 ) ( 0 , 1 ) 2 .
For given values of x β j , j = 1 , 2 , we define:
m j x β j = inf 0 < y j < 1 ϕ j y j | x β j and M j x β j = sup 0 < y j < 1 ϕ j y j | x β j , j = 1 , 2 .
Taking into account the condition defined in (3), bounds can be defined for the dependency parameter ω . However, as this parameter does not depend on the linear predictor, new extreme values are defined as: m j = max x β j m j ( x β j ) and M j = min x β j M j ( x β j ) , so that the bounds of the dependency parameter are:
max 1 m 1 m 2 , 1 M 1 M 2 ω min 1 m 1 M 2 , 1 M 1 m 2 .
The previous condition holds for every vector of covariates x, which implies that the dependency parameter must be located within the narrowest bounds. In practice, we will assume that the vectors observed in the sample dataset lead to the entire domain of values of linear predictors x β j , j = 1 , 2 . In the insurance context, where we will discuss our illustration, we assume that all possible risk profiles that can be insured by the company are already present in the portfolio.
For each vector of covariate observations x, we can also obtain the covariance between the dependent variables as:
c o v ( Y 1 , Y 2 ) = ω v 1 ( x ) v 2 ( x ) ,
where v j ( x ) = 0 1 y j ϕ j ( y j | x β j ) f j ( y j | x β j ) d y j , j = 1 , 2 . The correlation is obtained by dividing by the product of standard deviations.
There exist many possible specifications for the kernel functions ϕ j , j = 1 , 2 (see [18] [for a description of kernel functions proposed in the literature). When fitting the bivariate Beta distribution without covariates, Gupta and Wong [3] propose a kernel function such as ϕ j = 2 F j 1 , where F j is the cumulative distribution function (cdf). This specification has the advantage that the bounds for the dependency parameter are given by 1 ω 1 for any vector x. However, the previous model does not allow obtaining closed expressions for some magnitudes of interest, such as the conditioned moments. In this work, we propose to use kernels ϕ j = y j r E ( Y j r ) , where r is a value to be determined by the analyst. Next, some results obtained for the particular case of the Sarmanov distribution with marginal B e t a ( α , β ) distribution with r = 1 are analyzed. These cases intuitively correspond to a situation of linear dependency, controlled by the dependence parameter ω .

2.1. The Bivariate Beta GLM Model

The pdf of a random variable Y with B e t a ( α , β ) distribution, with parameters α , β > 0 , is:
f Y ( y ; α , β ) = Γ ( α + β ) Γ ( α ) Γ ( β ) y α 1 ( 1 y ) β 1 = 1 B ( α , β ) y α 1 ( 1 y ) β 1
and its cdf is:
F Y ( y ; α , β ) = B ( y , α , β ) B ( α , β ) ,
where Γ ( · ) and B ( · , · ) are the Gamma and Beta functions, respectively, and B ( y , · , · ) is the incomplete Beta function.
The Beta regression was proposed by Ferrari and Cribari-Neto [1], with the reparametrization μ = α α + β and ψ = α + β , so that:
f ( y ; μ , ψ ) = 1 B μ ψ , ( 1 μ ) ψ y μ ψ 1 ( 1 y ) ( 1 μ ) ψ 1 ,
where E ( Y ) = μ , with 0 < μ < 1 , and V ( Y ) = μ ( 1 μ ) ( 1 + ψ ) , with ψ > 0 , where ψ 1 is the scale parameter. We note that, given the values of μ and ψ , it holds that V ( Y ) < 0.25 . The specification as GLM is defined as (note that we use μ ( x ) to emphasize that μ depends on the linear predictor):
g μ ( x ) = x β ,
where g [ · ] is a link function that can be defined in different ways, in this work, we use the logit link, g μ ( x ) = log μ ( x ) 1 μ ( x ) .
To simplify the notation from now on, we eliminate the linear predictors in the conditioned part. The pdf associated with the bivariate random vector ( Y 1 , Y 2 ) with a Sarmanov distribution and Beta GLM marginals that will be called the Sarmanov-Beta-GLM is ():
f Y 1 , Y 2 ( y 1 , y 2 ) = 1 B μ 1 ( x ) ψ 1 , ( 1 μ 1 ( x ) ) ψ 1 y 1 μ 1 ( x ) ψ 1 1 ( 1 y 1 ) ( 1 μ 1 ( x ) ) ψ 1 1 × 1 B μ 2 ( x ) ψ 2 , ( 1 μ 2 ( x ) ) ψ 2 y 2 μ 2 ( x ) ψ 2 1 ( 1 y 2 ) ( 1 μ 2 ( x ) ) ψ 2 1 × 1 + ω ( y 1 μ 1 ( x ) ) ( y 2 μ 2 ( x ) ) , y 1 , y 2 ( 0 , 1 ) .
For the previous expression to be a pdf, the dependency parameter must be located within the bounds defined in (4), which, for the kernel functions that we propose, are:
max 1 max x β 1 μ 1 ( x ) max x β 2 μ 2 ( x ) , 1 min x β 1 1 μ 1 ( x ) min x β 2 1 μ 2 ( x ) ω min 1 min x β 1 1 μ 1 ( x ) max x β 2 μ 2 ( x ) , 1 max x β 1 μ 1 ( x ) min x β 2 1 μ 2 ( x ) .
The bivariate cdf associated with a Sarmanov-Beta-GLM is obtained directly from the double integral of the bivariate pdf defined in (6):
F Y 1 , Y 2 ( y 1 , y 2 ) = B y 1 , ψ 1 μ 1 ( x ) , ( 1 μ 1 ( x ) ) ψ 1 B ψ 1 μ 1 ( x ) , ( 1 μ 1 ( x ) ) ψ 1 × B y 2 , ψ 2 μ 2 ( x ) , ( 1 μ 2 ( x ) ) ψ 2 B ψ 2 μ 2 ( x ) , ( 1 μ 2 ( x ) ) ψ 2 × 1 + ω B y 1 , ψ 1 μ 1 ( x ) + 1 , ( 1 μ 1 ( x ) ) ψ 1 B ψ 1 μ 1 ( x ) , ( 1 μ 1 ( x ) ) ψ 1 μ 1 ( x ) B y 1 , ψ 1 μ 1 ( x ) , ( 1 μ 1 ( x ) ) ψ 1 B ψ 1 μ 1 ( x ) , ( 1 μ 1 ( x ) ) ψ 1 × B y 2 , ψ 2 μ 2 ( x ) + 1 , ( 1 μ 2 ( x ) ) ψ 2 B ψ 2 μ 2 ( x ) , ( 1 μ 2 ( x ) ) ψ 2 μ 2 ( x ) B y 2 , ψ 2 μ 2 ( x ) , ( 1 μ 2 ( x ) ) ψ 2 B ψ 2 μ 2 ( x ) , ( 1 μ 2 ( x ) ) ψ 2 ,
where y 1 , y 2 ( 0 , 1 ) .
Proposition 1.
The conditioned pdf is:
f Y 1 | Y 2 ( y 1 | Y 2 = y 2 ) = 1 B μ 1 ( x ) ψ 1 , ( 1 μ 1 ( x ) ) ψ 1 ) y 1 μ 1 ( x ) ψ 1 1 ( 1 y 1 ) ( 1 μ 1 ( x ) ) ψ 1 1 × 1 + ω ( y 1 μ 1 ( x ) ) ( y 2 μ 2 ( x ) ) , y 1 , y 2 ( 0 , 1 )
and similarly for f Y 2 | Y 1 ( y 2 | Y 1 = y 1 ) . Integrating the previous expression, the conditional cdf is obtained as
F Y 1 | Y 2 ( y 1 | Y 2 = y 2 ) = F 1 ( y 1 ) × 1 + ω y 2 μ 2 ( x ) 1 μ 1 ( x ) ω y 2 μ 2 ( x ) y 1 ( 1 y 1 ) ψ 1 μ 1 ( x ) f 1 ( y 1 ) , y 1 , y 2 ( 0 , 1 ) .
Proof. 
The conditioned pdf is obtained directly as
f Y 1 | Y 2 ( y 1 | Y 2 = y 2 ) = f Y 1 , Y 2 ( y 1 , y 2 ) f Y 2 ( y 2 ) .
Integrating the result of f Y 1 | Y 2 ( y 1 | Y 2 = y 2 ) in (9), we obtain:
F Y 1 | Y 2 ( y 1 | Y 2 = y 2 ) = 0 y 1 f 1 ( t ) d t + ω y 2 μ 2 ( x ) 0 y 1 f 1 ( t ) t μ 1 ( x ) d t = F 1 ( y 1 ) + ω y 2 μ 2 ( x ) B y 1 , ψ 1 μ 1 ( x ) + 1 , ( 1 μ 1 ( x ) ) ψ 1 B ψ 1 μ 1 ( x ) , ( 1 μ 1 ( x ) ) ψ 1 μ 1 ( x ) F 1 ( y 1 ) .
In addition, since
B y 1 , ψ 1 μ 1 ( x ) + 1 , ( 1 μ 1 ( x ) ) ψ 1 B ψ 1 μ 1 ( x ) , ( 1 μ 1 ( x ) ) ψ 1 = B y 1 , ψ 1 μ 1 ( x ) , ( 1 μ 1 ( x ) ) ψ 1 B ψ 1 μ 1 ( x ) , ( 1 μ 1 ( x ) ) ψ 1 y 1 μ 1 ( x ) ψ 1 ( 1 y 1 ) ( 1 μ 1 ( x ) ) ψ 1 ψ 1 μ 1 ( x ) B ψ 1 μ 1 ( x ) , ( 1 μ 1 ( x ) ) ψ 1 = F 1 ( y 1 ) y 1 ( 1 y 1 ) ψ 1 μ 1 ( x ) f 1 ( y 1 ) ,
then, by substituting the previous expression in (11), then (10) follows directly. □
The conditioned quantile is obtained from the inverse of expression (10), for which a numerical method (such as Newton’s method) can be used.
Proposition 2.
The conditional expectation is:
E Y 1 | Y 2 = y 2 = μ 1 ( x ) + ω ( y 2 μ 2 ( x ) ) V Y 1 | x ,
where V Y 1 | x = μ 1 ( x ) ( 1 μ 1 ( x ) ) ( ψ 1 + 1 ) is the variance, which also depends on the vector of covariates. Similarly, E Y 2 | Y 1 = y 1 can be found.
Proof. 
The conditional expectation is obtained directly by solving the integral:
E Y 1 | Y 2 = y 2 = 0 1 y 1 f Y 1 | Y 2 ( y 1 | Y 2 = y 2 ) d y 1 = 0 1 y 1 f Y 1 ( y 1 ) d y 1 × 1 + ω ( y 1 μ 1 ( x ) ) ( y 2 μ 2 ( x ) ) = 0 1 y 1 f Y 1 ( y 1 ) d y 1 + ω ( y 2 μ 2 ( x ) ) 0 1 y 1 2 f Y 1 ( y 1 ) d y 1 μ 1 ( x ) 0 1 y 1 f Y 1 ( y 1 ) d y 1 = μ 1 ( x ) + ω y 2 μ 2 ( x ) E ( Y 1 2 | x ) μ 1 ( x ) 2 = μ 1 ( x ) + ω y 2 μ 2 ( x ) V Y 1 | x .
Likewise, the corresponding result is obtained for E Y 2 | Y 1 = y 1 . □
Proposition 3.
From (5), the conditional covariance which depends on the vector of covariates x is:
c o v ( Y 1 , Y 2 ) = ω V ( Y 1 ) V ( Y 2 ) = ω μ 1 ( x ) ( 1 μ 1 ( x ) ) ( ψ 1 + 1 ) μ 2 ( x ) ( 1 μ 2 ( x ) ) ( ψ 2 + 1 )
and the correlation is:
c o r r ( Y 1 , Y 2 ) = ω μ 1 ( x ) ( 1 μ 1 ( x ) ) ( ψ 1 + 1 ) μ 2 ( x ) ( 1 μ 2 ( x ) ) ( ψ 2 + 1 ) .
Proof. 
Note that the covariance and the correlation are calculated directly if, in expression (5), we see that:
v j ( x ) = 0 1 y j ϕ j ( y j | x β j ) f j ( y j | x β j ) d y j = 0 1 y j ( y j μ j ( x ) ) f j ( y j | x β j ) d y j = E ( Y j 2 | x ) μ j ( x ) 2 , j = 1 , 2
 □
The dependence parameter of the model proposed in Gupta and Wong [3], which uses kernel functions ϕ j = 2 F j 1 , j = 1 , 2 , is located in the interval 1 ω 1 and is the same for all x. Our proposal bounds the dependence parameter to the narrowest interval among those obtained from all x. However, the advantage of our proposal is that our model allows for obtaining closed expressions for some magnitudes of interest such as bivariate moments (covariance) and conditional moments. In the numerical analysis section, we also compare the correlations estimated from our model and that of Gupta and Wong [3].

2.2. Estimation

In practice, we start from a bivariate sample of n observations. Let us denote the sample information as ( Y i 1 , Y i 2 ) , i = 1 , , n , where for each i we know the values of the covariates X i = ( X i 1 , , X i k ) . Our objective is to estimate the parameter vectors β j , the scale parameters, ψ j and j = 1 , 2 , and the dependency parameter ω , from the maximization of the logarithm of the likelihood function associated with the Sarmanov distribution:
l β 1 , β 2 , ψ 1 , ψ 2 , ω = i = 1 n log f 1 ( Y i 1 | X i β 1 ) + i = 1 n log f 2 ( Y i 2 | X i β 2 ) + i = 1 n log 1 + ω ϕ 1 ( Y i 1 | X i β 1 ) ϕ 2 ( Y i 2 | X i β 2 ) = l 1 ( β 1 , ψ 1 ) + l 2 ( β 2 , ψ 2 ) + l 12 ( ω , β 1 , β 2 , ψ 1 , ψ 2 ) ,
The maximization of (15) cannot be carried out directly without considering that the parametric space is restricted and, in addition, as it was shown in expression (4), the bounds of the dependence parameter are closely related to the parameters of the marginals. Thus, in the maximization process, infeasible solutions will often be reached unless a careful numerical procedure is specifically designed. One way to address these difficulties is to rely on the IFM (Inference from Margin) method that has been widely used in the estimation of copulas see [19] [for a review]. For the estimation of the Sarmanov distribution, the IFM was already used by Bolancé and Vernic [10] for the case of GLM marginals with Negative Binomial distributions.
The IFM method is implemented as follows:
Inicialization. The parameters for the marginals are estimated as:
β ^ 1 ( 0 ) , ψ ^ 1 ( 0 ) = max β 1 , ψ 1 l 1 ( β 1 , ψ 1 )
β ^ 2 ( 0 ) , ψ ^ 2 ( 0 ) = max β 2 , ψ 2 l 2 ( β 2 , ψ 2 ) .
For the initial estimation, function betareg() of betareg R package is used. With these parameters of the marginals, we start the iterative process in the two steps described below.
Step 1. Given the estimated marginal parameters in iteration m 1 and taking into account the limits of the dependence parameter ω defined in (4), with function optim() and the L-BFGS-B method using R, we estimate ω from the maximization of the likelihood function given fixed values of the marginal parameters, which is:
ω ^ ( m ) = max ω l ω | 12 ω | β ^ 1 ( m 1 ) , β ^ 2 ( m 1 ) , ψ ^ 1 ( m 1 ) , ψ ^ 2 ( m 1 ) ,
where l ω | 12 is the likelihood as a function of ω given the estimated parameters for the marginals in iteration m 1 .
Step 2. Given the estimated dependency parameter ω ^ ( m ) in step 1, the marginal parameters are re-estimated in iteration m as:
β ^ 1 ( m ) , ψ ^ 1 ( m ) , β ^ 2 ( m ) , ψ ^ 2 ( m ) = max β 1 , ψ 1 , β 2 , ψ 2 l 12 | ω ( β 1 , ψ 1 , β 2 , ψ 2 | ω ^ ( m ) ) ,
where l 12 | ω is the likelihood as a function of the marginal parameters given the dependence parameter estimated in step 1. The above maximization is also performed with function optim() and the L-BFGS-B method of R.
Steps 1 and 2 described above are repeated until reaching the convergence criterion based on the differences between parameter estimates obtained in two consecutive iterations.
Remark 1.
In the initialization process, if the dependent variables contain zeros or ones, the following correction Y ˜ j = ( Y j * ( n 1 ) + 0.5 ) / n , j = 1 , 2 was proposed by Smithson and Verkuilen [20].
In practice, the algorithm described above is based on the optimization of conditional likelihood functions and not on the likelihood function defined in (15). However, in the last stage, the parameters estimated with the IFM method can be used as initial parameters in the process of maximizing the full likelihood function defined in (15). For this purpose, function optim() and method L-BFGS-B of R are used again.
Remark 2.
To estimate the Sarmanov model proposed by Gupta and Wong [3], it is not necessary to use the two-step process, since the bounds of the dependence parameter do not depend on the parameters of the marginal distributions.

3. Numerical Analysis

We analyse a database corresponding to a car insurance portfolio, in which part of the variables have been measured via a telematic system. The objective of our analysis is to model the joint behaviour of the percentage of kilometres driven above the posted speed limits ( Y 1 ) and percentage of kilometres driven at night ( Y 2 ). It is well known that both variables are related to the risk of having an accident. In Table 1, we show the main descriptive statistics of the dependent variables and the covariates used in the modelling process. For the estimation of the Sarmanov-Beta-GLM, the dependent variables have been transformed as indicated in Remark 1 in Section 2.2. Furthermore, to avoid very low coefficient values due to the scale of some covariates, variables age ( X 1 ), age of driving license ( X 2 ) and age of the vehicle ( X 5 ) have been divided by 10; the vehicle power variable ( X 6 ) is divided by 100 and the total annual distance driven in kilometres ( X 7 ) is divided by 1000. In addition, note that, in this study, we have included a variable denoting the driver’s gender ( X 3 ) and an indicator of private garage ( X 4 ) as covariates.
The last row of Table 1 shows the Pearson correlation between the two dependent variables. This correlation is compared with the corresponding parameter estimate obtained from the Sarmanov model with marginal Beta proposed here and with the one proposed by Gupta and Wong [3], from now on the GW model. With this objective, Table 2 shows the dependence parameters estimated with both models, and the AIC and BIC statistics without including the covariates and including them. Using expression (14) and without covariates, from the dependence parameters ω ^ = 14.883 , it can be deduced that the estimated correlation is 0.0601 , which is within the confidence interval of the Pearson correlation as shown in the last row of Table 1. On the contrary, if we use the five parameter Beta distribution, the (residual) correlation that is obtained from the numerical calculation of expression (5) is practically zero. This means that the association is captured by the bivariate model. Comparing both models, with and without covariates, using the AIC and BIC statistics, the results of Table 2 show that the fit is better for the model proposed here than it is for the GW model.
Table 3 shows the results of our Sarmanov-Beta-GLM using different vectors of covariates. Model I includes all the explanatory variables, among which we have the age ( X 1 ), the age of the driving license ( X 2 ) and the total distance driven annually ( X 7 ), these three variables are associated with driving experience. To analyze the robustness of the results, in Model II, age ( X 1 ) is eliminated, and, in addition, in Model III, the age of a driver’s license ( X 2 ) is also eliminated. The results of Model I show that the effect of age is negative on both Y 1 and Y 2 that the effect of the driver’s license age is positive on Y 1 and negative on Y 2 and the effect of total distance, X 7 , is positive on both dependent variables. By eliminating age ( X 1 ) in Model II, the signs of the parameters associated with X 2 and X 7 are maintained, although the value is smaller in the case of X 2 and remains practically the same for X 7 . After eliminating variables X 1 and X 2 , we see that the effect of the total annual distance driven remains practically the same. If we observe the effects of the rest of covariates, these are practically the same in models I, II, and III. A man driver ( X 3 ) with a powerful vehicle ( X 6 ) would have larger Y 1 and Y 2 than the rest, all other characteristics being the same. However, using parking at night ( X 4 ) has a positive effect on the percentage of speeding distance ( Y 1 ) and a negative effect on the percentage of night-time driving ( Y 2 ); the opposite happens with the age of the vehicle ( X 5 ). The effect of X 5 indicates that, when the vehicle is older, drivers tends to diminish the percent of speed driving, while night-time driving is larger.
To visualize the dependence between Y 1 and Y 2 in different quantiles, the following three examples of insured drivers are graphically analysed:
  • Profile 1 corresponds to a 27-year-old man, who drives about 7000 kilometres per year, with a 7-year-old driving license, with parking, with a vehicle of about 8 years and 100 HP.
  • Profile 2 corresponds to a 20-year-old man, who drives about 4000 kilometres per year, with a 2-year-old driving license, with parking, with a vehicle of about 2 years and 75 HP.
  • Profile 3 corresponds to a 36-year-old man, who drives about 10,000 kilometres per year, with a 15-year-old driving license, without parking, with a vehicle of about 15 years and 200 HP.
Profile 1 represents the average insured individual of the portfolio; Profile 2 is a younger man driver, less experienced than Profile 1 and with a newer and less powerful vehicle; finally, Profile 3 is an older man driver, more experienced than Profile 1 and an older and more powerful vehicle. Figure 1 represents different quantiles of the variable kilometres driven above the speed limit ( Y 1 ) in the y-axis given the values of the percentage of kilometres driven at night ( Y 2 ) for Profile 1 in the x-axis. Quantiles have been obtained from the expression (10). Note that, if the dependence parameter was zero, all the curves would remain constant. The adjusted dependence structure results in the represented conditional quantiles having a negative nonlinear relationship and, furthermore, the curves for the different quantile levels are non-parallel. Figure 1 indicates that, for Profile 1, the higher the percentage of kilometres driven at night ( Y 2 ), the greater the caution in driving and, therefore, the lower the percentage of distance driven above the speed limits ( Y 1 ). The same quantiles at 75 % (plot on the left) and 95 % (plot on the right) confidence levels are represented in Figure 2. These plots show that the curves are non-parallel and that Profile 3 is the most risky, followed by Profiles 1 and 2.

4. Conclusions

We have developed a bivariate model based on the Sarmanov distribution with marginal Beta GLM which has allowed us to model two important variables in modern motor insurance telematics databases. Our model is an alternative to a proposal previously made by Gupta and Wong [3] based on what is known as Morgenstern’s distribution, which is a particular case of the Sarmanov distribution. Our proposal allows for obtaining closed expressions for some magnitudes of interest, such as the bivariate cdf and conditioned moments, covariance and correlation, which are fundamental in risk analysis. We have shown that our Sarmanov-Beta-GLM model presents better fits than previous proposals also based on the Sarmanov distribution.
The results of our case study have shown that, for a specific example, although the dependence parameter is positive, which directly implies that, in the mean, the relationship between the conditioned mean and the values of the variable that conditions is positive, the conditional quantiles show that the relationship between the conditioned quantile, and the value of the conditioning variable may be negative for high quantile levels, a result that is consistent with the expected behaviour of drivers.

Author Contributions

Conceptualization, C.B. and M.G.; methodology, C.B.; software, A.P.; validation, C.B., M.G., and A.P.; formal analysis, A.P.; investigation, A.P.; resources, M.G.; data curation, M.G.; writing—original draft preparation, C.B.; writing—review and editing, M.G. and A.P.; visualization, C.B.; supervision, C.B.; project administration, M.G.; funding acquisition, M.G. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by the Spanish Ministry of Science and Innovation grant PID2019–105986GB-C21, Fundación BBVA Research on Big Data and ICREA Academia.

Acknowledgments

We thank seminar participants and members of the Riskcenter, Universitat de Barcelona.

Conflicts of Interest

The authors declare no conflict of interest.The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript, or in the decision to publish the results.

References

  1. Ferrari, S.; Cribari-Neto, F. Beta regression for modelling rates and proportions. J. Appl. Stat. 2004, 31, 799–815. [Google Scholar] [CrossRef]
  2. Arnold, B.; Ng, H. Flexible bivariate Beta distributions. J. Multivar. Anal. 2011, 102, 1194–1202. [Google Scholar] [CrossRef] [Green Version]
  3. Gupta, A.; Wong, C. On three and five parameter bivariate beta distributions. Metrika 1985, 32, 85–91. [Google Scholar] [CrossRef]
  4. Olkin, I.; Liu, R. A bivariate beta distribution. Stat. Probab. Lett. 2003, 62, 407–412. [Google Scholar] [CrossRef]
  5. Olkin, I.; Trikalinos, T. Constructions for a bivariate beta distribution. Stat. Probab. Lett. 2015, 96, 54–60. [Google Scholar] [CrossRef] [Green Version]
  6. Sarmanov, O. Generalized normal correlation and two-dimensional frechet classes. Doclady Soviet Math. 1966, 168, 596–599. [Google Scholar]
  7. Lee, M. Properties and applications of the sarmanov family of bivariate distributions. Commun. Stat. Theory Methods 1996, 25, 1207–1222. [Google Scholar]
  8. Bairamov, I.; Altinsoy, B.; Kerns, G. On generalized Sarmanov bivariate distributions. TWMS J. Appl. Eng. Math. 2011, 1, 86–97. [Google Scholar]
  9. Morgenstern, D. Einfache beispiele zweidimen-sionaler verteilungen. Mitteilingsblatt Math. Stat. 1956, 8, 234–235. [Google Scholar]
  10. Bolancé, C.; Vernic, R. Multivariate count data generalized linear models: Three approaches based on the sarmanovdistribution. Insur. Math. Econ. 2019, 85, 89–103. [Google Scholar] [CrossRef] [Green Version]
  11. Guillen, M.; Nielsen, J.; Ayuso, M.; Pérez-Marín, A. The use of telematics devices to improve automobile insurance rates. Risk Anal. 2019, 39, 662–672. [Google Scholar] [CrossRef] [PubMed]
  12. Ayuso, M.; Guillen, M.; Nielsen, J. Improving automobile insurance ratemaking using telematics: Incorporating mileage and driver behaviour data. Transportation 2019, 46, 735–752. [Google Scholar] [CrossRef] [Green Version]
  13. Pérez-Marín, A.; Guillen, M. Semi-autonomous vehicles: Usage-based data evidences of what could be expected from eliminating speed limit violations. Accid. Anal. Prev. 2019, 123, 99–106. [Google Scholar] [CrossRef] [PubMed]
  14. Pérez-Marín, A.; Ayuso, M.; Guillen, M. Do young insured drivers slow down after suffering an accident? Transp. Res. Part F Psychol. Behav. 2019, 62, 690–699. [Google Scholar] [CrossRef]
  15. Pérez-Marin, A.; Guillen, M.; Alcañiz, M.; Bermúdez, L. Quantile regression with telematics information to assess the risk of driving above the posted speed limit. Risks 2019, 7, 80. [Google Scholar] [CrossRef] [Green Version]
  16. Pesantez-Narvaez, J.; Guillen, M.; Alcañiz, M. Predicting motor insurance claims using telematics data-xgboost versus logistic regression. Risks 2019, 7, 70. [Google Scholar] [CrossRef] [Green Version]
  17. Sun, S.; Bi, J.; Guillen, M.; Pérez-Marín, A. Assessing driving risk using internet of vehicles data: An analysis based on generalized linear models. Sensors 2020, 20, 2712. [Google Scholar] [CrossRef] [PubMed]
  18. Bahraoui, Z.; Bolancé, C.; Pelican, E.; Vernic, R. On the bivariate Sarmanov distribution and copula. An application on insurance data using truncated marginal distributions. Stat. Oper. Res. Trans. SORT 2015, 39, 209–230. [Google Scholar]
  19. Joe, H.; Xu, J. The estimation method of inference functions for margins for multivariate models. Open Collect. 1996. [Google Scholar] [CrossRef]
  20. Smithson, M.; Verkuilen, J. A better lemon squeezer? maximum-likelihood regression with beta-distributed dependent variables. Psychol. Methods 2006, 11, 54–71. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Quantiles of percentage of kilometres driven over the speed limit ( Y 1 ) in the y-axis for Profile 1 given the values of percentage of kilometres driven at night ( Y 2 ) in the x-axis.
Figure 1. Quantiles of percentage of kilometres driven over the speed limit ( Y 1 ) in the y-axis for Profile 1 given the values of percentage of kilometres driven at night ( Y 2 ) in the x-axis.
Mathematics 08 02020 g001
Figure 2. Quantiles of percentage of kilometres driven over the speed limit ( Y 1 ) for each driver profile given the values of percentage of kilometres driven at night ( Y 2 ), (left) 75 % level and (right) 95 % level.
Figure 2. Quantiles of percentage of kilometres driven over the speed limit ( Y 1 ) for each driver profile given the values of percentage of kilometres driven at night ( Y 2 ), (left) 75 % level and (right) 95 % level.
Mathematics 08 02020 g002
Table 1. Definition of variables and descriptive statistics: mean, standard deviation (STD), minimum (Min) and Maximum (Max). The last row shows the linear correlation between dependent variables and a confidence interval at the 95 % level.
Table 1. Definition of variables and descriptive statistics: mean, standard deviation (STD), minimum (Min) and Maximum (Max). The last row shows the linear correlation between dependent variables and a confidence interval at the 95 % level.
VariableDescriptionMeanSTDMinMax
Y 1 Percentage of kilometres driven above the speed limit0.0630.0680.0000.704
Y 2 Percentage of kilometres driven at night0.0690.0640.0001.000
X 1 Age of the driver27.5653.09419.84936.904
X 2 Age if driver License7.1743.0531.81015.910
X 3 Gender (=1 Men, =0 Women)0.4890.5000.0001.000
X 4 Night parking (=1 yes, 0=no)0.7740.4180.0001.000
X 5 Age of the vehicle8.7494.1741.93820.468
X 6 Power of the vehicle in Horse Power (HP)97.22627.77212.000500.000
X 7 Total Km7159.5104191.7531.59050,035.560
ρ Pearson correlation between dependent variables (CI)0.070 (0.057,0.082)
Table 2. Estimated dependence from Sarmanov-Beta models and goodness of fit criteria.
Table 2. Estimated dependence from Sarmanov-Beta models and goodness of fit criteria.
ω ^ (p-Value)AICBIC
Proposed ModelNo covariates14.883 (<0.001)−171,282.2−171,241.5
With all covariates2.388      (0.055)−177,508.8−177,354.4
GW ModelNo Covariates0.002      (0.346)−171,165.4−171,124.8
With all covariates0.002      (0.356)−177,497.2−177,342.8
Table 3. Parameter estimates (p-values) for the Sarmanov-Beta models and goodness of fit statistics.
Table 3. Parameter estimates (p-values) for the Sarmanov-Beta models and goodness of fit statistics.
Model IModel IIModel III
Y1Y2Y1Y2Y1Y2
Cons.−3.055 (<0.001)−2.556 (<0.001)−3.819 (<0.001)−2.975 (<0.001)−3.796 (<0.001)−3.061 (<0.001)
X 1 −0.339 (<0.001)−0.185 (<0.001)--
X 2 0.294 (<0.001)−0.052      (0.018)0.048      (0.002)−0.187 (<0.001)-
X 3 0.097 (<0.001)0.274 (<0.001)0.107 (<0.001)0.281 (<0.001)0.109 (<0.001)0.274 (<0.001)
X 4 0.108 (<0.001)−0.031      (0.007)0.107 (<0.001)−0.031      (0.007)0.107 (<0.001)−0.031      (0.007)
X 5 −0.043      (0.001)0.055 (<0.001)−0.043      (0.001)0.055 (<0.001)−0.043      (0.001)0.055 (<0.001)
X 6 0.653 (<0.001)0.077 (<0.001)0.654 (<0.001)0.079 (<0.001)0.664 (<0.001)0.038      (0.027)
X 7 0.045 (<0.001)0.035 (<0.001)0.046 (<0.001)0.035 (<0.001)0.046 (<0.001)0.035 (<0.001)
ϕ 1 18.480 (<0.001)18.300 (<0.001)18.294 (<0.001)
ϕ 2 14.823 (<0.001)14.782 (<0.001)14.703 (<0.001)
ω 2.388      (0.055)2.325      (0.059)2.214      (0.060)
AIC−177,508.8−177,238.5−177,113.5
BIC−177,354.4−177,100.3−176,991.6
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Bolancé, C.; Guillen, M.; Pitarque, A. A Sarmanov Distribution with Beta Marginals: An Application to Motor Insurance Pricing. Mathematics 2020, 8, 2020. https://doi.org/10.3390/math8112020

AMA Style

Bolancé C, Guillen M, Pitarque A. A Sarmanov Distribution with Beta Marginals: An Application to Motor Insurance Pricing. Mathematics. 2020; 8(11):2020. https://doi.org/10.3390/math8112020

Chicago/Turabian Style

Bolancé, Catalina, Montserrat Guillen, and Albert Pitarque. 2020. "A Sarmanov Distribution with Beta Marginals: An Application to Motor Insurance Pricing" Mathematics 8, no. 11: 2020. https://doi.org/10.3390/math8112020

APA Style

Bolancé, C., Guillen, M., & Pitarque, A. (2020). A Sarmanov Distribution with Beta Marginals: An Application to Motor Insurance Pricing. Mathematics, 8(11), 2020. https://doi.org/10.3390/math8112020

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