Next Article in Journal
Enhancing the Security and Privacy in the IoT Supply Chain Using Blockchain and Federated Learning with Trusted Execution Environment
Next Article in Special Issue
Extropy and Some of Its More Recent Related Measures for Concomitants of K-Record Values in an Extended FGM Family
Previous Article in Journal
Automatic Creation of 3D Documentation in CAD/BIM Based on Topology
Previous Article in Special Issue
The Log-Bimodal Asymmetric Generalized Gaussian Model with Application to Positive Data
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bivariate Unit-Weibull Distribution: Properties and Inference

by
Roger Tovar-Falón
*,†,
Guillermo Martínez-Flórez
*,† and
Luis Páez-Martínez
Departamento de Matemáticas y Estadística, Facultad de Ciencias Básicas, Universidad de Córdoba, Monteria 230002, Colombia
*
Authors to whom correspondence should be addressed.
These authors contributed equally to this work.
Mathematics 2023, 11(17), 3760; https://doi.org/10.3390/math11173760
Submission received: 13 July 2023 / Revised: 12 August 2023 / Accepted: 15 August 2023 / Published: 1 September 2023
(This article belongs to the Special Issue Probability, Statistics & Symmetry)

Abstract

:
In this article, we introduce a novel bivariate probability distribution that is absolutely continuous. Considering the Farlie–Gumbel–Morgenstern (FGM) copula and the unit-Weibull distribution, we can obtain a bivariate unit-Weibull distribution. We evaluate the main properties of the new proposal and use two estimation methods to estimate the parameter for the bivariate probability distribution. A brief Monte Carlo simulation study is conducted to assess the behavior of the employed estimation method and the characteristics of the estimators. Ultimately, as an illustration, a real-life application is presented, demonstrating the utility of the proposal.

1. Introduction

The probability distribution of proportion or rate data has been extensively studied as it is a prevalent finding in research across all domains of knowledge. The beta distribution is perhaps the most widely recognized. Furthermore, many extensions have been made, such as the beta regression model proposed by Ferrari and Cribari-Neto [1], studied later by Branscum et al. [2], and the beta distribution with an excess of zeros and/or ones of Ospina and Ferrari [3], which was extended to the class of inflated beta regression models, see [4], among others. Other outstanding works are those of Paolino [5], Cribari-Neto and Vasconcellos [6], Kieschnick and Mccullough [7] and Vasconcellos and Cribari-Neto [8].
Alternatives to the beta distribution for analysis of bounded data have been considered by other authors. Martínez-Flórez et al. [9] introduced a mixture model between a Bernoulli process and the power normal distribution to model proportions or rates with inflation at zero and/or one values. Given the flexibility of the power normal distribution to fit high degrees of kurtosis and skewness, this proposal is a viable alternative to the beta model. On the other hand, from the logarithmic transformation of the type X = exp ( Y ) , where Y is an absolutely continuous random variable with positive support, it is possible to define new families of adequate distributions to fit data in the interval ( 0 , 1 ) . Such obtained distributions are generically called unit-f, where f is the density to which the transformation is applied. Based on this idea, the unit-Gamma distribution, see [10], the log-Lindley distribution, see [11], and unit-Weibull, unit-Birnbaum–Saunders and unit-Lindley distributions, see work by Mazucheli et al. [12,13,14], have been developed. In addition, Martínez-Flórez et al. [15] introduced a class of distributions to model rates and proportions through an extension of the alpha-power extension of the skew normal distribution. The authors proposed a regression model based on this distribution for bounded responses.
An extension of the unit-Birnbaum–Saunders distribution and its associated regression model was introduced by Martínez-Flórez and Tovar-Falón [16], while another proposed method capable of fitting data between zeros and/or ones (inclusive) with inflation on the zeros/ones values was considered by Martínez-Flórez et al. [17].
Although there are many proposals developed to model proportion or rate data, there are few works that address the problem of this type of response from a multivariate perspective, that is, models capable of jointly modeling two or more responses in the interval ( 0 , 1 ) in such a way that the possible correlations existing in the responses can be considered. In relation to this, some authors have recently made some contributions, among which the works of Lemonte and Moreno-Arenas [18] and more recently Martínez-Flórez et al. [19] stand out.
This paper aims to introduce a novel multivariate distribution to handle responses in the region ( 0 , 1 ) × ( 0.1 ) jointly. The new distribution, which is called the bivariate unit-Weibull distribution is obtained from an extension of the univariate unit-Weibull distribution introduced by Mazucheli et al. [12] by using the Farlie–Gumbel–Morgenstern (FGM) copula, see [20].
In all areas of knowledge, the proposed model can be used to model proportion, rate, or index data. Particularly in the field of engineering, it is common to come across this type of information as a result of investigations, for example, to model the percentage of carbon dioxide and the concentration of ozone in the atmosphere, see [21,22], among others. The proposal can also prove to be advantageous for analyzing data sets from other domains of study, such as the social sciences and civic culture. For instance, it can be utilized to analyze data on laws pertaining to drunk driving and traffic fatalities in 48 states in the United States of America, see [23], or data on the Human Development Index and illiteracy rate, see [18,24].
The main benefits of using this proposal are related to the applicability of a usual regression model, which is the prediction or estimate of the mean response, given the values of the covariates considered in the model. Studying the relationship between the dependent variables and the covariates in such a way means that that this can be used in making decisions that help solve real-life problems.
The univariate unit-Weibull distribution has the follow cumulative distribution function (cdf) and probability density function (pdf), respectively.
F ( y ; α , β ) = exp α ( log y ) β ; 0 < y < 1 ,
where α , β > 0 , and
f ( y ; α , β ) = ( 1 / y ) × α β ( log y ) β 1 exp α ( log y ) β ; 0 < y < 1 ,
where α , β are the scale and shape parameters, respectively. More details and properties of the univariate unit-Weibull distribution can be found in [12].
The following describes the remainder of the article. Section 2 presents the bivariate UW distribution and some of its structural properties are also derived. In Section 3, the statistical inference process is carried out considering two estimation methods. The expected and observed information matrices are also derived. The outcomes of a simulation study and its corresponding discourse are presented in Section 4. Finally, application to real data is provided in Section 5.

2. Bivariate Unit-Weibull Distribution

In this section, a bivariate extension of the unit-Weibull (UW) distribution [12] is proposed. This extension is obtained from the Farlie–Gumbel–Morgenstern (FGM) copula discussed by Gumbel [20] and is called the bivariate unit-Weibull (BVUW) distribution. This new distribution is very attractive, since it allows modeling bivariate data whose response lies on the plane ( 0 , 1 ) 2 .
To introduce the BVUW distribution, we follow the idea of Almetwally et al. [25] using the FMG copula. Thus, according to the theorem of Sklar [26], Definition 1 follows.
Definition 1. 
A random vector X = ( X 1 , X 2 ) is said to have a BVUW distribution, if its joint pdf is given by,
f ( x 1 , x 2 ) = 1 x 1 α 1 β 1 ln x 1 β 1 1 e α 1 ln x 1 β 1 1 x 2 α 2 β 2 ln x 2 β 2 1 e α 2 ln x 2 β 2 × 1 + θ 1 2 e α 1 ln x 1 β 1 1 2 e α 2 ln x 2 β 2 , 0 < x 1 , x 2 < 1 ,
where 1 < θ < 1 and α 1 , α 2 , β 1 , β 2 > 0 .
For the joint pdf given in (3), the notation BVUW ( α 1 , α 2 , β 1 , β 2 , θ ) is used. The picture in Figure 1 shows how the BVUW distribution looks in three dimensions for certain values of α 1 , α 2 , β 1 , β 2 and θ . In the figure, one can observe the different shapes that the UW density can fit, for example, symmetrical and asymmetrical shapes in the shape of a bell or a bathtub.
The joint cdf of the BUW distribution can be written as
F ( x 1 , x 2 ) = e α 1 ln x 1 β 1 e α 2 ln x 2 β 2 1 + θ 1 e α 1 ln x 1 β 1 1 e α 2 ln x 2 β 2
where 1 < θ < 1 ; 0 < x 1 , x 2 < 1 ; α 1 , α 2 , β 1 , β 2 > 0 .
From Definition 1, we have the following theorem.
Theorem 1. 
Let ( X 1 , X 2 ) BVUW ( α 1 , α 2 , β 1 , β 2 , ρ ) , then
1. 
X j UW ( α j , β j ) for j = 1 , 2 .
2. 
The pdf of X 1 , given X 2 = x 2 , is
f ( x 1 x 2 ) = 1 x 1 α 1 β 1 ln x 1 β 1 1 u ( x 1 ) 1 + θ 1 2 u ( x 1 ) 1 2 u ( x 2 )
where u ( x j ) = e α j ln x j β j , for j = 1 , 2 ,
3. 
The pdf of X i , given X 2 = x 2 , is
f ( x i x j ) = 1 x i α i β i ln x i β i 1 e α i ln x i β i 1 + θ 1 2 u ( x i ) 1 2 u ( x j )
where u ( x i ) = e α i ln x i β i , for i = 1 , 2 .
4. 
The cdf X i , given X 2 = x 2 , is
F ( x i x j ) = ( 1 + θ ) u ( x i ) θ e 2 α i ln y i β i 2 θ u ( x i ) u ( x j ) + 2 θ e α j ln y j β j 2 α i ln y i β i .

2.1. Generating Random Variables

The conditional distribution method proposed by Nelsen [27] is employed to generate a random sample from a joint unit-Weibull distribution. It is noteworthy that the joint distribution function can be formulated as follows,
f ( x 1 , x 2 ) = f ( x 1 ) f ( x 2 x 1 )
By utilizing the subsequent procedures, we can generate a bivariate sample by utilizing the conditional approach.
  • From a uniform ( 0 , 1 ) distribution, generate U 1 and U 2 independently of each other.
  • Let X 1 = exp ln U 1 / α 1 1 / β 1 .
  • Allow F ( x 2 x 1 ) = U 2 to determine X 2 through numerical simulation.
  • Repeat Steps 1 to 3 n times to get the result ( x 1 i , x 2 i ) , i = 1 , 2 , , n .

2.2. Product Moments

Proposition 1. 
Let X = ( X 1 , X 2 ) be a random vector with a BVUW ( φ ) distribution, with φ = ( α 1 , α 2 , β 1 , β 1 , θ ) ; then, the rth and sth moments about zero can be computed in the following way
μ r s = E ( X 1 r X 2 s ) = 0 1 0 1 x 1 r x 2 s 1 x 1 α 1 β 1 ln x 1 β 1 1 e α 1 ln x 1 β 1 1 x 2 α 2 β 2 ln x 2 β 2 1 e α 2 ln x 2 β 2 × 1 + θ 1 2 e α 1 ln x 1 β 1 1 2 e α 2 ln x 2 β 2 d x 1 d x 2

2.3. Correlation Coefficient

Let X = ( X 1 , X 2 ) be a random vector with a BVUW ( φ ) distribution, then the correlation coefficient between X 1 and X 2 is given by:
ρ = c o r ( Y 1 , Y 2 ) = E ( Y 1 Y 2 ) E ( Y 1 ) E ( Y 2 ) σ 1 σ 2
where
σ j = E Y j E ( Y j ) 2 for j = 1 , 2 .
By means of a numerical simulation carried out in the R Development Core Team [28] software, it was found that the range of possible values for the correlation coefficient is:
1 3 ρ 1 3
Table 1 shows some values of the correlation coefficient for certain parameters.

2.4. Reliability Function

The joint survival function can be expressed as a copula of its marginal survival functions using the reliability function, where X 1 and X 2 are random variables with survival functions F ¯ ( x 1 ) and F ¯ ( x 2 ) , see Almetwally et al. [25].
Proposition 2. 
The reliability function of the marginal unit-Weibull distributions for the univariate case is given as follows:
R ( x j ) = 1 F ( x j ) = 1 e α j ln x j β j ; 0 < x j < 1 , α , β > 0 , j = 1 , 2 .
Proposition 3. 
The joint survival function for the FGM copula Nelsen [27] is
H ( x 1 , x 2 ) = 1 F ( x 1 ) F ( x 2 ) + C F ( x 1 ) , F ( x 2 ) .
It then follows that the reliability function of the joint FGM distribution for the bivariate UW distribution is
H ( x 1 , x 2 ) = 1 e α 1 ln x 1 β 1 e α 2 ln x 2 β 2 + e α 1 ln x 1 β 1 e α 2 ln x 2 β 2 1 + θ 1 e α 1 ln x 1 β 1 1 e α 2 ln x 2 β 2 .
Proposition 4. 
Basu [29] defines the hazard function for the bivariate case as
h ( x 1 , x 2 ) = f ( x 1 , x 2 ) H ( x 1 , x 2 ) .
Then, the hazard function of the FGM joint BUW distribution is,
h ( x 1 , x 2 ) = 1 x 1 α 1 β 1 ln x 1 β 1 1 e α 1 ln x 1 β 1 1 y 2 α 2 β 2 ln x 2 β 2 1 e α 2 ln x 2 β 2 1 F ( x 1 ) F ( x 2 ) + e α 1 ln x 1 β 1 e α 2 ln x 2 β 2 1 + θ 1 e α 1 ln x 1 β 1 1 e α 2 ln x 2 β 2 × 1 + θ 1 2 e α 1 ln x 1 β 1 1 2 e α 2 ln x 2 β 2 .

3. Statistical Inference for the BUW Distribution

In this section, two distinct estimation techniques are presented that are employed to estimate the parameters of the BVUW distribution, namely, maximum likelihood estimation (MLE) and estimation by inference functions for margins (IFM). These methods were also considered by Almetwally et al. [25].

3.1. Maximum Likelihood Estimation (MLE)

Under this parametric method, the MLE is found to estimate the parameters of the model jointly and in a single step. The log-likelihood function is given by
ln L = i = 1 n [ ln ( f 1 ( x 1 i ) f 2 ( x 2 i ) c ( F 1 ( x 1 i ) , F 2 ( x 2 i ) ; θ ) ) ]
To obtain estimates for each parameter, we use the log-likelihood function. We define
a ( x j ; α j , β j ) = 1 2 F ( x j ; α j , β j ) = 1 2 e α j ln x j β j , j = 1 , 2 .
A bivariate unit-Weibull distribution has a likelihood function defined by
L = α 1 α 2 β 1 β 2 n i = 1 n 1 x 1 i 1 x 2 i ln x 1 i β 1 1 ln x 2 i β 2 1 e i = 1 n α 1 ln x 1 β 1 i = 1 n α 2 ln x 2 β 2 × i = 1 n 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) .
The corresponding log-likelihood function of a BVUW distribution is subsequently provided by
ln L = n ( ln α 1 + ln β 1 + ln α 2 + ln β 2 ) + i = 1 n ln 1 x 1 i + i = 1 n ln 1 x 2 i + ( β 1 1 ) i = 1 n ln ln x 1 i + ( β 2 1 ) i = 1 n ln ln x 2 i i = 1 n α 1 ln x 1 β 1 i = 1 n α 2 ln x 2 β 2 + i = 1 n ln 1 + θ ( a ( y 1 i ; α 1 , β 1 ) ) ( a ( y 2 i ; α 2 , β 2 ) )
To obtain the estimations, we equate the derivatives of the log-likelihood function with respect to each parameter to zero. This is performed by determining the equation score in the following manner:
ln L α 1 = n α 1 i = 1 n ln x 1 i β 1 + i = 1 n 2 θ ( a ( x 2 i ; α 2 , β 2 ) ) ln x 1 i β 1 e α 1 ln x 1 i β 1 ( 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) )
ln L α 2 = n α 2 i = 1 n ln x 2 i β 2 + i = 1 n 2 θ ( a ( x 1 i ; α 1 , β 1 ) ) ln x 2 i β 2 e α 2 ln x 2 i β 2 ( 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) )
ln L β 1 = n β 1 + i = 1 n ln ( ln x 1 i ) i = 1 n α 1 ln x 1 i β 1 ln ( ln x 1 i ) i = 1 n 2 α 1 θ ( a ( x 2 i ; α 2 , β 2 ) ) ln ( ln x 1 i ) ln x 1 i β 1 e α 1 ln x 1 i β 1 ( 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) )
ln L β 2 = n β 2 + i = 1 n ln ( ln x 2 i ) i = 1 n α 1 ln x 2 i β 2 ln ( ln x 2 i ) i = 1 n 2 α 2 θ ( a ( x 1 i ; α 1 , β 1 ) ) ln ( ln x 2 i ) ln x 2 i β 2 e α 2 ln x 2 i β 2 ( 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) )
and
ln L θ = i = 1 n ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) ( 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) )
The MLE φ ^ = ( α ^ 1 , β ^ 1 , α ^ 2 , β ^ 2 , θ ^ ) can be obtained by simultaneously solving the score equations.
It is possible to obtain the MLE φ ^ = ( α ^ 1 , β ^ 1 , α ^ 2 , β ^ 2 , θ ^ ) by solving the score equation,
ln L θ = 0 , ln L β j = 0 , ln L α j = 0 , j = 1 , 2
The estimates are obtained by performing iterative numerical methods such as Newton–Raphson or quasi-Newton methods.

3.2. Information Matrix

In this section, we shall explicate the components of the information matrix for the BUW distribution that have been observed and expected. According to the definition, the elements of the observed information matrix are found using the expression
j φ p , φ q = 2 ln L ( φ ) φ p φ q , φ p , φ q ( α 1 , α 2 , β 1 , β 2 , θ )
where L ( φ ) is the log-likelihood function associated with the parameter vector φ = ( α 1 , α 2 , β 1 , β 2 , θ ) . These elements are presented in detail in Appendix A. The expected information matrix is determined by taking the expected value of the observed information matrix and using the expression.
i φ p , φ q = 1 n E 2 ln L ( φ ) φ p φ q
The above values are obtained numerically by using iterative methods, and therefore the information matrix can also be obtained numerically. Using the theory of large samples, the maximum likelihood estimator φ ^ asymptotically follows a normal distribution with
φ ^ d N 5 φ , Σ φ
where Σ φ = Var ( φ ^ ) = ( i φ p , φ q ) 1 . The approximation to N 5 φ , ( i φ p , φ q ) 1 can be used to construct confidence intervals for α 1 , α 2 , β 1 , β 2 and θ ; these are given by φ p z 1 δ / 2 σ ^ p p for φ p ( α 1 , α 2 , β 1 , β 2 , θ ) , where σ ^ p p is on the diagonal of the matrix Σ φ for each parameter and is the quantile 100 ( δ / 2 ) % of the standard normal distribution.

3.3. Estimation by Inference Functions for Margins (IFM)

This two-step estimation parametric method was proposed by Joe [30]. Each marginal distribution parameter is evaluated separately in the initial stage.
ln L 1 = j = 1 n ln f 1 ( y 1 j , δ 1 ) ; ln L 2 = j = 1 n ln f 2 ( y 2 j , δ 2 )
The subsequent step entails maximizing the log-likelihood function of the copula density by utilizing the maximum likelihood estimates of the marginals F ^ 1 ( x 1 j , δ ) and F ^ 2 ( x 2 j , δ ) to estimate the copula parameter. For a UW marginal distribution, the log-likelihood function is defined as follows.
ln L j = i = 1 n ln 1 x j i + n ln ( α j ) + ln ( β j ) + ( β j 1 ) i = 1 n ln ( ln x j i ) α j i = 1 n ln x j i β j
The maximum likelihood estimators ( α ^ 1 , β ^ 1 , α ^ 2 , β ^ 2 ) can be obtained by simultaneously solving the likelihood equations,
ln L β j = 0 , ln L α j = 0 , j = 1 , 2
That is,
ln L j β j = n β j + i = 1 n ( ln x j i ) α j i = 1 n ln x j i β j ln ( ln x j i ) = 0
ln L j α j = n α j i = 1 n ln x j i β j = 0
then,
F ^ j ( x j ) = e α j ln x j β j .
As per the preceding procedure, the IFM estimation of a bivariate unit-Weibull distribution is defined as follows,
ln L I F M = i = 1 n ln 1 + θ ( 1 2 F ^ 1 ( x 1 i ) ) ( 1 2 F ^ 2 ( x 2 i ) ) .
Based on this, the differentiation of the log-likelihood function with respect to θ thus is
ln L I F M θ = i = 1 n ( a ( x 1 i , α ^ 1 , β ^ 1 ) ) ( a ( x 2 i , α ^ 2 , β ^ 2 ) ) ( 1 + θ ( a ( x 1 i , α ^ 1 , β ^ 1 ) ) ( a ( x 2 i , α ^ 2 , β ^ 2 ) ) ) .
The estimate of the θ parameter is computed numerically by letting
ln L I F M θ = 0
The closed form expression for the maximum likelihood estimator θ ^ is not available, and its computation necessitates numerical execution using a nonlinear optimization algorithm such as optimize or uniroot in R Development Core Team [28].
For the IFM method, the asymptotic variance matrix does not have a closed form. Joe [30] defines the inference functions as,
i = 1 n g ( Y ; η )
where g = ( g 1 , , g m , g d ) , g j = l j α j for j = 1 , , m , corresponding to the equation,
g d = l δ
Thus, the asymptotic variance matrix is given by,
V = ( D g 1 ) M g ( D g 1 ) ,
where M g = Cov ( g ( Y ; η ) ) = E ( gg ) and D g = E g ( Y ; η ) η .
Other approaches can be used to perform the inference on the model in practical situations, for example, Bayesian methods. Nevertheless, in this work, the classical estimation of the parameters is considered using the methods described above. Based on the research conducted by Kundu and Gupta [31], Zhou et al. [32] and Peralta et al. [33], Bayesian methods can be proposed.

4. Monte Carlo Simulation Study

In this particular section, a Monte Carlo simulation study was conducted with the objective of comparing two copula-based estimation techniques, namely maximum likelihood and estimation by inference of marginal functions, in order to determine the BVUW distribution parameters.
For this simulation, the maxLik function of the statistical software [28] was used and data from the BVUW distribution were generated considering the following values of the parameters α 1 , β 1 , α 2 , β 2 and θ :
  • Scenario 1: α 1 = 1.5 , β 1 = 3.5 , α 2 = 0.5 , β 2 = 2.0 , θ = 0.25 .
  • Scenario 2: α 1 = 1.5 , β 1 = 3.5 , α 2 = 0.5 , β 2 = 2.0 , θ = 0.75 .
  • Scenario 3: α 1 = 1.5 , β 1 = 3.5 , α 2 = 0.5 , β 2 = 2.0 , θ = 0.25 .
  • Scenario 4: α 1 = 1.5 , β 1 = 3.5 , α 2 = 0.5 , β 2 = 2.0 , θ = 0.75 .
  • Scenario 5: α 1 = 2.5 , β 1 = 1.5 , α 2 = 1.5 , β 2 = 1 , θ = 0.25 .
  • Scenario 6: α 1 = 2.5 , β 1 = 1.5 , α 2 = 1.5 , β 2 = 1 , θ = 0.75 .
For each scenario, 5000 random samples were generated using the algorithm presented in Section 2.1 for the sample sizes n = 30 , 50 , 75 , 100 , 200 and 500. As quality measures to evaluate the benefits of the MLE and IFM methods, the bias, mean square error (MSE), confidence interval length (CIL) and the interval coverage rate (ICR) were used. These measurements were calculated as
B i a s ( δ ^ ( j ) ) = 1 5000 i = 1 5000 ( δ ^ i ( j ) δ ( j ) ) , M S E ( δ ^ ( j ) ) = 1 5000 i = 1 5000 ( δ ^ i ( j ) δ ) 2
and
C I L ( δ ^ ( j ) ) = 1 5000 i = 1 5000 2 z ( 0.975 ) × e e ( δ ^ ( j ) )
where δ ^ i ( j ) is the estimate of δ ( j ) for the ith sample, and e e ( δ ^ ( j ) ) is the standard error of δ ^ i ( j ) . The ICR is calculated as the proportion of the confidence interval that contain the true value of the parameter below the nominal 95% level. The outcomes of the simulation study are shown in Table A1, Table A2 and Table A3 in Appendix B.
Based on the data presented in the tables, it can be observed that, in general, as the sample size increases, the bias, MSE and CIL of all parameter estimators tend to decrease with the implementation of MLE and IFM methods. In regard to the ICR, the values obtained are very similar to the reference value of 95%. On the other hand, for large n sample sizes, the MLE and IFM methods produce very similar results in terms of bias, MSE and CIL values. Furthermore, both estimation methods have better results in terms of estimations when the copula parameter is not high, that is, for values of the parameter θ in the interval ( 0.6 , 0.6 ) .
The IFM method produces better results when estimating the copula parameter theta. This phenomenon can be attributed in part to the fact that the IFM methodology comprises two estimation steps: the first step involves the estimation of the marginal distribution parameters, followed by the estimation of the copula parameter by incorporating the estimates of the previous marginal distribution parameters.

5. Illustration

This application demonstrates the relevance of the WUB model in a practical context. To accomplish this goal, a collection of facts is considered, focusing on the legislation governing drunk driving and fatalities from traffic collisions in states of the United States of America during the period from 1980 to 2004, as revealed by a thorough investigation conducted by Freeman [23].
The database is accessible within the woolridge library of R Development Core Team [28] software, where it is designated as driving. It contains information related to current legislation, accident records and demographic characteristics. For this particular application, the variables unemployment rate ( x 1 ) and the percent population aged 14 to 24 ( x 2 ) were employed, and the BUW distribution was subsequently adjusted. Likewise, to compare with our proposal, the bivariate Johnson SB (BVSJB) distribution of Lemonte and Moreno-Arenas [18], the bivariate normal (BVN) distribution and the conditional bivariate skew normal (BVSN) distribution of Arnold et al. [34] were studied.
For the purpose of comparison, the Akaike information criterion (AIC) [35] and the Bayesian information criterion (BIC) [36] were employed. These measures are defined in the following:
AIC = 2 ( θ ^ ) + 2 p and BIC = 2 ( θ ^ ) + n log ( n )
The parameter p represents the number of parameters and the log-likelihood function ^ ( · ) was evaluated at the MLEs of the parameters. It is best to choose the model with the smallest AIC or BIC. We used the maxLik function of the statistical package [28] to fit the bivariate model.
The parameters of these models were estimated using the maximum likelihood method. The standard errors are shown in parentheses in Table 2. According to the AIC and BIC criteria, the BVUW distribution provides the best fit.
The graphs in Figure 2 show the contours of the BVN, BVSN and BVUW distributions. We carried out a bivariate Kolmogorov–Smirnov test. For the special case of the BVUW model,
d n ( WUB ) = max 3.1061671 × 10 21 , 5.1055311 × 10 20 = 5.105531 × 10 20 ;
therefore, it is possible to conclude that the BVUW model is well suited to the present data set.

6. Concluding Remarks

Although there are several proposals for fitting proportion, rate or index data, very few focus on the problem of fitting multivariate responses, for example, in the region ( 0 , 1 ) × ( 0 , 1 ) .
In this work, a new bivariate distribution was proposed to model data whose response falls within the region ( 0 , 1 ) × ( 0 , 1 ) . It is possible to obtain the new distribution by extending the univariate unit-Weibull distribution introduced by Mazucheli et al. [12]. For the proposed distribution, the main properties such as the joint density function and the joint distribution function were studied. The moments were also studied and two methods of estimation of the parameters involved in the model were considered. The results of a Monte Carlo-type simulation study showed good asymptotic properties of the parameter estimators and that the two estimation methods are equally efficient. Finally, application to a real data set showed that the proposed model is a viable alternative to analyze proportion, rate or index data in a multivariate context.
Future work could use the results of this proposal to consider how to extend the case of regression models and make inferences about the models from a Bayesian perspective.

Author Contributions

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

Funding

The research of R. Tovar-Falón and G. Martinez-Flórez was supported partially by project: Resolución de Problemas de Situaciones Reales Usando Análisis Estadístico a través del Modelamiento Multidimensional de Tasas y Proporciones; Esquemas de Monitoreamiento para Datos Asimétricos no Normales y una Estrategia Didáctica para el Desarrollo del Pensamiento Lógico-Matemático. Universidad de Córdoba, Colombia, Acta de Compromiso Número FCB-05-19.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Details about available data are given in Section 5.

Acknowledgments

G. Martínez-Flórez and R. Tovar-Falón acknowledge the support given by Universidad de Córdoba, Montería, Colombia.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A. Elements of the Observed Information Matrix

This section presents the elements of the observed information matrix presented in Section 3.
The following elements are obtained using the expression in (7) and after some algebraic simplifications:
j α 1 , α 1 = n α 1 2 + i = 1 n 2 θ ( a ( x 2 ; α 2 , β 2 ) ) ln x 1 i 2 β 1 e α 1 ln x 1 i β 1 ( 1 + θ ( a ( x 1 ; α 1 , β 1 ) ) ( a ( x 2 ; α 2 , β 2 ) ) ( 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) ) 2 + i = 1 n 2 θ ( a ( x 2 ; α 2 , β 2 ) ) ln x 1 i 2 β 1 e α 1 ln x 1 i β 1 2 θ ( a ( x 2 ; α 2 , β 2 ) ) e α 1 ln x 1 i β 1 ( 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) ) 2
j α 1 , α 2 = i = 1 n 4 θ ln x 1 i β 1 e α 1 ln x 1 i β 1 ln x 2 i β 2 e α 2 ln x 2 i β 2 ( 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) ) 2
j α 1 , β 1 = 2 θ ( a ( x 2 i ; α 2 , β 2 ) ) ln x 1 i β 1 ln ( ln x 1 i ) e α 1 ln x 1 i β 1 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) × ( 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) ) 1 α 1 ln x 1 i β 1 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) + 2 θ ( a ( x 2 i ; α 2 , β 2 ) ) ln x 1 i β 1 ln ( ln x 1 i ) e α 1 ln x 1 i β 1 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) × 2 θ α 1 ( a ( x 2 i ; α 2 , β 2 ) ) ln x 1 i β 1 e α 1 ln x 1 i β 1 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) + i = 1 n ln x 1 i β 1 ln ( ln x 1 i )
j α 1 , β 2 = i = 1 n 4 θ α 2 ln x 1 i β 1 e α 1 ln x 1 i β 1 ln ( ln x 2 i ) ln x 2 i β 2 e α 2 ln x 2 i β 2 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) 2
j α 1 , θ = i = 1 n 2 ( a ( x 2 i ; α 2 , β 2 ) ) ln x 1 i β 1 e α 1 ln x 1 i β 1 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) 2
j α 2 , α 2 = n α 2 2 + i = 1 n 2 θ ( a ( x 1 ; α 1 , β 1 ) ) ln x 2 i 2 β 2 e α 2 ln x 2 i β 2 ( 1 + θ ( a ( x 1 ; α 1 , β 1 ) ) ( a ( x 2 ; α 2 , β 2 ) ) 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) 2 + i = 1 n 2 θ ( a ( x 1 ; α 1 , β 1 ) ) ln x 2 i 2 β 2 e α 2 ln x 2 i β 2 2 θ ( a ( x 2 ; α 2 , β 2 ) ) e α 1 ln x 1 i β 1 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) 2
j α 2 , β 1 = i = 1 n 4 θ α 1 ln x 2 i β 2 e α 2 ln x 2 i β 2 ln ( ln x 1 i ) ln x 1 i β 1 e α 1 ln x 1 i β 1 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) 2
j α 2 , β 2 = i = 1 n 2 θ ( a ( x 1 i ; α 1 , β 1 ) ) ln x 2 i β 2 ln ( ln x 2 i ) e α 2 ln x 2 i β 2 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) × ( 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) ) 1 α 2 ln x 2 i β 2 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) + i = 1 n 2 θ ( a ( x 1 i ; α 1 , β 1 ) ) ln x 2 i β 2 ln ( ln x 2 i ) e α 2 ln x 2 i β 2 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) × 2 θ α 2 ( a ( x 1 i ; α 1 , β 1 ) ) ln x 2 i β 2 e α 2 ln x 2 i β 2 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) + i = 1 n ln x 2 i β 2 ln ( ln x 2 i )
j α 2 , θ = i = 1 n 2 ( a ( x 1 i ; α 1 , β 1 ) ) ln x 2 i β 2 e α 2 ln x 2 i β 2 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) 2
j β 1 , β 1 = i = 1 n 2 α 1 θ ( a ( x 2 i ; α 2 , β 2 ) ) ln ( ln x 1 i ) 2 ln x 1 i β 1 e α 1 ln x 1 i β 1 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) × 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) 1 α 1 ln x 1 i β 1 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) + i = 1 n 4 α 1 2 θ 2 ( a ( x 2 i ; α 2 , β 2 ) 2 ) ln ( ln x 1 i ) 2 ln x 1 i 2 β 1 e 2 α 1 ln x 1 i β 1 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) 2 + i = 1 n α 1 ln ( ln x 1 i ) 2 ln x 1 i β 1 + n β 1 2
j β 1 , β 2 = i = 1 n 2 α 1 θ ln ( ln x 1 i ) ln x 1 i β 1 e α 1 ln x 1 i β 1 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) 2
j β 1 , θ = i = 1 n 2 α 1 ( a ( x 2 i ; α 2 , β 2 ) ) ln ( ln x 1 i ) ln x 1 i β 1 e α 1 ln x 1 i β 1 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) 2
j β 2 , β 2 = i = 1 n 2 α 2 θ ( a ( x 1 i ; α 1 , β 1 ) ) ln ( ln x 2 i ) 2 ln x 2 i β 2 e α 2 ln x 2 i β 2 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) × 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) 1 α 2 ln x 2 i β 2 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) + i = 1 n 4 α 2 2 θ 2 ( a ( x 1 i ; α 1 , β 1 ) 2 ) ln ( ln x 2 i ) 2 ln x 2 i 2 β 2 e 2 α 2 ln x 2 i β 2 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) 2 + i = 1 n α 2 ln ( ln x 2 i ) 2 ln x 2 i β 2 + n β 2 2
j β 2 , θ = i = 1 n 2 α 2 ( a ( x 1 i ; α 1 , β 1 ) ) ln ( ln x 2 i ) ln x 2 i β 2 e α 2 ln x 2 i β 2 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) 2
j θ , θ = i = 1 n ( a ( x 1 i ; α 1 , β 1 ) 2 ) ( a ( x 2 i ; α 2 , β 2 ) 2 ) 1 + θ ( a ( x 1 i ; α 1 , β 1 ) ) ( a ( x 2 i ; α 2 , β 2 ) ) 2

Appendix B. Simulation Tables

Table A1. Quality measure (QM) of the estimates obtained under the MLE and IMF methods: Scenarios 1 and 2.
Table A1. Quality measure (QM) of the estimates obtained under the MLE and IMF methods: Scenarios 1 and 2.
Scenario 1MLE MethodIFM Method
n QM α ^ 1 β ^ 1 α ^ 2 β ^ 2 θ ^ α ^ 1 β ^ 1 α ^ 2 β ^ 2 θ ^
30Bias0.04050.16800.00480.0845−0.06470.04080.16760.00440.0854−0.0701
MSE0.09070.32650.01680.12110.19560.09100.32290.01660.11950.2005
CIL1.11732.06570.47601.17422.07771.11572.06450.47491.17362.0586
ICR0.95490.95420.92460.93860.96890.95730.95420.92450.93900.9741
50Bias0.02350.09060.00290.0543−0.03200.02480.09260.00290.0543−0.0342
MSE0.04820.16720.00970.06170.15350.04890.16820.00970.06140.1558
CIL0.85011.56140.36780.89161.61330.85081.56240.36780.89171.6055
ICR0.95500.95430.93310.93380.94510.95330.95330.93310.93930.9428
75Bias0.01890.05510.00480.0296−0.01810.01970.05520.00480.0296−0.0206
MSE0.03230.10620.00630.03750.11210.03260.10670.00620.03730.1113
CIL0.69041.25910.30110.71761.32450.69091.25960.30130.71771.3186
ICR0.94480.95500.93730.93660.94410.94550.95440.94080.93740.9415
100Bias0.01190.04050.00260.0222−0.01540.01230.04000.00290.0219−0.0166
MSE0.02390.08080.00470.02600.08690.02390.08080.00470.02590.0859
CIL0.59441.08510.26010.61891.15320.59481.08530.26040.61901.1485
ICR0.94770.95100.94030.95440.94770.94500.94900.93890.95370.9477
200Bias0.00790.02580.00080.0083−0.00070.00800.02580.00090.0081−0.0013
MSE0.01190.03830.00220.01250.04390.01190.03830.00220.01250.0437
CIL0.41840.76300.18360.43420.81820.41870.76330.18380.43430.8160
ICR0.94670.95330.94530.95400.93930.94730.95270.94670.95270.9387
500Bias0.00250.0063−0.00080.0062−0.00100.00250.0062−0.00080.0061−0.0012
MSE0.00440.01460.00090.00520.01720.00440.01460.00090.00520.0172
CIL0.26340.47940.11600.27430.51920.26360.47960.11600.27440.5183
ICR0.95070.95470.94670.94270.95470.95130.95270.94800.94200.9540
Scenario 2
30Bias0.05930.18010.00290.0908−0.28920.06250.19200.00340.0916−0.2993
MSE0.09970.31860.01590.11330.22480.10490.33750.01650.11850.2341
CIL1.13012.07030.47401.17711.97241.13292.07810.47441.17782.0189
ICR0.95680.95680.92900.94240.93020.95570.95470.92410.93890.9505
50Bias0.02910.10070.00120.0637−0.18230.02630.10210.00060.0637−0.1826
MSE0.04950.16780.00930.06440.12460.04950.17030.00940.06510.1251
CIL0.85101.56510.36580.89481.52550.85171.56820.36690.89621.5299
ICR0.95450.95160.93700.93020.94860.95360.95360.93660.92900.9612
75Bias0.02290.06790.00260.0417−0.12630.02380.06730.00300.0404−0.1264
MSE0.03210.10700.00600.03780.07770.03260.10770.00610.03810.0783
CIL0.69001.26220.29900.72031.22940.69271.26440.30060.72151.2356
ICR0.94860.95930.93980.94860.95750.94770.95820.94250.94420.9747
100Bias0.01720.05220.00150.0302−0.10070.01720.05230.00150.0297−0.1016
MSE0.02400.08270.00470.02660.05760.02430.08300.00470.02660.0575
CIL0.59461.08750.25860.62001.07060.59671.08970.25990.62111.0661
TCR0.94520.94860.93510.95030.96460.94320.94910.93070.95240.9766
200Bias0.00950.03250.00040.0101−0.03680.01020.03130.00050.0099−0.0363
MSE0.01200.03910.00220.01260.02770.01220.03900.00220.01260.0278
CIL0.41690.76260.18230.43340.73980.41930.76470.18370.43470.7382
ICR0.94630.95530.94330.95300.96050.94530.95420.94240.95490.9749
500Bias0.00230.0078−0.00090.0063−0.00680.00260.0070−0.00080.0061−0.0052
MSE0.00440.01450.00090.00520.01250.00440.01470.00090.00520.0129
CIL0.26190.47820.11510.27350.46670.26360.47970.11600.27440.4621
TCR0.95630.95490.94940.94600.97060.95060.95260.94380.94450.9675
Table A2. Quality measure of the estimates obtained under the MLE and IMF methods: Scenario 3 and 4.
Table A2. Quality measure of the estimates obtained under the MLE and IMF methods: Scenario 3 and 4.
Scenario 3MLE MethodIFM Method
n QM α ^ 1 β ^ 1 α ^ 2 β ^ 2 θ ^ α ^ 1 β ^ 1 α ^ 2 β ^ 2 θ ^
30Bias0.03740.16990.00680.08260.08780.04290.17240.00740.08210.0927
MSE0.08680.32880.01670.11870.20610.09400.32990.01680.11820.2156
CIL1.11482.06700.47771.17322.08151.11762.06670.47671.17152.0661
ICR0.95270.95660.93100.93800.97440.95330.95630.93370.93520.9721
50Bias0.02250.09000.00430.05170.04000.02140.09030.00470.05110.0450
MSE0.04870.16880.00990.06020.15090.04860.16820.00980.05980.1514
CIL0.84961.56050.36850.89061.61700.84891.56060.36860.89041.6097
ICR0.95320.95460.93620.93760.94890.95190.95340.93850.93500.9519
75Bias0.01750.05190.00450.03000.01580.01760.05180.00440.03030.0138
MSE0.03280.10630.00630.03700.11330.03260.10620.00620.03700.1136
CIL0.68961.25780.30090.71811.32170.68991.25820.30110.71841.3177
ICR0.94200.95560.94130.93650.94270.94500.95520.94090.93820.9402
100Bias0.01140.03930.00310.0214−0.00510.01130.03920.00310.0212−0.0044
MSE0.02390.08070.00470.02570.08500.02380.08090.00470.02580.0844
CIL0.59401.08470.26030.61891.14870.59431.08510.26050.61911.1441
ICR0.94890.94820.94080.95090.95090.94560.94890.94150.94890.9469
200Bias0.00790.02580.00100.0082−0.00050.00800.02580.00090.00830.0001
MSE0.01190.03830.00220.01240.04460.01190.03830.00220.01240.0444
CIL0.41840.76300.18370.43420.81820.41870.76330.18380.43440.8158
ICR0.94930.95130.95070.95070.94070.94730.95270.95400.95270.9400
500Bias0.00250.0062−0.00080.0060−0.00130.00250.0062−0.00080.0061−0.0011
MSE0.00440.01460.00090.00510.01730.00440.01460.00090.00510.0172
CIL0.26340.47940.11600.27430.51930.26360.47960.11600.27440.5184
ICR0.95130.95400.95070.94470.95730.95130.95270.95130.94470.9573
Scenario 4
30Bias0.04020.17420.00830.08990.29110.04630.18640.00930.09260.2995
MSE0.08440.32410.01660.11440.21600.09880.34810.01720.11690.2258
CIL1.11482.07070.47771.17621.99781.12072.07540.47781.17742.0209
ICR0.95880.96310.92620.93600.93930.95080.95590.92310.93640.9579
50Bias0.02350.10470.00440.05440.19410.02710.10170.00490.05500.1966
MSE0.04960.17120.00990.05910.12920.05150.17080.01030.06030.1295
CIL0.84881.56590.36770.89181.52090.85231.56610.36870.89261.5376
ICR0.95530.95530.93290.93770.94360.94780.95730.93070.94020.9658
75Bias0.01670.06580.00270.03650.13150.01980.05980.00210.03730.1324
MSE0.03250.11360.00600.03790.08130.03290.11300.00610.03810.0809
CIL0.68721.26060.29920.71981.21480.69101.26100.30030.72151.2337
ICR0.94000.94720.94630.94180.93460.94350.94960.94260.94170.9664
100Bias0.01140.04950.00250.02450.09050.01040.04950.00160.02610.0895
MSE0.02380.08420.00450.02500.05440.02380.08500.00450.02530.0543
CIL0.59181.08560.25880.61921.04810.59401.08850.26010.62131.0537
ICR0.95030.94350.94090.95720.95630.94520.94270.94350.95360.9755
200Bias0.00840.02970.00130.01000.03870.00890.02870.00110.01000.0399
MSE0.01200.03860.00220.01250.02850.01200.03830.00220.01250.0285
CIL0.41660.76170.18250.43370.73750.41890.76380.18390.43500.7376
ICR0.94820.95120.94970.94740.96840.94610.95280.95060.95130.9813
500Bias0.00240.0071−0.00080.00630.00400.00230.0069−0.00070.00610.0036
MSE0.00440.01470.00090.00510.01270.00440.01480.00090.00510.0129
CIL0.26190.47810.11510.27350.46500.26350.47970.11610.27440.4621
ICR0.95260.95600.94640.94440.95470.95170.95230.95100.94620.9598
Table A3. Quality measure of the estimates obtained under the MLE and IMF methods: Scenario 5 and 6.
Table A3. Quality measure of the estimates obtained under the MLE and IMF methods: Scenario 5 and 6.
Scenario 5MLE MethodIFM Method
n QM α ^ 1 β ^ 1 α ^ 2 β ^ 2 θ ^ α ^ 1 β ^ 1 α ^ 2 β ^ 2 θ ^
30Bias0.14080.06870.06650.0435−0.06830.14190.07160.06530.0424−0.0720
MSE0.35770.05720.11520.03020.19450.36230.05920.11250.02990.1997
CIL2.09700.88351.13900.58802.08342.09620.88481.13620.58662.0627
ICR0.96640.95780.94450.93910.97660.96330.95410.94420.93880.9755
50Bias0.07800.03910.04200.0274−0.03100.08040.03970.04290.0271−0.0353
MSE0.17160.03090.06130.01550.15400.17400.03090.06110.01540.1553
CIL1.56120.66930.86160.44601.61601.56250.66960.86200.44581.6061
ICR0.95920.95290.94240.93390.94870.95880.95320.94550.93920.9441
75Bias0.05520.02390.03130.0149−0.01760.05590.02370.03190.0148−0.0206
MSE0.10940.01970.03700.00940.11230.11010.01960.03700.00930.1113
CIL1.25550.53970.69630.35881.32531.25580.53980.69670.35891.3186
ICR0.96320.95370.94690.93670.94420.95920.95440.94620.93740.9415
100Bias0.03670.01740.02060.0110−0.01490.03740.01720.02120.0109−0.0166
MSE0.08000.01480.02680.00650.08720.08030.01480.02670.00650.0859
CIL1.07530.46510.59800.30941.15321.07570.46510.59840.30951.1484
ICR0.96310.95100.93700.95440.94840.96240.94900.93620.95370.9477
200Bias0.02370.01110.00650.0041−0.00120.02380.01100.00650.0041−0.0013
MSE0.03930.00700.01180.00310.04360.03930.00700.01180.00310.0437
CIL0.75320.32700.41800.21710.81830.75350.32710.41830.21720.8160
ICR0.94530.95330.94600.95400.94000.94800.95270.94870.95270.9387
500Bias0.00690.00270.00150.0031−0.00100.00690.00270.00150.0031−0.0012
MSE0.01470.00270.00460.00130.01720.01470.00270.00460.00130.0172
CIL0.47160.20550.26330.13720.51920.47170.20550.26340.13720.5183
ICR0.95330.95470.95070.94270.95470.95330.95270.94930.94200.9533
Scenario 6
30Bias0.17750.07470.06460.0428−0.28930.19010.08340.06640.0470−0.2994
MSE0.38540.05790.11210.02810.22210.42100.06150.11140.02980.2341
CIL2.13270.88581.13640.58722.00012.14760.89121.13690.58952.0130
ICR0.97010.95900.93690.94350.93910.96840.95680.93990.93890.9515
50Bias0.08540.04330.04420.0308−0.17910.08580.04350.04370.0320−0.1838
MSE0.16870.03060.05520.01570.12000.17260.03120.05590.01620.1253
CIL1.56430.67080.86000.44711.53511.56690.67210.86190.44821.5338
ICR0.96390.95130.95610.93180.95810.96300.95260.95540.92980.9649
75Bias0.06690.02970.03560.0200−0.12690.06720.02880.03590.0202−0.1264
MSE0.10820.01970.03650.00950.07760.10910.01980.03660.00950.0783
CIL1.26070.54120.69590.36001.24391.26280.54190.69840.36071.2354
ICR0.96540.95920.95040.94590.96900.96160.95820.95380.94420.9738
100Bias0.04860.02230.02440.0150−0.09830.05000.02240.02340.0148−0.1014
MSE0.08070.01530.02660.00670.05760.08230.01520.02680.00660.0575
CIL1.08010.46600.59720.31001.06761.08270.46700.59920.31061.0655
ICR0.95820.94650.94060.94980.96320.96080.94910.93910.95250.9750
200Bias0.03070.01390.00680.0049−0.03610.02940.01340.00710.0048−0.0359
MSE0.04020.00720.01180.00320.02770.04020.00720.01190.00310.0278
CIL0.75410.32680.41600.21670.74130.75570.32770.41840.21730.7375
ICR0.94800.95540.95460.95240.96730.94540.95430.95500.95500.9742
500Bias0.00700.00330.00110.0030−0.00590.00730.00290.00140.0030−0.0051
MSE0.01470.00270.00450.00130.01260.01470.00270.00460.00130.0129
CIL0.47020.20490.26170.13670.46660.47180.20560.26340.13720.4620
ICR0.95440.95580.94830.94550.97410.95270.95200.94660.94460.9655

References

  1. Ferrari, S.; Cribari-Neto, F. Beta regression for modelling rates and proportions. Appl. Stat. 2004, 31, 799–815. [Google Scholar] [CrossRef]
  2. Branscum, A.J.; Johnson, W.O.; Thurmond, M.C. Bayesian beta regression: Applications to household expenditure data and genetic distance between foot-and-mouth diseases viruses. Aust. N. Z. J. Stat. 2007, 49, 287–301. [Google Scholar] [CrossRef]
  3. Ospina, R.; Ferrari, S.L.P. Inflated beta distribution. Stat. Pap. 2010, 51, 111–126. [Google Scholar] [CrossRef]
  4. Ospina, R.; Ferrari, S.L.P. A general class of zero-or-one inflated beta regression models. Comput. Stat. Data Anal. 2012, 56, 1609–1623. [Google Scholar] [CrossRef]
  5. Paolino, P. Maximum likelihood estimation of models with beta-distributed dependent variables. Political Anal. 2001, 9, 325–346. [Google Scholar] [CrossRef]
  6. Cribari-Neto, F.; Vasconcellos, K.L.P. Nearly unbiased maximum likelhood estimation for the beta distribution. J. Stat. Comput. Simul. 2002, 72, 107–118. [Google Scholar] [CrossRef]
  7. Kieschnick, R.; Mccullough, B.D. Regression analysis of variates observed on (0, 1). Stat. Model. 2003, 3, 193–213. [Google Scholar] [CrossRef]
  8. Vasconcellos, K.L.P.; Cribari-Neto, F. Improved maximum likelihood estimation in a new class of beta regression models. Braz. J. Probab. Stat. 2005, 19, 13–31. [Google Scholar]
  9. Martínez-Flórez, G.; Bolfarine, H.; Gómez, H.W. Power-models for proportions with zero/one excess. Appl. Math. Inf. Sci. 2018, 12, 293–303. [Google Scholar] [CrossRef]
  10. Grassia, A. On a family of distributions with argument between 0 and 1 obtained by transformation of the Gamma distribution and derived compound distributions. Aust. J. Stat. 1977, 19, 108–114. [Google Scholar] [CrossRef]
  11. Gómez-Déniz, E.; Sordo, M.A.; Calderín-Ojeda, E. The Log-Lindley distribution as an alternative to the Beta regression model with applications in insurance. Insur. Math. Econ. 2013, 54, 49–57. [Google Scholar] [CrossRef]
  12. Mazucheli, J.; Menezes, A.F.B.; Ghitany, M.E. The unit-Weibull distribution and associated inference. J. Appl. Probab. Stat. 2018, 13, 1–22. [Google Scholar]
  13. Mazucheli, J.; Menezes, A.; Dey, S. The unit-Birnbaum-Saunders distribution with applications. Chil. J. Stat. 2018, 9, 47–57. [Google Scholar]
  14. Mazucheli, J.; Menezes, A.F.B.; Chakraborty, S. On the one parameter unit-Lindley distribution and its associated regression model for proportion data. J. Appl. Stat. 2019, 49, 700–714. [Google Scholar] [CrossRef]
  15. Martínez-Flórez, G.; Azevedo-Farias, R.B.; Tovar-Falón, R. New Class of Unit-Power-Skew-Normal Distribution and Its Associated Regression Model for Bounded Responses. Mathematics 2022, 10, 3035. [Google Scholar] [CrossRef]
  16. Martínez-Flórez, G.; Tovar-Falón, R. New Regression Models Based on the Unit-Sinh-Normal Distribution: Properties, Inference, and Applications. Mathematics 2021, 9, 1231. [Google Scholar] [CrossRef]
  17. Martínez-Flórez, G.; Gómez, H.W.; Tovar-Falón, R. Modeling Proportion Data with Inflation by Using a Power-Skew-Normal/Logit Mixture Model. Mathematics 2021, 9, 1989. [Google Scholar] [CrossRef]
  18. Lemonte, A.J.; Moreno-Arenas, G. On a multivariate regression model for rates and proportions. J. Appl. Stat. 2019, 46, 1084–1106. [Google Scholar] [CrossRef]
  19. Martínez-Flórez, G.; Lemonte, A.J.; Moreno-Arenas, G.; Tovar-Falón, R. The Bivariate Unit-Sinh-Normal Distribution and Its Related Regression Model. Mathematics 2022, 10, 3125. [Google Scholar] [CrossRef]
  20. Gumbel, E.G. Bivariate exponential distributions. J. Am. Stat. Assoc. 1960, 292, 698–707. [Google Scholar] [CrossRef]
  21. Gokhale, S.; Khare, M. Vehicle wake factor for heterogeneous traffic in urban environments. Int. J. Environ. Pollut. 2007, 30, 97–105. [Google Scholar] [CrossRef]
  22. Martínez-Flórez, G.; Azevedo-Farias, R.B.; Tovar-Falón, R. An Exponentiated Multivariate Extension for the Birnbaum-Saunders Log-Linear Model. Mathematics 2022, 10, 1299. [Google Scholar] [CrossRef]
  23. Freeman, D.G. Drunk Driving Legislation and Traffic Fatalities: New Evidence on BAC 08 Laws. Contemp. Econ. Policy 2007, 25, 293–308. [Google Scholar] [CrossRef]
  24. Martínez-Flórez, G.; Vergara-Cardozo, S.; Tovar-Falón, R.; Rodriguez-Quevedo, L. The Multivariate Skewed Log-Birnbaum–Saunders Distribution and Its Associated Regression Model. Mathematics 2023, 11, 1095. [Google Scholar] [CrossRef]
  25. Almetwally, E.M.; Muhammed, H.Z.; El-Sherpieny, E.S.A. Bivariate Weibull Distribution: Properties and Different Methods of Estimation. Ann. Data. Sci. 2020, 7, 163–193. [Google Scholar] [CrossRef]
  26. Sklar, A. Random variables, joint distributions, and copulas. Kybernetica 1973, 9, 449–460. [Google Scholar]
  27. Nelsen, R.B. An Introduction to Copulas, 2nd ed.; Springer: New York, NY, USA, 2010. [Google Scholar]
  28. R Development Core Team. R: A Language and Environment for Statistical Computing; R Foundation for Statistical Computing: Vienna, Austria, 2021; Available online: http://www.R-project.org (accessed on 31 July 2021).
  29. Basu, A. Bivariate failure rate. J. Am. Stat. Assoc. 1971, 66, 103–104. [Google Scholar] [CrossRef]
  30. Joe, H. Asymptotic efficiency of the two-stage estimation method for copula-based models. J. Multivar. Anal. 2005, 94, 401–419. [Google Scholar] [CrossRef]
  31. Kundu, D.; Gupta, A.K. Bayes estimation for the Marshall-Olkin bivariate Weibull distribution. Comput. Stat. Data Anal. 2013, 57, 271–281. [Google Scholar] [CrossRef]
  32. Zhou, S.; Xu, A.; Tang, Y.; Shen, L. Fast Bayesian Inference of Reparameterized Gamma Process With Random Effects. IEEE Trans. Reliab. 2023, 1–14. [Google Scholar] [CrossRef]
  33. Peralta, D.; DeOliveira, R.P.; Achcar, J.A. A hierarchical Bayesian analysis for bivariate Weibull distribution under left-censoring scheme. J. Appl. Stat. 2023, 1–20. [Google Scholar] [CrossRef]
  34. Arnold, B.C.; Castillo, E.; Sarabia, J.M. Conditionally specified multivariate skewed distributions. Sankhya Indian J. Stat. Ser. A 2002, 64, 206–226. [Google Scholar]
  35. Akaike, H. A new look at statistical model identification. IEEE Trans. Autom. Control 1974, 19, 716–722. [Google Scholar] [CrossRef]
  36. Schwarz, G. Estimating the dimension of a model. Ann. Stat. 1978, 6, 461–464. [Google Scholar] [CrossRef]
Figure 1. Joint probability density function of the BVUW distribution with parameter vector φ = ( α 1 , α 2 , β 1 , β 2 , θ ) for (a) φ = ( 1.0 , 1.5 , 0.5 , 2.5 , 0.75 ) (b) φ = ( 1.5 , 1.5 , 2.0 , 2.0 , 0.75 ) , (c) φ = ( 2.0 , 2.0 , 3.5 , 3.5 , 0.25 ) , (d) φ = ( 1.0 , 1.0 , 0.5 , 0.5 , 0.25 ) , (e) φ = ( 2.5 , 2.5 , 1.0 , 1.0 , 0.75 ) , (f) φ = ( 2.5 , 2.5 , 1.0 , 2.5 , 0.75 ) .
Figure 1. Joint probability density function of the BVUW distribution with parameter vector φ = ( α 1 , α 2 , β 1 , β 2 , θ ) for (a) φ = ( 1.0 , 1.5 , 0.5 , 2.5 , 0.75 ) (b) φ = ( 1.5 , 1.5 , 2.0 , 2.0 , 0.75 ) , (c) φ = ( 2.0 , 2.0 , 3.5 , 3.5 , 0.25 ) , (d) φ = ( 1.0 , 1.0 , 0.5 , 0.5 , 0.25 ) , (e) φ = ( 2.5 , 2.5 , 1.0 , 1.0 , 0.75 ) , (f) φ = ( 2.5 , 2.5 , 1.0 , 2.5 , 0.75 ) .
Mathematics 11 03760 g001
Figure 2. Contour plots for the fitted models. (a) BVN, (b) BVSN and (c) BVUW.
Figure 2. Contour plots for the fitted models. (a) BVN, (b) BVSN and (c) BVUW.
Mathematics 11 03760 g002
Table 1. Correlation coefficient values for some selected parameter values.
Table 1. Correlation coefficient values for some selected parameter values.
α 1 α 2 β 1 β 2 θ ρ
3.02.01.01.0−0.75−0.235
3.02.01.01.0−0.25−0.078
3.02.01.01.00.250.078
3.02.01.01.00.750.235
1.52.51.01.0−0.99−0.317
1.52.51.01.00.500.160
4.51.51.01.0−0.99−0.307
4.51.51.01.00.990.307
7.02.01.01.0−0.90−0.269
7.02.01.01.0−0.50−0.150
7.02.01.01.0−0.10−0.030
3.51.81.01.00.900.281
2.0121.02.5−0.80−0.258
3.51.81.81.50.990.325
0.93.01.01.0−0.990−0.315
Table 2. MLE (SE) for the fitted models.
Table 2. MLE (SE) for the fitted models.
ParametersBVSJBBVNBVSNBVUW
α 1 0.1554(0.0010)0.1532(0.0005)0.1546(0.0005)18.0497(0.3906)
α 2 0.0624(0.0018)0.0595(0.0005)0.0602(0.0005)9.3641(0.1949)
β 1 5.0188(0.9010) 1.9351(0.0032)
β 2 2.7763(0.1020) 3.0170(0.0096)
σ 1 1.7204(0.6972)0.0003(0.00001)0.0188(0.0003)
σ 2 8.0839(0.2439)0.0004(0.00001)0.0207(0.0004)
λ 3.2156(0.5746)0.0001(0.00001)0.6973(0.0594)0.9990(0.0038)
AIC−10,980.07−12,250.64−12,245.92−12,447.81
BIC−10,944.44−12,225.19−12,220.47−12,422.36
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

Tovar-Falón, R.; Martínez-Flórez, G.; Páez-Martínez, L. Bivariate Unit-Weibull Distribution: Properties and Inference. Mathematics 2023, 11, 3760. https://doi.org/10.3390/math11173760

AMA Style

Tovar-Falón R, Martínez-Flórez G, Páez-Martínez L. Bivariate Unit-Weibull Distribution: Properties and Inference. Mathematics. 2023; 11(17):3760. https://doi.org/10.3390/math11173760

Chicago/Turabian Style

Tovar-Falón, Roger, Guillermo Martínez-Flórez, and Luis Páez-Martínez. 2023. "Bivariate Unit-Weibull Distribution: Properties and Inference" Mathematics 11, no. 17: 3760. https://doi.org/10.3390/math11173760

APA Style

Tovar-Falón, R., Martínez-Flórez, G., & Páez-Martínez, L. (2023). Bivariate Unit-Weibull Distribution: Properties and Inference. Mathematics, 11(17), 3760. https://doi.org/10.3390/math11173760

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