Next Article in Journal
Estimating General Parameters from Non-Probability Surveys Using Propensity Score Adjustment
Next Article in Special Issue
Bivariate Mixed Poisson and Normal Generalised Linear Models with Sarmanov Dependence—An Application to Model Claim Frequency and Optimal Transformed Average Severity
Previous Article in Journal
Heating a 2D Thermoelastic Half-Space Induced by Volumetric Absorption of a Laser Radiation
Previous Article in Special Issue
A Sarmanov Distribution with Beta Marginals: An Application to Motor Insurance Pricing
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Lorenz Surfaces Based on the Sarmanov–Lee Distribution with Applications to Multidimensional Inequality in Well-Being

by
José María Sarabia
1,*,† and
Vanesa Jorda
2,*,†
1
Department of Quantitative Methods, CUNEF University, Leonardo Prieto Castro 2, 28040 Madrid, Spain
2
Department of Economics, University of Cantabria, Avda de los Castros s/n, 39005 Santander, Spain
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2020, 8(11), 2095; https://doi.org/10.3390/math8112095
Submission received: 14 October 2020 / Revised: 16 November 2020 / Accepted: 17 November 2020 / Published: 23 November 2020
(This article belongs to the Special Issue Multivariate Sarmanov Distributions and Applications)

Abstract

:
The purpose of this paper is to derive analytic expressions for the multivariate Lorenz surface for a relevant type of models based on the class of distributions with given marginals described by Sarmanov and Lee. The expression of the bivariate Lorenz surface can be conveniently interpreted as the convex linear combination of products of classical and concentrated univariate Lorenz curves. Thus, the generalized Gini index associated with this surface is expressed as a function of marginal Gini indices and concentration indices. This measure is additively decomposable in two factors, corresponding to inequality within and between variables. We present different parametric models using several marginal distributions including the classical Beta, the GB1, the Gamma, the lognormal distributions and others. We illustrate the use of these models to measure multidimensional inequality using data on two dimensions of well-being, wealth and health, in five developing countries.

1. Introduction

Even though a consensus has slowly emerged among scholars that economic inequality is a multidimensional construct, most of the early work in this area relies almost exclusively on income variables. In the last decades, however, academics’ interest has shifted from income alone to a more comprehensive conception of well-being, which has an intrinsic multidimensional nature (see, e.g., [1,2,3]). Academics have repeatedly argued that the same monetary amount may yield different standards of quality of life for people with different needs [4,5]. Hence, the level of income is a crude proxy of well-being that an individual enjoys.
The renewed interest in well-being inequality has motivated the development of different approaches to measuring this phenomenon. Ideally, the evaluation of multidimensional inequality should consider both inequality within each dimension and the degree of association between dimensions. Multidimensional measures that consider both types of inequalities have been proposed by [6,7,8,9,10,11,12,13]. However, depending on their properties, different measures might reveal completely different results. When this happens, graphical methods can be a powerful tool to explore where the redistributive movements take place at different parts of the distribution. The extension of the univariate Lorenz curve (LC) to higher dimensions is therefore essential to further explore the distributional patterns of well-being. The three existing definitions of the multivatiate LC were proposed by Taguchi [14,15], Arnold [16] and Koshevoy and Mosler [17], who introduced the concepts of Lorenz zonoid and Gini zonoid index.
In this paper, using the definition proposed by Arnold [16], we derive analytic expressions for the bivariate Lorenz surface for different formulations of the underlying bivariate distribution. A preliminary account of this formulation is provided in [18]. In the present paper, we generalize these results and apply them to the study of multidimensional inequality in well-being. We study a relevant type of models based on the class of bivariate distributions with given marginals described by Sarmanov and Lee [19,20]. This model presents several advantages. The expression of the bivariate Lorenz surface can be easily expressed as a convex linear combination of classical and generalized LCs. The generalized Gini index is also available in a closed form and can be decomposed into two factors, corresponding to the equality within and between variables. We explore several models with different types of marginals. We also discuss the extension to higher dimensions.
In the following section, we present preliminary results including the definition of univariate Lorenz and concentration curves and a short review of the three different definitions of bivariate Lorenz surfaces proposed in the literature. In Section 3, we obtain a closed formula of the bivariate Sarmanov–Lee Lorenz surface and its corresponding generalized Gini index. The decomposition of this index into two factors is also derived. Thereafter, we explore different models of bivariate Lorenz surfaces in Section 4. We study the extension to higher dimensions in Section 5. We provide an application to measurement of multidimensional inequality in well-being in Section 6. Finally, some conclusions are sketched in Section 7.

2. Preliminary Results

2.1. Univariate Lorenz and Concentration Curves

We denote the class of univariate distributions functions with positive finite expectations by L and denote by L + the class of all distributions in L with F ( 0 ) = 0 corresponding to non-negative random variables. We use the following definition by Gastwirth [21].
Definition 1.
The Lorenz curve L of a random variable X with cumulative distribution function F L is,
L ( u ; F ) = 0 u F 1 ( y ) d y 0 1 F 1 ( y ) d y = 0 u F 1 ( y ) d y E ( X ) , 0 u 1 ,
where
F 1 ( y ) = sup { x : F ( x ) y } , 0 y < 1 , = sup { x : F ( x ) < 1 } , y = 1 ,
is the right continuous inverse distribution function or quantile function corresponding to F.
Now, we present the concept of concentration curve introduced in [22]. Let g ( x ) be a continuous function of x such that its first derivative exists and g ( x ) 0 . If the mean E F [ g ( X ) ] exits, then one can define
L g ( y ; F ) = 0 x g ( x ) d F ( x ) E F [ g ( X ) ]
where y = g ( x ) and f ( x ) and F ( x ) are, respectively, the probability density function (pdf) and the cumulative distribution function (cdf) of the random variable X. The implicit relation between L g ( g ( x ) ; F ) and F ( x ) will be called the concentration curve of the function g ( X ) . If η g ( x ) and η g * ( x ) denote the elasticities of g ( x ) and g * ( x ) , the concentration curve for the function g ( x ) will lie above (below) the concentration curve for the function g * ( x ) if η g ( x ) is less (greater) than η g * ( x ) for all x 0 . The concentration curve admits the simple implicit representation,
L g ( u ; F ) = 1 E F [ g ( X ) ] 0 u g [ F 1 ( t ) ] d t ,
which will be used in the following sections.

2.2. The Three Altermative Definitions of Bivariate Lorenz Surface

In this section, we discuss briefly the three definitions proposed in the literature for constructing a bivariate Lorenz surface. These definitions were suggested by Taguchi [14,15], Arnold [16] and Koshevoy and Mosler [17]. We consider the bivariate definition, which can be extended to higher dimensions.
First, we present the Taguchi [14] definition of bivariate Lorenz surface, which was the first proposal.
Definition 2.
The Lorenz surface of F 12 is the set of points ( s , t , L ( s , t ) ) in R + 3 defined by,
s = 0 u 0 v f 12 ( x 1 , x 2 ) d x 1 d x 2 , t = 0 u 0 v x 1 f 12 ( x 1 , x 2 ) d x 1 d x 2 E ( X 1 ) ,
L ( s , t ; F 12 ) = 0 u 0 v x 2 f 12 ( x 1 , x 2 ) d x 1 d x 2 E ( X 2 ) .
This definition was investigated by Taguchi [14,15,23]. The definition is not symmetric in ( s , t ) , and its extension to higher dimensions does not look simple.
The following definition was proposed by Arnold [16,24], which is an extension of (1) to higher dimensions. Let X = ( X 1 , X 2 ) be a bivariate random variable with bivariate probability distribution function F 12 on R + 2 having finite second and positive first moments. We denote by F i , i = 1 , 2 the marginal cdf corresponding to X i , i = 1 , 2 , respectively.
Definition 3.
The Lorenz surface of F 12 is the graph of the function,
L ( u 1 , u 2 ; F 12 ) = 0 s 1 0 s 2 x 1 x 2 d F 12 ( x 1 , x 2 ) 0 0 x 1 x 2 d F 12 ( x 1 , x 2 ) ,
where
u 1 = 0 s 1 d F 1 ( x 1 ) , u 2 = 0 s 2 d F 2 ( x 2 ) , 0 u 1 , u 2 1 .
Now, we denote by F a the one-point distribution at a R + 2 , which corresponds to the bivariate egalitarian distribution at a. In consequence, the egalitarian distribution has bivariate Lorenz surface L ( u 1 , u 2 ; F a ) = u 1 u 2 .
The two-attribute Gini index G ( F 12 ) is defined as,
G ( F 12 ) = 4 0 1 0 1 [ u 1 u 2 L ( u 1 , u 2 ; F 12 ) ] d u 1 d u 2
where the egalitarian surface is given by L ( u 1 , u 2 ; F a ) = u 1 u 2 . This index evaluates well-being inequality using a composite index of the form:
W ( x 1 , x 2 ) = x 1 · x 2 .
Given the multiplicative nature of the index, well-being dimensions are neither perfect substitutes nor complements. The marginal rate of substitution is x 2 / x 1 .
Previous definition has not been explored in detail in the literature. We highlight some of its properties:
  • The marginal LCs can be obtained as L ( u 1 ; F 1 ) = lim u 2 L ( u 1 , u 2 ; F 12 ) and L ( u 2 ; F 2 ) = lim u 1 L ( u 1 , u 2 ; F 12 ) .
  • The bivariate Lorenz surface does not depend on changes of scale in the marginals.
  • If F 12 is a product distribution function, then
    L ( u 1 , u 2 ; F 12 ) = L ( u 1 ; F 1 ) L ( u 2 ; F 2 ) ,
    which is just the product of the marginal LCs.
  • In the case of a product distribution, the two-attribute Gini defined in (4) can be written as,
    1 G ( F 12 ) = [ 1 G ( F 1 ) ] [ 1 G ( F 2 ) ] .
The following definition of Lorenz surface was initially proposed by Koshevoy [25] and further results were derived by Koshevoy and Mosler [17,26] and Mosler [27]. Denote by L + 2 de set of all 2-dimensional non-negative random vectors X with finite positive marginal expectations. Let Ψ ( 2 ) denote the class of all measurable functions from R + 2 to [ 0 , 1 ] .
Definition 4.
Let X L + 2 . The Lorenz zonoid L ( X ) of the random vector X = ( X 1 , X 2 ) with distribution F 12 is defined as,
L ( X ) = ψ ( x ) d F 12 ( x ) , x 1 ψ ( x ) E ( X 1 ) d F 12 ( x ) , x 2 ψ ( x ) E ( X 2 ) d F 12 ( x ) : ψ Ψ ( 2 ) , = E [ ψ ( X ) ] , E [ X 1 ψ ( X ) ] E ( X 1 ) , E [ X 2 ψ ( X ) ] E ( X 2 ) : ψ Ψ ( 2 ) .
The Lorenz zonoid is a convex American football-subset of the 3-dimensional unit cube that includes the points ( 0 , 0 , 0 ) and ( 1 , 1 , 1 ) . Extension to higher dimensions is quite direct. If we consider some bivariate parametric families of distributions, for example the bivariate lognormal distribution (see [28] Section 4.8) used for modelling the joint distribution of household size and income, the computation of previous Lorenz zonoid is not simple. The computation of the formulas involved in L ( X ) are not explicit in many situations.

2.3. Bivariate Lorenz Surface Based on Bivariate Beta-Generated Distributions

The Arnold’s Lorenz surface can be evaluated implicitly in some relevant bivariate families of income distributions. We present a broad class of Lorenz surfaces based on the class of bivariate beta-generated distributions introduced in [29]. The class of bivariate beta-generated distributions is defined in terms of the joint density by,
g F ( x 1 , x 2 ) = k F 1 ( x 1 ) a 1 F 2 ( x 2 ) b 1 ( 1 F 1 ( x 1 ) ) b + c 1 ( 1 F 2 ( x 2 ) ) a + c 1 f 1 ( x 1 ) f 2 ( x 2 ) ( 1 F 1 ( x 1 ) F 2 ( x 2 ) ) a + b + c ,
where k = Γ ( a + b + c ) Γ ( a ) Γ ( b ) Γ ( c ) , F i ( x i ) and f i ( x i ) , i = 1 , 2 , are genuiene cdfs and pdfs, respectively, and a , b , c > 0 . The marginal distributions in (6) belong to the class of beta-generated distributions [30] and are given by X 1 B G ( a , c ; F 1 ) and X 2 B G ( b , c ; F 2 ) , which belong to the class L , that is, the support of F i , i = 1 , 2 is [ 0 , ) and the first moment exits for both marginals. According to Arnold and Sarabia [31], we define
G X 1 , X 2 ( 1 , 1 ) ( x 1 , x 2 ; a , b , c ) = 1 E ( X 1 X 2 ) 0 x 1 0 x 2 t 1 t 2 g F ( t 1 , t 2 ) d t 1 d t 2
where g F is defined in (6), provided that E ( X 1 X 2 ) < . Then using Equation 3.5 in Arnold and Sarabia [31] we get the bivariate Lorenz surface,
L ( u 1 , u 2 ; a , b , c ) = G X 1 , X 2 ( 1 , 1 ) ( G X 1 1 ( x 1 ) , G X 2 1 ( x 2 ) ; a , b , c )
Expression (8) includes and generalized several relevant models. For example, if we take F i ( x i ) = 1 1 / ( 1 + x i ) , i = 1 , 2 in (6) we obtain the class of bivariate beta distributions of the second kind defined by
f X 1 , X 2 ( x 1 , x 2 ; α ̲ ) = Γ ( α 1 + α 2 + α 3 ) Γ ( α 1 ) Γ ( α 2 ) Γ ( α 3 ) x 1 α 1 1 x 2 α 2 1 ( 1 + x 1 + x 2 ) α 1 + α 2 + α 3 ,
where x 1 , x 2 > 0 , α ̲ = ( α 1 , α 2 , α 3 ) with α i > 0 , i = 1 , 2 , 3 , which was considered by Arnold and Sarabia [31]. The Lorenz surface is given by,
L ( u 1 , u 2 ; α 1 , α 2 , α 3 ) = F X 1 , X 2 ( F X 1 1 ( u 1 ) , F X 2 1 ( u 2 ) ; α ˜ ̲ )
where α ˜ ̲ = ( α 1 + 1 , α 2 + 1 , α 3 2 ) , with α 3 > 2 and F X 1 , X 2 ( x 1 , x 2 ) is the joint cdf of the bivariate beta distributions of the second kind distribution, which can be found in [32].

2.4. Bivariate Generalized Gini Index

In this section we introduce a bivariate version of the generalized Gini index defined by,
G ( ν ) ( X ) = 1 ν ( ν 1 ) 0 1 ( 1 u ) ν 2 L X ( u ) d u ,
where ν > 1 . If ν = 2 , we obtain the usual Gini index. When ν increases, higher weights are attached to small incomes. The limiting case, when ν goes to infinity, depends solely on the lowest income. This definition would be appropriate if we accepted the judgement, introduced by Rawls, that social welfare should be evaluated from the needs the poorest society member. This index was proposed by Donaldson and Weymark [33], Kakwani [34] and was studied by Yitzhaki [35].
Let ( X 1 , X 2 ) a bivariate random variable with E ( X i ) < , i = 1 , 2 , E ( X 1 X 2 ) < and bivariate Lorenz surface L 12 ( u 1 , u 2 ) . The bivariate generalized Gini index is defined by,
G 12 ( ν 1 , ν 2 ) ( X 1 , X 2 ) = 1 k ν 1 k ν 2 0 1 0 1 ( 1 u 1 ) ν 1 2 ( 1 u 2 ) ν 2 2 L 12 ( u 1 , u 2 ) d u 1 d u 2 ,
where k ν = ν ( ν 1 ) and ν 1 , ν 2 > 1 . We write G ¯ 12 ( ν 1 , ν 2 ) ( X 1 , X 2 ) = 1 G 12 ( ν 1 , ν 2 ) ( X 1 , X 2 ) . For the special case ν 1 = ν = 2 , we obtain the bivariate Gini index defined in (4). If X 1 and X 2 are independent random variables we have,
G ¯ 12 ( ν 1 , ν 2 ) ( X 1 , X 2 ) = G ¯ 1 ( ν 1 ) ( X 1 ) G ¯ 2 ( ν 2 ) ( X 2 ) .

3. The Bivariate Sarmanov–Lee Lorenz Surface

The bivariate Lorenz surface (3) admits the following simple explicit representation,
L ( u 1 , u 2 ; F 12 ) = 1 E ( X 1 X 2 ) 0 u 1 0 u 2 A ( x 1 , x 2 ) d x 1 d x 2 , 0 u 1 , u 2 1 ,
where
A ( x 1 , x 2 ) = F 1 1 ( x 1 ) F 2 1 ( x 2 ) f 12 ( F 1 1 ( x 1 ) , F 2 1 ( x 2 ) ) f 1 ( F 1 1 ( x 1 ) ) f 2 ( F 2 1 ( x 2 ) ) .
Now, we continue with the so-called bivariate Sarmanov–Lee Lorenz surface and as a previous step, we present the bivariate Sarmanov–Lee distribution.

3.1. The Bivariate Sarmanov–Lee Distribution

Let X = ( X 1 , X 2 ) be a bivariate Sarmanov–Lee (SL) distribution with joint pdf,
f ( x 1 , x 2 ) = f 1 ( x 1 ) f 2 ( x 2 ) 1 + w ϕ 1 ( x 1 ) ϕ 2 ( x 2 ) ,
where f 1 ( x ) and f 2 ( y ) are the univariate pdf marginals, ϕ i ( t ) , i = 1 , 2 are bounded nonconstant functions such that,
ϕ i ( t ) f i ( t ) d t = 0 , i = 1 , 2 ,
and w is a real number which satisfies the condition 1 + w ϕ 1 ( x 1 ) ϕ 2 ( x 2 ) 0 for all x 1 and x 2 . We denote μ i = E ( X i ) = t f i ( t ) d t , i = 1 , 2 , σ i 2 = v a r ( X i ) = ( t μ i ) 2 f i ( t ) d t , i = 1 , 2 and ν i = E [ X i ϕ i ( X i ) ] = t ϕ i ( t ) f i ( t ) d t , i = 1 , 2 . Properties of this family has been explored by Lee [19]. Moments and regressions of this family can be easily obtained. The product moment is E ( X 1 X 2 ) = μ 1 μ 2 + w ν 1 ν 2 , and the regression of X 2 on X 1 is given by
E ( X 2 | X 1 = x 1 ) = μ 2 + w ν 2 ϕ 1 ( x 1 ) .
The copula associated to (13) is,
C ( u 1 , u 2 ; w , ϕ ) = u 1 u 2 + 0 u 1 0 u 2 ϕ ˜ 1 ( s 1 ) ϕ ˜ 2 ( s 2 ) d s 1 d s 2 ,
where ϕ ˜ i ( s i ) = ϕ i ( F i 1 ( s i ) ) , i = 1 , 2 , and F i ( x i ) are the cdf of X . The pdf of the copula associated to (13) is,
c ( u 1 , u 2 ; w , ϕ ) = C ( u 1 , u 2 ; w , ϕ ) u 1 u 2 = 1 + w ϕ 1 ( F 1 1 ( u 1 ) ) ϕ 2 ( F 2 1 ( u 2 ) ) .
Note that (13) and its associated copula has two components: a first component corresponding to the marginal distributions and the second component defining the structure of dependence, given by the parameter w and the functions ϕ i ( u ) , i = 1 , 2 . These two components will be translated to the structure of the associated bivariate Lorenz surface, and the corresponding bivariate Gini index.
In comparison with other copulas, the Sarmanov–Lee copula has several advantages. Its joint pdf and cdf are quite simple and the covariance structure in general is not limited. On the other hand, many of its different probabilistic features can be obtained in a explicit form, which is very useful from a computational point of view. Finally, the SL distribution includes several relevant special cases including the classical Farlie–Gumbel–Morgenstern (FGM) distribution, and the variations proposed in [36,37].

3.2. The Bivariate Sl Lorenz Surface

The bivariate SL Lorenz surface [18] is obtained using (13) in Equation (3).
Theorem 1.
Let X = ( X 1 , X 2 ) a bivariate Sarmanov–Lee distribution with joint pdf (13) with non-negative marginals satisfying E [ X 1 ] < , E [ X 2 ] < and E [ X 1 X 2 ] < . Then, the bivariate Lorenz surfece is given by,
L S L ( u 1 , u 2 ; F 12 ) = π L ( u 1 ; F 1 ) L ( u 2 ; F 2 ) + ( 1 π ) L g 1 ( u 1 ; F 1 ) L g 2 ( u 2 ; F 2 ) ,
where
π = μ 1 μ 2 E ( X 1 X 2 ) = μ 1 μ 2 μ 1 μ 2 + w ν 1 ν 2 ,
and L ( u i ; F i ) , i = 1 , 2 are the LCs of the marginal distribution X i , i = 1 , 2 respectively, and L g i ( u i ; F i ) , i = 1 , 2 represent the concentration curves of the random variables g i ( X i ) = X i ϕ i ( X i ) , i = 1 , 2 , respectively.
The interpretation of (14) is the following: the bivariate Lorenz surface can be written as a convex linear combination of two components: a first component corresponding to the product of the marginal LCs, that is the marginal components, and a second component corresponding to the product of the concentration LCs, which corresponds with the structure dependence component.

3.3. Bivariate Generalized Gini Index

The following result provides a convenient write of the two-attribute bivariate generalized Gini defined in (4). This expression generalized previous result by Sarabia and Jordá [18] and permits a simple decomposition of the equality in two factors: a first factor which represent the equality within variables and a second factor corresponding to the equality between variables.
Theorem 2.
Let ( X 1 , X 2 ) have a bivariate Sarmanov–Lee distribution with bivariate Lorenz surface L 12 ( u 1 , u 2 ) . The two-attribute bivariate generalized Gini index defined in (10) admits the following representation,
G ¯ 12 ( ν 1 , ν 2 ) ( X 1 , X 2 ) = π G ¯ 1 ( ν 1 ) ( X 1 ) G ¯ 2 ( ν 2 ) ( X 2 ) + ( 1 π ) G ¯ g 1 ( ν 1 ) ( X 1 ) G ¯ g 2 ( ν 2 ) ( X 2 ) ,
where G i ( ν i ) ( X i ) , i = 1 , 2 are the generalized Gini indices of the marginal LCs, and G g i ( ν i ) ( X i ) , i = 1 , 2 represent the concentration indices of the concentration if the curves L ( u i , g i ) , i = 1 , 2 in the sense of Definition 2.
Proof. 
The proof is direct using expression (14) and taking into account the definition of bivariate generalized Gini index given in (10). □
Then the overall equality OE ( ν 1 , ν 2 ) can be decomposed in two factors,
OE ( ν 1 , ν 2 ) = EW ( ν 1 , ν 2 ) + EB ( ν 1 , ν 2 ) ,
where
OE ( ν 1 , ν 2 ) = G ¯ 12 ( ν 1 , ν 2 ) ( X 1 , X 2 ) , EW ( ν 1 , ν 2 ) = π G ¯ 1 ( ν 1 ) ( X 1 ) G ¯ 2 ( ν 2 ) ( X 2 ) , EB ( ν 1 , ν 2 ) = ( 1 π ) G ¯ g 1 ( ν 1 ) ( X 1 ) G ¯ g 2 ( ν 2 ) ( X 2 ) ,
and EW ( ν 1 , ν 2 ) represents the equality within variables and the second factor EB ( ν 1 , ν 2 ) represent the equality between variables. Note that the decomposition the structure of dependence of the underlying bivariate income distribution through the functions g i , i = 1 , 2 .

4. Bivariate Lorenz Surface Models

In this section we consider several relevant models based on the previous methodology.

4.1. Bivariate Power Lorenz Surfaces Based on the Fgm Family

Let X = ( X 1 , X 2 ) be a bivariate FGM distribution with Power marginals and joint pdf,
f ( x 1 , x 2 ; a , w ) = a 1 a 2 x 1 a 1 1 x 2 a 2 1 1 + w ( 1 2 x 1 a 1 ) ( 1 2 x 2 a 2 ) , 0 x 1 , x 2 1 ,
with a i > 0 , i = 1 , 2 and 1 w 1 . Note that (17) belongs to the Sarmanov–Lee family with ϕ i ( x i ) = 1 2 x i a i , i = 1 , 2 and marginal distributions F i ( x i ) = x i a i , i = 1 , 2 . The product moment is,
E ( X 1 X 2 ) = a 1 a 2 [ ( 1 + 2 a 1 ) ( 1 + 2 a 2 ) + w ] ( 1 + a 1 ) ( 1 + a 2 ) ( 1 + 2 a 1 ) ( 1 + 2 a 2 ) ,
and the marginal LC, L ( u i ; F i ) = u i 1 / a i + 1 , 0 u i 1 , i = 1 , 2 . After computing the concentration curves and using (14) we obtain,
L ( u 1 , u 2 ; a , w ) = ( 1 + 2 a 1 ) ( 1 + 2 a 2 ) ( 1 + 2 a 1 ) ( 1 + 2 a 2 ) + w u 1 1 / a 1 + 1 u 2 1 / a 2 + 1 H ( u 1 , u 2 ) ,
where
H ( u 1 , u 2 ) = 1 + w [ 1 + 2 a 1 ( u 1 1 ) + 2 u 1 ] [ 1 + 2 a 2 ( u 2 1 ) + 2 u 2 ] ( 1 + 2 a 1 ) ( 1 + 2 a 2 ) ,
with 0 u 1 , u 2 1 . The function H ( u , v ) represents the dependence structure of the bivariate Lorenz surface. Note that (18) can also be written of the form,
L ( u 1 , u 2 ; a , w ) = u 1 1 / a 1 + 1 u 2 1 / a 2 + 1 ( A + B u 1 + C u 2 + D u 1 u 2 ) ( 1 + 2 a 1 ) ( 1 + 2 a 2 ) + w , 0 u 1 , u 2 1 ,
where A , B , C , D are constants which depend on a 1 , a 2 and w. Using Equation (4), the bivariate Gini index is given by
G ( a 1 , a 2 , w ) = ( 1 + 3 a 1 ) ( 1 + 3 a 2 ) ( 1 + 2 a 1 + 2 a 2 ) + w h ( a 1 , a 2 ) ( 1 + 3 a 1 ) ( 1 + 3 a 2 ) [ ( 1 + 2 a 1 + 2 a 2 ) + w ] ,
where
h ( a 1 , a 2 ) = 1 + 3 ( a 1 + a 2 ) + [ 13 4 ( a 1 + a 2 ) + 4 a 1 a 2 ] a 1 a 2 .
If we set a 1 = a 2 = 1 we get the bivariate Lorenz surface [38],
L 12 ( u 1 , u 2 , w ) = u 1 2 u 2 2 [ 9 + w ( 4 u 1 3 ) ( 4 u 2 3 ) ] 9 + w , 0 u 1 , u 2 1 .
For this last surface, the bivariate generalized Gini index (10) is given by,
G ¯ ( ν 1 , ν 2 ) ( X 1 , X 2 ) = k ν 1 k ν 2 36 [ w ( ν 1 1 ) ( ν 2 1 ) + ( 3 + ν 1 ) ( 3 + ν 2 ) ] Γ ( ν 1 ) Γ ( ν 2 ) ( 9 + w ) Γ ( 4 + ν 1 ) Γ ( 4 + ν 2 ) ,
where k ν = ν ( ν 1 ) , with ν > 1 .

4.2. Bivariate Sarmanov–Lee Lorenz Surfaces with Beta and Gb1 Marginals

Let X i B e ( a i , b i ) , i = 1 , 2 be two classical beta distributions with pdf,
f i ( x i ; a i , b i ) = x i a i 1 ( 1 x i ) b i 1 B ( a i , b i ) , 0 x i 1 , i = 1 , 2
where B ( a i , b i ) = Γ ( a i ) Γ ( b i ) / Γ ( a i + b i ) , for i = 1 , 2 . This distribution has been proposed as a model of income distribution by McDonald [39], and is commonly used in the literature. If we take the mixing functions ϕ i ( x i ) = x i μ i , where μ i = E ( X i ) = a i / ( a i + b i ) , i = 1 , 2 , the bivariate SL distribution is,
f 12 ( x 1 , x 2 ) = f 1 ( x 1 ; a 1 , b 1 ) f 2 ( x 2 ; a 2 , b 2 ) 1 + w x 1 a 1 a 1 + b 1 x 2 a 2 a 2 + b 2 ,
where w satisfies,
( a 1 + b 1 ) ( a 2 + b 2 ) max { a 1 a 2 , b 1 b 2 } w ( a 1 + b 1 ) ( a 2 + b 2 ) max { a 1 b 2 , a 2 b 1 } .
A relevant property of this family is that it can be expressed as a linear combination of products of univariate beta densities.
The univariate LC of the univariate classical beta distribution if given by [40],
L ( u i ; F i ) = G B e ( a i + 1 , b i ) [ G B e ( a i , b i ) 1 ( u i ) ] , i = 1 , 2 ,
where G B e ( a , b ) ( z ) = ( 1 / B ( a , b ) ) 0 z t a 1 ( 1 t ) b 1 d t represents the CDF of a classical Beta distribution. In consequence, the concentration curve can be written as,
L g i ( u i ; F i ) = E ( X i 2 ) G B e ( a i + 2 , b i ) ( G B e ( a i , b i ) 1 ( u i ) ) E ( X i ) 2 L ( u i ; F i ) v a r ( X i ) , i = 1 , 2 .
Note that ν i = E F i [ X i ϕ i ( X i ) ] = v a r ( X i ) , i = 1 , 2 . Finally, combining (20) with (21) in (14), we obtain the bivariate beta Lorenz surface.
This model can be extended easily to the SL distribution with generalized beta of the first type (GB1) marginals, with pdf,
f ( x i ; a i , b i , p i ) = p i x i a i p i 1 ( 1 x i p i ) b i 1 B ( a i , b i ) , 0 x i 1 , i = 1 , 2 ,
and mixing function,
ϕ i ( x i ) = x i μ i = x i Γ ( a i + 1 / p i ) Γ ( a i + b i ) Γ ( a i + b i + 1 / p i ) Γ ( a i ) , i = 1 , 2 .

4.3. Bivariate Sl Lorenz Surfaces with Gamma Marginals

Let X i G ( α i , λ i ) , i = 1 , 2 be classical gamma distributions with pdf f i ( x i ) = λ i α i x i α i 1 e λ i x i / Γ ( α i ) , where x i > 0 and α i , λ i > 0 , i = 1 , 2 . If we consider the mixing function ϕ i ( x i ) = e x i L ( 1 ; α i , λ i ) , where L ( s , α , λ ) = ( 1 + s / λ ) α is the Laplace transform of a gamma distributions [19], we can construct the SL distribution with gamma marginals defined as,
f ( x 1 , x 2 ) = f 1 ( x 1 ) f 2 ( x 2 ) 1 + w e x 1 L ( 1 ; α 1 , λ 1 ) e x 2 L ( 1 ; α 2 , λ 2 ) ,
with x 1 , x 2 > 0 . Note that (22) can be written as a linear combination of products of univariate gamma distributions. Then, if we define,
g i ( x i ) = x i e x i x i L ( 1 ; α i , λ i ) , i = 1 , 2 ,
the bivariate SL Lorenz surfaces with gamma marginals is,
L G ( u 1 , u 2 ; F 12 ) = π L ( u 1 ; F 1 ) L ( u 2 ; F 2 ) + ( 1 π ) L g 1 ( u 1 ; F 1 ) L g 2 ( u 2 ; F 2 ) ,
where
L ( u i ; F i ) = H α i + 1 [ H α i 1 ( u i ) ] , i = 1 , 2 ,
and
L g i ( u i ; F i ) = ( 1 + λ i ) 1 H α i + 1 [ ( 1 + 1 λ i ) H α i 1 ( u i ) ] λ i 1 H α i + 1 [ H α i 1 ( u i ) ] ( 1 + λ i ) 1 λ i 1 , i = 1 , 2 ,
being H α ( x ) = 0 x t α 1 e t d t / Γ ( α ) the cdf of a classical gamma distribution.

4.4. Bivariate Sl Lorenz Surfaces with Lognormal Marginals

The lognormal distribution is one of the most relevant distributions in the analysis of income and wealth data. There are several classes of multivariate distributions with lognormal marginals. The classical multivariate lognormal distribution is defined in terms of a monotone marginal transformation of the classical multivariate normal distribution. Other possible alternatives have been proposed by Sarabia et al. [41] based on conditional specification. In both previous models, we can obtain a bivariate Lorenz surface. If we consider X i LN ( μ X i , σ X i 2 ) , i = 1 , 2 two lognormal distributions, we define a SL distribution, where the mixing functions are given by,
ϕ i ( x i ) = f i ( x i ) 0 f i 2 ( x i ) d x i = f i ( x i ) 1 2 σ X i π exp μ X i + σ X i 2 4 , i = 1 , 2 .
If we define g i ( x i ) = x i ϕ i ( x i ) = x i f i ( x i ) x i 0 f i 2 ( x i ) d x i , i = 1 , 2 we can use (14) for obtaining the bivariate Lorenz surface.

4.5. Other Classes of Bivariate Lorenz Surfaces

The class of Lorenz surfaces based on the classical FGM bivariate distribution with Pareto marginals has been considered in [18]. Other alternative families of bivariate Lorenz surfaces can be constructed, including models with marginals specified in terms of univariate LCs [42] and models based on mixture of distributions [43], which permits to incorporate heterogeneity factors in the inequality analysis.
Income variables typically behave differently at the tails of the distribution [44]. In particular, while the generalised beta distribution of the second kind seems to be a suitable model for the bulk of the income distribution [45], incomes above the 95th percentile seem to be better represented by a Pareto distribution [46]. In this regard, the models proposed by Arnold and Sarabia [38], based on mixtures of Lorenz surfaces and bivariate mixture distributions, might be a suitable alternative to model multivariate income distributions.

5. Extensions to Higher Dimensions and Stochastic Dominance

In this section we discuss briefly how to extend the results in previous sections to dimensions higher than two. As well, we include some comments about dominance in higher dimensions.
First, we consider the general definition of Lorenz surface (3). Let X = ( X 1 , , X m ) be a random vector in L + m with joint CDF F 12 m ( x 1 , , x m ) . The multivariate Arnold’s Lorenz surface can be defined as,
L ( u ; F 12 m ) = 0 s 1 0 s m i = 1 m x i d F 12 m ( x 1 , , x m ) E i = 1 m X i ,
where u i = 0 s i d F i ( x i ) , i = 1 , 2 , , m and 0 u 1 , , u m 1 .
The m-dimensional version of the Sarmanov–Lee distribution is defined as,
f ( x 1 , , x m ) = i = 1 m f i ( x i ) 1 + R ϕ 1 ϕ m Ω m ( x 1 , , x m ) ,
where
R ϕ 1 ϕ m Ω m ( x 1 , , x m ) = 1 i 1 i 2 m w i 1 i 2 ϕ i 1 ( x i 1 ) ϕ i 2 ( x i 2 ) = 1 i 1 i 2 i 3 m w i 1 i 2 i 3 ϕ i 1 ( x i 1 ) ϕ i 2 ( x i 2 ) ϕ i 3 ( x i 3 ) + + w 12 m i = 1 m ϕ i ( x i ) ,
and ϕ i ( x i ) , i = 1 , , k are bounded nonconstant functions such as ϕ i ( t ) f i ( t ) d t = 0 , i = 1 , , k , Ω m = { w i 1 i 2 , w i 1 i 2 i 3 , , w 12 m } , w i 1 i 2 , w i 1 i 2 i 3 , are real numbers and k 2 and m 3 . The set of real numbers Ω m is such that 1 + R ϕ 1 ϕ m Ω m ( x 1 , , x m ) 0 , ( x 1 , , x m ) R m .
Using expression (23), the m-variate Lorenz surface takes de form,
L ( u ; F 12 m ) = w 0 i = 1 m L ( u i ; F i ) + i 1 , , i k w i 1 i k j = i 1 i k L g j ( u j ; F j ) j = i k + 1 i m L ( u j ; F j ) + w 12 m i = 1 m L g i ( u i ; F i ) ,
where g i ( x i ) = x i ϕ i ( x i ) , i = 1 , 2 , , m .
We conclude this section with some references on multidimensional inequality based on zonoid inclusions. Koshevoy [25] and Koshevoy and Mosler [17,26] have studied the partial order of multivariate distributions of transferable attributes induced by the inclusion of either lift zonoids or Lorenz zonoids, defined in (5). In this sense, Koshevoy and Mosler [26] provide the normative rationale for these orderings, using the concept of price majorization (PM) originally proposed by Kolm [9]. See also Andreoli and Zoli [47,48].
Making use of PM order, the distribution of attributes X 1 is always preferred to X 2 by an inequality averse policy maker that is concerned in maximizing the welfare of a distribution of budgets obtained by aggregating the attributes with a given set of prices (see [48]). Mosler [27] provides an interesting result by considering the normative foundations using the extended lift zonoid inclusion order as a robust multidimensional inequality criterion.
From a normative point of view, the interest of the concept of positive PM is related to stochastic orderings that are monotone in attributes and sensitive to the effects of correlation increasing transformations (see Tsui [12,13]). Andreoli and Zoli [48] conceive PM with non-negatives prices as an extension of some multivariate stochastic dominance definition.

6. Application: Multidimensional Inequality in Well-Being

6.1. Data and Estimation Methods

To illustrate the implementation of multidimensional Lorenz surfaces and their associated Gini indices, we collected data from the Demographic and Health Surveys (DHS) of five sub-Saharan African countries: Angola (2015), Uganda (2011), Ethiopia (2011), Haiti (2012) and Nigeria (2012). The DHS are nationally representative household surveys, regularly conducted in more than 85 low- and middle-income countries since 1984. Standard DHS are the largest surveys, designed to cover the national population with sample sizes of between 5000 and 30,000 households.
The DHS were designed to expand the data collected by the World Fertility Surveys and Contraceptive Prevalence Surveys. Hence, these surveys provide rich data on fertility, family planning, health, and nutrition of women and children in developing nations. Although the Household and the Individual modules in the DHS cover the whole population, the data used in this illustration is retrieved from the Household Members Recode. In this module, the age range of the samples are limited to women aged between 15 and 49, but for males the maximum age is not the same in all surveys. Three of the countries considered in this application (Ethiopia, Haiti and Nigeria) include males aged 15–59. Surveys in Angola and Uganda include only men up to the age of 54. The limitation in the age range might introduce a source of bias in our estimates because older population tend to have lower income and worse health, although potentially higher wealth. We are well aware that different age ranges for males might weaken the comparability of the results across countries. In this case, however, this is not a serious issue given the illustrative purpose of the application.
For this illustration, we assume that the level of well-being of the individual i depends solely on two dimensions: wealth ( Y 1 ) and health ( Y 2 ). These two dimensions are included in the set of most prominent aspects of quality of life suggested by Stiglitz et al. [49]. Wealth is captured through a composite index of asset ownership. Compare with income and consumption, household assets seem to provide a better approximation of the economic status of the household. The accurate measurement of household income is extremely difficult due to undereporting and undercoverage issues. This bias is more prevalent in developing countries with large informal sectors. Expenditures are a more reliable measure of household economic status. However, the assessment of consumption expenditures also presents some drawbacks related to the type of expenditures included and the period of time that should be covered. The DHS consider wealth as an underlying unobserved variable, which can be approximated as a linear combination of 30 household assets and utility services associated with a household’s economic position. The weight assigned to each variable is defined using principal component analysis. A potential limitation of this index is the no consideration of financial wealth, which, as a result, underestimates the dispersion of wealth distributions. See [50] for a thorough review on the limitations of the wealth index and their potential solutions.
Health variables in the DHS fit into two categories: those collected objectively and those based on self- reporting. Standard DHS collect data on several nutrition-related biomarkers since the late 90 s. Among them, the hemoglobin level is available in most surveys. Low hemoglobin levels seem to be associated with more physician visits and low body mass index [51]. The prior literature also suggests that there is a close link between low hemoglobin level, malnutrition, and limited access to medical care [52]. Hence, we believe this indicator can be a reliable proxy for the health condition of the population.
We move now to the estimation of the theoretical bivariate distribution of wealth and health. The first step is to normalise wealth and health variables to be bounded between 0 and 1. To do so, we follow the standard procedure used by the United Nations Development Program (UNDP) in the construction of the Human Development Index. The following formula is used to normalise the wealth index and the hemoglobin level:
X i = Y i min ( Y i ) max ( Y i ) min ( Y i ) , i = 1 , 2 .
The normalised variables are used to estimate the bidimensional distribution of well-being and the level of inequality. Since the variables have support on [ 0 , 1 ] , the classical beta distribution is a suitable model for the marginal distributions. The dependence between dimensions is modelled using the Sarmanov–Lee distribution, studied in Section 4.2.
Let X = ( X 1 , X 2 ) be a bivariate distribution with joint pdf (19) and let ( x 11 , x 21 ) , , ( x 1 n , x 2 n ) a sample of n individuals from (19). For the estimation of the parameters ( a 1 , b 1 , a 2 , b 2 , w ) we proceed in two steps:
  • Moments estimation of the marginal distributions. We define,
    m i = 1 n j = 1 n x i j , s i 2 = 1 n j = 1 n ( x i j m i ) 2 , i = 1 , 2
    so, point estimates of the couples ( a i , b i ) , i = 1 , 2 are,
    a ^ i = m i ( m i m i 2 s i 2 ) s i 2 ,
    b ^ i = ( 1 m i ) ( m i m i 2 s i 2 ) s i 2 .
  • Moment estimation of the dependence parameter. The estimate of w is based on the simple relation ρ = w σ 1 σ 2 . Then, if r denotes the sample linear correlation coefficient, and s i , i = 1 , 2 the sample standard deviation of the marginal distributions X i , i = 1 , 2 , the point estimate of w is,
    w ^ = r s 1 · s 2 .
Standard errors of ( a ^ i , b ^ i ) , i = 1 , 2 and w can be easily computed using parametric bootstrap. The Gini index of the marginals distributions is given by,
G ( B e ( a i , b i ) ) = Γ ( a i + b i ) Γ ( a i + 1 2 ) Γ ( b i + 1 2 ) Γ ( a i + b i + 1 2 ) Γ ( a i + 1 ) Γ ( b i ) π , i = 1 , 2 .

6.2. Results

In this section, we present the results of applying the concepts developed in this paper to wealth and health data in five developing countries. Using the methodology described in the previous section we have estimated the parameters of the Sarmanov–Lee distribution considering the beta distribution for modelling marginal distributions (Table 1). Standard errors have been simulated by Monte Carlo with 1000 replications.
Table 2 presents information about the level of well-being inequality in each of the analysed nations. Looking first at unidimensional Gini indices, our estimates suggest that Ethiopia is the most egalitarian country in terms of wealth, with a Gini index of 0.257. The most unequal nation, Nigeria, presents a Gini index of 0.6094. Moving now to health inequality, Angola is by far the country with the lowest Gini index, with a value of 0.1048. Ethiopia and Haiti present a similar level of inequality, around 0.62. Relatively close to this value, we find Nigeria and Uganda, the countries with the highest level of health inequality.
Turning our attention to the bivariate Gini index, our estimates reveal that Angola presents the most equal distribution of well-being if health and wealth are jointly considered. It is followed by Ethiopia, with a Gini index of 0.6921. Then, we find Haiti, where multivariate inequality raises to 0.7273. Uganda and Nigeria are again the most unequal nations, with Gini indices of 0.7954 and 0.864, respectively.
Although Angola presents the lowest bivariate Gini index, this country is not the most equal in terms of wealth. This indicates that some sectors of the population (principally the poorest) might be more equally distributed in other countries. To provide a complete picture of the distribution of well-being, it is therefore essential to examine the bidimensional Lorenz surface. Figure 1 shows the bidimensional Lorenz surfaces for wealth and health in Ethiopia and Angola (panel a) and Uganda and Angola (panel b). The comparison of well-being distributions in Ethiopia and Angola reveals that there is not Lorenz dominance, since the Lorenz surface of Ethiopia lies above the Lorenz surface of Angola in the upper-right corner. This suggests that well-being among the poorest individuals with reasonably good health is more equally distributed in Ethiopia. If we now compare the Lorenz surface of Angola with the Lorenz surface of Uganda, the country with the highest bivariate Gini index, we observe a clear pattern of dominance.
We conclude this illustration with the decomposition of the bivariate Gini index (Equations (15) and (16)). Overall equality (expressed as one minus the bivariate Gini index) is decomposed in terms of equality within each dimension and equality between dimensions. The last two columns of Table 2 show the estimates for these two components for the five countries involved in the analysis. These results suggest that multidimensional well-being equality is mainly determined by the within-dimensional distribution, while equality between dimensions represents a residual proportion.

7. Conclusions

In this paper, using the definition proposed by Arnold [16], we presented closed expressions for the bivariate Lorenz surface and its associated Gini index, assuming different formulations of the underlying bivariate income distribution. We discussed a relevant type of models based on the class of bivariate distributions with given marginals described by Sarmanov and Lee [19,20]. For the marginal distributions, We considered several candidates including the Beta, the GB1, the Gamma and the lognormal distributions. Extensions of the bivariate Lorenz surface to higher dimensions were also included. The implementation of multivariate Lorenz surfaces introduced in this paper was illustrated with an application to inequality in well-being in five developing countries.

Author Contributions

Conceptualization, J.M.S.; Data curation, V.J.; Investigation, J.M.S.; Methodology, J.M.S.; Software, and V.J.; Writing—original draft, V.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was partly funded by the Ministerio de Economía y Competitividad, grant number PID2019-105986GB-C22.

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.

Abbreviations

The following abbreviations are used in this manuscript:
LCLorenz curve
GB1Generalised beta of the first kind
FGMFarlie–Gumbel–Morgenstern
SLSarmanov–Lee
DHSDemographic and Health Surveys
UNDPUnited Nations Development Program

References

  1. Sen, A. The Concept of Development. In Handbook of Development Economics; Chenery, H., Srinivasan, T.N., Eds.; Elsevier: Amsterdam, The Netherlands, 1988; pp. 9–26. [Google Scholar]
  2. Sen, A. Development as Capabilities Expansion. J. Dev. Plan. 1989, 19, 41–58. [Google Scholar]
  3. Sen, A. Development as Freedom; Oxford University Press: Oxford, UK, 1999. [Google Scholar]
  4. Jorda, V.; López-Noval, B.; Sarabia, J.M. Distributional dynamics of life satisfaction in Europe. J. Happiness Stud. 2019, 20, 1015–1039. [Google Scholar] [CrossRef] [Green Version]
  5. Jorda, V.; Sarabia, J.M. International convergence in well-being indicators. Soc. Indic. Res. 2015, 120, 1–27. [Google Scholar] [CrossRef]
  6. Atkinson, A.B. Multidimensional deprivation: Contrasting social welfare and counting approaches. J. Econ. Inequal. 2003, 1, 51–65. [Google Scholar] [CrossRef]
  7. Atkinson, A.B.; Bourguignon, F. The comparison of multi-dimensioned distributions of economic status. Rev. Econ. Stud. 1982, 49, 183–201. [Google Scholar] [CrossRef]
  8. Decancq, K.; Lugo, M.A. Inequality of Wellbeing: A Multidimensional Approach. Economica 2012, 79, 721–746. [Google Scholar] [CrossRef]
  9. Kolm, S.C. Multidimensional Equalitarianisms. Q. J. Econ. 1977, 91, 1–13. [Google Scholar] [CrossRef]
  10. Maasoumi, E. The measurement and decomposition of multi-dimensional inequality. Econometrica 1986, 54, 991–997. [Google Scholar] [CrossRef] [Green Version]
  11. Slottje, D.J. Relative Price Changes and Inequality in the Size Distribution of Various Components. J. Bus. Econ. Stat. 1987, 5, 19–26. [Google Scholar]
  12. Tsui, K.Y. Multidimensional generalizations of the relative and absolute inequality indices: The Atkinson-Kolm-Sen approach. J. Econ. Theory 1995, 67, 251–265. [Google Scholar] [CrossRef]
  13. Tsui, K.Y. Multidimensional inequality and multidimensional generalized entropy measures: An axiomatic derivation. Soc. Choice Welf. 1999, 16, 145–157. [Google Scholar] [CrossRef] [Green Version]
  14. Taguchi, T. On the two-dimensional concentration surface and extensions of concentration coefficient and Pareto distribution to the two-dimensional case-I. Ann. Inst. Stat. Math. 1972, 24, 355–382. [Google Scholar] [CrossRef]
  15. Taguchi, T. On the two-dimensional concentration surface and extensions of concentration coefficient and Pareto distribution to the two-dimensional case-II. Ann. Inst. Stat. Math. 1972, 24, 599–619. [Google Scholar] [CrossRef]
  16. Arnold, B.C. Pareto Distributions; International Co-Operative Publishing House: Fairland, MD, USA, 1983. [Google Scholar]
  17. Koshevoy, G.; Mosler, K. The Lorenz zonoid of a multivariate distribution. J. Am. Stat. Assoc. 1996, 91, 873–882. [Google Scholar] [CrossRef]
  18. Sarabia, J.M.; Jordá, V. Bivariate Lorenz Curves based on the Sarmanov–Lee Distribution. In Topics in Statistical Simulation, Springer Proceedings in Mathematics & Statistics; Springer: New York, NY, USA, 2014; Volume 114, pp. 447–455. [Google Scholar]
  19. Lee, M.-L.T. Properties of the Sarmanov Family of Bivariate Distributions. Commun. Stat. Theory Methods 1996, 25, 1207–1222. [Google Scholar]
  20. Sarmanov, O.V. Generalized Normal Correlation and Two-Dimensional Frechet Classes. Doklady 1966, 168, 596–599. [Google Scholar]
  21. Gastwirth, J.L. A general definition of the Lorenz curve. Econometrica 1971, 39, 1037–1039. [Google Scholar] [CrossRef]
  22. Kakwani, N.C. Applications of Lorenz Curves in Economic Analysis. Econometrica 1977, 45, 719–728. [Google Scholar] [CrossRef]
  23. Taguchi, T. On the structure of multivariate concentration - some relationships among the concentration surface and two variate mean difference and regressions. Comput. Stat. Data Anal. 1988, 6, 307–334. [Google Scholar] [CrossRef]
  24. Arnold, B.C. Majorization and the Lorenz Curve; Lecture Notes in Statistics 43; Springer: New York, NY, USA, 1987. [Google Scholar]
  25. Koshevoy, G. Multivariate Lorenz majorization. Soc. Choice Welf. 1995, 12, 93–102. [Google Scholar] [CrossRef]
  26. Koshevoy, G.; Mosler, K. Multivariate Gini indices. J. Multivar. Anal. 1997, 60, 252–276. [Google Scholar] [CrossRef]
  27. Mosler, K. Multivariate Dispersion, Central Regions and Depth: The Lift Zonoid Approach; Springer: Berlin, Germany, 2002. [Google Scholar]
  28. Kleiber, C.; Kotz, S. Statistical Size Distributions in Economics and Actuarial Sciences; John Wiley: Hoboken, NJ, USA, 2003. [Google Scholar]
  29. Sarabia, J.M.; Prieto, F.; Jordá, V. Bivariate beta-generated distributions with applications to well-being data. J. Stat. Distrib. Appl. 2014, 1, 1–15. [Google Scholar] [CrossRef] [Green Version]
  30. Alexander, C.; Cordeiro, G.M.; Ortega, E.M.M.; Sarabia, J.M. Generalized beta-generated distributions. Comput. Stat. Data Anal. 2012, 56, 1880–1897. [Google Scholar] [CrossRef]
  31. Arnold, B.C.; Sarabia, J.M. Analytic Expressions for Multivariate Lorenz Surfaces. Sankhya A Indian J. Stat. 2018, 80, 84–111. [Google Scholar] [CrossRef]
  32. Balakrishnan, N.; Lai, C.-D. Continuous Bivariate Distributions; Springer: New York, NY, USA, 2009. [Google Scholar]
  33. Donaldson, D.; Weymark, J.A. A single parameter generalization of the Gin index of inequality. J. Econ. Theory 1980, 22, 67–86. [Google Scholar] [CrossRef]
  34. Kakwani, N.C. Income Inequality and Poverty, Methods and Estimation and Policy Applications; Oxford University Press: New York, NY, USA, 1980. [Google Scholar]
  35. Yitzhaki, S. On an extension of the Gini inequality index. Int. Econ. Rev. 1983, 24, 617–628. [Google Scholar] [CrossRef]
  36. Huang, J.S.; Kotz, S. Modifications of the Farlie-Gumbel-Morgenstern distributions. A tough hill to climb. Metrika 1999, 49, 135–145. [Google Scholar] [CrossRef]
  37. Bairamov, I.; Kotz, S. On a new family of positive quadrant dependent bivariate distributions. Int. Math. J. 2003, 3, 1247–1254. [Google Scholar]
  38. Arnold, B.C.; Sarabia, J.M. Majorization and the Lorenz Order with Applications in Applied Mathematics and Economics; Springer: New York, NY, USA, 2018. [Google Scholar]
  39. McDonald, J.B. Some generalized functions for the size distribution of income. Econometrica 1984, 52, 647–663. [Google Scholar] [CrossRef]
  40. Sarabia, J.M. Parametric Lorenz Curves: Models and Applications. In Modeling Income Distributions and Lorenz Curves; Chotikapanich, D., Ed.; Springer: New York, NY, USA, 2008; pp. 167–190. [Google Scholar]
  41. Sarabia, J.M.; Castillo, E.; Pascual, M.; Sarabia, M. Bivariate Income Distributions with Lognormal Conditionals. J. Econ. Inequal. 2007, 5, 371–383. [Google Scholar] [CrossRef] [Green Version]
  42. Sarabia, J.M.; Castillo, E.; Slottje, D. An Ordered Family of Lorenz Curves. J. Econ. 1999, 91, 43–60. [Google Scholar] [CrossRef]
  43. Sarabia, J.M.; Castillo, E.; Pascual, M.; Sarabia, M. Mixture Lorenz Curves. Econ. Lett. 2005, 89, 89–94. [Google Scholar] [CrossRef]
  44. Hlasny, V. Nonresponse Bias in Inequality Measurement: Cross-Country Analysis Using Luxembourg Income Study Surveys. Soc. Sci. Q. 2020, 101, 712–731. [Google Scholar] [CrossRef]
  45. Jorda, V.; Sarabia, J.M.; Jäntti, M. Estimation of income inequality from grouped data. arXiv 2018, arXiv:1808.09831. [Google Scholar]
  46. Jenkins, S.P. Pareto models, top incomes and recent trends in UK income inequality. Economica 2017, 84, 261–289. [Google Scholar] [CrossRef] [Green Version]
  47. Andreoli, F.; Zoli, C. Measuring Dissimilarity; Working Paper Series; Department of Economics University of Verona: Verona, Italy, 2014; Volume 23. [Google Scholar]
  48. Andreoli, F.; Zoli, C. From unidimensional to multidimensional inequality: A review. Metron 2020, 78, 5–42. [Google Scholar] [CrossRef]
  49. Stiglitz, J.E.; Sen, A.; Fitoussi, J.P. Report by the Commission on the Measurement of Economic Performance and Social Progress; Citeseer: University Park, PA, USA, 2009. [Google Scholar]
  50. Hlasny, V.; AlAzzawi, S. Asset inequality in the MENA: The missing dimension? Q. Rev. Econ. Financ. 2019, 73, 44–55. [Google Scholar] [CrossRef]
  51. Balarajan, Y.; Ramakrishnan, U.; Ozaltin, E.; Shankar, A.H.; Subramanian, S.V. Anaemia in low-income and middle-income countries. Lancet 2011, 378, 2123–2135. [Google Scholar] [CrossRef]
  52. Scanlon, K.S.; Yip, R.; Schieve, L.A.; Cogswell, M.E. High and low hemoglobin levels during pregnancy: Differential risks for preterm birth and small for gestational age. Obstet. Gynecol. 2000, 96, 741–748. [Google Scholar] [CrossRef]
Figure 1. Bivariate Lorenz surfaces: evolution of multidimensional inequality in wealth and health. The curves are defined by Equation (14) with components defined in Equations (20) and (21).
Figure 1. Bivariate Lorenz surfaces: evolution of multidimensional inequality in wealth and health. The curves are defined by Equation (14) with components defined in Equations (20) and (21).
Mathematics 08 02095 g001aMathematics 08 02095 g001b
Table 1. Parameter estimates for the Sarmanov–Lee distribution with classical beta marginals (Equations (24)–(26)). Bootstrap standard errors in parenthesis.
Table 1. Parameter estimates for the Sarmanov–Lee distribution with classical beta marginals (Equations (24)–(26)). Bootstrap standard errors in parenthesis.
IncomeHealthDependence
a ^ 1 b ^ 1 a ^ 2 b ^ 2 w
Angola1.26602.095217.831829.82693.6167
(0.0188)(0.0438)(1.2273)(4.2834)(0.7054)
Ethiopia2.82695.27410.23380.60672.3398
(0.0758)(0.335)(0.0147)(0.027)(0.114)
Haiti1.57943.87060.38141.89322.8026
(0.0284)(0.0855)(0.0231)(0.0832)(0.1965)
Nigeria0.45323.10330.13240.33061.1249
(0.012)(0.0684)(0.0068)(0.0144)(0.1341)
Uganda0.94732.20200.19940.63141.0116
(0.0185)(0.0539)(0.0211)(0.0534)(0.2055)
Table 2. Bivariate Gini index (Equation (15)) and marginal Gini indices (27) and Decomposition of the equality (Equation (16)) for the Sarmanov–Lee distribution with classical beta marginals.
Table 2. Bivariate Gini index (Equation (15)) and marginal Gini indices (27) and Decomposition of the equality (Equation (16)) for the Sarmanov–Lee distribution with classical beta marginals.
InequalityEquality
WealthHealthBidimensionalTotalWithinBetween
Gini VariablesVariables
Angola0.35160.10480.39430.60570.57660.0291
Ethiopia0.25700.62600.69210.30790.26080.0471
Haiti0.34670.62140.72730.27270.22940.0434
Nigeria0.60940.66780.86400.13600.11740.0186
Uganda0.41980.66800.79540.20460.17990.0247
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sarabia, J.M.; Jorda, V. Lorenz Surfaces Based on the Sarmanov–Lee Distribution with Applications to Multidimensional Inequality in Well-Being. Mathematics 2020, 8, 2095. https://doi.org/10.3390/math8112095

AMA Style

Sarabia JM, Jorda V. Lorenz Surfaces Based on the Sarmanov–Lee Distribution with Applications to Multidimensional Inequality in Well-Being. Mathematics. 2020; 8(11):2095. https://doi.org/10.3390/math8112095

Chicago/Turabian Style

Sarabia, José María, and Vanesa Jorda. 2020. "Lorenz Surfaces Based on the Sarmanov–Lee Distribution with Applications to Multidimensional Inequality in Well-Being" Mathematics 8, no. 11: 2095. https://doi.org/10.3390/math8112095

APA Style

Sarabia, J. M., & Jorda, V. (2020). Lorenz Surfaces Based on the Sarmanov–Lee Distribution with Applications to Multidimensional Inequality in Well-Being. Mathematics, 8(11), 2095. https://doi.org/10.3390/math8112095

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