Next Article in Journal
Non-Fragile Sampled Control Design for an Interconnected Large-Scale System via Wirtinger Inequality
Next Article in Special Issue
Maximum Penalized Likelihood Estimation of the Skew–t Link Model for Binomial Response Data
Previous Article in Journal
Computational Study of Shocked V-Shaped N2/SF6 Interface across Varying Mach Numbers
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bivariate Pareto–Feller Distribution Based on Appell Hypergeometric Function

by
Christian Caamaño-Carrillo
1,
Moreno Bevilacqua
2,3,
Michael Zamudio-Monserratt
4 and
Javier E. Contreras-Reyes
5,*
1
Departamento de Estadística, Facultad de Ciencias, Universidad del Bío-Bío, Concepción 4081112, Chile
2
Facultad de Ingeniería y Ciencia, Universidad Adolfo Ibañez, Viña del Mar 2580335, Chile
3
Dipartimento di Scienze Ambientali, Informatica e Statistica, Ca’ Foscari University of Venice, 30172 Venice, Italy
4
NCT/Odisseia-Observatory of Socio-Environmental Dynamics (UnB), Universidade de Brasilia, Brasilia 70910-900, Brazil
5
Instituto de Matemática, Física y Estadística, Facultad de Ingeniería y Negocios, Universidad de Las Américas, Sede Viña del Mar, 7 Norte 1348, Viña del Mar 2531098, Chile
*
Author to whom correspondence should be addressed.
Submission received: 25 July 2024 / Revised: 25 September 2024 / Accepted: 1 October 2024 / Published: 9 October 2024
(This article belongs to the Special Issue Advances in Statistical Simulation and Computing)

Abstract

:
The Pareto–Feller distribution has been widely used across various disciplines to model “heavy-tailed” phenomena, where extreme events such as high incomes or large losses are of interest. In this paper, we present a new bivariate distribution based on the Appell hypergeometric function with marginal Pareto–Feller distributions obtained from two independent gamma random variables. The proposed distribution has the beta prime marginal distributions as special case, which were obtained using a Kibble-type bivariate gamma distribution, and the stochastic representation was obtained by the quotient of a scale mixture of two gamma random variables. This result can be viewed as a generalization of the standard bivariate beta I (or inverted bivariate beta distribution). Moreover, the obtained bivariate density is based on two confluent hypergeometric functions. Then, we derive the probability distribution function, the cumulative distribution function, the moment-generating function, the characteristic function, the approximated differential entropy, and the approximated mutual information index. Based on numerical examples, the exact and approximated expressions are shown.

1. Introduction

The Pareto distribution has been a key tool in economics and for modeling wealth distribution, as well as in other areas like insurance and finance, where capturing extreme events is crucial [1]. However, real-world data often exhibit more complex structures than those described by the classical Pareto distribution [2]. In this context, Feller (1971) [3] proposed an extension called the Pareto–Feller distribution which includes an additional shape parameter, allowing the distribution to better fit a wider range of phenomena with thicker or heavy tails as needed [1].
The Pareto–Feller distribution has been widely used across various disciplines to model “heavy-tailed” phenomena, where extreme events such as high incomes or large losses are of interest. The Pareto–Feller distribution emerged from the need for a distribution that offers greater flexibility in data modeling by introducing an additional parameter to control the tail shape and skewness, thus providing a more accurate description of empirical data compared with the standard Pareto distribution [2]. This distribution has found applications in various fields such as risk theory, river flow modeling, and natural disaster analysis due to its ability to represent both heavy tails and asymmetric distributions. The additional flexibility provided by this distribution is especially valuable in situations where conventional distributions, like the classical Pareto distribution, fail to capture the observed variability and extreme behavior in the data. Notable variants include Pareto type I, used in wealth analysis [1], the generalized Pareto distribution for modeling extreme events [4], and the Lomax or Pareto type II distribution, which is applied in survival analysis [5]. Other important variants are Pareto type IV, which offers greater flexibility in tail shapes [6], and the truncated Pareto distribution, which is used in scenarios with physical upper limits [7]. Additionally, the log-logistic distribution, which is employed in contexts similar to the Pareto distribution but with greater flexibility in the tails, is used in survival analysis and system failure studies [8].
The Pareto–Feller distribution can be constructed as a location-scale transformation of the ratio of two independent gamma-distributed random variables. This method allows the distribution to capture a wide range of tail behaviors and offers flexibility in modeling heavy-tailed phenomena. The use of gamma distributions for generating such models is well documented in the statistical literature. Specifically, more flexible models can be obtained through an appropriate transformation of a bivariate gamma distribution or independent copies of it. This approach is appealing because the correlation of the transformed bivariate gamma distribution is directly tied to the correlation of the original bivariate gamma distribution. Such constructions are commonly used in spatial data modeling. Examples include t distributed spatial models [9], Weibull spatial models [10], and Poisson spatial models [11]. On the other hand, Kotz et al. (2004) [12] described similar methods for transforming distributions via ratios of gamma variables which are widely applicable in reliability and survival analysis contexts [12].
Since the construction of the Pareto–Feller distribution is related to the gamma distribution, it is necessary to first define the bivariate gamma distribution to develop the bivariate case of the Pareto–Feller distribution. Thus, we start by defining a sequence of independent normal random variables and show how these lead to a gamma distribution as follows. Let Z i k , i = 1 , , ν with ν > 2 be a sequence of independent standardized normal random variables whose correlation function is given by C o r r ( Z i k , Z j k ) = ρ , i j , and let
W k = i = 1 ν Z i k 2 α , α > 0 , k = 1 , 2 .
Then, W k is a random variable with a gamma marginal distribution (i.e., W k G a m m a ( ν / 2 , α / 2 ) , where ν / 2 represents the shape parameter and α / 2 represents the rate parameter), with the probability density function (pdf) being given by
f W k = α ν / 2 2 ν / 2 Γ ( ν / 2 ) w k ν / 2 1 e α w k / 2 ,
where E [ W k ] = ν / α and V a r ( W k ) = 2 ν / α 2 [10].
The construction of a bivariate Pareto–Feller distribution is derived from the ratio of two bivariate gamma distributions as shown in Section 2. We consider the bivariate vector W = ( W 1 , W 2 ) , where the stochastic representation of W k , k = 1 , 2 is given in Equation (1). Vere-Jones [13] showed that the distribution of W has a correlated bivariate gamma distribution with the parameters ν > 0 and γ > 0 , while the pdf is given by
f W ( w ) = 2 ν γ ν ( w 1 w 2 ) ν / 2 1 e γ ( w 1 + w 2 ) 2 ( 1 ρ 2 ) Γ ν 2 ( 1 ρ 2 ) ν / 2 γ ρ 2 w 1 w 2 2 ( 1 ρ 2 ) 1 ν / 2 I ν / 2 1 γ ρ 2 w 1 w 2 ( 1 ρ 2 ) ,
where I α ( · ) is the usual modified Bessel function of the γ -order of the first kind. Gamma variables can be used as building blocks for the construction of flexible non-Gaussian variables. Henceforth, we will call W a gamma random vector with an underlying correlation ρ [14,15] such that the correlation of the gamma bivariate is ρ W = ρ 2 . Moreover, when ρ = 0 , Equation (2) can be written as the product of two independent gamma random variables (i.e., W k G a m m a ( ν / 2 , α / 2 ) , k = 1 , 2 ). Thus, zero pairwise correlation implies pairwise independence, as in the Gaussian case. The pdf f W was first discussed in [16], and its properties were studied in [17,18].
The bivariate Pareto–Feller distribution represents a less commonly discussed extension in the statistical literature. It builds upon the principles established by the univariate Pareto distribution. The authors of [19] addressed extreme value distributions and included discussions which may be relevant to bivariate generalizations. Additionally, the authors of [20] provided further insights into bivariate distributions, enriching the understanding of Pareto–Feller distributions. These references offer both theoretical and practical frameworks for researching and applying Pareto–Feller distributions in bivariate contexts. The latter works motivated this paper, which presented a bivariate Pareto–Feller distribution built from an Appell hypergeometric function.
This paper is organized as follows. Section 2 presents the bivariate Pareto–Feller distribution. In particular, the pdf, cumulative distribution function (cdf), joint moment-generating function (mgf), characteristic functions, cross-product moment function, mean, variance, covariance, and correlation function are presented. In Section 3, some approximations of the differential entropy and, consequently, the mutual information index are presented. Finally, some discussions and conclusions are presented in Section 4. All simulations included special hypergeometric functions such as the Appell hypergeometric one, and all were implemented in R 4.4.1. software [21] using the zipfR and hypergeo packages. All proofs of theorems and propositions can be found in the Appendix A.

2. Bivariate Pareto–Feller Distribution

Let us define a random variable V with support on the positive real line, defined as a scale mixture of two gamma random variables:
V = W R ,
where W G a m m a ( α / 2 , 1 ) , α > 0 and R G a m m a ( ν / 2 , 1 ) , ν > 0 . Then, V is a random variable with a marginal distribution of the beta I type or beta prime [12,22] and denoted by V B e ( ν / 2 , α / 2 ) , with the pdf given by
f V ( v ) = Γ ν + α 2 Γ ν 2 Γ α 2 v ν / 2 1 ( v + 1 ) ( ν + α ) / 2 .
The beta prime distribution is anchored to a shape parameter of the gamma distributions. This construction was previously proposed in [23,24,25].
Based on the stochastic representation in Equation (3), we consider the bivariate vector of V = ( V 1 , V 2 ) , where
V i = W i R i , i = 1 , 2 ;
Here, W = ( W 1 , W 2 ) and R = ( R 1 , R 2 ) are two correlated bivariate gamma distributions with correlations ρ W = ρ R = ρ 2 , where W i G a m m a ( α / 2 , 1 ) , α > 0 and R i G a m m a ( ν / 2 , 1 ) , ν > 0 , R i W j , i , j , i , j = 1 , 2 . Thus, V i B e ( ν / 2 , α / 2 ) , i = 1 , 2 .
A new bivariate distribution with a beta prime marginal distribution obtained from the Kibble-type bivariate gamma distribution given in Equation (2) is presented in the following theorem. This result can be viewed as a generalization of the standard bivariate beta I distribution (or inverted bivariate beta distribution) [26].
Theorem 1.
Let W and R be two independent gamma random variables, and let V = W R 1 . Then, the pdf of V = ( V 1 , V 2 ) R + × R + is given by
f V ( v ) = ( v 1 v 2 ) ν / 2 1 Γ 2 ν + α 2 [ ( v 1 + 1 ) ( v 2 + 1 ) ] ( ν + α ) / 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 × F 4 ν + α 2 , ν + α 2 ; ν 2 , α 2 ; ρ 2 v 1 v 2 ( v 1 + 1 ) ( v 2 + 1 ) , ρ 2 ( v 1 + 1 ) ( v 2 + 1 ) ,
where  F 4  is an Appell hypergeometric function of the fourth kind, defined as
F 4 ( a , b ; c , c ; w , z ) = k = 0 m = 0 ( a ) k + m ( b ) k + m w k z m k ! m ! ( c ) k ( c ) m , | w | + | z | < 1 .
The special functions F 4 and the Gaussian hypergeometric function F 1 2 are related through the identity
F 4 ( a , b ; c , c ; w , z ) = k = 0 ( a ) k ( b ) k z k k ! ( c ) k F 1 2 ( a + k , b + k ; c ; w ) , | w | + | z | < 1 .
where ( a ) k : = Γ ( a + k ) / Γ ( a ) , for which k I N { 0 } is the Pochhammer symbol. The Gaussian hypergeometric function is a special case of more general power series, where the generalized hypergeometric function is defined for p , q = 0 , 1 , 2 , as
F q p ( a 1 , a 2 , , a p ; b 1 , b 2 , , b q ; x ) : = k = 0 ( a 1 ) k , ( a 2 ) k , , ( a p ) k ( b 1 ) k , ( b 2 ) k , , ( b q ) k x k k !
When ρ = 0 , the pdf in Theorem 1 involves the product of two independent beta prime random variables B e ( ν / 2 , α / 2 ) .
We now consider a new random variable Y, defined as
Y : = μ + q c V 1 / p ,
where Y μ 0 , c = Γ ν 2 + 1 p Γ α 2 1 p Γ ν 2 Γ α 2 , q > 0 , and p > 0 . This random variable is a marginal Pareto–Feller distribution. Specifically, using the notation of [20], we marginally have Y P F ( μ , q / c , 1 / p , ν / 2 , α / 2 ) with a density defined by
f Y ( y ) = ( c p ) q c ( y μ ) q p ν / 2 1 Γ ν + α 2 Γ ν 2 Γ α 2 c ( y μ ) q p + 1 ( ν + α ) / 2 , y > μ ,
and a mean and variance given by
E [ Y ] = μ + q , V a r ( Y ) = q 2 Γ ν 2 + 2 p Γ α 2 2 p c Γ ν 2 + 1 p Γ α 2 1 p 1 ,
respectively, with α p > 4 [3].
The Pareto–Feller distribution includes as special cases different types of Pareto random variables definitions (type I, II, III, and IV; see [20]) and the so-called beta prime one. If we consider μ = 0 , then Y P F ( 0 , q / c , 1 / p , ν / 2 , α / 2 ) , and Y = ( Y 1 , Y 2 ) is a random vector with a marginal Pareto–Feller distribution. A new bivariate distribution based on marginal Pareto–Feller distributions is presented in the following theorem.
Theorem 2.
Let Y = ( Y 1 , Y 2 ) , where Y i : = q i c V i 1 / p , i = 1 , 2 . The pdf of Y is given by
f Y ( y ) = ( c p ) 2 q 1 q 2 c 2 y 1 y 2 q 1 q 2 p ν / 2 1 Γ 2 ν + α 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 c y 1 q 1 p + 1 c y 2 q 2 p + 1 ( ν + α ) / 2 × F 4 ν + α 2 , ν + α 2 ; ν 2 , α 2 ; ρ 2 ( c 2 y 1 y 2 ) p ( q 1 q 2 ) p c y 1 q 1 p + 1 c y 2 q 2 p + 1 , ρ 2 c y 1 q 1 p + 1 c y 2 q 2 p + 1 .
The pdf in Theorem 2 considers an Appell hypergeometric function F 4 . We can write this as a series of hypergeometric functions F 1 2 .
Figure 1 illustrates the pdf of Equation (7) for some parameters. When ρ increases, the largest values of y 1 and y 2 in the pdf are produced. However, these values depend on the other parameters. Independent of the ρ value, the pdf is close at the origin ( y 1 , y 2 ) ( 0 , 0 ) with a positive bias, and it decays exponentially for the smallest values ( ν = 4 , α = 4 , q 1 = 4 , q 2 = 4 , and p = 4 ). When ρ = 0.9 and ν , α = 8 or 12, for example, the pdf has more symmetry and variability but less bias. Note that the pdf of Y , given in Equation (5), is symmetric for negative values of ρ . Specifically, the same representations hold for ρ = 0.25 , 0.5 , 0.9 from Figure 1 while keeping the other parameters fixed.
Theorem 3.
The joint cdf of Y = ( Y 1 , Y 2 ) in Equation (7) can be expressed as
F Y ( Y 1 t 1 , Y 2 t 2 ) = c 2 t 1 t 2 q 1 q 2 p ν / 2 Γ 2 ν + α 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 × k = 0 m = 0 ν + α 2 k + m 2 ρ 2 k + 2 m k ! m ! ν 2 k α 2 m ( ν / 2 + k ) 2 c 2 t 1 t 2 q 1 q 2 p k × F 1 2 ν + α 2 + k + m , ν 2 + k ; ν 2 + k + 1 ; c t 1 q 1 p × F 1 2 ν + α 2 + k + m , ν 2 + k ; ν 2 + k + 1 ; c t 2 q 2 p .
Proposition 1.
The joint mgf and characteristic functions of Y = ( Y 1 , Y 2 ) given in Equation (7) are
( a ) M Y ( t 1 , t 2 ) = ( c p ) 2 q 1 q 2 Γ 2 ν + α 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 × k = 0 m = 0 ν + α 2 k + m ν + α 2 k + m ν 2 k α 2 m k ! m ! ρ 4 k m ( c 2 / q 1 q 2 ) ( p ν / 2 + p k 1 ) × 1 p q c p ( ν / 2 + k ) Γ ( ν / 2 + k ) Γ ( k + m + ( ( ν + α ) / 2 ) Γ ( m + ( ( ν + α ) / 2 ) + ν / 2 ) 2 ( t 1 ) 2 1 ( t 1 ) 4 ( t 2 ) 2 1 ( t 2 ) 4 ,
with t < 0 and
( b ) ϕ Y ( t 1 , t 2 ) = ( c p ) 2 q 1 q 2 Γ 2 ν + α 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 × k = 0 m = 0 ν + α 2 k + m ν + α 2 k + m ν 2 k α 2 m k ! m ! ρ 4 k m ( c 2 / q 1 q 2 ) ( p ν / 2 + p k 1 ) × 1 p q c p ( ν / 2 + k ) Γ ( ν / 2 + k ) Γ ( k + m + ( ( ν + α ) / 2 ) Γ ( m + ( ( ν + α ) / 2 ) + ν / 2 ) 2 × ( i t 1 ) 2 1 ( i t 1 ) 4 ( i t 2 ) 2 1 ( i t 2 ) 4 .
Proposition 2.
The cross-product moment of Y = ( Y 1 , Y 2 ) in Equation (7) can be expressed as
E [ Y 1 a Y 2 b ] = q 1 a q 2 b Γ α 2 a p Γ α 2 b p Γ ν 2 + a p Γ ν 2 + b p c a + b Γ 2 ν 2 Γ 2 α 2 × F 1 2 a p , b p ; ν 2 ; ρ 2 F 1 2 a p , b p ; α 2 ; ρ 2 .
Proposition 2 illustrates that the cross-product moment is the product of two Gaussian hypergeometric functions. Corollary 1 is straightforward from Proposition 2, where the expected value and variance of a marginal Pareto–Feller random variable Y k are presented as well as the covariance and correlation between two marginal Pareto–Feller random variables ( Y 1 and Y 2 ).
Corollary 1.
If Y = ( Y 1 , Y 2 ) , then it has a pdf according to Equation (7). According to Proposition 2, we have the following:
1. 
E [ Y i ] = q i , i = 1 , 2 .
2. 
V a r ( Y i ) = q i 2 Γ ν 2 + 2 p Γ α 2 2 p c Γ ν 2 + 1 p Γ α 2 1 p 1 if α p > 4 for i = 1 , 2 .
3. 
C o v ( Y 1 , Y 2 ) = q i 2 Γ α 2 1 p Γ α 2 1 p Γ ν 2 + 1 p Γ ν 2 + 1 p c 2 Γ 2 ν 2 Γ 2 α 2 × F 1 2 1 p , 1 p ; ν 2 ; ρ 2 F 1 2 1 p , 1 p ; α 2 ; ρ 2 .
4. 
Let ρ Y C o r r ( Y 1 , Y 2 ) . Thus, we have
ρ Y = Γ 2 ν 2 + 1 p Γ 2 α 2 1 p Γ ν 2 Γ α 2 Γ ν 2 + 2 p Γ α 2 2 p Γ 2 ν 2 + 1 p Γ 2 α 2 1 p × F 1 2 1 p , 1 p ; ν 2 ; ρ 2 F 1 2 1 p , 1 p ; α 2 ; ρ 2 1 .
Figure 2 shows the correlation ρ Y for some density parameters. When ν and α increase, ρ Y increases. More specifically, when p increases, the correlation ρ Y increases with small-to-large values of ν and α . When ρ increases, ρ Y increases, as does its maximum value (from 0.05 to 0.80). Nevertheless, Corollary 1(4) illustrates that the correlation ρ Y does not depend on either q or c.

3. Differential Entropy and Mutual Information Index

The differential entropy of Y is an information uncertainty measure [27]. The differential entropy of Y = ( Y 1 , Y 2 ) with a pdf f Y ( y ) is
H ( Y ) = E Y [ log { f Y ( Y ) } ] = 0 0 f Y ( y ) log f Y ( y ) d y 1 d y 2 ,
and measures the contained information in Y based on its pdf’s parameters.
On the other hand, the mutual information index (MII) [28] between Y 1 and Y 2 under a dependence assumption ( ρ 0 ) is
M ( Y ) = E log f Y ( y 1 , y 2 ) f Y 1 ( y 1 ) f Y 2 ( y 2 ) = 0 0 f Y ( y 1 , y 2 ) log f Y ( y 1 , y 2 ) f Y 1 ( y 1 ) f Y 2 ( y 2 ) d y 1 d y 2 .
Proposition 3
([29]). Let Y i P F ( μ , q i / c , 1 / p , ν / 2 , α / 2 ) . The entropy of Y i ( i = 1 , 2 ) is
H ( Y i ) = log c p q i Γ ( ν + α 2 ) Γ ( α 2 ) Γ ( ν 2 ) + 1 p α 2 ψ α 2 ψ ν 2 + ν + α 2 ψ ν + α 2 ψ ν 2 ,
where ψ ( X ) = d log Γ ( x ) d x is the digamma function.
Proposition 4
([30]). Let x = u / n . We have that
log ( 1 + x ) = x + O ( n 2 ) ,
as n .
Proposition 5.
If Y = ( Y 1 , Y 2 ) has the pdf given in Equation (7), then the following are true:
(a) 
E Y [ Y i p ] = q i p Γ ν 2 + 1 Γ α 2 1 c p Γ ν 2 Γ α 2 .
(b) 
E Y [ log ( Y i ) ] = 1 p ( 1 ρ 2 ) ( ν + α ) / 2 k = 0 m = 0 ν 2 k α 2 m ρ 2 k + 2 m k ! m ! ψ ν 2 + k ψ α 2 + m log c q i p .
Proposition 6.
An approximation of the differential entropy of Y = ( Y 1 , Y 2 ) is
H ( Y ) log ( c p ) 2 q 1 q 2 Γ 2 ν + α 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 + 1 p ν 2 log c 2 q 1 q 2 + E Y [ log Y 1 ] + E Y [ log Y 2 ] + c p ν + α 2 E Y [ Y 1 p ] q 1 p + E Y [ Y 2 p ] q 2 p ,
where E Y [ Y i p ] and E Y [ log Y i ] , i = 1 , 2 are obtained from parts (a) and (b) of Proposition 5, respectively.
From Equation (15), an MII between Y 1 and Y 2 is expressed in terms of the marginal and joint differential entropies [28,31]. Then, using Propositions 3 and 6, the MII between Y 1 and Y 2 can be approximated as follows:
M ( Y ) = H ( Y 1 ) + H ( Y 2 ) H ( Y ) 2 log c p Γ ( ν + α 2 ) Γ ( α 2 ) Γ ( ν 2 ) log q 1 q 2 + log ( c p ) 2 q 1 q 2 Γ 2 ν + α 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 1 p ν 2 log c 2 q 1 q 2 + E Y [ log Y 1 ] + E Y [ log Y 2 ] c p ν + α 2 E [ Y 1 p ] q 1 p + E [ Y 2 p ] q 2 p + 2 p α ψ α 2 ψ ν 2 + ( ν + α ) ψ ν + α 2 ψ ν 2 .
One particular case is when p = 1 . Thus, using Corollary 1(1), we have
M ( Y ) 1 + ν 2 log c 2 q 1 q 2 + 2 log Γ 2 ( ν + α 2 ) Γ 2 ( ν 2 ) ν + α 2 log ( 1 ρ 2 ) + ( 2 α ) ψ α 2 ψ ν 2 + ( ν + α ) ψ ν + α 2 ψ ν 2 c 1 ν 2 E Y [ log Y 1 ] + E Y [ log Y 2 ] .
Figure 3 illustrates the behavior of the approximated MII obtained in Equation (16) while assuming several values for Y and p = 1 . We observed that M ( Y ) > 0 and increased for q i 0 . As in the correlation function case (Figure 2), the MII increased when the correlation parameter ρ increased.

4. Concluding Remarks

We presented a representation of the Pareto–Feller distribution with a scale mixture of two gamma random variables. The respective stochastic representation was obtained by the quotient of a scale mixture of two gamma random variables. Then, the resulting bivariate density considered the products of two confluent hypergeometric functions. In particular, the probability distribution function, cumulative distribution function, moment generation function, covariance function, correlation function, cross-product moments, and approximations for the differential entropy and, as a consequence, the mutual information index were derived. Some numerical examples illustrated the behavior of the provided expressions.
Some inferential aspects can be addressed in a future work, such as (1) a numerical approach for optimization of the log-likelihood function; (2) the pseudo-likelihood method, considering the optimization of an objective function which depends on a bivariate pdf; (3) the model’s identifiability; (4) a Bayesian approach; and (5) an extention to the multivariate case. We also encourage researching the consideration of a Pareto–Feller distribution in modeling nonnegative bivariate data.

Author Contributions

C.C.-C., M.B., M.Z.-M., and J.E.C.-R. wrote the paper and contributed the reagents, analysis, and materials; M.Z.-M. and J.E.C.-R. conceived, designed, and performed the experiments. All authors have read and agreed to the published version of the manuscript.

Funding

C. Caamaño-Carrillo was partially supported by grant FONDECYT 11220066 from the Chilean government and DIUBB 2120538 IF/R from the University of Bío-Bío. M. Bevilacqua acknowledges financial support from grants FONDECYT 1200068 and ANID/PIA/ANILLOS ACT210096 from the Chilean government. J. Contreras-Reyes’s research was supported by FONDECYT (Chile) grant No. 11190116.

Data Availability Statement

No new data were created or analyzed in this study. Data sharing is not applicable.

Acknowledgments

The authors thank the editor and three anonymous referees for their helpful comments and suggestions.

Conflicts of Interest

The authors declare that there are no conflicts of interest in the publication of this paper.

Appendix A

Proof of Theorem 1.
Following [10], note that the bivariate distribution of the gamma random vector can be rewritten in terms of the hypergeometric function using the identity F 1 0 ( ; b ; x ) = Γ ( b ) x ( 1 b ) / 2 I b 1 ( 2 x ) , given by
f H ( h ) = ( h 1 h 2 ) δ / 2 1 e ( h 1 + h 2 ) ( 1 ρ 2 ) Γ 2 δ 2 ( 1 ρ 2 ) δ / 2 F 1 0 δ 2 ; ρ 2 h 1 h 2 ( 1 ρ 2 ) 2
Under the transformations w 1 = v 1 r 1 and w 2 = v 2 r 2 in Equation (A1) with the Jacobian J ( ( w 1 , w 2 ) ( v 1 , v 2 ) ) = r 1 r 2 , and by using series expansion of the hypergeometric function F 1 0 , we have
f V ( v ) = R + 2 f W | R ( w | r ) f R ( r ) J d r = R + 2 ( v 1 v 2 ) ν / 2 1 ( r 1 r 2 ) ( ν + α ) / 2 1 e ( v 1 r 1 + v 2 r 2 ) ( 1 ρ 2 ) e ( r 1 + r 2 ) ( 1 ρ 2 ) Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 × F 1 0 ν 2 ; ρ 2 v 1 v 2 r 1 r 2 ( 1 ρ 2 ) 2 F 1 0 α 2 ; ρ 2 r 1 r 2 ( 1 ρ 2 ) 2 d r = k = 0 m = 0 R + 2 ( v 1 v 2 ) ν / 2 1 ( r 1 r 2 ) ( ν + α ) / 2 1 e ( v 1 + 1 ) r 1 ( 1 ρ 2 ) e ( v 2 + 1 ) r 2 ( 1 ρ 2 ) Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 k ! m ! ν 2 k α 2 m × ρ 2 ( m h ) v 1 v 2 r 1 r 2 ( 1 ρ 2 ) 2 k ρ 2 r 1 r 2 ( 1 ρ 2 ) 2 m d r = ( v 1 v 2 ) ν / 2 1 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 × k = 0 m = 0 I ( k , m ) k ! m ! ν 2 k α 2 m ρ 2 v 1 v 2 ( 1 ρ 2 ) 2 k ρ 2 ( 1 ρ 2 ) 2 m ,
where, using Fubini’s theorem and Equation (3.381.4) in [32], we obtain
I ( k , m ) = R + r 1 ( ν + α ) / 2 + k + m 1 e ( v 1 + 1 ) r 1 ( 1 ρ 2 ) d r i R + r 2 ( ν + α ) / 2 + k + m 1 e ( v 2 + 1 ) r 2 ( 1 ρ 2 ) d r 2 = Γ 2 ν + α 2 + k + m ( 1 ρ 2 ) v 1 + 1 ν + α 2 + k + m ( 1 ρ 2 ) v 2 + 1 ν + α 2 + k + m
In addition, by combining Equations (A3) and (A2), we obtain
f V 12 ( v 12 ) = ( v 1 v 2 ) ν / 2 1 Γ 2 ν + α 2 [ ( v 1 + 1 ) ( v 2 + 1 ) ] ( ν + α ) / 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 × k = 0 m = 0 ν + α 2 k + m 2 k ! m ! ν 2 k α 2 m ρ 2 v 1 v 2 ( v 1 + 1 ) ( v 2 + 1 ) k ρ 2 ( v 1 + 1 ) ( v 2 + 1 ) m = ( v 1 v 2 ) ν / 2 1 Γ 2 ν + α 2 [ ( v 1 + 1 ) ( v 2 + 1 ) ] ( ν + α ) / 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 × F 4 ν + α 2 , ν + α 2 ; ν 2 , α 2 ; ρ 2 v 1 v 2 ( v 1 + 1 ) ( v 2 + 1 ) , ρ 2 ( v 1 + 1 ) ( v 2 + 1 )
Proof of Theorem 2.
Under the transformations v 1 = ( y 1 / q 1 ) p and v 2 = ( y 2 / q 2 ) p in Equation (5) with a Jacobian J ( ( v 1 ; v 2 ) ( y 1 ; y 2 ) ) = ( p / q ) 2 ( y 1 y 2 / ( q 1 q 2 ) ) p 1 , the pdf of Y is given by Equation (7). □
Proof of Theorem 3.
Using series expansion of the Appell hypergeometric function, we obtain
F Y ( Y 1 t 1 , Y 2 t 2 ) = p 2 c p ν ( q 1 q 2 ) p ν / 2 Γ 2 ν + α 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 0 t 1 0 t 2 ( y 1 y 2 ) p ν / 2 1 c y 1 q 1 p + 1 c y 2 q 2 p + 1 ( ν + α ) 2 × F 4 ν + α 2 , ν + α 2 ; ν 2 , α 2 ; ρ 2 ( c 2 y 1 y 2 ) p ( q 1 q 2 ) p c y 1 q 1 p + 1 c y 2 q 2 p + 1 , ρ 2 c y 1 q 1 p + 1 c y 2 q 2 p + 1 d y = p 2 c p ν ( q 1 q 2 ) p ν / 2 Γ 2 ν + α 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 0 t 1 0 t 2 ( y 1 y 2 ) p ν / 2 1 c y 1 q 1 p + 1 c y 2 q 2 p + 1 ( ν + α ) 2 × k = 0 m = 0 ν + α 2 k + m 2 k ! m ! ν 2 k α 2 m × ρ 2 ( c 2 y 1 y 2 ) p ( q 1 q 2 ) p c y 1 q 1 p + 1 c y 2 q 2 p + 1 k ρ 2 c y 1 q 1 p + 1 c y 2 q 2 p + 1 m d y = p 2 c p ν ( q 1 q 2 ) p ν / 2 Γ 2 ν + α 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 k = 0 m = 0 ν + α 2 k + m 2 ρ 2 m I ( k , m ) k ! m ! ν 2 k α 2 m c 2 p ρ 2 ( q 1 q 2 ) p k
Let z 1 = y 1 p and z 2 = y 2 p be transformations in Equation (A4), where 0 < z i < t i p , i = 1 , 2 . Using Fubini’s theorem and Equation (3.1941.1) in [32], we obtain
I ( k , m ) = 1 p 2 0 t 1 p z 1 ν / 2 + k 1 c q 1 p z 1 + 1 ν + α 2 + k + m d z 1 × 0 t 2 p z 2 ν / 2 + k 1 c q 2 p z 2 + 1 ν + α 2 + k + m d z 2 = ( t 1 t 2 ) p ν / 2 + p k p 2 ( ν / 2 + k ) 2 F 1 2 ν + α 2 + k + m , ν 2 + k ; ν 2 + k + 1 ; c t 1 q 1 p F 1 2 ν + α 2 + k + m , ν 2 + k ; ν 2 + k + 1 ; c t 2 q 2 p
In addition, by combining Equations (A4) and (A5), we obtain Equation (8). □
Proof of Proposition 1.
For part (a), by using the definition of the mgf and series expansion of the Appell hypergeometric function, we obtain
M Y ( t 1 , t 2 ) = p 2 c p ν ( q 1 q 2 ) p ν / 2 Γ 2 ν + α 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 × R + 2 e t 1 y 1 + t 2 y 2 y 1 p ν / 2 1 y 2 p ν / 2 1 c y 1 q 1 p + 1 c y 2 q 2 p + 1 ( ν + α ) 2 × k = 0 m = 0 ν + α 2 k + m 2 k ! m ! ν 2 k α 2 m × ρ 2 ( c 2 y 1 y 2 ) p ( q 1 q 2 ) p c y 1 q 1 p + 1 c y 2 q 2 p + 1 k ρ 2 c y 1 q 1 p + 1 c y 2 q 2 p + 1 m d y = p 2 c p ν ( q 1 q 2 ) p ν / 2 Γ 2 ν + α 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 × k = 0 m = 0 ν + α 2 k + m 2 ρ 2 m I ( k , m ) k ! m ! ν 2 k α 2 m c 2 p ρ 2 ( q 1 q 2 ) p k .
Using Fubini’s theorem in Equation (A6), we obtain
I ( k , m ) = 0 e t 1 y 1 y 1 ( p ν / 2 ) + p k 1 c y 1 q 1 p + 1 ( ν + α ) 2 + k + m d y 1 × 0 e t 2 y 2 y 2 ( p ν / 2 ) + p k 1 c y 2 q 2 p + 1 ( ν + α ) 2 + k + m d y 2 = 1 p q c p ( ν / 2 + k ) Γ ( ν / 2 + k ) Γ ( k + m + ( ( ν + α ) / 2 ) Γ ( m + ( ( ν + α ) / 2 ) + ν / 2 ) e ( t 1 y 1 ) 1 ( t 1 ) 2 | y 1 = 0 y 1 = × 1 p q c p ( ν / 2 + k ) Γ ( ν / 2 + k ) Γ ( k + m + ( ( ν + α ) / 2 ) Γ ( m + ( ( ν + α ) / 2 ) + ν / 2 ) e ( t 2 y 2 ) 1 ( t 2 ) 2 | y 2 = 0 y 2 = = 1 p q c p ( ν / 2 + k ) Γ ( ν / 2 + k ) Γ ( k + m + ( ( ν + α ) / 2 ) Γ ( m + ( ( ν + α ) / 2 ) + ν / 2 ) 2 1 ( t i ) 2 1 1 ( t i ) 2 .
Note that if f ( y ) = e t i y i and
g ( y ) = y i ( p ν / 2 ) + p k 1 c q i p ( y i ) p + 1 ( k + m + ( ν + α ) / 2 ) d y i ,
then we can apply integration by substitution for 0 f ( y ) g ( y ) d y . With the change in variables u i = e t i y i d u i = t i e t i y i , and by letting
d v i = y i ( p ν / 2 ) + p k 1 c q i p ( y i ) p + 1 ( k + m + ( ν + α ) / 2 ) d y i ,
then v i = 0 d v i , where Equation (3.241.4.11) in [32] was used with t i < 0 , i = 1 , 2 . By combining Equations (A6) and (A7), we obtain
M Y ( t 1 , t 2 ) = ( c p ) 2 q 1 q 2 Γ 2 ν + α 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 × k = 0 m = 0 ν + α 2 k + m ν + α 2 k + m ν 2 k α 2 m k ! m ! ρ 4 k m ( c 2 / q 1 q 2 ) ( p ν / 2 + p k 1 ) × 1 p q c p ( ν / 2 + k ) Γ ( ν / 2 + k ) Γ ( k + m + ( ( ν + α ) / 2 ) Γ ( m + ( ( ν + α ) / 2 ) + ν / 2 ) 2 × ( t 1 ) 2 1 ( t 1 ) 4 ( t 2 ) 2 1 ( t 2 ) 4 .
For part (b), and by following the proof of part (a), the proof of the characteristic function of Y is straightforward. □
Proof of Proposition 2.
By the definition of cross-product moment, and using series expansion of the Appell hypergeometric function of the fourth kind, we obtain
E ( Y 1 a Y 2 b ) = p 2 c p ν ( q 1 q 2 ) p ν / 2 Γ 2 ν + α 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 R + 2 y 1 p ν / 2 + a 1 y 2 p ν / 2 + b 1 c y 1 q 1 p + 1 c y 2 q 2 p + 1 ( ν + α ) 2 × F 4 ν + α 2 , ν + α 2 ; ν 2 , α 2 ; ρ 2 ( c 2 y 1 y 2 ) p ( q 1 q 2 ) p c y 1 q 1 p + 1 c y 2 q 2 p + 1 , ρ 2 c y 1 q 1 p + 1 c y 2 q 2 p + 1 d y = p 2 c p ν ( q 1 q 2 ) p ν / 2 Γ 2 ν + α 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 × R + 2 y 1 p ν / 2 + a 1 y 2 p ν / 2 + b 1 c y 1 q 1 p + 1 c y 2 q 2 p + 1 ( ν + α ) 2 × k = 0 m = 0 ν + α 2 k + m 2 k ! m ! ν 2 k α 2 m × ρ 2 ( c 2 y 1 y 2 ) p ( q 1 q 2 ) p c y 1 q 1 p + 1 c y 2 q 2 p + 1 k ρ 2 c y 1 q 1 p + 1 c y 2 q 2 p + 1 m d y = p 2 c p ν ( q 1 q 2 ) p ν / 2 Γ 2 ν + α 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 k = 0 m = 0 ν + α 2 k + m 2 ρ 2 m I ( k , m ) k ! m ! ν 2 k α 2 m c 2 p ρ 2 ( q 1 q 2 ) p k ,
Using Fubini’s theorem and Equation (3.241.411) from [32] in Equation (A8), we obtain
I ( k , m ) = R + y 1 p ν / 2 + a + p k 1 c y 1 q 1 p + 1 ν + α 2 + k + m d y 1 × R + y 2 p ν / 2 + b + p k 1 c y 2 q 2 p + 1 ν + α 2 + k + m d y 2 = Γ ν 2 + a p + k Γ α 2 a p + m Γ ν 2 + b p + k Γ α 2 b p + m p 2 Γ 2 ν + α 2 + k + m × q 1 c p ν 2 + a + p k q 2 c p ν 2 + b + p k .
By combining Equations (A8) and (A9), we obtain
E ( Y 1 a Y 2 b ) = q 1 a q 2 b Γ ν 2 + a p Γ ν 2 + b p Γ α 2 a p Γ α 2 b p c a + b Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ν + α 2 × k = 0 m = 0 ν 2 + a p k ν 2 + b p k α 2 a p m α 2 b p m k ! m ! ν 2 k α 2 k ρ 2 ( k + m ) = q 1 a q 2 b Γ ν 2 + a p Γ ν 2 + b p Γ α 2 a p Γ α 2 b p c a + b Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ν + α 2 × F 1 2 ν 2 + a p , ν 2 + b p ; ν 2 ; ρ 2 F 1 2 α 2 a p , α 2 b p ; α 2 ; ρ 2
Finally, using a Euler transformation, we have Equation (11). □
Proof of Proposition 5.
For part (a), we have
E Y ( Y i p ) = p 2 c p ν ( q i q j ) p ν / 2 Γ 2 ν + α 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 × k = 0 m = 0 ν + α 2 k + m 2 ρ 2 m I ( k , m ) k ! m ! ν 2 k α 2 m c 2 p ρ 2 ( q i q j ) p k .
By using Fubini’s theorem and Equation (3.241.411) from [32] in Equation (A10), we obtain
I ( k , m ) = R + y i p ν / 2 + p k + p 1 c y i q i p + 1 ν + α 2 + k + m d y i × R + y j p ν / 2 + p k 1 c y j q j p + 1 ν + α 2 + k + m d y j = Γ ν 2 + k + 1 Γ α 2 + m 1 Γ ν 2 + k Γ α 2 + m p 2 Γ 2 ν + α 2 + k + m q i c p ν 2 + p k + p q j c p ν 2 + p k .
By combining Equations (A10) and (A11), we obtain
E Y ( Y i p ) = q i p Γ ν 2 + 1 Γ α 2 1 c p Γ ν 2 Γ α 2 ( 1 ρ 2 ) ν + α 2 × k = 0 m = 0 ν 2 + 1 k ν 2 k α 2 1 m α 2 m k ! m ! ν 2 k α 2 k ρ 2 ( k + m ) = q i p Γ ν 2 + 1 Γ α 2 1 c p Γ ν 2 Γ α 2 ( 1 ρ 2 ) ν + α 2 F 1 2 ν 2 + 1 , ν 2 ; ν 2 ; ρ 2 F 1 2 α 2 1 , α 2 ; α 2 ; ρ 2 .
Using the identity F 1 2 ( a , b ; b ; x ) = ( 1 x ) a in the last equality, we obtain
E Y ( Y i p ) = q i p Γ ν 2 + 1 Γ α 2 1 c p Γ ν 2 Γ α 2 .
For part (b), we have
E Y ( log Y i ) = p 2 c p ν ( q i q j ) p ν / 2 Γ 2 ν + α 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 × k = 0 m = 0 ν + α 2 k + m 2 ρ 2 m I ( k , m ) k ! m ! ν 2 k α 2 m c 2 p ρ 2 ( q i q j ) p k .
Using Fubini’s theorem and Equation (3.241.411) from [32] in Equation (A12), we obtain
I ( k , m ) = R + log ( y i ) y i p ν / 2 + p k 1 c y i q i p + 1 ν + α 2 + k + m d y i × R + y j p ν / 2 + p k 1 c y j q j p + 1 ν + α 2 + k + m d y j = Γ 2 ν 2 + k Γ 2 α 2 + m p 3 Γ 2 ν + α 2 + k + m q i c p ν 2 + p k q j c p ν 2 + p k × ψ ν 2 + k ψ α 2 + m log c q i p .
By combining Equations (A12) and (A13), we obtain
E Y ( log Y i ) = 1 p ( 1 ρ 2 ) ( ν + α ) / 2 × k = 0 m = 0 ν 2 k α 2 m ρ 2 k + 2 m k ! m ! ψ ν 2 + k ψ α 2 + m log c q i p .
Proof of Proposition 6.
In evaluating the density of Theorem 2 in Equation (14), we have
H ( Y ) = log ( c p ) 2 q 1 q 2 Γ 2 ν + α 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 p ν 2 1 log c 2 q 1 q 2 + E Y [ log Y 1 ] + E Y [ log Y 2 ] + ν + α 2 E Y log c y 1 q 1 p + 1 + E Y log c y 2 q 2 p + 1 E Y log F 4 ν + α 2 , ν + α 2 ; ν 2 , α 2 ; ρ 2 ( c 2 y 1 y 2 ) p ( q 1 q 2 ) p c y 1 q 1 p + 1 c y 2 q 2 p + 1 , ρ 2 c y 1 q 1 p + 1 c y 2 q 2 p + 1 .
By assuming in the last Appell hypergeometric function that its sum converges at the first term (i.e., when k = m = 0 ) [14], we have that E Y [ log F 4 ( · ) ] 0 . Aside from this, considering Proposition 4, we have that
E Y log c Y i q i p + 1 c q i p E Y [ Y i p ] , i = 1 , 2 .
Then, the differential entropy of Y can be approximated by
H ( Y ) log ( c p ) 2 q 1 q 2 Γ 2 ν + α 2 Γ 2 ν 2 Γ 2 α 2 ( 1 ρ 2 ) ( ν + α ) / 2 p ν 2 1 log c 2 q 1 q 2 + E Y [ log Y 1 ] + E Y [ log Y 2 ] + c p ν + α 2 E Y [ Y 1 p ] q 1 p + E Y [ Y 2 p ] q 2 p .
This concludes the proof. □

References

  1. Johnson, N.; Kotz, S.; Balakrishnan, N. Continuous Univariate Distributions, 2nd ed.; John Wiley & Sons Ltd.: New York, NY, USA, 1994; Volume 1. [Google Scholar]
  2. Arnold, B. Pareto Distributions, 2nd ed.; Chapman and Hall/CRC: Boca Raton, FL, USA, 2015. [Google Scholar]
  3. Feller, W. An Introduction to Probability Theory and Its Applications, 2nd ed.; John Wiley & Sons: Hoboken, NJ, USA, 1971; Volume 2. [Google Scholar]
  4. Pickands, J. Statistical Inference Using Extreme Order Statistics. Ann. Stat. 1975, 3, 119–131. [Google Scholar]
  5. Lomax, K.S. Business failures: Another example of the analysis of failure data. J. Am. Stat. Assoc. 1954, 49, 847–852. [Google Scholar] [CrossRef]
  6. Arnold, B. Pareto Distributions; International Cooperative Publishing House: Fairland, MD, USA, 1983. [Google Scholar]
  7. Cohen, C.A. Truncated and Censored Samples: Theory and Applications, 1st ed.; CRC Press: Boca Raton, FL, USA, 1991. [Google Scholar]
  8. Gupta, R.C.; Akman, O.; Lvin, S. A study of log—Logistic model in survival analysis. Biometrical J. 1999, 41, 431–443. [Google Scholar] [CrossRef]
  9. Bevilacqua, M.; Caamaño-Carrillo, C.; Arellano-Valle, R.B.; Morales-Oñate, V. Non-Gaussian geostatistical modeling using (skew) T processes. Scand. J. Stat. 2021, 48, 212–245. [Google Scholar] [CrossRef]
  10. Bevilacqua, M.; Caamaño-Carrillo, C.; Gaetan, C. On modeling positive continuous data with spatiotemporal dependence. Environmetrics 2020, 31, e2632. [Google Scholar] [CrossRef]
  11. Morales-Navarrete, D.; Bevilacqua, M.; Caamaño-Carrillo, C.; Castro, L.M. Modeling point referenced spatial count data: A Poisson process approach. J. Am. Stat. Assoc. 2022, 119, 1–14. [Google Scholar] [CrossRef]
  12. Kotz, S.; Balakrishnan, N.; Johnson, N.L. Continuous Multivariate Distributions, Volume 1: Models and Applications; John Wiley & Sons: Hoboken, NJ, USA, 2004. [Google Scholar]
  13. Vere-Jones, D. Ergodic properties of nonnegative matrices. I. Pac. J. Math. 1967, 22, 361–386. [Google Scholar] [CrossRef]
  14. Caamaño-Carrillo, C.; Contreras-Reyes, J.E. A generalization of the bivariate gamma distribution based on generalized hypergeometric functions. Mathematics 2022, 10, 1502. [Google Scholar] [CrossRef]
  15. Caamano-Carrillo, C.; Contreras-Reyes, J.E.; González-Navarrete, M.; Sánchez, E. Bivariate superstatistics based on generalized gamma distribution. Eur. Phys. J. B 2020, 93, 43. [Google Scholar] [CrossRef]
  16. Krishnamoorthy, A.S.; Parthasarthy, M.D. A multivariate gamma-type distribution. Ann. Math. Stat. 1951, 22, 549–557. [Google Scholar] [CrossRef]
  17. Krishnaiah, P.R.; Rao, M.M. Remarks on a multivariate gamma distribution. Math. Assoc. Am. 1961, 68, 342–346. [Google Scholar] [CrossRef]
  18. Royen, T. Integral representations and approximations for multivariate gamma distributions. Ann. Inst. Stat. Math. 2004, 59, 499–513. [Google Scholar] [CrossRef]
  19. Kotz, S.; Nadarajah, S. Extreme Value Distributions: Theory and Applications; Imperial College Press: London, UK, 2004. [Google Scholar]
  20. Arnold, B. Univariate and multivariate pareto models. J. Stat. Distrib. Appl. 2014, 1, 1–16. [Google Scholar] [CrossRef]
  21. R Core Team. A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2023; ISBN 3-900051-07-0. Available online: http://www.R-project.org (accessed on 24 July 2024).
  22. Dubey, S.D. Compound gamma, beta and F distributions. Metrika 1970, 16, 27–31. [Google Scholar] [CrossRef]
  23. Bowman, K.O.; Shenton, L.R.; Gailey, P.C. Distribution of the ratio of gamma variates. Commun. Stat. Simul. Comput. 1998, 27, 1–19. [Google Scholar] [CrossRef]
  24. Pham-Gia, T.; Duong, Q.P. The Generalized Beta- and F-distributions in statistical modelling. Math. Comput. Model. 1989, 12, 1613–1625. [Google Scholar] [CrossRef]
  25. Crooks, G. Field Guide to Continuous Probability Distributions. V1.0.0; Berkeley Institute for Theoretical Sciences: Berkeley, CA, USA, 2006. [Google Scholar]
  26. Tiao, G.; Guttman, I. The multivariate inverted beta distribution with applications. J. Am. Stat. Assoc. 1965, 60, 793–805. [Google Scholar] [CrossRef]
  27. Cover, T.M.; Thomas, J.A. Elements of Information Theory; Wiley & Son, Inc.: New York, NY, USA, 2006. [Google Scholar]
  28. Contreras-Reyes, J.E. Mutual information matrix based on Rényi entropy and application. Nonlinear Dyn. 2022, 110, 623–633. [Google Scholar] [CrossRef]
  29. Tahmasebi, S.; Behboodian, J. Shannon entropy for the Feller-Pareto (FP) family and order statistics of FP subfamilies. Appl. Math. Sci. 2010, 4, 495–504. [Google Scholar]
  30. Contreras-Reyes, J.E. Asymptotic form of the Kullback—Leibler divergence for multivariate asymmetric heavy-tailed distributions. Physica A 2014, 395, 200–208. [Google Scholar] [CrossRef]
  31. Contreras-Reyes, J.E. Mutual information matrix based on asymmetric Shannon entropy for nonlinear interactions of time series. Nonlinear Dyn. 2021, 104, 3913–3924. [Google Scholar] [CrossRef]
  32. Gradshteyn, I.; Ryzhik, I. Table of Integrals, Series, and Products, 7th ed.; Academic Press: New York, NY, USA, 2007. [Google Scholar]
Figure 1. Bivariate pdf of Equation (7) for some parameter combinations.
Figure 1. Bivariate pdf of Equation (7) for some parameter combinations.
Axioms 13 00701 g001
Figure 2. Correlation ρ Y of Corollary 1(4) for some parameter combinations.
Figure 2. Correlation ρ Y of Corollary 1(4) for some parameter combinations.
Axioms 13 00701 g002
Figure 3. Approximated mutual information index of Pareto–Feller distribution assuming p = 1 , α = 5 , ν = 4 , (a) ρ = 0.25 , and (b) ρ = 0.75 .
Figure 3. Approximated mutual information index of Pareto–Feller distribution assuming p = 1 , α = 5 , ν = 4 , (a) ρ = 0.25 , and (b) ρ = 0.75 .
Axioms 13 00701 g003
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

Caamaño-Carrillo, C.; Bevilacqua, M.; Zamudio-Monserratt, M.; Contreras-Reyes, J.E. Bivariate Pareto–Feller Distribution Based on Appell Hypergeometric Function. Axioms 2024, 13, 701. https://doi.org/10.3390/axioms13100701

AMA Style

Caamaño-Carrillo C, Bevilacqua M, Zamudio-Monserratt M, Contreras-Reyes JE. Bivariate Pareto–Feller Distribution Based on Appell Hypergeometric Function. Axioms. 2024; 13(10):701. https://doi.org/10.3390/axioms13100701

Chicago/Turabian Style

Caamaño-Carrillo, Christian, Moreno Bevilacqua, Michael Zamudio-Monserratt, and Javier E. Contreras-Reyes. 2024. "Bivariate Pareto–Feller Distribution Based on Appell Hypergeometric Function" Axioms 13, no. 10: 701. https://doi.org/10.3390/axioms13100701

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