Next Article in Journal
Hamiltonian Formulation for Continuous Systems with Second-Order Derivatives: A Study of Podolsky Generalized Electrodynamics
Next Article in Special Issue
Constrained Bayesian Method for Testing Equi-Correlation Coefficient
Previous Article in Journal
Exact Solutions to Fractional Schrödinger–Hirota Equation Using Auxiliary Equation Method
Previous Article in Special Issue
Research on Change Point Detection during Periods of Sharp Fluctuations in Stock Prices–Based on Bayes Method β-ARCH Models
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Bivariate Length-Biased Exponential Distribution under Progressive Type-II Censoring: Incorporating Random Removal and Applications to Industrial and Computer Science Data

1
Department of Statistics, Faculty of Science, King Abdulaziz University, Jeddah 21589, Saudi Arabia
2
Department of Mathematics and Statistics, Faculty of Science, Imam Mohammad Ibn Saud Islamic University (IMSIU), Riyadh 11432, Saudi Arabia
3
Department of Statistics, Faculty of Business Administration, Delta University for Science and Technology, Gamasa 11153, Egypt
4
Department of Statistics, Mathematics and Insurance, Benha University, Benha 13511, Egypt
*
Authors to whom correspondence should be addressed.
Axioms 2024, 13(10), 664; https://doi.org/10.3390/axioms13100664
Submission received: 10 July 2024 / Revised: 15 September 2024 / Accepted: 22 September 2024 / Published: 26 September 2024
(This article belongs to the Special Issue Applications of Bayesian Methods in Statistical Analysis)

Abstract

:
In this paper, we address the analysis of bivariate lifetime data from a length-biased exponential distribution observed under Type II progressive censoring with random removals, where the number of units removed at each failure time follows a binomial distribution. We derive the likelihood function for the progressive Type II censoring scheme with random removals and apply it to the bivariate length-biased exponential distribution. The parameters of the proposed model are estimated using both likelihood and Bayesian methods for point and interval estimators, including asymptotic confidence intervals and bootstrap confidence intervals. We also employ different loss functions to construct Bayesian estimators. Additionally, a simulation study is conducted to compare the performance of censoring schemes. The effectiveness of the proposed methodology is demonstrated through the analysis of two real datasets from the industrial and computer science domains, providing valuable insights for illustrative purposes.

1. Introduction

In numerous life test studies, accurately recording the lifetimes of test units can be challenging. An experimenter may choose to end the life test prematurely to save time or reduce costs. Consequently, the test is deemed censored, meaning the collected data includes precise failure times for the failed units and running times for the non-failed units. For example, accidental breakage of some experimental units or subjects dropping out for personal reasons during clinical trials can occur midway through the trial. In certain life-testing studies involving expensive units, it may be advantageous to remove some units early from the test, allowing their use in other tests. Additionally, in some clinical trials, the survival duration after treatment may span several years, prompting the experimenter to conclude the study before observing the survival durations for all individuals. Censored data emerges in these scenarios when the experimenter does not acquire complete information for all units or individuals under study.
Various types of censoring arise depending on the data collection methodology in a life-testing experiment. Let’s examine a life-testing experiment involving n items under examination. Assume that the experiment must be concluded at a predetermined time, denoted as T. In this scenario, one can only acquire failure times that are less than or equal to T, resulting in what is known as Type-I censored data. Alternatively, instead of predetermining the total time of the experiment, the experimenter may decide to conclude the experiment after observing the first r failures. In such cases, the data are referred to as Type-II censored. The two most common censoring schemes, Type-I and Type-II, have a limitation in that systems can only be removed from the test at the final termination time. However, in many practical applications, systems need to be removed at different points during the test, which is known as progressive censoring. Progressive censoring allows flexibility in the removal of systems throughout the testing process, beyond just the final termination time. For further details, refer to Balakrishnan and Aggarwala [1], Balakrishnan [2], Balakrishnan and Cramer [3], and El-Sherpieny et al. [4].
The copula function has garnered considerable attention as a contemporary area of study in the realm of statistics. Its significance lies in its ability to serve as a powerful and efficient tool for modeling bivariate and multivariate distributions. This versatility finds extensive practical applications across diverse disciplines, including economics, finance, hydrology, computer science, sports, psychology, social sciences, climatology, meteorology, environmental studies, and medical data analysis. A copula is a type of multivariate distribution function in which the individual margins are uniformly distributed on the interval [0, 1]. In this context, we focus solely on the bivariate copula distribution. To explore copulas and their properties further, comprehensive information can be found in the works of Nelsen [5]. Recent contributions to bivariate models incorporating copula functions and their applications can be found in Flores [6], Suzuki et al. [7], Verril et al. [8], Peres et al. [9], Peres et al. [10], Samanthi et al. [11], Elgawad et al. [12], Barakat et al. [13], Qura et al. [14], Fayomi et al. [15,16], Chesneau et al. [17], and Almetwally et al. [18], among others.
In the domain of bivariate censoring, Balakrishnan and Kim [19] laid the foundation by introducing the likelihood function for a bivariate distribution, specifically focusing on Type-II censoring. Subsequently, Balakrishnan and Kim [20,21] extended their work by presenting the likelihood function based on a Progressive Type-II censoring scheme for a bivariate distribution. Kim et al. [22] examined the use of a Type-II censoring scheme in a bivariate context for the estimation of unknown parameters in the bivariate generalized exponential distribution. Bai et al. [23] developed Type-I progressive interval censoring schemes for the Marshall-Olkin bivariate Weibull (MOBW) distribution, focusing on how percentile bootstrap methods impact the construction of confidence intervals for unknown parameters. Aly et al. [24] derived the likelihood function for progressive Type-I censoring and applied it to the Marshall-Olkin bivariate Kumaraswamy distribution. Muhammed [25] derived the likelihood function for a bivariate Weibull distribution, which is part of the bivariate hazard power parameter family of distributions, based on complete, Type-I, and Type-II censored samples.
In the context of progressive Type-II censored samples, EL-Morshedy et al. [26] have explored both maximum likelihood and Bayesian methods for estimating parameters of the bivariate Burr X-G family with various distributions. Their discussion is specifically centered around Type-II censored data.El-Sherpieny et al. [27] acquired Bayesian and non-Bayesian estimations for the parameters of the bivariate generalized Rayleigh distribution, utilizing the Clayton copula function within the framework of progressive Type-II censoring schemes. El-Sherpieny et al. [4] investigated the point and interval estimation of unknown parameters in the Farlie-Gumbel-Morgenstern (FGM) bivariate Weibull distribution based on progressive Type-II censored samples. Additionally, the study introduced two bootstrap confidence intervals. Muhamed [28] developed a generalized likelihood function for the Marshal-Olkin bivariate class of distributions under progressive Type II censoring, applied to the bivariate Dagum distribution, and focused on maximum likelihood estimation for the unknown parameters.
In this paper, we extend the work of Arun et al. [29] by considering the estimation of the bivariate length-biased exponential distribution using the Farlie-Gumbel-Morgenstern copula (FGM-BLBE) parameters under a progressive Type-II censoring scheme with random removal, based on the maximum likelihood and Bayesian estimation methods. Furthermore, the interval estimation of the FGM-BLBE parameters is obtained using asymptotic and bootstrap confidence intervals under the progressive Type-II censoring scheme. The performance of the proposed estimators is evaluated through extensive simulations. The optimal censoring scheme is proposed based on three distinct optimality criteria: mean squared error (MSE), bias, and length of the confidence interval (L.CI). Ultimately, this effort aims to establish a comprehensive guideline for identifying the most effective censoring scheme, providing valuable insights for statisticians. Applications of two real datasets are introduced to confirm the validity of the model.
The rest of this paper is organized as follows. The description and formulation of the model are provided in Section 2. The maximum likelihood estimation (MLE), along with approximate confidence intervals, is presented in Section 3. Section 4 covers Bayesian estimates and the corresponding highest posterior density (HPD) credible intervals. In Section 5, the construction of confidence intervals for unknown parameters uses both asymptotic and bootstrap methods. In Section 6, the Monte Carlo simulation results and data analysis are presented. Section 7 introduces and analyzes two real datasets to evaluate the model. Finally, concluding remarks are given in Section 8.

2. Description and Formulation of the Model

Consider random variables Y 1 and Y 2 with their joint cumulative distribution function (CDF) denoted as F Y 1 , Y 2 ( y 1 , y 2 ) and their joint probability density function (PDF) as f Y 1 , Y 2 ( y 1 , y 2 ) . Additionally, let F Y 1 ( y 1 ) and F Y 2 ( y 2 ) represent the marginal CDFs of Y 1 and Y 2 , respectively, and their marginal PDFs denoted as f Y 1 ( y 1 ) and f Y 2 ( y 2 ) , respectively. Sklar [30] introduced the concept of a copula as a linking function that connects the joint distribution function to its marginal distribution functions using the relationship F Y 1 , Y 2 ( y 1 , y 2 ) = P ( Y 1 y 1 , Y 2 y 2 ) = C ( F Y 1 ( y 1 ) , F Y 2 ( y 2 ) ) , and the associated joint density is f Y 1 , Y 2 ( y 1 , y 2 ) = f Y 1 ( y 1 ) f Y 2 ( y 2 ) c ( F Y 1 ( y 1 ) , F Y 2 ( y 2 ) ) , where c is the copula density.
In the literature, various types of copulas have been extensively discussed, each possessing distinct properties and characteristics. One prominent example is the Farlie-Gumbel-Morgenstern copula, which belongs to a widely recognized parametric copula family introduced by Gumbel [31]. The joint CDF and PDF of the Farlie-Gumbel-Morgenstern (FGM) copula are provided as follows:
C ( u , v ) = u v 1 + θ ( 1 u ) ( 1 v ) , 1 θ 1 ,
and
c ( u , v ) = 1 + θ ( 1 2 u ) ( 1 2 v ) , 1 θ 1 ,
where u = F Y 1 ( y 1 , Ψ 1 ) , v = F Y 2 ( y 2 , Ψ 2 ) , and Ψ 1 and Ψ 2 are vectors of the parameters for the variables Y 1 and Y 2 , respectively. The parameter θ serves as a dependence parameter that relies on the underlying random variables, and when set to θ = 0 , it represents the independent case, signifying complete absence of dependence between the variables. The FGM copula is thus simple, flexible, and possesses tail dependence, allowing for parametric control over the level of dependence. Furthermore, its analytical properties are advantageous, and its historical significance as one of the earliest introduced copula families contributes to the development and understanding of copula theory.
Dara and Ahmad [32] introduced the length-biased exponential (LBE) model, also referred to as the moment exponential (ME) model, by assigning weights to the exponential (E) model. Their research demonstrated that the LBE distribution offers greater adaptability compared to the E model. The CDF, probability density function (PDF) and hazard rate function (HRF) of the LBE model with scale parameter λ are as follows:
F Y ( y ; λ ) = 1 1 + y λ e y λ , y > 0 , λ > 0 ,
f Y ( y ; λ ) = y λ 2 e y λ , y > 0 , λ > 0 ,
and
h Y ( y ; λ ) = y λ 2 1 + y λ .
Arun et al. [29] introduced the bivariate length-biased exponential distribution using the Farlie-Gumbel-Morgenstern approach. This distribution, denoted as the FGM-BLBE distribution, serves as a bivariate counterpart to the length-biased exponential distribution.
The CDF, PDF, and HRF of an FGM-BLBE distribution are obtained as follows:
F Y 1 , Y 2 ( y 1 , y 2 ) = 1 1 + y 1 λ 1 e y 1 λ 1 1 1 + y 2 λ 2 e y 2 λ 2 1 + θ 1 + y 1 λ 1 e y 1 λ 1 1 + y 2 λ 2 e y 2 λ 2 , y 1 , y 2 > 0 ,   λ 1 , λ 2 > 0 ,   1 θ 1 ,
f Y 1 , Y 2 ( y 1 , y 2 ) = y 1 λ 1 2 e y 1 λ 1 y 2 λ 2 2 e y 2 λ 2 1 + θ 2 1 + y 1 λ 1 e y 1 λ 1 1 2 1 + y 2 λ 2 e y 2 λ 2 1 ,
and
h Y 1 , Y 2 ( y 1 , y 2 ) = f Y 1 , Y 2 ( y 1 , y 2 ) S Y 1 , Y 2 ( y 1 , y 2 ) ,
where
S Y 1 , Y 2 ( y 1 , y 2 ) = 1 + y 1 λ 1 e y 1 λ 1 1 + y 2 λ 2 e y 2 λ 2 1 + θ 1 1 + y 1 λ 1 e y 1 λ 1 1 1 + y 2 λ 2 e y 2 λ 2 ,
S is the survival function (SF), which is expressed as S Y 1 , Y 2 ( y 1 , y 2 ) = C S Y 1 ( y 1 ) , S Y 2 ( y 2 ) .
The hazard rate function quantifies the probability of an event (such as failure) occurring per unit time, given that the system has survived up until that moment. By analyzing variations in the hazard rate over time or across different values of x and y, we can evaluate the dynamic behavior and risk profile of the FGM-BLBE distribution, as illustrated in Figure 1 and Figure 2.
Consider a random sample y 1 1 : n , y 2 1 : n , y 1 2 : n , y 2 2 : n , , y 1 n : n , y 2 n : n drawn from a bivariate distribution with cumulative distribution function F Y 1 , Y 2 y 1 , y 2 and probability density function f Y 1 , Y 2 y 1 , y 2 . When arranging the values y 1 i , where i = 1 , 2 , , n , in ascending order, we obtain y 1 1 : n < y 1 2 : n < < y 1 n : n . Here, the term y 2 [ i : n ] represents the concomitant of the i t h order statistic. The consideration of these concomitants becomes crucial when choosing bivariate and prediction techniques that rely on the rankings of Y 1 . To delve deeper into the topic of concomitants of order statistics, you can refer to the work of David and Nagaraja [33] and Ke Wang [34].
Under the progressive type-II censoring scheme, the following assumptions will be employed and can be described as follows:
  • Suppose there are n independent units undergoing life testing, and these units have the FGM-BLBE distribution.
  • Generate progressive samples as ( y 1 1 : m : n , y 2 [ 1 : m : n ] ) < ( y 1 2 : m : n , y 2 [ 2 : m : n ] ) < < ( y 1 m : m : n , y 2 [ m : m : n ] ) , subject to a progressive censoring scheme denoted as R i , where i = 1 , 2 , , m (which is less than or equal to n) and m is predetermined as the number of complete failures to be observed.
  • Identify the m completely observed failure times as Y 1 i : m : n ( R i ) , where i = 1 , , m . At the time of the first failure, ( y 1 1 : m : n , y 2 [ 1 : m : n ] ) , R 1 follows a binomial ( n m , p ) distribution, representing the number of units randomly removed from the remaining ( n 1 ) surviving items. At the time of the second failure, ( y 1 2 : m : n , y 2 [ 2 : m : n ] ) , R 2 follows a binomial ( n m R 1 , p ) distribution, indicating the number of units randomly removed from the remaining ( n 2 R 1 ) units, and so on. The process continues until the m-th failure occurs at any time, at which point all the remaining ( n m R 1 R 2 R m 1 ) units are removed.
  • Assume that removing an individual unit from Y 1 during the test is independent of the others, and all units have the same removal probability, denoted as p. In such a scenario, the number of units removed at each failure time follows a binomial distribution. Then, the distribution of units removed is denoted as R i binomial ( n m j = 1 m 1 R j , p ) , and R m = n m j = 1 m 1 R j for j = 2 , 3 , , m 1 .
The full likelihood function, represented as L ( y 1 i : m : n , y 2 [ i : m : n ] , ) , based on a progressive type-II censored sample with binomial removals, is defined for any vector of parameters ∇, as:
L ( y 1 i : m : n , y 2 [ i : m : n ] , ) = L 1 ( y 1 i : m : n , y 2 [ i : m : n ] , ) P R = r .
Balakrishnan and Kim [20] define the likelihood function L 1 ( y 1 i : m : n , y 2 [ i : m : n ] ) for the case of progressive type II censored samples with R removal as follows:
L 1 ( y 1 i : m : n , y 2 [ i : m : n ] , ) = C i = 1 m f Y 1 , Y 2 ( y 1 i : m : n , y 2 [ i : m : n ] ) 1 F Y 1 ( y 1 i : m : n ) R i ,
where C is a constant which does not depend on vector parameter ∇ and is determined by C = n n R 1 1 n R 1 R 2 2 n R 1 R 2 R m 1 m + 1 . The likelihood function of the binomial distribution is given by:
P ( R = r ) = P R 1 = r 1 , R 2 = r 2 , , R m 1 = r m 1 , = n m ! n m j = 1 m 1 r j ! j = 1 m 1 r j p j = 1 m 1 r j 1 p ( m 1 ) ( n m ) j = 1 m 1 ( m j ) r j ,
where p denotes a fixed constant probability, it is determined by the experimenter’s level of experience regarding the number of removals.
We observe that, as L 1 ( y 1 i : m : n , y 2 [ i : m : n ] , ) does not rely on the binomial parameter p, the Maximum Likelihood Estimator (MLE) of p can be derived directly by maximizing Equation (12). Hence, the MLE of p follows as
p ^ = i = 1 m 1 r i ( m 1 ) ( n m ) i = 1 m 1 ( m i 1 ) r i .
Similarly, P ( R = r ) does not depend on the parameter vector ∇, then the MLE of ∇ can be readily obtained by maximizing L 1 ( y 1 i : m : n , y 2 [ i : m : n ] , ) .

3. MLE of the FGM-BLBE Distribution under Progressive Type-II Censoring

In this section, we focus on the estimation of FGM-BLBE parameters in the presence of progressive Type-II censored samples, employing the Maximum Likelihood Estimators (MLEs).
Let Y 1 1 : m : n , Y 2 1 : m : n , Y 1 2 : m : n , Y 2 2 : m : n , , Y 1 m : m : n , Y 2 m : m : n represent a progressively Type-II censored sample with random removal from the FGM-BLBE distribution, characterized by the PDF in Equation (7), and defined by the parameter vector = ( λ 1 , λ 2 , θ ) . The log-likelihood function is structured as follows:
( ) = 2 m ln ( λ 1 ) + ln ( λ 2 ) + i = 1 m ln ( y 1 i : m : n ) + i = 1 m ln ( y 2 [ i : m : n ] ) 1 λ 1 i = 1 m y 1 i : m : n 1 λ 2 i = 1 m y 2 [ i : m : n ] + i = 1 m ln 1 + θ 2 ϕ ( y 1 i : m : n , λ 1 ) 1 2 η ( y 2 [ i : m : n ] , λ 2 ) 1 + i = 1 m R i ln ϕ ( y 1 i : m : n , λ 1 ) ,
where
ϕ ( y 1 i : m : n , λ 1 ) = 1 + y 1 i : m : n λ 1 e y 1 i : m : n λ 1 ,
and
η ( y 2 [ i : m : n ] , λ 2 ) = 1 + y 2 [ i : m : n ] λ 2 e y 2 [ i : m : n ] λ 2 .
The partial derivatives of Equation (14) with respect to the unknown parameters can be computed as follows:
( ) λ 1 = 2 m λ 1 + 1 λ 1 2 i = 1 m y 1 i : m : n + i = 1 m 2 θ 2 η ( y 2 [ i : m : n ] , λ 2 ) 1 B 1 i : m : n ( y 1 i : m : n , λ 1 ) 1 + θ 2 ϕ ( y 1 i : m : n , λ 1 ) 1 2 η ( y 2 [ i : m : n ] , λ 2 ) 1 + i = 1 m R i B 1 i : m : n ( y 1 i : m : n , λ 1 ) ϕ ( y 1 i : m : n , λ 1 ) ,
( ) λ 2 = 2 m λ 2 + 1 λ 2 2 i = 1 m y 2 [ i : m : n ] + i = 1 m 2 θ 2 ϕ ( y 1 i : m : n , λ 1 ) 1 B 2 [ i : m : n ] ( y 2 [ i : m : n ] , λ 2 ) 1 + θ 2 ϕ ( y 1 i : m : n , λ 1 ) 1 2 η ( y 2 [ i : m : n ] , λ 2 ) 1 ,
and
( ) θ = i = 1 m 2 ϕ ( y 1 i : m : n , λ 1 ) 1 2 η ( y 2 [ i : m : n ] , λ 2 ) 1 1 + θ 2 ϕ ( y 1 i : m : n , λ 1 ) 1 2 η ( y 2 [ i : m : n ] , λ 2 ) 1 ,
where
B 1 i : m : n ( y 1 i : m : n , λ 1 ) = λ 1 ϕ ( y 1 i : m : n , λ 1 ) = y 1 i : m : n 2 λ 1 3 e y 1 i : m : n λ 1 ,
and
B 2 [ i : m : n ] ( y 2 [ i : m : n ] , λ 2 ) = λ 2 η ( y 2 [ i : m : n ] , λ 2 ) = y 2 [ i : m : n ] 2 λ 2 3 e y 2 [ i : m : n ] λ 2 .
The MLE of the parameters in the FGM-BLBE model, denoted by ^ , under the progressive Type-II censoring scheme is obtained by solving a set of nonlinear equations derived from the likelihood function. Due to the complexity of these equations, finding an analytical solution is often infeasible. As a result, numerical techniques, such as the Newton-Raphson method, are employed to approximate the optimal values of the parameters.

4. Bayesian Estimation

In this segment, we explore the Bayesian inference of undisclosed parameters in the context of latent failure times following the FGM-BLBE distribution, utilizing a progressive Type-II censoring scheme. Bayesian estimates (BEs) and Highest Posterior Density (HPD) credible intervals were computed for the undisclosed parameters.

4.1. Prior Distributions

In this subsection, we made the assumption that the prior distributions of λ 1 , and λ 2 , are independent and are characterized by Gamma distributions with parameters ( a 1 , b 1 ) , and ( a 2 , b 2 ) , respectively. The copula parameter θ is assumed to follow a uniform distribution over the interval ( 1 , 1 ) . Consequently, the prior distributions for λ 1 , λ 2 , and θ take the following forms:
π ( λ 1 ) = b 1 a 1 Γ a 1 λ 1 a 1 1 e b 1 λ 1 λ 1 > 0 a n d a 1 , b 1 > 0 ,
π ( λ 2 ) = b 2 a 2 Γ a 2 λ 2 a 2 1 e b 2 λ 2 λ 2 > 0 a n d a 2 , b 2 > 0 ,
and
π ( θ ) = 1 2 1 < θ < 1 ,
then, the joint prior is given by
π ( λ 1 , λ 2 , θ ) = π ( λ 1 ) π ( λ 2 ) π ( θ ) = 1 2 b 1 a 1 Γ a 1 λ 1 a 1 1 e b 1 λ 1 b 2 a 2 Γ a 2 λ 2 a 2 1 e b 2 λ 2 .

4.2. Posterior Distribution

Using the combined equations of prior and likelihood, the posterior density function for the parameters λ 1 , λ 2 , and θ in the context of the FGM-BLBE model, employing a progressive Type-II censoring scheme, can be expressed in the following manner:
π ( λ 1 , λ 2 , θ | x ) = L ( y 1 , y 2 | λ 1 , λ 2 , θ ) π ( λ 1 , λ 2 , θ ) 0 0 0 L ( y 1 , y 2 | λ 1 , λ 2 , θ ) π ( λ 1 , λ 2 , θ ) d λ 1 d λ 2 d θ λ 1 2 m + a 1 λ 2 2 m + a 2 e λ 1 i = 1 m y 1 i : m : n λ 1 2 + b 1 e λ 2 i = 1 m y 2 i : m : n λ 2 2 + b 2 i = 1 m ϕ ( y 1 i : m : n , λ 1 ) R i i = 1 m 1 + θ 2 1 + y 1 i : m : n λ 1 e y 1 i : m : n λ 1 1 2 1 + y 2 i : m : n λ 2 e y 2 i : m : n λ 2 1 .
The Bayesian estimator for any function involving variables λ 1 , λ 2 , and θ , denoted as ζ ( λ 1 , λ 2 , θ ) within the context of the square error loss (SEL) function, is represented by the posterior mean. This is expressed as ζ ˜ ( λ 1 , λ 2 , θ ) and can be derived in the following manner:
ζ ˜ ( λ 1 , λ 2 , θ ) = 0 0 0 ζ ( λ 1 , λ 2 , θ ) L ( y 1 , y 2 | λ 1 , λ 2 , θ ) π ( λ 1 , λ 2 , θ ) d λ 1 d λ 2 d θ 0 0 0 L ( y 1 , y 2 | λ 1 , λ 2 , θ ) π ( λ 1 , λ 2 , θ ) d λ 1 d λ 2 d θ .
The Bayesian estimator for ζ ( λ 1 , λ 2 , θ ) lacks a closed-form expression, necessitating the utilization of an approximation method to compute the estimate presented in Equation (20). We suggest employing the Markov Chain Monte Carlo (MCMC) method to obtain Bayesian estimators (BEs) and Highest Posterior Density (HPD) credible intervals for the unknown parameters. This method proves particularly advantageous in Bayesian inference due to its focus on posterior distributions that often pose challenges for mathematical analysis. The Metropolis-Hastings (MH) algorithm is employed, initiating with the generation of a candidate sample ζ * from the proposal distribution q ( . ) . It’s crucial to note that samples from the proposal distribution aren’t automatically accepted as posterior samples; rather, acceptance is determined probabilistically based on the acceptance probability.
The MH sampling algorithm, introduced by Metropolis et al. [35] and further developed by Hastings (1970), can be elucidated in the following manner:
  • Step 1: Commence by initializing any initial estimate ∇, denoted as ζ ( 0 ) = ( λ ^ 1 , λ ^ 2 , θ ^ ) .
  • Step 2: Set the iteration index as t = 1 .
  • Step 3: Generate a candidate point ζ * using a normal proposal distribution q ( ζ ) = N ( ζ ^ , var ( ζ ^ ) ) .
  • Step 4: For the given candidate point ζ * , compute the acceptance probability.
    A ( ζ ( t 1 ) , ζ * ) = min 1 , π ( ζ * | x ̲ ) q ( ζ ( t 1 ) | ζ * ) π ( ζ ( t 1 ) | x ̲ ) q ( ζ * | ζ ( t 1 ) ) .
  • Step 5: Generate a sample from a uniform distribution, i.e., u U ( 0 , 1 ) .
    If u A ( ζ ( t 1 ) , ζ * ) A c c e p t ζ * = ζ ( t ) , u A ( ζ ( t 1 ) , ζ * ) A c c e p t ζ ( t ) = ζ ( t 1 ) .
  • Step 6: Increment the iteration index: t = t + 1 , and repeat steps 2–5 M times until obtaining M samples, resulting in ( λ 1 ( t ) , λ 2 ( t ) , θ ( t ) ) for t = 1 , 2 , M .
From the random samples of size M drawn from the posterior density, a subset of initial samples may be discarded (burn-in), and the remaining samples can be utilized to calculate Bayesian estimators (BEs). The BEs for ζ concerning the SEL function are given by:
ζ ˜ = t = + 1 M ζ ( t ) ( M ) .
Here, ⋎ represents the number of burn-in samples, and the summation is performed over t = + 1 , , M .
  • Step 7: Under the LINEX loss function, the BEs of ζ are obtained as:
ζ = 1 c ln t = + 1 M e c ζ ( t ) ( M ) .

5. Confidence Intervals

In this section, two different approaches for constructing confidence intervals (CIs) for the unknown parameters of the FGM-BLBE model under progressive Type-II censoring are introduced. These techniques encompass the asymptotic confidence interval (ACI) and the bootstrap confidence interval, which is subdivided into percentile bootstrap and bootstrap-t methods for λ j , where j = 1 , 2 , and θ .

5.1. Asymptotic Confidence Intervals

One of the most widely used approaches for establishing confidence intervals for parameters involves leveraging the asymptotic normality of the maximum likelihood estimator (MLE). Specifically, this approach relies on the asymptotic variance-covariance matrix of the MLE of the parameters, denoted as the Fisher information matrix I ( ) . This matrix is computed from the negative second derivatives of the natural logarithm of the likelihood function, evaluated at the estimated parameter values ^ = ( λ ^ 1 , λ ^ 2 , θ ^ ) . The asymptotic variance–covariance matrix of the parameter vector ∇ can be represented as
I ( ^ ) = I λ ^ 1 λ ^ 1 I λ ^ 2 λ ^ 1 I λ ^ 2 λ ^ 2 I θ ^ λ ^ 1 I θ ^ λ ^ 2 I θ ^ θ ^ ,
where V ( ^ ) = I 1 ( ^ ) . Confidence intervals for parameter ∇ can be constructed based on the asymptotic normality of the MLE. Specifically, a 100 ( 1 α ) % confidence interval for parameter ∇ can be calculated as λ ^ l ± Z 0.025 I λ ^ l λ ^ l and θ ^ ± Z 0.025 I θ ^ θ ^ , where l = 1 , 2 and Z 0.025 represents the percentile of the standard normal distribution with a right tail probability of α 2 .

5.2. Highest Posterior Density

To obtain the Highest Posterior Density (HPD) credible intervals for ζ , arrange the samples ζ t = ( θ t , λ 1 t , λ 2 t ) for t = 1 , 2 , M in ascending order, as ζ [ 1 ] , ζ [ 2 ] , ζ [ M ] after burn-in, specifically as ζ [ + 1 ] , ζ [ + 2 ] , ζ [ M ] . For any arbitrary 0 < α < 1 , the 100 ( 1 α ) % two-sided HPD credible interval of ζ can be computed as:
ζ [ ( α / 2 ) ( M ) ] , ζ [ ( 1 ( α / 2 ) ) ( M ) ] .
Subsequently, the HPD credible interval of ζ t is determined by selecting the interval with the smallest width.

6. Simulation

The simulation section in statistical analysis is a powerful tool for estimating the parameters of a distribution. Simulations involve generating artificial data based on the FGM-BLBE distribution, allowing researchers to explore the behavior of these models and estimate the associated parameters. This method is particularly useful when analytical solutions are challenging or unavailable. The primary goals of the simulation section are to assess the performance of statistical methods (MLE and Bayesian with different loss functions), evaluate the robustness of parameter estimation techniques, and gain insights into the potential variability of results under different scenarios. By mimicking the underlying data generation process, simulations provide a controlled environment for testing the statistical properties of estimation procedures. The objectives of this section are as follows:
  • Assess the accuracy and precision of parameter estimates of FGM-BLBE distribution based on progressive Type-II Censoring.
  • Compare the performance of various estimation methods and conclude the best estimation method for parameter of FGM-BLBE distribution based on progressive Type-II Censoring.
  • Investigate the impact of sample size, size of censored sample and distributional assumptions on estimation results.

Simulation Design

The simulation design of this study assesses the effectiveness of theoretical outcomes, encompassing both point and interval estimators, through different estimation methods such as maximum likelihood and Bayesian methodologies based on symmetric and asymmetric loss functions. A Monte Carlo simulation study is conducted, employing initial parameter values of the FGM-BLBE distribution, specifically ( λ 1 , λ 2 , θ ) = (2.4, 2.8, 0.4), (0.6, 2.8, 0.4), (0.6, 1.4, 0.4), and ( 1.5 , 2 , 0.8 ) . By varying combinations of n (sample sizes), m (effective sample size), and p (parameter of the binomial removal pattern), a comprehensive set of 5000 progressively Type-II censored samples is generated, incorporating binomial removal from the FGM-BLBE distribution.
Moreover, varying values of ( n , m , p ) are considered. Within each configuration, the MLEs, Bayesian estimators (BEs), and their associated ACIs or HPD intervals are examined and assessed. Additionally, various combinations of ( n , m ) such as (50, 35), (50, 45), (100, 75), and (100, 90) are considered, with fixed values of p = 0.3 , 0.5 , and 0.8 . For each configuration, the MLEs, Bayesian estimates, and their associated asymptotic confidence intervals (ACIs) or HPD intervals are assessed with a confidence level of 95%.
In the simulation study, it’s important to highlight that the effectiveness of the suggested point estimates is evaluated by assessing biases, and mean squared errors (MSEs). As for the proposed interval estimates, their performance is evaluated by considering average interval lengths (AILs) (LACI and LCCI, for MLE and Bayesian, respectively) and coverage probability (CP).
In the simulation study, Bayesian estimators (BEs) are computed using the SELF and LINEX functions, specifically for c values of 0.5 and 2. It is essential to note that the hyperparameter values are determined through an elicitation method, where the choice of the hyperparameter method relies on prior knowledge about the unknown parameters. These informative priors are derived from the MLEs for ( λ 1 , λ 2 ) , obtained by equating the mean and variance of ( λ ^ 1 , λ ^ 2 ) to the mean and variance of the considered priors (Gamma priors). Here, t ranges from 1 to M, where M is the number of samples available from the FGM-BLBE distribution. By equating the mean and variance of ( λ 1 t , λ 2 t ) with those of the Gamma priors, the hyperparameters can be determined, as described by Dey et al. [36].
Based on the generated data, Maximum Likelihood Estimates (MLEs) and their associated 95% ACIs or HPD are computed. It’s noteworthy that, when obtaining MLEs, the initial guess values are assumed to be identical to the true parameter values. Subsequently, the hyperparameter values are derived. These hyperparameter values are then incorporated to calculate the desired estimates. Finally, utilizing the MH algorithm for BEs, 2000 burn-in samples are excluded from the total 12,000 MCMC samples generated from the posterior density. Table 1, Table 2, Table 3 and Table 4 provide insights into these observations as follows:
  • BEs for parameters are calculated using two distinct loss functions, namely, the SEL and LINEX functions. Estimates associated with the LINEX function are obtained for c values of 0.5 and 2.
  • All average estimates and their corresponding MSEs for both methods are presented.
  • In the case of fixed censoring schemes (with identical parameters for binomial removal), as the effective size increases (i.e., n or m, or their combinations), the Bias, MSEs, and AILs of both MLEs and BEs for parameters of FGM-BLBE decrease.
  • For fxed value of n, when effective sample size m increases, the simulated MSEs decreases for most cases.
  • Based on Bias and MSEs, the Bayes estimates under SEL and LINEX provide better results than other estimates for MLEs.
  • As the sample size (n) and effective sample size (m) increase, the AILs for intervals related to the parameters of FGM-BLBE decrease.
  • The HPD credible intervals outperform the confidence intervals of the MLEs based on AIL and Coverage Probability (CP).
  • In terms of CP, the confidence intervals of the MLEs exhibit an average coverage probability of 95The performance of both classical and Bayesian estimates is deemed quite satisfactory.
  • Notably, the performance of BEs is better under the LINEX loss function compared to the SEL function.
Table 1. MLE and Bayesian: λ 1 = 2.4 ,   λ 2 = 2.8 ,   θ = 0.4 .
Table 1. MLE and Bayesian: λ 1 = 2.4 ,   λ 2 = 2.8 ,   θ = 0.4 .
λ 1 = 2.4 ,   λ 2 = 2.8 ,   θ = 0.4 MLESELFLINEX 1 (c = 0.5)LINEX 2 (c = 2)
npm BiasMSELACICPBiasMSELCCIBiasMSELCCIBiasMSELCCI
500.335 λ 1 −0.03280.08591.142495.00%0.02990.02300.5356−0.02730.01990.49700.09280.03490.5771
λ 2 0.01240.12921.409096.10%0.06530.03930.6593−0.02610.02960.60330.15660.06760.7453
θ 0.08490.46792.661995.40%0.06310.15781.5195−0.07650.17161.52990.07760.15861.4898
45 λ 1 −0.00610.06230.978595.50%0.02800.01930.4932−0.00860.01600.46240.08830.02830.5255
λ 2 −0.00980.10651.279397.70%0.04070.03040.5656−0.02210.02620.53190.10920.04450.6137
θ 0.01970.25381.974395.60%−0.02790.13841.4122−0.07000.15521.42930.04320.13421.3645
0.535 λ 1 −0.00990.09391.201595.20%0.04440.02570.5542−0.01380.02080.52390.10850.03980.5948
λ 2 0.01790.12891.406295.80%0.06890.03980.7032−0.01930.02970.63770.16090.06920.7777
θ 0.06501.23084.343595.10%−0.03660.16601.5987−0.11810.18641.61720.04420.16151.5679
45 λ 1 −0.00810.05670.931995.70%0.03300.01690.4774−0.01290.01420.45060.08270.02480.5081
λ 2 −0.00540.08301.129595.90%0.04430.02590.5930−0.01840.02130.56120.11280.04070.6421
θ 0.00860.26122.004395.20%−0.03480.14411.4432−0.10910.16031.48880.03830.14061.4236
0.835 λ 1 −0.01920.10101.244495.20%0.03960.02240.5497−0.01880.01820.51200.10390.03590.5991
λ 2 −0.00520.15981.567995.20%0.05590.03660.6698−0.02390.02900.61350.14490.06180.7362
θ 0.14752.12143.642396.30%−0.03730.16541.5142−0.12070.18561.53420.04440.16181.4888
45 λ 1 −0.01220.06661.011197.30%0.03530.01880.5002−0.01080.01580.47250.08520.02740.5342
λ 2 −0.00150.09331.198397.70%0.04670.02680.6009−0.01610.02190.56270.11520.04180.6486
θ 0.03170.29822.138299.90%−0.01980.13731.4085−0.09280.15241.45120.04250.13521.3700
1000.375 λ 1 −0.01480.03810.762995.20%0.01840.01060.3849−0.00940.00960.37330.04760.01330.3999
λ 2 −0.00330.05530.921894.60%0.02770.01580.4686−0.01470.01400.45520.06800.02110.4926
θ −0.00510.12401.380995.30%−0.02220.09641.1995−0.07530.10561.21450.03010.09361.1848
90 λ 1 −0.00960.03010.679095.90%0.01690.00810.3380−0.00630.00740.32680.04100.01000.3473
λ 2 −0.00270.04110.794595.60%0.02070.01120.4087−0.01110.01020.39620.05410.01460.4212
θ −0.00450.09991.239196.10%−0.02030.08201.1113−0.07360.09061.14680.02600.07961.0764
0.575 λ 1 −0.01750.03790.760994.80%0.01620.01020.3824−0.01140.00940.37090.04510.01280.3982
λ 2 −0.00500.04970.874294.50%0.02580.01360.4494−0.01210.01210.43130.06590.01850.4711
θ 0.01210.16161.575996.00%−0.01790.10411.2251−0.07320.11371.25360.03710.10161.2036
90 λ 1 0.00170.03310.713195.40%0.01520.00970.3500−0.00100.00880.34470.04170.01200.3572
λ 2 −0.00030.04380.820895.80%0.02500.01270.4004−0.00690.01130.38790.05830.01640.4152
θ 0.00170.10491.270096.60%−0.01530.08261.0660−0.07070.09141.09880.01930.07951.0361
0.875 λ 1 −0.01150.04000.783695.80%0.02110.01100.3793−0.00670.00990.36940.05020.01390.3950
λ 2 −0.00650.04930.870694.60%0.02480.01410.4378−0.01320.01260.41860.06480.01900.4604
θ 0.01680.13161.421394.60%−0.02690.10361.2091−0.08120.11361.21110.02700.10011.2003
90 λ 1 −0.00930.03140.693696.70%0.01740.00880.3454−0.00600.00800.33790.04180.01080.3530
λ 2 −0.00300.04040.788494.90%0.02240.01150.3998−0.00930.01040.38560.05540.01490.4137
θ 0.01640.10661.279194.70%−0.00650.08691.1233−0.05550.09371.14340.01940.08571.1097
Table 2. MLE and Bayesian: λ 1 = 0.6 ,   λ 2 = 2.8 ,   θ = 0.4 .
Table 2. MLE and Bayesian: λ 1 = 0.6 ,   λ 2 = 2.8 ,   θ = 0.4 .
λ 1 = 0.6 ,   λ 2 = 2.8 ,   θ = 0.4 MLESELFLINEX 1 (c = 0.5)LINEX 2 (c = 2)
npm BiasMSELACICPBiasMSELCCIBiasMSELCCIBiasMSELCCI
500.335 λ 1 −0.00470.00540.287095.00%0.01140.00170.14610.00770.00160.14290.01530.00190.1495
λ 2 0.02770.10331.255795.00%0.06240.03660.6941−0.01960.02810.64280.15490.06420.7735
θ 0.02970.33602.270494.20%−0.04150.15941.5305−0.12100.18121.56690.03790.15281.4824
45 λ 1 0.00340.00400.248795.40%0.01040.00130.12580.00740.00120.12420.01340.00140.1275
λ 2 0.02410.09121.180695.20%0.05070.02960.6108−0.01310.02400.57080.12050.04580.6623
θ 0.02740.27962.053195.70%−0.00860.13721.4309−0.07950.14861.48110.02610.13811.4265
0.535 λ 1 0.00100.00470.268995.10%0.01200.00150.13460.00820.00140.13260.01590.00170.1364
λ 2 0.02540.11341.316994.20%0.05630.03730.6965−0.02510.03000.65140.14780.06300.7502
θ 0.07300.40352.474794.40%−0.03340.16101.5234−0.11140.17651.52830.04390.16001.4947
45 λ 1 −0.00070.00400.247295.40%0.00960.00120.12620.00660.00110.12450.01260.00130.1278
λ 2 0.00770.09171.187294.40%0.04490.02700.5973−0.01800.02230.56260.11370.04180.6484
θ 0.06470.28872.091995.10%−0.01590.13901.4569−0.08570.15261.47260.04150.13701.4280
0.835 λ 1 0.00210.00620.308994.70%0.01300.00160.13900.00920.00150.13770.01690.00180.1412
λ 2 0.04400.15581.538394.40%0.06840.04040.7062−0.01380.03060.65670.16120.07020.7839
θ 0.10840.95643.811994.70%−0.00720.15351.5017−0.08510.16801.54140.07030.15431.4590
45 λ 1 −0.00130.00380.243294.80%0.00910.00110.12200.00610.00110.12080.01210.00120.1236
λ 2 0.02160.09001.173394.70%0.05250.02940.6045−0.01140.02350.57670.12270.04620.6564
θ 0.06080.24441.924395.50%−0.00710.13341.3996−0.07710.14581.42960.06070.13361.3723
1000.375 λ 1 −0.00140.00240.193695.10%0.00540.00070.09860.00370.00070.09790.00720.00070.0990
λ 2 0.01890.04940.871294.00%0.02730.01480.4558−0.01110.01310.43880.06790.01990.4729
θ 0.01040.13281.428895.30%−0.02520.10191.2475−0.08050.11201.23960.02920.09871.2122
90 λ 1 −0.00050.00200.176595.70%0.00450.00060.08950.00300.00050.08870.00600.00060.0901
λ 2 0.01370.04620.841395.20%0.02480.01350.4145−0.00720.01220.40390.05820.01720.4269
θ 0.01040.10201.244195.70%−0.01500.08561.1287−0.06420.09351.17310.02830.08331.0873
0.575 λ 1 0.00380.00260.198595.80%0.00570.00070.10060.00390.00070.09920.00750.00080.1014
λ 2 0.01970.05170.888295.30%0.03200.01530.4358−0.00640.01320.42200.07270.02090.4541
θ 0.01480.12161.366594.90%−0.02310.10411.2150−0.07770.11531.25610.03110.09911.1716
90 λ 1 −0.00290.00200.173996.80%0.00390.00050.08630.00240.00050.08560.00540.00060.0867
λ 2 0.01910.04030.783395.50%0.02860.01180.4000−0.00350.01040.39080.06210.01570.4125
θ 0.01100.09851.229895.10%−0.02250.08251.0821−0.07460.09111.10930.02400.07951.0663
0.875 λ 1 0.00190.00230.188094.10%0.00700.00060.09360.00530.00060.09300.00890.00070.0945
λ 2 0.02030.04510.829595.00%0.03290.01460.4401−0.00530.01250.42320.07330.02010.4602
θ 0.03910.12341.368995.10%−0.02120.09421.1688−0.06550.10311.19800.04190.09241.1256
90 λ 1 −0.00120.00190.169096.10%0.00490.00050.08350.00340.00050.08250.00640.00050.0842
λ 2 0.01810.04020.783596.50%0.02380.01170.4021−0.00480.01050.38730.05740.01530.4194
θ 0.02530.10521.268295.70%−0.02010.08411.0994−0.06100.09131.11560.02660.08221.0762
Table 3. MLE and Bayesian: λ 1 = 0.6 ,   λ 2 = 1.4 ,   θ = 0.4 .
Table 3. MLE and Bayesian: λ 1 = 0.6 ,   λ 2 = 1.4 ,   θ = 0.4 .
λ 1 = 0.6 ,   λ 2 = 1.4 ,   θ = 0.4 MLESELFLINEX 1 (c = 0.5)LINEX 2 (c = 2)
npm BiasMSELACICPBiasMSELCCIBiasMSELCCIBiasMSELCCI
500.335 λ 1 −0.01320.00510.279694.80%0.01170.00160.14150.00800.00140.13900.01560.00170.1436
λ 2 0.00730.02780.653395.40%0.02990.00880.33720.00900.00730.32420.05220.01150.3513
θ 0.09470.39872.448494.90%−0.03080.17631.5871−0.10970.19091.57190.05750.17591.5793
45 λ 1 −0.00820.00390.243294.90%0.00900.00120.12900.00600.00110.12763.00000.00130.1312
λ 2 0.00670.02150.571895.70%0.02750.00720.29590.00810.00610.28580.04470.00890.3084
θ 0.05580.23931.906295.30%−0.01660.13511.3935−0.08470.14981.42070.05140.13071.3536
0.535 λ 1 0.00450.00600.304594.50%0.01250.00150.14030.00870.00140.13800.01640.00170.1430
λ 2 0.01120.03080.687294.30%0.03010.00900.34300.00930.00750.32610.05230.01170.3608
θ 0.08471.02153.949994.20%−0.02640.15981.5356−0.10440.17811.61910.05190.15541.4894
45 λ 1 −0.00110.00410.252694.80%0.01050.00130.12560.00750.00120.12380.01360.00140.1276
λ 2 0.00270.02290.593395.20%0.02080.00680.29580.00480.00600.28820.03750.00820.3069
θ 0.05270.23481.889294.80%−0.02600.12971.4052−0.09960.14461.44220.03970.12561.3729
0.835 λ 1 −0.00340.00550.290494.50%0.01180.00170.14600.00800.00160.14390.01570.00190.1487
λ 2 −0.00130.02970.675495.60%0.02760.00950.34730.00700.00800.33020.04950.01210.3645
θ 0.11880.36472.322294.60%−0.00950.16131.5605−0.08930.17331.57900.06950.16471.5288
45 λ 1 −0.00280.00370.237695.20%0.00840.00110.12400.00540.00110.12250.01150.00120.1259
λ 2 0.00100.02190.579095.83%0.02610.00710.30060.00610.00600.28780.04320.00880.3129
θ 0.07300.27032.018995.50%−0.00810.13671.4128−0.08280.14851.44420.05770.13681.4020
1000.375 λ 1 −0.00070.00210.180894.70%0.00550.00060.09530.00370.00060.09430.00730.00070.0963
λ 2 0.01110.01280.441394.60%0.01700.00380.22200.00720.00350.21840.02710.00450.2295
θ 0.03140.12181.362995.00%−0.02090.09901.2019−0.07590.10851.22190.03330.09621.1645
90 λ 1 0.00060.00200.173595.40%0.00540.00060.08690.00300.00050.08620.00710.00060.0874
λ 2 0.00100.01020.396995.30%0.01010.00300.20490.00210.00280.20100.01830.00330.2077
θ 0.03050.10681.266495.20%−0.01780.08851.1292−0.06590.09561.14230.02990.08621.1189
0.575 λ 1 −0.00180.00230.186595.10%0.00550.00060.09340.00370.00060.09290.00730.00070.0939
λ 2 0.00730.01310.448195.20%0.01500.00380.22750.00520.00350.22160.02510.00440.2321
θ 0.03700.12441.375893.70%−0.02690.09681.1932−0.07970.10741.21780.02570.09281.1561
90 λ 1 0.00130.00210.178896.50%0.00510.00060.08940.00340.00060.08880.00710.00060.0903
λ 2 0.00120.01070.402495.80%0.01460.00330.21450.00470.00300.21180.02290.00370.2182
θ 0.02440.10771.283794.50%−0.02310.08571.1155−0.07480.09661.14660.01790.08081.0828
0.875 λ 1 −0.00140.00240.191395.90%0.00580.00070.09470.00400.00060.09370.00760.00070.0958
λ 2 0.00340.01150.421095.70%0.01310.00340.21710.00340.00320.21440.02300.00390.2224
θ 0.00990.12231.371195.40%−0.04000.09761.1740−0.09840.11211.22060.01760.09161.1566
90 λ 1 −0.00120.00200.177596.70%0.00460.00060.08810.00310.00050.08720.00610.00060.0886
λ 2 0.00190.01050.401495.90%0.01140.00300.20680.00300.00280.20290.01970.00330.2108
θ 0.00610.10111.247096.30%−0.03510.08441.1097−0.09700.09591.1398−0.00240.07881.0932
Table 4. MLE and Bayesian: λ 1 = 1.5 ,   λ 2 = 2 ,   θ = 0.8 .
Table 4. MLE and Bayesian: λ 1 = 1.5 ,   λ 2 = 2 ,   θ = 0.8 .
MLESELFLINEX 1 (c = 0.5)LINEX 2 (c = 2)
npm BiasMSELACICPBiasMSELCCIBiasMSELCCIBiasMSELCCI
500.335 λ 1 −0.00700.03430.726094.20%0.03120.00930.34030.00850.00760.32630.05520.01240.3597
λ 2 0.00640.07221.053594.70%0.03910.01760.4679−0.00150.01430.43770.08300.02530.5004
θ 0.24611.18244.154094.10%−0.07420.24691.9682−0.23960.30212.01280.08180.25291.9908
45 λ 1 −0.00590.02760.651495.70%0.02300.00780.31780.00520.00670.31040.04170.00950.3326
λ 2 −0.00100.04720.851795.40%0.03330.01360.40880.00130.01130.39120.06690.01840.4356
θ 0.14900.47812.648195.50%−0.07080.20221.7733−0.23100.25461.78410.05210.19971.7729
0.535 λ 1 −0.00390.03300.712394.40%0.02920.00960.35210.00650.00790.33270.05320.01260.3718
λ 2 0.00920.06611.007694.60%0.04870.01980.48060.00750.01540.44950.09350.02900.5141
θ 0.18210.80073.436194.80%−0.09440.26292.0018−0.26600.32351.99560.06600.27301.9872
45 λ 1 −0.00350.02560.626795.00%0.02470.00770.31670.00610.00660.30630.04350.00960.3271
λ 2 −0.00290.04780.857295.10%0.03380.01390.42870.00220.01160.40700.06740.01880.4513
θ 0.12310.44192.562195.00%−0.08810.21491.7729−0.23580.26461.80310.04860.21461.7874
0.835 λ 1 0.00620.04160.799994.60%0.03170.01040.36270.00890.00860.34610.05580.01350.3863
λ 2 0.01730.07301.057694.50%0.05060.01990.47890.00950.01540.45200.09500.02890.5167
θ 0.20361.83355.250294.80%−0.09300.26931.9472−0.26370.33011.96010.06640.27051.9056
45 λ 1 −0.00510.02580.629194.70%0.02330.00710.29630.00560.00620.28900.04180.00890.3123
λ 2 0.01070.05050.880095.60%0.04080.01460.41830.00880.01180.39690.07490.02010.4439
θ 0.13350.48122.669796.50%−0.08880.21041.7752−0.23190.25251.73740.04360.21331.7664
1000.375 λ 1 −0.00440.01540.486195.50%0.01370.00430.24120.00300.00400.23570.02460.00490.2473
λ 2 −0.00410.02910.668995.70%0.01880.00790.3236−0.00080.00720.31400.03850.00950.3345
θ 0.06900.15871.538895.00%−0.06460.12211.3426−0.16120.15111.38900.02350.11381.3015
90 λ 1 0.00350.01170.424596.20%0.01250.00340.21160.00260.00310.20850.02440.00390.2154
λ 2 −0.00080.02120.570796.10%0.01520.00580.2821−0.00070.00530.27690.03150.00680.2913
θ 0.04210.11031.292296.40%−0.05720.09151.1658−0.15390.11781.20330.00320.08111.1233
0.575 λ 1 −0.00370.01500.481094.80%0.01470.00430.23910.00400.00390.23360.02570.00490.2447
λ 2 0.00340.02620.634595.30%0.01910.00720.31630.00300.00650.30900.03870.00880.3248
θ 0.07010.16011.545095.10%−0.07230.11091.1736−0.16700.14101.24430.01420.10131.1556
90 λ 1 −0.00250.01190.428095.40%0.01230.00330.21120.00330.00300.20770.02140.00370.2152
λ 2 0.00310.02100.567795.70%0.01870.00590.28500.00280.00530.27910.03510.00700.2924
θ 0.04900.11921.340595.60%−0.06970.09591.1814−0.15060.12131.21890.00490.08611.1466
0.875 λ 1 0.00180.01420.467194.70%0.01660.00430.25010.00570.00390.24370.02790.00500.2547
λ 2 −0.00790.04910.868395.10%0.03140.16650.3152−0.00530.00670.30570.07860.02610.3280
θ 0.03350.14021.462895.00%−0.09680.12171.3807−0.19220.15571.4195−0.00850.10871.3358
90 λ 1 −0.00070.01120.414795.20%0.01220.00320.20460.00320.00290.20110.02140.00360.2089
λ 2 0.00610.02360.601895.70%0.02170.00680.30320.00510.00610.29590.03820.00810.3118
θ 0.03040.10951.287496.00%−0.08550.09241.1046−0.16850.12131.1597−0.00820.07981.0759

7. Applications

In this section, we discussed two bivariate real data sets. The FGM-BLBE distribution has been used based on progressive Type-II censoring with random removal.

7.1. Iron Material Jobs Data

The Iron material jobs data set has been obtained from Dasgupta [37]. In the context of tasks involving iron sheets, we employ a perforation procedure that involves drilling a total of four holes, two on each arm, in an L-shaped rectangular sheet measuring 100 mm by 150 mm. This efficient process is conducted on a 100-ton press operating at a speed of 250 strokes per hour, known as piercing. Simultaneously creating two holes at a time, the perforated L-shaped iron sheet is intended for use in mini- or light-duty truck chassis.
Following the piercing, a burr is formed in a circular shape around the hole on the opposite side of the metal sheet. The intense pressure applied during the piercing causes an irregular shift in the contact surface, leading to uneven and aggressive elevation of metal granules along the rim of the circular hole. The size of the burr is influenced by the characteristics of the metal grains and the piercing load. Skrotzki et al. [38] noted that factors such as composition, melting point, cooling pace, thermal and constitutional under-cooling, and convection all impact the grain structure and texture of metals.
Subsequently, the burr is removed through chamfering using a drill. The dial gauge, employed to measure burr in the provided data sets, has a minimum count of 20 microns (m), equivalent to 0.02 mm. For the first set of 50 observations on burr, the hole diameter and sheet thickness are 12 mm and 3.15 mm, respectively. In the second data set of 50 observations, the hole diameter and sheet thickness are 9 mm and 2 mm, respectively.
During the actual works, one hole is selected and positioned based on a predetermined orientation, and readings of the hole diameter are taken in relation to that specific hole. It is important to note that the two data sets are associated with the comparison of two computers. For additional details regarding data on iron material jobs, refer to Dasgupta [37]. Now, we used this data.
The MLE estimates and different measures for each bivariate models for iron material jobs data and MLE and Bayesian estimation under progressive type II with random removal for Iron material jobs data are presented in Table 5 and Table 6, respectively, depicting various instances of the progressive type II scheme. Specifically, when m = 35 and p = 0.3, and 0.8, the corresponding value of R is shown as follows:
  • Data 1: m = 35, p = 0.3:
  • x: 0.02, 0.02, 0.04, 0.04, 0.06, 0.08, 0.08, 0.08, 0.12, 0.12, 0.12, 0.14, 0.14, 0.14, 0.14, 0.14, 0.14, 0.16, 0.16, 0.16, 0.16, 0.16, 0.18, 0.18, 0.18, 0.22, 0.22, 0.24, 0.24, 0.24, 0.24, 0.26, 0.26, 0.26, 0.32
  • y: 0.12, 0.16, 0.06, 0.04, 0.24, 0.16, 0.32, 0.18, 0.26, 0.06, 0.04, 0.14, 0.24, 0.22, 0.02, 0.24, 0.14, 0.06, 0.16, 0.22, 0.06, 0.16, 0.02, 0.14, 0.32, 0.24, 0.18, 0.16, 0.18, 0.14, 0.16, 0.22, 0.18, 0.26, 0.04
  • R: 4 3 2 1 1 1 0 3 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
  • m = 35, p = 0.8:
  • x: 0.02, 0.02, 0.04, 0.06, 0.08, 0.08, 0.08, 0.08, 0.12, 0.12, 0.12, 0.14, 0.14, 0.14, 0.14, 0.14, 0.16, 0.16, 0.16, 0.16, 0.16, 0.16, 0.18, 0.22, 0.22, 0.22, 0.24, 0.24, 0.24, 0.24, 0.24, 0.26, 0.26, 0.32, 0.32
  • y: 0.12, 0.16, 0.06, 0.24, 0.16, 0.32, 0.02, 0.18, 0.26, 0.06, 0.04, 0.14, 0.24, 0.02, 0.24, 0.14, 0.06, 0.16, 0.14, 0.06, 0.06, 0.16, 0.32, 0.08, 0.24, 0.14, 0.16, 0.12, 0.18, 0.14, 0.16, 0.22, 0.18, 0.22, 0.04
  • R: 3 10 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0
Given the absence of information on the unknown parameters, we employ non-informative priors with a1 = b1 = a2 = b2 = 0 in this illustrative example. Additionally, we presume informative gamma priors for the parameters of FGM-BLBE. MLEs and BEs, along with their standard errors, are computed for both case 1 and case 2, and the results are presented in Table 6. This table reveals the point estimates for the unknown parameters λ 1 , λ 2 , and θ obtained through maximum likelihood and Bayesian methods. The findings in Table 6 highlight that the HPD intervals are marginally shorter than the other confidence intervals.
The FGM-BLBE distribution is compared to the following bivariate models in Table 5 as FGM Kavya-Manoharan Pareto (FGMKMP) (Fayomi et al. [16]), Frank bivariate X Lindely (FBXL), Gumbel bivariate X Lindely (GBXL), Clayton bivariate X Lindely (CBXL), Clayton Kavya-Manoharan exponential (CKME) (Fayomi et al. [16]), FGM Fréchet (FGMF) (Almetwally and Muhammed [39]), AMH Fréchet (AMHF) (Almetwally and Muhammed [39]), and Bivariate FGM Generalized Exponential (FGMGE) (Abd Elaal and Jarwan [40]). Table 5 describes the MLEs of the distribution’s parameters and presents the goodness-of-fit metrics: Akaike Information Criterion (AIC), Corrected Akaike Information Criterion (CAIC), Bayesian Information Criterion (BIC), Hannan-Quinn Information Criterion (HQIC), Cramér-von Mises (CVM), and Anderson-Darling (AD). According to Table 5 findings, the FGM-BLBE distribution performs better than the FGMKMP, FBXL, CBXL, GBXL, FGMF, AMHF, and UEPFGMGE distributions for the iron material jobs data set.
Figure 3 presents histograms for each variable of the iron material jobs data, which are displayed in the diagonal elements, where the X-set, ranging from 0.02 to 0.32, exhibits variability with mean around 0.16. Distribution skewness is apparent, showcasing a diverse spread of values with some concentration around the median. While for Y-set, spanning values from 0.02 to 0.32, displays variability with a mean around 0.15. It exhibits some concentration around the median, indicating a skewed distribution with instances of recurring values, notably around 0.14 and 0.22. In the upper triangle of the matrix, numerical values represent the correlation coefficients between pairs of variables of iron material jobs data as 0.16. While in the lower triangle of the matrix, pairs panels function is particularly useful for exploring relationships between multiple variables simultaneously, providing a visual summary of the iron material jobs data to identify patterns, trends, and potential correlations.
Figure 4 presents scatter, contour, and density plots for the FGM (Farlie-Gumbel-Morgenstern) distribution applied to iron material jobs data under different censoring schemes. The scatter plots provide a visual representation of the relationship between variables under different censoring scenarios, showing how the data points distribute across the given parameters. The density plots offer insights into the probability density function of the data, highlighting areas with higher concentrations of probability mass, while the contour plots delineate the levels of the density function, helping to visualize the gradient and spread of the data across the parameter space.
The scatter plot in R visually depicts the relationship between variables X and Y of iron material jobs data in a three-dimensional space with different binomial removal as (p = 0, p = 0.3, and p = 0.8). The red dots represent data points, revealing patterns and trends in the distribution. This 3D scatter plot offers an insightful view of the dataset’s structure and inter-variable interactions of iron material jobs data. The density plot of FGM copula, showcasing the FGM copula for variables X and Y, offers a captivating three-dimensional representation. This visualization vividly communicates the dependence structure, allowing for a nuanced understanding of the joint behavior of X and Y within the FGM copula framework with different copula parameter under different binomial removal as (p = 0, p = 0.3, and p = 0.8). The contour plot effectively visualizes the density of the FGM copula for variables X and Y. Contour lines delineate regions of similar density, revealing intricate patterns in the joint distribution. This plot provides valuable insights into the dependence structure between X and Y within the FGM copula framework for iron material jobs data.

7.2. Computers Times

The data was obtained from Oliveira et al. [41]. The dataset comprises 50 simulated primitive computer series systems, each equipped with a processor and memory. The operational state of the computer relies on the proper functioning of both components. Assuming that the system undergoes a latent deterioration process, the deterioration advances rapidly over a short time frame, rendering the system more vulnerable to shocks. This heightened susceptibility allows a lethal shock to randomly damage either the first, second, or both components. Due to the possibility of a fatal shock affecting both components simultaneously, the reliability of the independence premise is questionable. Consequently, we employed the FGM copula to investigate this matter.
The FGM-BLBE distribution is compared to the following bivariate models in Table 7 as FGM Kavya-Manoharan Pareto (FGMKMP) (Fayomi et al. [16]), FGM Fréchet (FGMF) (Almetwally and Muhammed [39]), and AMH Fréchet (AMHF) (Almetwally and Muhammed [39]). Table 7 described the MLEs of the distribution’s parameters and showed the goodness of fit metrics, AIC, CAIC, BIC, HQIC, CVM, and AD. According to Table 7 findings, the FGM-BLBE distribution performs better than the FGMKMP, FGMP, FGMF, and AMHF distributions for the computer times data set.
To identify outliers in the processor and memory data, we generated Figure 5 and employed matrix histogram, scatter plot with a boxplot to distinguish between different groups. Conversely, Figure 5 delved into the dependent plot of processor and memory data, highlighting dependencies and the density distribution of numerical data. The examination of Figure 5 reveals that the data exhibits right-skewed shapes, asymmetry, and a weak negative correlation.
The MLE estimates and various metrics for individual bivariate models in computer times data, as well as the MLE and Bayesian estimates under the progressive Type II with random removal, are displayed in Table 7 and Table 8. These tables represent different scenarios of the progressive Type II scheme, particularly when m = 40 and p = 0.3, and 0.8. The associated value of R is provided for each case.
  • Data set 2: Case 1: m = 40 p = 0.3:
  • x: 0.1058, 0.1058, 0.1115, 0.1955, 0.3309, 0.4614, 0.5520, 0.6270, 0.7947, 0.8584, 0.9096, 0.9867, 0.9867, 0.9938, 1.0051, 1.0833, 1.0833, 1.0986, 1.1739, 1.1917, 1.3288, 1.3482, 1.5254, 1.6465, 1.7494, 1.9292, 2.1000, 2.1396, 2.3364, 2.3757, 2.5913, 3.3887, 3.5202, 3.6621, 3.6621, 3.6621, 4.3435, 5.0533, 5.7561, 5.7561.
  • y: 0.1058, 0.1058, 0.1115, 0.1955, 0.3309, 0.8584, 0.5520, 1.7289, 0.7947, 1.9556, 0.6214, 0.9867, 0.9867, 1.7689, 1.0051, 3.3059, 3.3059, 1.0986, 3.3857, 0.0801, 0.9689, 1.9705, 4.4088, 2.0197, 2.3643, 3.9291, 2.0513, 2.1548, 0.1624, 2.7953, 2.5913, 1.9796, 1.4095, 0.0026, 0.0026, 0.0026, 1.0001, 2.3238, 0.3212, 0.3212.
  • R: 3 2 1 1 0 1 0 2 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.
  • Case 2: m = 40 p = 0.8:
  • x: 0.1058, 0.1058, 0.1115, 0.1181, 0.3309, 0.5079, 0.5520, 0.5784, 0.6270, 0.8503, 0.8584, 0.9096, 0.9867, 0.9867, 0.9938, 1.0051, 1.0833, 1.0833, 1.0986, 1.1739, 1.1917, 1.3482, 1.3640, 1.3989, 1.5254, 1.6465, 1.7494, 1.9386, 2.1000, 2.1396, 2.3364, 2.3757, 2.5913, 3.3887, 3.5202, 3.6621, 3.6621, 5.0533, 5.7561, 5.7561.
  • y: 0.1058, 0.1058, 0.1115, 0.0884, 0.3309, 5.3535, 0.5520, 1.8795, 1.7289, 2.8578, 1.9556, 0.6214, 0.9867, 0.9867, 1.7689 1.0051, 3.3059, 3.3059, 1.0986, 3.3857, 0.0801, 1.9705, 1.3640, 4.1268, 4.4088, 2.0197, 2.3643, 4.0043, 2.0513, 2.1548, 0.1624, 2.7953, 2.5913, 1.9796, 1.4095, 0.0026, 0.0026, 2.3238, 0.3212, 0.3212.
  • R: 2 7 1 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0 0.
Figure 5 discussed histograms for each variable of computers times data are displayed in diagonal elements, where the X-set, ranging from 0.1058 to 5.7561. Mean-centered around 1.5662, the distribution displays variability with instances of higher values, notably around 3.6621. While for Y-set, spanning values from 0.0026 to 5.3535, displaying variability with peaks around 3.3059 and 4.0043. Mean-centered at 1.3983, the distribution suggests diverse patterns. In the upper triangle of the matrix, numerical values represent the correlation coefficients between pairs of variables of computers times data as −0.03 (negative week correlation). While in the lower triangle of the matrix, pairs panels function is particularly useful for exploring relationships between multiple variables simultaneously, providing a visual summary of the computers times data to identify patterns, trends, and potential correlations. Figure 6 discussed Scatter, contour, and density plots For FGM copula for computers times data under different censoring schemes.

8. Conclusions

In this paper, we explored both Bayesian and non-Bayesian methods for estimating the parameters of the FGM-BLBE distribution when data is subject to progressive Type-II censoring with random removal. The binomial removal method was employed to eliminate failure times and generate censored samples. Our findings are relevant to various real-world applications. Beyond traditional confidence intervals, we examined bootstrap confidence intervals using both the percentile and bootstrap-t methods. Our simulations revealed that Bayesian estimation using MCMC—particularly under the LINEX loss function—yielded accurate parameter estimates for the FGM-BLBE distribution. This conclusion is supported by low bias, MSE, and narrow confidence intervals. Importantly, increasing the number of censored samples (m) further improved the accuracy of all estimators. Notably, Bayesian estimation outperformed non-Bayesian methods, and Bayes credible intervals proved more reliable than asymptotic confidence intervals. To illustrate practical applications, we presented two real-world examples from industry and computer science, demonstrating successful parameter estimation within the FGM-BLBE distribution.

Author Contributions

Conceptualization, M.E.Q. and E.M.A.; methodology, A.F., E.M.A. and M.E.Q.; software, E.M.A.; formal analysis, M.E.Q.; writing—original draft, M.E.Q.; writing—review and editing, A.F.; visualization, E.M.A.; project administration, A.F. and M.E.Q.; funding acquisition, A.F. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

All data are available in the manuscript.

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Balakrishnan, N.; Aggarwala, R. Progressive Censoring: Theory, Methods, and Applications; Springer Science & Business Media: Berlin, Germany, 2000. [Google Scholar]
  2. Balakrishnan, N. Progressive censoring methodology: An appraisal. Test 2007, 16, 211–259. [Google Scholar] [CrossRef]
  3. Balakrishnan, N.; Cramer, E. The Art of Progressive Censoring; Statistics for Industry and Technology; Springer: New York, NY, USA, 2014. [Google Scholar]
  4. El-Sherpieny, E.S.A.; Muhammed, H.Z.; Almetwally, E.M. Progressive Type-II censored samples for bivariate Weibull distribution with economic and medical applications. Ann. Data Sci. 2022, 11, 51–85. [Google Scholar] [CrossRef]
  5. Nelsen, R.B. An Introduction to Copulas; Springer Science & Business Media: Berlin/Heidelberg, Germany, 2007; ISBN 978-0387286594. [Google Scholar]
  6. Flores, A.Q. Testing Copula Functions as a Method to Derive Bivariate Weibull Distributions. In Proceedings of the American Political Science Association (APSA), Annual Meeting 2009, Toronto, ON, Canada, 3–6 September 2009; Volume 4. [Google Scholar]
  7. Suzuki, A.K.; Louzada-Neto, F.; Cancho, V.G.; Barriga, G.D. The FGM bivariate lifetime copula model: A bayesian approach. Adv. Appl. Stat. 2011, 21, 55–76. [Google Scholar]
  8. Verrill, S.P.; Evans, J.W.; Kretschmann, D.E.; Hatfield, C.A. Asymptotically efficient estimation of a bivariate Gaussian–Weibull distribution and an introduction to the associated pseudo-truncated Weibull. Commun. Stat.-Theory Methods 2015, 44, 2957–2975. [Google Scholar] [CrossRef]
  9. Peres, M.V.D.O.; Achcar, J.A.; Martinez, E.Z. Bivariate modified Weibull distribution derived from Farlie-Gumbel-Morgenstern copula: A simulation study. Electron. J. Appl. Stat. Anal. 2018, 11, 463–488. [Google Scholar]
  10. Peres, M.V.D.O.; de Oliveira, R.P.; Achcar, J.A.; Martinez, E.Z. The Bivariate Defective Gompertz Distribution Based on Clayton Copula with Applications to Medical Data. Austrian J. Stat. 2020, 51, 144–168. [Google Scholar] [CrossRef]
  11. Samanthi, R.G.M.; Sepanski, J. On bivariate Kumaraswamy-distorted copulas. Commun. Stat.-Theory Methods 2022, 51, 2477–2495. [Google Scholar] [CrossRef]
  12. Elgawad, M.A.A.; Barakat, H.M.; Alawady, M.A.; El-Rahman, D.A.A.; Husseiny, I.A.; Hashem, A.F.; Alotaibi, N. Extropy and Some of Its More Recent Related Measures for Concomitants of K-Record Values in an Extended FGM Family. Mathematics 2023, 11, 4934. [Google Scholar] [CrossRef]
  13. Barakat, H.M.; Nigm, E.M.; Syam, A.H. Concomitants of order statistics and record values from Bairamov-Kotz-Becki-FGM bivariate-generalized exponential distribution. Filomat 2018, 32, 3313–3324. [Google Scholar] [CrossRef]
  14. Qura, M.E.; Fayomi, A.; Kilai, M.; Almetwally, E.M. Bivariate power Lomax distribution with medical applications. PloS ONE 2023, 18, e0282581. [Google Scholar] [CrossRef]
  15. Fayomi, A.; Almetwally, E.M.; Qura, M.E. A novel bivariate Lomax-G family of distributions: Properties, inference, and applications to environmental, medical, and computer science data. AIMS Math. 2023, 8, 17539–17584. [Google Scholar] [CrossRef]
  16. Fayomi, A.; Almetwally, E.M.; Qura, M.E. Exploring New Horizons: Advancing Data Analysis in Kidney Patient Infection Rates and UEFA Champions League Scores Using Bivariate Kavya–Manoharan Transformation Family of Distributions. Mathematics 2023, 11, 2986. [Google Scholar] [CrossRef]
  17. Chesneau, C. Three-Dimensional Copulas: A Generalized Convex Mixture Copulas Strategy. Comput. J. Math. Stat. Sci. 2024, 3, 240–257. [Google Scholar] [CrossRef]
  18. Almetwally, E.M.; Fayomi, A.; Qura, M.E. Advanced Copula-Based Models for Type II Censored Data: Applications in Industrial and Medical Settings. Mathematics 2024, 12, 1774. [Google Scholar] [CrossRef]
  19. Balakrishnan, N.; Kim, J.A. EM algorithm for Type-II right censored bivariate normal data. In Parametric and Semiparametric Models with Applications to Reliability, Survival Analysis, and Quality of Life; Birkhäuser: Boston, MA, USA, 2004; pp. 177–210. [Google Scholar]
  20. Balakrishnan, N.; Kim, J.A. EM algorithm and optimal censoring schemes for progressively Type-II censored bivariate normal data. In Advances in Ranking and Selection, Multiple Comparisons, and Reliability: Methodology and Applications; Birkhäuser: Boston, MA, USA, 2005; pp. 21–45. [Google Scholar]
  21. Balakrishnan, N.; Kim, J.A. Point and interval estimation for bivariate normal distribution based on progressively Type-II censored data. Commun. Stat. Methods 2005, 34, 1297–1347. [Google Scholar] [CrossRef]
  22. Kim, S.W.; Ng, H.K.T.; Jang, H. Estimation of parameters in a bivariate generalized exponential distribution based on Type-II censored samples. Commun. Stat.-Simul. Comput. 2016, 45, 3776–3797. [Google Scholar] [CrossRef]
  23. Bai, X.; Shi, Y.; Liu, B.; Fu, Q. Statistical inference of Marshall-Olkin bivariate Weibull distribution with three shocks based on progressive interval censored data. Commun. Stat.-Simul. Comput. 2019, 48, 637–654. [Google Scholar] [CrossRef]
  24. Aly, H.M.; Muhammed, H.Z.; Abuelamayem, O.A. Estimation of the bivariate Kumaraswamy lifetime distribution under progressive type-I censoring. J. Data Sci. 2020, 18, 739–749. [Google Scholar] [CrossRef]
  25. Muhammed, H.Z. Estimating the shape parameters of a new version of bivariate weibull distribution based on complete, Type I and Type II censored samples. Int. J. Reliab. Appl. 2020, 21, 53–78. [Google Scholar]
  26. EL-Morshedy, M.; Alhussain, Z.A.; ATTA, D.; Almetwlly, E.M.; Eliwa, M.S. Bivariate Burr X generator of distributions: Properties and estimation methods with applications to complete and type-II censored samples. Mathematics 2020, 8, 264. [Google Scholar] [CrossRef]
  27. El-Sherpieny, E.S.A.; Almetwally, E.M.; Muhammed, H.Z. Bayesian and non-bayesian estimation for the parameter of bivariate generalized Rayleigh distribution based on clayton copula under progressive type-II censoring with random removal. Sankhya A 2021, 85, 1205–1242. [Google Scholar] [CrossRef]
  28. Muhammed, H. Analysis of Dependent Variables Following Marshal-Olkin Bivariate Distributions in the Presence of Progressive Type II Censoring. Stat. Optim. Inf. Comput. 2023, 11, 694–708. [Google Scholar] [CrossRef]
  29. Arun, S.P.; Chesneau, C.; Maya, R.; Irshad, M.R. Farlie–Gumbel–Morgenstern Bivariate Moment Exponential Distribution and Its Inferences Based on Concomitants of Order Statistics. Stats 2023, 6, 253–267. [Google Scholar] [CrossRef]
  30. Sklar, M. Fonctions de répartition à n dimensions et leurs marges. Ann. L’Isup 1959, 8, 229–231. [Google Scholar]
  31. Gumbel, E.J. Bivariate exponential distributions. J. Am. Stat. Assoc. 1960, 55, 698–707. [Google Scholar] [CrossRef]
  32. Dara, T.; Ahmad, M. Recent Advances in Moment Distributions and Their Hazard Rate. Ph.D. Thesis, National College of Business Administration and Economics, Lahore, Pakistan, 2012. [Google Scholar]
  33. David, H.A.; Nagaraja, H.N. Concomitants of order statistics. Handb. Stat. 1998, 16, 487–513. [Google Scholar]
  34. Wang, K. On Concomitants of Order Statistics. Ph.D. Thesis, The Ohio State University, Columbus, OH, USA, 2008. [Google Scholar]
  35. Metropolis, N.; Rosenbluth, A.W.; Rosenbluth, M.N.; Teller, A.H.; Teller, E. Equations of State Calculations by Fast Computing Machines. J. Chem. Phys. 1953, 21, 1087–1091. [Google Scholar] [CrossRef]
  36. Dey, S.; Dey, T.; Luckett, D.J. Statistical Inference for The Generalized Inverted Exponential Distribution based on Upper Record Values. Math. Comput. Simul. 2016, 120, 64–78. [Google Scholar] [CrossRef]
  37. Dasgupta, R. On the distribution of burr with applications. Sankhya B 2011, 73, 1–19. [Google Scholar] [CrossRef]
  38. Skrotzki, W.; Kegler, K.; Tamm, R.; Oertel, C.G. Grain structure and texture of cast iron aluminides. Cryst. Res. Technol. J. Exp. Ind. Crystallogr. 2005, 40, 90–94. [Google Scholar] [CrossRef]
  39. Almetwally, E.M.; Muhammed, H.Z. On a bivariate Fréchet distribution. J. Stat. Appl. Probab. 2020, 9, 1–21. [Google Scholar]
  40. Abd Elaal, M.K.; Jarwan, R.S. Inference of bivariate generalized exponential distribution based on copula functions. Appl. Math. Sci. 2017, 11, 1155–1186. [Google Scholar]
  41. Oliveira, R.P.; Achcar, J.A.; Mazucheli, J.; Bertoli, W. A new class of bivariate Lindley distributions based on stress and shock models and some of their reliability properties. Reliab. Eng. Syst. Saf. 2021, 211, 107528. [Google Scholar] [CrossRef]
Figure 1. Plots of joint density of FGM-BLBE distribution for some parameter values.
Figure 1. Plots of joint density of FGM-BLBE distribution for some parameter values.
Axioms 13 00664 g001aAxioms 13 00664 g001b
Figure 2. Plots of joint hazard rate of FGM-BLBE distribution for some parameter values.
Figure 2. Plots of joint hazard rate of FGM-BLBE distribution for some parameter values.
Axioms 13 00664 g002
Figure 3. Matrix for iron material jobs data: Histograms, Correlation, and pairs panels.
Figure 3. Matrix for iron material jobs data: Histograms, Correlation, and pairs panels.
Axioms 13 00664 g003
Figure 4. Scatter, contour, and density plots For FGM with iron material jobs data under different censoring schemes.
Figure 4. Scatter, contour, and density plots For FGM with iron material jobs data under different censoring schemes.
Axioms 13 00664 g004
Figure 5. Correlation matrix for computers times data II.
Figure 5. Correlation matrix for computers times data II.
Axioms 13 00664 g005
Figure 6. Scatter, countour, and density plots For FGM with computer times data II under different censoring schemes.
Figure 6. Scatter, countour, and density plots For FGM with computer times data II under different censoring schemes.
Axioms 13 00664 g006
Table 5. MLE estimates and different measures for each bivariate models: iron material jobs data.
Table 5. MLE estimates and different measures for each bivariate models: iron material jobs data.
λ 1 β 1 λ 2 β 2 θ AICCAICBICHQICCVMAD
FGM-BLBEEstimtes0.07959-0.07398-0.79705−205.98594 −205.46420−200.24987−203.801617.3000142.32165
St.Er0.00780 0.00728 0.46327
FGMKMPEstimtes6.01178-6.59510-2.34381−147.11384−146.59210−141.37777−144.929517.3093442.69029
St.Er0.84166 0.90998 0.48686
FBXLEstimtes6.66253-7.22629-3.52154−170.55050−170.02876−164.81443−168.366177.3448742.41380
St.Er0.91581 0.99686 1.16191
CBXLEstimtes6.82193-7.43917-1.00952−171.76838−171.24664−166.03231−169.584057.3410442.45320
St.Er0.88092 0.96401 0.37817
GBXLEstimtes9.00178-11.74399-1.55245−84.14827−83.62654−78.41221−81.963957.3161842.37532
St.Er0.97309 1.01774 0.16223
CKMEEstimtes4.79516-5.28292-1.17668−158.81737−158.29563−153.08130−156.633047.3580242.49229
St.Er0.73600 0.81263 0.41577
FGMFEstimtes0.092751.248450.082171.197090.89482−144.68363−143.31999−135.12351−141.043087.3256242.49152
St.Er0.011310.118680.010490.116980.62744
AMHFEstimtes0.093911.242020.082911.190880.44762−144.12488−142.76125−134.56477−140.484347.3151742.40317
St.Er0.011390.118370.010500.116560.28883
FGMGEEstimtes3.1475011.462582.7212911.505830.64172−203.87218−202.50854−194.31207−200.231637.3288942.31873
St.Er0.707561.563050.591011.614050.42020
Table 6. MLE and Bayesian estimation under progressive type II with random removal: iron material jobs data.
Table 6. MLE and Bayesian estimation under progressive type II with random removal: iron material jobs data.
MLEBayesian Gamma PriorBayesian Non-Informative Prior
p EstimtesSt.ErLACIEstimtesSt.ErLCCICPEstimtesSt.ErLCCICP
0.3 λ 1 0.07980.00940.03670.08170.00930.03570.94970.08530.01460.05590.9500
λ 2 0.07790.00930.03630.07950.00910.03570.94950.08450.01550.05900.9498
θ 0.58980.59322.32550.18530.47791.65330.94960.38570.62062.31370.9497
0.8 λ 1 0.08120.00970.03800.08270.00930.03680.94980.08630.01530.05770.9499
λ 2 0.07360.00890.03480.07480.00870.03140.94980.08060.01550.05880.9498
θ 0.44290.64542.53010.23170.48341.66120.94930.19510.72452.63550.9499
Table 7. MLE estimates and different measures for each bivariate models: Data II.
Table 7. MLE estimates and different measures for each bivariate models: Data II.
FGM-BLBEFGMKMPFGMPFGMFAMHF
EstimtesSt.ErEstimtesSt.ErEstimtesSt.ErEstimtesSt.ErEstimtesSt.Er
λ 1 0.856140.085620.845940.129991.190890.161040.635280.110990.623080.10579
β 1 ---0.888580.088870.887870.08895
λ 2 0.815780.081000.973250.154441.314750.179780.245550.089250.280820.10113
β 2 ---0.420120.038890.416040.03737
θ 0.355290.265280.942480.342170.855560.324450.727270.458710.677550.18285
AIC350.69989354.05397352.84962386.88554385.38092
CAIC351.22163354.57571363.37136388.24918386.74456
BIC356.43596359.79004358.58569396.44566394.94104
HQIC352.88422356.23830355.03395390.52609389.02147
CVM7.676247.714197.726387.915177.81652
AD42.4936242.7382442.5670843.0515742.91217
Table 8. MLE and Bayesian estimation under progressive type II with random removal: Data 2.
Table 8. MLE and Bayesian estimation under progressive type II with random removal: Data 2.
MLEBayesian Gamma PriorBayesian Non-informative Prior
p EstimtesSt.ErLACIEstimtesSt.ErLCCICPEstimtesSt.ErLCCICP
0.3 λ 1 0.92800.10390.40740.94130.10530.40920.94980.98990.16390.62950.9500
λ 2 0.70370.07840.30750.71220.07720.30040.94990.74030.12240.47130.9499
θ 0.27550.31581.23780.33490.33821.27850.94960.11890.33141.25190.9499
0.8 λ 1 0.84690.09550.37430.86180.09550.37540.94990.90420.15290.57270.9499
λ 2 0.84990.09480.37170.86160.09530.37670.94990.88860.14290.54050.9500
θ 0.23160.31801.24640.23180.02011.40300.94990.23160.02031.23660.9499
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

Fayomi, A.; Almetwally, E.M.; Qura, M.E. Bivariate Length-Biased Exponential Distribution under Progressive Type-II Censoring: Incorporating Random Removal and Applications to Industrial and Computer Science Data. Axioms 2024, 13, 664. https://doi.org/10.3390/axioms13100664

AMA Style

Fayomi A, Almetwally EM, Qura ME. Bivariate Length-Biased Exponential Distribution under Progressive Type-II Censoring: Incorporating Random Removal and Applications to Industrial and Computer Science Data. Axioms. 2024; 13(10):664. https://doi.org/10.3390/axioms13100664

Chicago/Turabian Style

Fayomi, Aisha, Ehab M. Almetwally, and Maha E. Qura. 2024. "Bivariate Length-Biased Exponential Distribution under Progressive Type-II Censoring: Incorporating Random Removal and Applications to Industrial and Computer Science Data" Axioms 13, no. 10: 664. https://doi.org/10.3390/axioms13100664

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