Next Article in Journal
Torricelli’s Law in Fractal Space–Time Continuum
Previous Article in Journal
A Matrix Approach to Vertex-Degree-Based Topological Indices
Previous Article in Special Issue
Non-Parametric Estimation of the Renewal Function for Multidimensional Random Fields
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Applying the Laplace Transform Procedure, Testing Exponentiality against the NBRUmgf Class

1
Department of Mathematics, Faculty of Science, Zagazig University, Zagazig 44519, Egypt
2
Department of Management Information Systems, College of Business and Economics, Qassim University, Buraydah 51452, Saudi Arabia
3
Department of Mathematics and Statistics, Faculty of Management Technology and Information Systems, Port Said University, Port Said 42526, Egypt
*
Author to whom correspondence should be addressed.
Mathematics 2024, 12(13), 2045; https://doi.org/10.3390/math12132045
Submission received: 9 May 2024 / Revised: 22 June 2024 / Accepted: 27 June 2024 / Published: 30 June 2024
(This article belongs to the Special Issue Parametric and Nonparametric Statistics: From Theory to Applications)

Abstract

:
This paper addresses a hypothesis testing problem for comparing exponentially distributed data against a new class termed “New Better than Renewal Used in Moment Generating Function” ( N B R U m g f ). A measure of departure from exponentiality is constructed using the Laplace transform, followed by the development of a U-statistic-based test for the hypothesis. Additionally, a test based on the goodness of fit approach is examined as a special case. The asymptotic normality of the proposed statistic is introduced, and Pitman’s asymptotic efficiency of the two tests is computed and compared with other tests. Percentiles of the test statistics are computed for certain sample sizes in the case of complete data, and the powers of the tests are computed for popular reliability distributions. Finally, practical applications of the proposed tests are demonstrated in multiple cases.
MSC:
62Gxx; 62N01; 62N02; 62N03; 62N86

1. Introduction

The investigation of several life distribution classes was introduced in reliability by Barlow and Proschan (1963) [1], Siddiqui and Bryson (1969) [2], Barlow and Proschan (1981) [3], Lee (1990) [4], and Serfling (1980) [5]. These kinds of life distributions are currently used in biometrics, engineering, maintenance, and the social and biological sciences. Consequently, there has been an increasing interest from statisticians and reliability analysts in modeling survival data with life distributions classified according to certain features of aging.
The NBRUmgf class is particularly significant in this context because it generalizes and encompasses a variety of other well-known classes of life distributions, such as new better than used (NBU), New Better than Used in Expectation (NBUE), and Increasing Failure Rate (IFR). This makes the NBRUmgf class a versatile and powerful tool for reliability modeling and life data analysis.
Everestus et al. (2022) [6] examined 91 distinct exponentiality tests. It was claimed that, while some tests were against specific kinds of alternative distributions, others were universally consistent. Their simulation study’s outcomes revealed that some tests that performed well under one set of shape alternatives would not do as well under another set of shape alternatives. Conversely, some tests that seemed less potent at small sample numbers also seemed more potent at large sample sizes. Ultimately, there was total concurrence between the outcomes of the simulated studies and the real-life implementations.
In more recent times, several classes of aging distributions of life have been further established within documented texts. For instance, using the goodness of fit and U-statistic approaches, respectively, Abu-Yousef et al. (2020, 2022) [7,8] introduced two novel test statistics for comparing the new is preferred over the used convex ordering function which generates moments (NBUCmgf) life class distribution to the exponential distribution. In regard to a goodness of fit methodology, Abu Yousef and El-toony (2022) [9] suggested a new test statistic for comparing the better average increasing concave in Laplace transform order (UBAC (2)L) class of life distribution to exponentiality.
A new nonparametric class of life distribution known as new better than renewal used in moment-generating function (NBRUmgf) was proposed by Hassan and Said (2021) [10,11]. The NBRUmgf class extends the concepts of NBU and NBUE by incorporating moment-generating functions, thus providing a more comprehensive framework for analyzing life distributions with various aging properties. Subsequently, two new test statistics were developed based on moment inequalities and the U-statistic, respectively, for testing exponentiality versus (NBRUmgf). El-Arishy et al. (2019) have identified several classes, among which are renewal new better than used in moment-generating function (RNBUmgf) [12]; renewal new superior to that utilized in the Laplace transform order (RNBUL), as noted by Mahmoud et al. (2019) [13]; new superior to that utilized in the Laplace transform ( N B R U L ), as noted by Mahmoud et al. (2018) [14]; El-Arishy et al. (2017) discuss the application of “new better than renewal” ( N B R U r p ) in the RP order [15]; Al-Zahrani and Stoyanov (2010) [16] discuss harmonic new better than renewal used in expectation ( H N B R U E ); Elbatal (2009) discusses renewal new better than renewal used ( R N B R U ) [17]; Mahmoud et al. (2008) [18] discuss the new better-than-average failure rate ( N B A F R ) and the new better-than-used failure rate ( N B U F R ); Abdel-Aziz (2007) [19] discusses renewal new better than renewal used in expectation ( R N B R U E ); Mahmoud et al. (2004) [20] discuss renewal new better than used; and Abouammoh et al. (1994) [21] discusses renewal new better than renewal used classes ( N B R U ).
The following is the paper’s composition: A novel category of life distributions known as new better than renewal is defined, and its relationships are used in the moment-generating function (NBRUmgf) in Section 2. In Section 2, a test statistic computed using the Laplace transform technique is provided for evaluating exponentiality against NBRUmgf. As a specific example, a test based on the goodness of fit approach is also suggested, and its properties are examined. Section 3 calculates PAE and PARE for a few widely used reliability distributions. In Section 4, we simulate the critical points of the Monte Carlo null distribution and tabulate the critical values of our two-test statistic.
For sample sizes n = 5 , 10 , 11 , 15 ( 5 )   40 , 43 , 45 , 50 . The power of two tests is also computed for a few widely used reliability distributions in Section 5. In Section 6, applications for full data are covered. In Section 7, we finally provide a conclusion for our work.

2. Moment-Generating Function Class Uses New Better than Renewal

Consider a lifetime component X with distribution function F ( x ) . A series of mutually identical components that are independent of the initial component will replace the component in the event of a failure. Over time, the remaining life distribution of the component operating at time t is referred to as the stationary renewal distribution if the system is renewed indefinitely. This can be defined as follows:
W F t = μ 1 0 t F ¯ u d u , t 0
The corresponding density function of the stationary renewal distribution w F t and the renewal survival function is given by
W ¯ F t = μ 1 t F ¯ u d u
The failure rate of the renewal distribution W F t is given by
λ w t = w F t W ¯ F t = F ¯ ( t ) t F ¯ u d u = (   μ t   ) 1 , 0 t < ,
where the mean time to failure ( μ ) is defined by
μ = 0 F ¯ t d t
The random residual at age t life is defined by X t = [ T t | T > t ] with survival function F ¯ t   x = F ¯ ( x + t ) F ¯ ( t )   ,   x ,   t 0 . The average residual life of X t is given by
μ t = t F ¯ ( x ) d x F ¯ ( t ) ,   t ,   x 0   ;   F ¯ ( t ) > 0
For more details, see Barlow and Proschan [3].
Definition 1.
If X is a random variable with a non-negative survival function F ¯ x , then X is said to have new better (worse) than renewal used distribution, denoted by F N B R U ( N W R U ), if
W ¯ F x + t W ¯ F t F ¯ x ;     x 0 , t 0
Definition 2.
A non-negative random variable X is said to have new better than renewal used in moment-generating function (denoted by F N B R U m g f ) if, and only if,
0 e s x W ¯ F x + t d x   W ¯ F t 0 e s x F ¯ x d x   ;   x ,   t   ,   s 0
Definition 3.
Put s = 0 in Equation (7); a non-negative random variable X is said to have new better than renewal used in the convex ordering (denoted by F N B R U C ) if, and only if,
x W ¯ F x + t d x   W ¯ F t x F ¯ x d x     ,   x ,   t 0
Then, we have
N B R U   N B R U C N B R U m g f

3. Testing Exponentiality Versus NBRUmgf Class for Complete Data

This section tests the null hypothesis. Hypothesis H 0   : F is exponential with mean μ F against an alternative hypothesis that H 1 : F belongs to N B R U m g f and is not exponential. We draw a sample X 1 , X 2 , X 3 , . X n     from a population with distribution F and use the symbol δ ( s , β ) as a measure of departure from exponentiality. According to Equation (7), the measure of departure from exponentiality based on the Laplace transform technique is defined by the following equation:
δ s , β = 0 e β t   [   W ¯ F t 0 e s x F ¯ x d x 0 e s x W ¯ F x + t d x   ]   d t ;   s ,   β 0
In the theorem below, the measure of departure from exponentiality for the NBRUmgf class is derived and assumed that every moment is temporary and exists.
Theorem 1.
If X is a random variable with distribution function F ( x ) and F ( x ) relates to the NBRUmgf class, then
δ s , β = φ ( s ) ( β ) s β 2 + μ s β 1 s 2 β + s 1 s β 2 φ s 1 s β 2 1 β 2 β + s β + 1 s 2 β
where
φ s = 0 e s x d F ( x )   ;   β = 0 e β x d F ( x )
Proof. 
Since F is N B R U m g f , then by multiplying Equation (7) by e β t   and integrating both sides with respect to t over [0, ∞), we find
0 e β t 0 e s x W ¯ F x + t d x d t     0 e β t   W ¯ F t d t 0 e s x F ¯ x d x
Setting
I 1 = 0 e s x F ¯ x d x = E [ 0 I X > x e s x d x ] = E [ 0 x e s x d x ]   = 1 s φ s 1
I 2 = 0 e β t W ¯ F t d t = E 0 e β t X t I X > t d t = E 0 x e β t X t d t   = 1 μ ( μ   β + β β 2 1 β 2   )
I 3 = 0 e s v W ¯ F v d v = E 0 x e s v X v I X > v d v = 1 μ   μ s + φ s s 2 1 s 2
R . H . S .   o f   Equation 11 = I 1 I 2 = 1 μ [   φ s μ s β + φ s β s β 2 φ s s β 2 μ s β β s β 2 + 1 s β 2   ] ,
where
I X > x = { 0 i f   x X 1 i f   x < X
Let u = x + t , then d u = d x , and v = t , and then d v = d t ; J = 1 ;
L . H . S   o f Equation 11 = 0   e β v v   e s ( u v )   W ¯ F ( u ) d u d v = 0   e β u 0 v   e s ( v u )   W ¯ F ( v )   d u d v = 0   e s v   W ¯ F ( v ) 0 v   e u ( β + s )   d u d v = 1 ( β + s ) 0   e s v   W ¯ F ( v ) d v 0   e β v   W ¯ F ( v ) d v = 1 ( β + s ) I 3 I 2 = 1 μ μ s β + φ ( s ) s 2 ( β + s ) ( β ) β 2 ( β + s ) + ( s β ) β 2 s 2
Substituting from Equations (15) and (16) into Equation (9), we obtain Equation (10).
Note that under H 0 : δ s , β = 0 , while under H 1 : δ s , β > 0 . □
Corollary 1.
If we set β = 1 in Equation (10), then the following equation becomes the measure of departure from exponentiality based on the goodness of fit approach:
δ s , 1 = φ s ( 1 ) ) s + μ s 1 s 2 1 + s 1 s φ s 1 s 1 1 + s ( 1 ) + 1 s 2 ,
and note that under H 0 : δ s , 1 = 0 , while under   H 1 : δ s , 1 > 0 .

3.1. Test Statistic Empirical for NBRUmgf

To estimate δ s , β , let X 1 , X 2 , X 3 , . X n   be a sample selected at random from a distributed population F. To make the test statistic δ ^ s , β scale invariant, we set
s , β = δ s , β μ ¯
Suppose that F ¯ n x , a n d   ^ n s , β denote the empirical distribution of the survival function F ¯ ( x ) and the empirical estimate of s , β , respectively. We find
^ n s , β = 1 n 2 X ¯   i = 1 n j = 1 n ( X i , X j ) ;
X i , X j = e s X i   ×   e β X j s β 2 + X i s β 1 s 2 β + s 1 s β 2 e s X j 1 s β 2 1 β 2 β + s e β X i + 1 s 2 β
F ¯ n x = 1 n i = 1 n I ( X i > x ) ,                   d F x = 1 n
Note that X i , X j is not symmetric, and we define the symmetric kernel as
Ψ X 1 , X 2 = 1 2 ! R X i , X j
where the summation is over all arrangements of X i , X j ; this leads to ^ s , β corresponding to the U-statistic provided by the next equation (see Lee [4]).
U n = 1 n 2 i < j Ψ X i , X j
The asymptotic normality of ^ n s , β may be summed up by the upcoming theorem.
Theorem 2.
(1) 
n ,   n ( ^ n s , β s , β ) is asymptotically normal with zero mean and variance σ 2 where σ 2 is given by
σ 2 = V a r   [   e s X β + e β X   φ s s β 2 1 s 2 β + s + 1 s β 2 φ s + e s X + 2 s 2 β 1 s β 2 1 β 2 β + s e β X + ( β )   + X φ ( s ) + μ e s X s β ]
(2) 
Under H 0 , the variance σ 0 2 is
σ 0 2 = a 8 1 + β a 9 1 s + a 10 1 s 2 + 2 a 5     s β 1 s 3 + a 5 2 1 2 s + a 11 + a 12
where
a 1 = 1 s 2 ( β + s ) + 1 s β 2 ;   a 2 = 1 s β 2 1 β 2 ( β + s ) ;   a 3 = 1 ( 1 s ) ( 1 + β ) ; a 4 = 1 ( 1 + β ) ( 1 2 s ) ;   a 5 = 1 s β a 1 ;   a 6 = 1 ( 1 + β s ) ;   a 7 = 1 ( 1 + 2 β ) ; a 8 = a 4 s 2 β 4 2 a 2 a 6 + a 3 s β 2 4 a 2 a 3 s β 8 a 2 s 2 β ;   a 9 = 2 2 a 1 a 3 + a 2 a 7 a 5 a 6 s β 2 8 a 5 s 2 β a 10 = 2 a 3 s 2 β 3 + 2 s 2 β 2 2 a 1 s β + + a 7 s 2 β 4 + 3 a 2 2 + a 1 2 2 a 1 a 5 ; a 11 = 2 a 3 a 6 s 2 β 4 + 2 a 4 a 5 s β 2 + 6 a 3 s 3 β 3 + 2 a 3 2 s 2 β 3 + 4 s 4 β 2 ; a 12 = 2 a 1 a 2 a 3 3 a 2 a 3 a 5 + a 7 a 2 2 2 a 2 a 5 a 6
Proof. 
First, as n , we have X ¯ P μ and ^ n s , β P s , β . Using the results of Serfling [5], we note that as n ;   ^ n s , β has an asymptotically normal with a mean of zero and a variance of σ 2 which is given by
σ 2 = V a r ( ϕ 1 X + ϕ 2 ( X ) )
ϕ 1 X = E X 1 , X 2 X 1 ;   ϕ 2 X = E X 1 , X 2 X 2
Recall the definition of X 1 , X 2 in Equation (19); then, it is easy shown that
ϕ 1 X = e s X       E   e β X s β 2 + X s β 1 s 2 β + s 1 s β 2 E e s X 1 s β 2 1 β 2 β + s e β X + 1 s 2 β
ϕ 2 X = E e s X       e β X s β 2 + μ s β 1 s 2 β + s 1 s β 2 e s X 1 s β 2 1 β 2 β + s E e β X + 1 s 2 β
From Equations (26) and (27), Equation (22) can be proved.
Under H 0 , where F 0 x = 1 e x , x 0 , it is easy to prove that
μ 0 = E ϕ 1 X + ϕ 2 X = 2 β + 2 s 2 β + 2 s s β 2 ( 1 + β ) ( 1 s ) ( β + s ) + 2 ( s β 2 + s 2 β β 2 s β s 2 + β 2 β 2 s + β s β s 2 ) s 2 β 2 β + s 1 s = 0
and the variance σ 0 2 = E ϕ 1 X + ϕ 2 X 2 reduces to Equation (23). □
Corollary 2.
If we option to β = 1 in Equation (23), we obtain the variance in the case of goodness of fit approach as
σ 0 2 ( s , 1 ) = a 8 2 a 9 1 s + a 10 1 s 2 + 2 a 5     s 1 s 3 + a 5 2 1 2 s + a 11 + a 12

3.2. Pitman’s Asymptotic Efficiency

In order to assess this procedure’s quality, Pitman asymptotic efficiencies (PAEs) of our test are computed for the alternative distributions below in the class (since they are in the N B R U m g f ). Then, we perform comparison testing of various famous classes of life distribution with our suggested test based on PAEs.
(i)
Linear failure rate distribution (LFR):
F ¯ 1   x = e X θ 2 X 2 ,   x 0 ,   θ 0
(ii)
Makeham distribution:
F ¯ 2   x = e x θ ( X + e X 1 ) ,   x 0 ,   θ 0
(iii)
Weibull distribution:
F ¯ 3   x = e X θ   ,   x 0 ,   θ 0
Note that, for θ = 1, F ¯ 3   x   decreases the exponential distribution, but in the case of θ = 0 ,   F ¯ 1   x and F ¯ 2   x reduce to exponential distribution. The P A E δ θ s , β is defined by
PAE δ θ s , β = 1 σ 0 d δ θ s , β d θ θ θ 0
where
δ θ s , β = φ θ ( s ) θ ( β ) s β 2 + μ θ s β 1 s 2 β + s 1 s β 2 φ θ s 1 s β 2 1 β 2 β + s θ β + 1 s 2 β
The P A E δ θ s , β can be written as
PAE δ θ s , β = 1 σ 0 1 s β 2 φ ´ θ s 0 e β x d F θ x + 0 e s x d F θ x   ´ θ β + μ θ s β 1 s 2 β + s 1 s β 2 φ θ ´ s + μ ´ θ s β φ θ s 1 s β 2 1 β 2 β + s ´ θ β θ θ 0
Upon using Equations (28) and (29), we obtain the P A E δ θ s , β of these above three alternative distributions of our procedure, which are, respectively,
(i)
Linear failure rate distribution ( θ 0 = 0 ) :
P A E δ θ s , β , F ¯ 1   x = 1 σ 0 a 3 ( 1 s ) + s β 2 a 5 β 2 1 s 3 1 s β 1 s + a 3 ( 1 + β ) s β 2 a 2 s p 1 + β 3
(ii)
Makeham distribution ( θ 0 = 0 ) :
P A E δ θ s , β , F ¯ 2   x   = 1 σ 0 a 3 ( 1 + β ) 2 s β 2 a 2 s β ( 2 + β ) 1 + β 3 ( 1 s ) a 3 + 2 β 2 a 5 β 2 2 s 1 s 2 1 2 s β 1 s
(iii)
Weibull distribution ( θ 0 = 1 ) :
P A E δ θ s , β , F ¯ 3   x = 1 σ 0     1 s + s l o g 1 s 1 s a 3 + s β 2 a 5 s β 2 1 s 2   0.422784 s β 1 s     + 1 + β + β l o g 1 + β 1 + β a 3 s β 2 a 2 s β 2 1 + β 2    
By using Equations (30)–(32) and the Mathematica program 12.3, we obtain Table 1 below, which shows our suggested test’s Pitman asymptotic efficiencies δ θ s , β for LFR, Makeham, and Weibull distributions at different values of s ,   β .
From Table 1, it is shown that these alternative distributions have maximum values of δ θ s , β at s = 0.9   ,   β = 2 and minimum values at s = 0.58 ,   β = 5 .
In Table 2, we display Pittman’s asymptotic efficiencies for the proposed test statistics   δ θ s , β and other test statistics.
Using Equation (33), we obtain Table 3 below, which illustrates our test statistics’ Pitman’s asymptotic efficiency ( δ ^ n ) with respect to other test statistics (T).
P A R E = e   δ θ s , β ) ,   T =   P A E ( δ θ ) P A E ( T )  
From Table 2 and Table 3, it is obvious that our test statistic is more efficient than other tests because e δ θ ,   T   > 1 .

4. Critical Points of the Monte Carlo Null Distribution

The present part presents the Monte Carlo critical points for the null distribution for ^ n s , β in Equation (18), which are simulated based on 10,000 samples of sizes n = 5,10,11,15 ( 5 )   40,43,46,50 from the standard exponential distributions using the Mathematica software 12.3. Table 4 and Table 5 provide the statistic’s upper percentile points. ^ n s , β with 10,000 replications at s = 0.9 ,   β = 2 and s = 0.58 , β = 1 ,   respectively.
Observations reveal that the critical values in Table 4 and Table 5 increase as the confidence level rises and falls with the growing sample size, with the exception of n = 40, 46 for a confidence level of 99% at s = 0.9 ,   β = 2 in Table 4.
It is clear that the critical values in Figure 1 and Figure 2 rise with increasing confidence levels and fall with increasing sample sizes, with the exception of n = 40, 46 for a confidence level of 99% at s = 0.9 ,   β = 2 in Figure 1.

5. The Estimated Power

The power of a test is defined as the probability of rejecting H 0 when it is false. In the calculation of the power from simulated samples, we calculate the proportion of times H 0 is rejected when it is false. The power of ^ n s , β is estimated at a confidence level of 1 α % ,   α = 0.05 with suitable parameters of θ = 1 ,   2 ,   3 ,   4 ,   5 at n = 10 ,   20 ,   30 . This is in regard to the three selections: failure rate in a linear approach and Makeham and Weibull distributions based on 10,000 replications.
From Table 6 and Table 7, it is clear that the power estimate of our test is good, and it increases when the value of the parameter θ increases and the sample size n increases for LFR, Makeham, and Weibull distributions. Note that the critical values in Table 4 and Table 5 increase as the sample size decreases. This means that the critical values increase when the value of the parameter θ decreases. Therefore, the power shown in Table 6 and Table 7 decreases as the N B R U m g f approaches the exponential distribution, because the proportion of times H 0 is rejected when it is false decreases with the decrease in sample size and the value of the parameter θ .

6. Applications for Complete Data

In this section, we apply the test on some real data sets to elucidate the applications of the NBRUmgf at a 95% confidence level.
Data Set 1
Consider the data set below in Abouammoh et al. [21], these statistics pertain to a group of forty Saudi Arabian Ministry of Health hospital patients with leukemia, a kind of blood cancer. In this case, we find ^ n s = 0.9 , β = 2 = 2.645587 and ^ n s = 0.58 , β = 1 = 1.160749   which is more than the critical values of 1.811924 of Table 4 and 0.213249 of Table 5, respectively, at the 95% upper percentile. Then, we accept H 1 , the alternative hypotheses, which states that the data set has NBRUmgfproperty and that it is not exponential.
0.3150.4960.6161.1351.2081.2631.4142.0252.036
2.1622.2112.3702.5322.6932.8052.9102.9123.192
3.2633.3483.3483.4273.4993.5343.7673.7513.858
3.9864.0494.2444.3234.3814.3924.3974.6474.753
4.9294.9735.0744.381
Data Set 2
Consider the following data set given in Grubbs [22]. This data set provides the times between the arrivals of 25 customers at a facility. In this case, we find ^ n s = 0.9 , β = 2 = 3.755778 and ^ n s = 0.58 , β = 1 = 5.131664 which are greater than the critical values of 2.661309 from Table 4 and 0.337393 from Table 5, respectively, at the 95% upper percentile. Hence, we accept H 1   ,   which states that the data set has NBRUmgf property, and it is not exponential.
1.802.892.933.033.153.433.483.57
3.853.923.984.064.114.134.164.23
4.344.374.534.624.654.844.914.99
5.17
Data Set 3
Consider the following data set which has been given in Lawless [23]. This data set shows the failure times, expressed in hours, for a particular kind of electrical insulation under continually rising voltage stress in an experiment. Here, we find ^ n s = 0.9 , β = 2 = 0.52761 and ^ n s = 0.58 , β = 1 = 0.038167 , which are less than the critical values of 5.054081 from Table 4 and 0.708224 from Table 5, respectively, at the 95% upper percentile. Hence, we accept the null hypothesis, which states that the data set has exponential property.
0.2050.3630.4070.7700.7200.7821.1781.2551.592
1.6352.310
Data Set 4
Consider the data set which has been given by Kots and Johnson [24], which represents the survival times (in years) after the diagnosis of 43 patients with a certain kind of leukemia. In this case, we find ^ n s = 0.9 , β = 2 = 4.251812 and ^ n s = 0.58 , β = 1 = 1.239598 , which are greater than the critical values of 1.700669 from Table 4 and 0.195653 from Table 5, respectively, at the 95% upper percentile. Hence, we accept H 1 , which states that the data set has NBRUmgf property and is not exponential.
0.0190.1290.1590.2030.4850.6360.7480.781
0.8691.1751.2061.2191.2191.2821.3561.362
1.4581.5641.5861.5921.7811.9231.9592.134
2.4132.4662.5482.6522.9513.0383.6003.655
3.7454.2034.6904.8885.1435.1675.6035.633
6.1926.6556.874
The present paper deals with a hypothesis testing problem for testing exponentiality against the N B R U m g f class. Two empirical measures of departure from exponentiality
^ n s , β = 1 n 2 X ¯   i = 1 n j = 1 n ( X i , X j )   ; and   ^ n s , 1 = 1 n 2 X ¯   i = 1 n j = 1 n X i , X j ;
where
X i , X j = e s X i   *   e β X j s β 2 + X i s β 1 s 2 β + s 1 s β 2 e s X j 1 s β 2 1 β 2 β + s e β X i + 1 s 2 β
has been obtained by applying the Laplace transform technique and goodness of fit approach, respectively.
Authors in [10,11] dealt with a hypothesis testing problem for testing exponentiality against the N B R U m g f class. Two empirical measures of departure from exponentiality
δ ^ n = 1 n 2 s 3 i n   j n   6 s X i + 3 s 2 X i 2 + 6 + s 3 X i 3 6 e s X j ,  
and
δ ^ n = 1 n 2 j n [ s 2 X i 2 e s X j 2 e s X i + 2 s X i + 2 ]
were obtained by applying moment inequalities and the U-statistic, respectively.
All authors in the present paper [10,11] obtained the same result (they reject the null hypothesis) which states that data sets 1, 2, and 4 have the NBRUmgf property and are not exponential (Large Data).
All authors in the present papers, refs. [8,10,11,13,25], obtained the same result (they accept the null hypothesis) which states that data set 3 has the exponential property (Small Data). Also, all authors in [7,8,9,13,25,26,27] obtained the same result (they reject the null hypothesis), which states that at least one of data set 1, data set 2, or data set 4 has the UBAC(2)L, NBUCmgf, RNBUL, NBUL, ODL, and NBURL properties, and are not exponential (Large Data), respectively.
All authors in [14] obtained a different result (they accept the null hypothesis), which states that data set 2 has exponential property (Large Data), and also in [28], they obtained a different result (they accept the null hypothesis), which states that data set 1 has the exponential property (Large Data.)
So, our test does not perform equally well in small and large samples.

7. Conclusions

Based on the Laplace transform technique and goodness of fit approach, new two-test statistics for assessing exponentiality versus a new class of life distributions termed new better than renewal used in the moment-generating function ( N B R U m g f ) is examined. The quality requirements of the two tests are validated using the widely used Pitman asymptotic efficiency criterion.
Pitman’s asymptotic efficiencies for the suggested tests are computed for a number of well-known alternative reliability distributions, including the Makeham, Weibull, and LFR distributions. It has been demonstrated that the test statistics we suggested are more effective than the others. The performance of the suggested tests is evaluated by computing and tabulating the upper percentiles and power. Lastly, the utility of our tests is demonstrated by applying them to some actual data.

Author Contributions

Conceptualization, N.A.H., M.M.S., R.A.E.-W.A. and T.R.; Methodology, N.A.H., M.M.S., R.A.E.-W.A. and T.R.; Software, N.A.H., M.M.S., R.A.E.-W.A. and T.R.; Validation, N.A.H., M.M.S., R.A.E.-W.A. and T.R.; Formal analysis, N.A.H., M.M.S., R.A.E.-W.A. and T.R.; Investigation, N.A.H., M.M.S., R.A.E.-W.A. and T.R.; Resources, N.A.H., M.M.S., R.A.E.-W.A. and T.R.; Data curation, N.A.H., M.M.S., R.A.E.-W.A. and T.R.; Writing—original draft, N.A.H., M.M.S., R.A.E.-W.A. and T.R.; Writing—review & editing, N.A.H., M.M.S., R.A.E.-W.A. and T.R.; Visualization, N.A.H., M.M.S., R.A.E.-W.A. and T.R.; Supervision, N.A.H., M.M.S., R.A.E.-W.A. and T.R.; Project administration, N.A.H., M.M.S., R.A.E.-W.A. and T.R.; Funding acquisition, T.R. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Data are contained within the article.

Acknowledgments

The Researchers would like to thank the Deanship of Graduate Studies and Scientific Research at Qassim University for financial support (QU-APC-2024-9/1).

Conflicts of Interest

The authors declare no conflicts of interest.

References

  1. Barlow, R.E.; Marshall, A.W.; Proschan, F. Properties of probability distributions with monotonic hazard rate. Ann. Math. Stat. 1963, 34, 348–350. [Google Scholar] [CrossRef]
  2. Bryson, M.C.; Siddiqui, M.M. Some criteria for ageing. J. Am. Stat. Assoc. 1969, 64, 1472–1483. [Google Scholar] [CrossRef]
  3. Barlow, R.E.; Proschan, F. Statistical Theory of Reliability and Life Testing; To Begin With: Silver Spring, MD, USA, 1981. [Google Scholar]
  4. Lee, A.J. U-Statistics; Marcel-Dekker: New York, NY, USA, 1990. [Google Scholar]
  5. Serfling, R.J. Approximation Theorems of Mathematical Statistics; John Wiley & Sons: New York, NY, USA, 1980. [Google Scholar]
  6. Everestus, O.O.; Mbanefo, S.M.; Madukaife, A.V.; Oladugba, V.A. A review of tests for exponentiality with Monte Carlo comparisons. Appl. Stat. 2022, 49, 1277–1304. [Google Scholar]
  7. Abu-Yousef, S.E.; Hassan, E.M.A.; Gerges, S.T. A new nonparametric class of life distributions based on ordering moment generating approach. J. Stoch. Appl. Probab. Lett. 2020, 7, 151–162. [Google Scholar]
  8. Abu-Yousef, S.E.; Gerges, S.T. Based on the goodness of fit approach. A new test statistic for testing NBUCmgf class of life distributions. Pak. J. Stat. 2022, 38, 129–144. [Google Scholar]
  9. Abu-Yousef, S.E.; EL-Toony, A.A. A new class of life distribution based on Laplace transform and its applications. Inf. Sci. Lett. 2022, 11, 355–362. [Google Scholar]
  10. Hassan, E.M.A.; Said, M.M. A new class of life distribution based on moment inequalities. Asian J. Probab. Stat. 2021, 13, 47–57. [Google Scholar] [CrossRef]
  11. Hassan, E.M.A.; Said, M.M. Testing exponentially against new better than renewal used in moment generating function class order based on U-statistics. J. Stat. Appl. Probab. 2019, 8, 7. [Google Scholar]
  12. El-Arishy, S.M.; Diab, L.S.; El-Atfy, E.S. Testing exponentiality against RNBUmgf based on Laplace transform technique. J. Stat. Appl. Probab. 2019, 8, 229–239. [Google Scholar] [CrossRef]
  13. Mahmoud, M.A.W.; Diab, L.S.; Radi, D.M. Testing exponentiality against renewal new better than used in Laplace transform order, applications in reliability. J. Egypt. Math. Soc. 2019, 27, 49. [Google Scholar] [CrossRef]
  14. Mahmoud, M.A.W.; El Sagheer, R.M.; Etman, W.B.H. Moment’s inequalities for NBRUL distributions with hypotheses testing applications. Austrian J. Stat. 2018, 47, 95–104. [Google Scholar]
  15. El-Arishy, S.M.; Diab, L.S.; El-Atfy, E.S. Moments Inequalities for NBRUrp Distributions with Hypotheses Testing Applications. Int. J. Comput. Appl. 2017, 159, 34–40. [Google Scholar]
  16. Al-Zahrani, B.; Stoyanov, J. Moment inequalities for HNBRUE with hypothesis testing application. Stat. Methodol. 2010, 7, 58–67. [Google Scholar] [CrossRef]
  17. EL-Battal, I.I. On testing statistics of renewal new better than renewal used. Int. J. Contemp. Math. Sci. 2009, 4, 17–29. [Google Scholar]
  18. Mahmoud, M.A.W.; Abdul Alim, A.N. A goodness of fit approach for testing NBUFR (NWUFR) and NBAFR (NWAFR) properties. Int. J. Reliab. Appl. 2008, 1, 125–140. [Google Scholar]
  19. Abdul Aziz, A.A. On testing exponentiality against RNBRUE alternatives. Appl. Math. Sci. 2007, 1, 1725–1736. [Google Scholar]
  20. Mahmoud, M.A.W.; EL-Arishy, S.M.; Diab, L.S. Testing renewal new better than used life distributions based on U-test. Appl. Math. Models 2004, 29, 784–796. [Google Scholar] [CrossRef]
  21. Abouammoh, A.M.; Abdulghani, S.A.; Qamber, I.S. On partial orderings and testing of new better than renewal used classes. Reliab. Eng. Syst. Saf. 1994, 43, 37–41. [Google Scholar] [CrossRef]
  22. Grubbs, F.E. Fiducial bounds on reliability for the two-parameter negative exponential distribution. Technometrics 1971, 13, 873–876. [Google Scholar] [CrossRef]
  23. Lawless, J.F. Statistical Models & Methods for Lifetime Data; John Wiley & Sons: New York, NY, USA, 1982. [Google Scholar]
  24. Stephens, M.; Kots, S.; Johnson, N.L. Encyclopedia of Statistical Sciences; Wiley: New York, NY, USA, 1982. [Google Scholar]
  25. Diab, L.S.; Kayid, M.; Mahmoud, M.A.W. Moments inequalities for NBUL distributions with hypotheses testing applications. Contemp. Eng. Sci. 2009, 2, 319–332. [Google Scholar]
  26. Diab, L.S.; El-Atfy, E.S. A moment inequality for overall decreasing life class of life distributions with hypothesis testing applications. Int. J. Math. Stat. Invent. 2017, 5, 62–71. [Google Scholar]
  27. Mahmoud, M.A.W.; EL-Sagheer, R.M.; Etman, W.B. Testing exponentiality against new better than renewal used in Laplace transform order. J. Stat. Appl. Pro. 2016, 5, 279–285. [Google Scholar] [CrossRef]
  28. Mahmoud, M.A.W.; Moshref, M.E.; Gadallah, A.M. On NBUL class at specific age. Int. J. Reliab. Appl. 2014, 15, 11–22. [Google Scholar]
Figure 1. Relation between sample size, critical values, and confidence intervals of statistic ^ n s , β at s = 0.9 ,   β = 2 .
Figure 1. Relation between sample size, critical values, and confidence intervals of statistic ^ n s , β at s = 0.9 ,   β = 2 .
Mathematics 12 02045 g001
Figure 2. Relation between sample size, critical values, and confidence intervals of statistic ^ n s , β at s = 0.58 ,   β = 1 .
Figure 2. Relation between sample size, critical values, and confidence intervals of statistic ^ n s , β at s = 0.58 ,   β = 1 .
Mathematics 12 02045 g002
Table 1. The Pitman asymptotic efficiencies of δ θ s , β at different values of s, β.
Table 1. The Pitman asymptotic efficiencies of δ θ s , β at different values of s, β.
Distributionβs = 0.58s = 0.59s = 0.6s = 0.7s = 0.8s = 0.9
LFR11.78171.44081.28971.07891.22681.6276
21.11151.17301.24252.51886.616930.933
30.96661.02731.09472.31206.282930.698
40.87930.93741.00172.15955.980729.925
50.81560.87120.93262.03865.725329.189
Makeham14.34953.49153.10332.46762.65353.3851
21.21621.24161.27181.79862.99986.8103
31.00701.03441.06521.55932.67996.3496
40.90550.93300.96321.43742.51606.1036
50.83640.86340.89291.35082.39715.9251
Weibull11.50971.21911.08870.85300.79480.6836
20.18150.20750.23560.68921.89696.8722
30.12420.14790.17340.59451.76716.7673
40.10280.12470.14840.54411.66036.5830
50.09020.11070.13290.50811.58386.4152
Table 2. The (PAEs) for LFR, Makeham, and Weibull distributions.
Table 2. The (PAEs) for LFR, Makeham, and Weibull distributions.
TestLFRMakehamWeibull
Abu-Yousef et al. (2022) [8]12.89860.74623.532
Abu-Yousef et al. (2022) [9]1.30000.5800---------
Hassan et al. (2021) [10]8.301592.904182.2092
Abu-Yousef et al. (2020) [7]5.6650.87652.1679
Mahmoud et al. (2019) [13]0.91500.17200.6180
El. Arishy et al. (2019) [12]0.956580.203950.73343
Our   test   δ θ s = 0.9 , β = 2 30.93306.81036.8722
Table 3. The Pitman’s asymptotic relative efficiencies (PAREs) of our test.
Table 3. The Pitman’s asymptotic relative efficiencies (PAREs) of our test.
PARELFRMakehamWeibull
e ( δ θ s = 0.58 , β = 1 ) ,   T   i n [8]) 2.3989.1273.532
e ( δ θ s = 0.58 , β = 1 ) ,   T   i n [9]) 23.79511.742-------
e ( δ θ s = 0.9 , β = 2 ,   T   i n [10]) 3.7262.3453.111
e ( δ θ s = 0.9 , β = 2 ,   T   i n [7]) 5.4607.7703.170
e ( δ θ s = 0.9 , β = 2 ,   T   i n [13]) 33.80739.59511.120
e ( δ θ s = 0.9 , β = 2 ,   T   i n [12]) 32.33733.3929.370
Table 4. The upper percentile points of the statistic ^ n s , β with 10,000 replications at s = 0.9, β = 2.
Table 4. The upper percentile points of the statistic ^ n s , β with 10,000 replications at s = 0.9, β = 2.
n90%95%99%
56.8324719.11140920.881602
103.7678765.49263916.612772
113.4877695.05408115.462071
152.7672894.16304012.836489
202.1643773.0871909.591994
251.8412002.6613098.859419
301.5766252.3017667.199750
351.3667941.9590735.227894
401.2206471.8119246.043173
431.1273701.7006695.076861
461.0570421.5150085.439105
500.9974641.4791344.638035
Table 5. The upper percentile points of the statistic ^ n s , β with 10,000 replications at s = 0.58, β = 1.
Table 5. The upper percentile points of the statistic ^ n s , β with 10,000 replications at s = 0.58, β = 1.
n90%95%99%
51.4255261.6548812.512848
100.6887240.8008921.272633
110.6224190.7082241.075039
150.4639810.5501871.034439
200.3521170.4152520.704349
250.2829360.3373930.573153
300.2354790.2769300.457279
350.2042160.2440690.409616
400.1789880.2132490.340378
430.1654150.1956530.328344
460.1547610.1812040.276341
500.1452530.1697420.279583
Table 6. The power estimates of ^ n s = 0.9 , β = 2 at α = 0.05.
Table 6. The power estimates of ^ n s = 0.9 , β = 2 at α = 0.05.
nθLFRMakehamWeibull
1010.10270.33160.4630
20.42020.97970.6214
30.68391.0000.7641
40.83951.0000.9673
50.97281.0000.9754
2010.10560.36870.5571
20.45900.99850.6431
30.76791.0000.7883
40.88931.0000.9780
50.97501.0000.9832
3010.10770.47720.5597
20.56160.99920.6782
30.80661.0000.8542
40.89051.0000.9812
50.98341.0000.9931
Table 7. The power estimates of ^ n s = 0.58 , β = 1 at α = 0.05.
Table 7. The power estimates of ^ n s = 0.58 , β = 1 at α = 0.05.
nθLFRMakehamWeibull
1010.19650.26210.5335
20.40990.99770.7452
30.62201.0000.8682
40.77751.0000.9786
50.99021.0000.9860
2010.31050.49320.6900
20.70651.0000.7743
30.92121.0000.9123
40.98271.0000.9892
50.99701.0000.9976
301041930.65970.7342
20.86531.0000.8547
30.98701.0000.9678
40.99961.0000.9901
51.0001.0000.9987
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

Hassan, N.A.; Said, M.M.; Attwa, R.A.E.-W.; Radwan, T. Applying the Laplace Transform Procedure, Testing Exponentiality against the NBRUmgf Class. Mathematics 2024, 12, 2045. https://doi.org/10.3390/math12132045

AMA Style

Hassan NA, Said MM, Attwa RAE-W, Radwan T. Applying the Laplace Transform Procedure, Testing Exponentiality against the NBRUmgf Class. Mathematics. 2024; 12(13):2045. https://doi.org/10.3390/math12132045

Chicago/Turabian Style

Hassan, Naglaa A., Mayar M. Said, Rasha Abd El-Wahab Attwa, and Taha Radwan. 2024. "Applying the Laplace Transform Procedure, Testing Exponentiality against the NBRUmgf Class" Mathematics 12, no. 13: 2045. https://doi.org/10.3390/math12132045

APA Style

Hassan, N. A., Said, M. M., Attwa, R. A. E. -W., & Radwan, T. (2024). Applying the Laplace Transform Procedure, Testing Exponentiality against the NBRUmgf Class. Mathematics, 12(13), 2045. https://doi.org/10.3390/math12132045

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