Stochastic comparisons of extreme order statistic from dependent and heterogeneous lower-truncated Weibull variables under Archimedean copula

: This article studies the stochastic comparisons of order statistics with dependent and heterogeneous lower-truncated Weibull samples under Archimedean copula. To begin, we obtain the usual stochastic and hazard rate orders of the largest and smallest order statistics from heterogeneous and dependent lower-truncated Weibull samples under Archimedean copula. Second, under Archimedean copula, we get the convex transform and the dispersive orders of the largest and smallest order statistics from dependent and heterogeneous lower-truncated Weibull samples. Finally, several numerical examples are given to demonstrate the theoretical conclusions.


Introduction
Order statistics are important in many probability areas, including reliability theory, auction theory, and operations research. Denote X k:n by the k-th order statistic of random variables X 1 , X 2 , ..., X n , k = 1, 2, ..., n. In reliability theory, X k:n characterizes the lifetime of a (n − k + 1)-out-of-n system. Specifically, X n:n and X 1:n can express as the lifetimes of parallel and series systems, respectively. In auction theory, X n:n and X 1:n represent the first-price sealed-bid auction and final price of the first-price procurement auction, respectively. Exploring the stochastic behavior of coherent systems is one of the most important subjects in reliability. Stochastic orders, which are valuable instruments for measuring the size and variability of random variables, have been frequently employed in stochastic comparisons of reliability systems.
Many scholars have spent decades studying stochastic comparisons of order statistics from heterogeneous and independent distributions. The usual stochastic order, the hazard rate order, and the likelihood ratio order are all stochastic orders that compare the "magnitude" of random variables. For example, [1] studied the largest order statistic in terms of the hazard rate and the dispersive order. [2] obtained ordering properties of the largest order statistic with two independent heterogeneous exponential samples with respect to the likelihood ratio order and the hazard rate order. [3] investigated the likelihood ratio order and the stochastic order between the smallest order statistic of independent and dependent samples. [4] focused on stochastic orders to compare the magnitudes of two parallel systems from Weibull distributions when one set of scale parameters majorizes the other. There have been lots of extensions works to the case of stochastic comparisons with independent heterogeneous Weibull distribution, for example, [5][6][7]. Besides, we study stochastic orders that compare the "variability" or the "dispersion" of random variables. The most important and common orders are the convex and dispersive orders. The variability of order statistics from heterogeneous and independent random variables has nice applications in the reliability theory and actuarial science, and has been studied by many researchers. For instance, the convex transform, star, and dispersive orders are used for comparing variability in probability distributions. [8] proved the star order between the largest order statistic from heterogeneous(homogeneous) and independent proportional hazard rates(PHR) models. [9] investigated the stochastic comparisons of the largest and the smallest order statistics with independent heterogeneous generalized exponential samples in terms of various stochastic orders. [10] gave some sufficient conditions for stochastic comparisons between the largest order statistic with exponentiated Weibull samples with respect to the usual stochastic, the likelihood ratio order, and the dispersive order. [11] compared two of the largest order statistics having heterogeneous exponentiated Weibull samples in terms of the reversed hazard rate and likelihood ratio orders. [12] considered stochastic comparisons of the smallest order statistic from the location-scale family of distributions with respect to different stochastic orders. [13] focused on some stochastic comparisons between the corresponding order statistics based on modified proportional hazard rates and modified proportional reversed hazard rates models. [14] presented some new ordering properties between two parallel systems comprising general independent heterogeneous samples in the sense of the usual stochastic and reversed hazard rate orders. [15] dealt with some stochastic comparisons of both the largest and the smallest order statistics comprising dependent Burr Type XII samples under the Archimedean copula with respect to the star order and the convex transform order. [16] considered the largest order statistic with Pareto samples and studied the effect of heterogeneity on the skewness of such systems. [17] established some properties of the new measures for various classes of symmetric and asymmetric distributions, and characterized the generalized Pareto distribution in terms of the convex transform order. For more investigations, one may refer to [18][19][20][21][22][23][24][25][26][27][28][29][30][31].
All of the above studies are restricted to the independent case. However, in most practical scenarios, due to the common environment or share the same workload, which results in the dependence among the samples. In recent years, the dependent samples of order statistics have attracted widespread attention in the academic world. [32] investigated the ordering properties of order statistics from random variables of Archimedean copulas. [33] studied order statistics from random variables following the scale model with respect to the usual stochastic order of the sample extremes and the second smallest order statistic, the dispersive order, and the star order of the sample extremes. [34] further obtained the stochastic properties of the largest and the smallest order statistics with heterogeneous and dependent PHR samples in the sense of the usual stochastic order. [35] discussed stochastic comparisons of the largest and the smallest order statistics with heterogeneous and dependent resilience-scaled samples. [36] provided sufficient conditions for the hazard rate order on the smallest and proportional hazard rates or scales, and the reversed hazard rate order on the largest of the sample with Archimedean copula. [37] obtained the variability of both the largest and the smallest order statistics of heterogeneous samples, sufficient conditions are established for the dispersive and the star orders between the smallest order statistic consisting of dependent samples having multiple-outlier proportional hazard rates and Archimedean copulas. [38] focused on the usual stochastic, star, and convex transform orders of extreme order statistic comparing heterogeneous and dependent extended exponential samples under Archimedean copulas. [39] carried out stochastic comparisons of the largest and the smallest order statistics with dependent heterogeneous Topp-Leone generated samples in terms of the usual stochastic order and the reversed hazard rate order. [40] investigated stochastic comparisons on extreme order statistic from heterogeneous and dependent samples following modified proportional reversed hazard rated and modified proportional hazard rates models. [41] provided distribution-free results to compare, in the usual stochastic order under some majorization conditions, coherent systems with heterogeneous and dependent components where the dependency structure can be defined by any copula. The stochastic comparisons with statistically dependent samples has attracted widespread attention among a lot of scholars. The interested readers may refer to [42][43][44][45][46][47][48][49][50][51][52][53][54].
Truncated Weibull distributions are comprehensively applied to product reliability, modeling product failures, quality control and life testing, so it is very meaningful that we research the ordering statistics results of truncated Weibull models. As discussed in [55,56], in many sampling settings, observations are limited to a subset of the population's potential values. As a result of the limiting process, only n of N prospective observations remain visible, while (N-n) are deleted. When N(or N-n) is known, the resulting collection of partial data is referred to as censored; otherwise, it is referred to as truncated. A defined number of objects are generally evaluated for a fixed amount of time in life testing, which is aimed to estimate the average life span of, for example, transistors from a certain production line. Because the lives of things surviving the life test are unknown, this results in a censored sample of lifetimes. In certain real-life testing circumstances, the total number of items on the exam is unknown. This happens, for example, when a certain unknown number of products are put through a life test and one of them has a specific fault that is only discovered after the item fails. If the lifespan of an item with this specific defect is the variable of interest, the sample is truncated in the sense that the number of missing observations with lifetime larger than the burn-in or testing period is unknown. Therefore, conducting sample research with Lower-truncated Weibull distribution has very important theoretical needs and practical significance. [57] discussed the stochastic comparisons of extreme order statistics with heterogeneous and independent lower-truncated Weibull samples. A random variable X is said to have doubly truncated Weibull distribution if its cumulative distribution function is . When T = 0 and L → ∞, it becomes two-parameter Weibull distribution. When T = 0 and L < ∞, it is upper-truncated Weibull distribution, and when T > 0 and L → ∞, it is lower-truncated Weibull distribution. For more details on Weibull distribution and its truncated forms, one may refer to [58,59].
In this paper, we will focus attention on the lower-truncated Weibull with cumulative distribution function and probability density function respectively. We denote X ∼ LT W(α, λ) if X has the distribution functions G(x). Motivated by the work of [57], this paper focus on studying the usual stochastic, the hazard rate, the convex transform and the dispersive orders of the extreme order statistic from dependent and heterogeneous lower-truncated Weibull samples. We derive the magnitude and variability with the heterogeneity considered in the shape and scale parameter. The difference from [57] is that we extend the independent situation to the dependent situation. The rest of this paper is organized as follows. In Section 2, we review some definitions of stochastic orders, majorizations, and some lemmas. In Section 3, we discuss the usual stochastic order and the hazard rate order of the smallest and the largest order statistics with dependent and heterogeneous lower-truncated Weibull samples. In Section 4, we get the convex transform and the dispersive orders of the largest and the smallest order statistics with dependent and heterogeneous lower-truncated Weibull samples. Finally, Section 5 concludes the paper with some remarks.

Preliminaries
Before proceeding to the main results, we briefly introduce some basic concepts about stochastic orders, majorization, and Archimedean copula, which are very useful tools to compare random variables arising from reliability theory, operations research, actuarial science, and so on. Furthermore, the notion 'a sgn = b' means that a and b have the same sign.

Stochastic orders
Let X be a non-negative random variable with distribution function F X (t), survival function F X (t) = 1 − F X (t), probability density function f X (t), the hazard rate function h X (t) = f X (t)/F X (t), and the right-continuous inverse F −1 X . Definition 1. A non-negative random variable X is said to be smaller than Y in the sense of It is known that but the reversed statement is not true in general. The convex transform is usually used to describe skewness. Skewness is an important quantity used for measuring the asymmetry in a distribution. For a unimodal distribution, nagative skewness indicates that the tail on the left side of the density function is longer than the right side, while positive skewness indicates that the tail on the right side of the density function is longer than the left side. The convex transform order in this regard implies that one distribution is more skewed to the right than the other distribution. It should be noticed that the convex transform suitable for a non-negative random variable, hence the implication X ≤ c Y ⇒ X ≤ * Y holds when X and Y are two non-negative random variables (see [60]). For more comprehensive discussions on various stochastic orders and their applications, one may refer to [61,62].

Majorization order
Majorization order is an important tool in establishing various inequaties arising from many research areas. Denote I n = 1, 2, ..., n.
Definition 2. For two real vectors a = (a 1 , a 2 , ..., a n ) and For more details on majorization, weak majorization and Schurconvex (Schur-convave) functions, one may refer to [63]. The sufficient and necessary conditions for the characterization of Schur-convex (Schur-concave) function is presented in the following lemma.
where ∂ψ(x)/∂x i denotes the partial derivative of ψ with respect to its i-th argument.

Archimedean copula
Archimedean copula has been used in reliability theory and many other areas due to its capability of capturing wide ranges of dependence and mathematical tractability. By definition, for a continuous and decreasing function ϕ : is called an Archimedean copula with the generator ψ if (−1) n−2 ψ (n−2) (x) is decreasing and convex and For more discussions on copulas and their properties, one may refer to [64,65].

Stochastic comparisons of the magnitude
In this section, we provide the ordering properties of the largest and the smallest order statistics with dependent and heterogeneous lower-truncated Weibull samples in terms of the usual stochastic order and the hazard rate order. We denote X ∼ LT W(α, λ, ψ) as the sample arising from non-negative random variables X 1 , X 2 , ..., X n assembled with an Archimedean copula having generator ψ, where X i ∼ LT W(α i , λ i ) for i = 1, 2, ..., n.

Orders of the largest order statistic
First, Theorem 1 gives sufficient conditions for the largest order statistic from dependent and heterogeneous lower-truncated Weibull samples in terms of the usual stochastic order.
Remark 2. It should be noted that for independent samples, ψ(x) = e −x , thus Theorem 1 and Theorem 2 generalizes Theorem 3.2 of [57] from independent samples to the case of statistically dependent components with Archimedean copulas. Remark 3. Proposition 3.7 in [41] is considerably more generic, because the random variables studied in such proposition have arbitrary distribution functions, however, in this paper, we explore a specific sort of distribution function, namely the lower-truncated Weibull distribution. In comparison to Proposition 3.7 in [41], the Theorem 1 and the Theorem 2 do not need the requirement ψ 1 is log-convex, hence our conditions are weaker for lower-truncated Weibull distribution in the largest order statistic.
We present the following Example 2 to illustrate the result of Theorem 2.

Orders of the smallest order statistic
First, we will give a stochastic comparison of the smallest order statistic from dependent and heterogeneous lower-truncated Weibull samples in terms of the usual stochastic order. Corollary 1. Let X ∼ LT W(α, λ, ψ 1 ) and Y ∼ LT W(β, λ, ψ 2 ). If ϕ 2 • ψ 1 be super-additive and ψ 1 is log-convex. Then, for λ > 0, x ≥ 1/λ. We have Proof. The distribution function of X is increasing in α, and the survival function of X is log-concave in α, which satisfies conditions in Proposition 3.16 of [41]. Hence, according to Proposition 3.16 of [41], we can easily establish the result of Corollary 1.  It is apparent to the ϕ 2 • ψ 1 be super-additive and ψ 1 is log-convex. The whole of survival function curves of X 1:3 and Y 1:3 are display in Figure 3, which confirmed that X 1:3 ≤ st Y 1:3 . Corollary 2. Let X ∼ LT W(α, λ, ψ 1 ) and Y ∼ LT W(α, µ, ψ 2 ). Suppose ϕ 2 • ψ 1 be super-additive.
Proof. (i) The distribution function of X is increasing in λ, and the survival function of X is logconcave in λ for 0 < α ≤ 1, which satisfies conditions in Proposition 3.16 of [41]. As a result, Corollary 2(i) can be obtained immediately from proposition 3.7 in [41]. (ii) The proof of Corollary 2(ii) can be obtained similarly that of Corollary 2(i), thus we omit the proof process of Corollary 2(ii). □ In the following, we will give a numerical example to illustrate the Corollary 2.
As is seen in Figure 4(ii), the survival function curve of X 1:3 is always beneath that of Y 1:3 , that is, X 1:3 ≤ st Y 1:3 , which verified Corollary 2 (ii).  The following Theorem 3 carries out stochastic comparison of the smallest order statistics from dependent and heterogeneous lower-truncated Weibull samples in the sense of the hazard rate order.
We present the following Example 5 to illustrate the result of Theorem 3.
As is seen in Figure 5, the hazard rate function curve of Y 1:3 is always beneath that of X 1:3 , that is, X 1:3 ≤ hr Y 1:3 , which coincides Theorem 3.   Figure 6, the difference between the hazard rate functions h X 1:n and h Y 1:n is not always non-negative and this means that Y 1:n ⩽̸ hr X 1:n . So the result in Theorem 3 couldn't hold if α ≥ 1/2.

Stochastic comparisons based on variability orders
In this section, we study the convex transform order and the dispersive order of the smallest and the largest order statistics from dependent and heterogeneous lower-truncated Weibull samples.

Stochastic comparisons based on the convex transform order
First, we compare the smallest and the largest order statistics between dependent and heterogeneous-homogeneous samples with respect to the convex transform order.
Proof. (i) The survival functions of X 1:n and Y 1:n are given bȳ respectively. It holds that The first and the second partial derivatives of F −1 Y 1:n (F X 1:n (x)) with respect to x are and it is sufficient to have 0 < α ≤ 1 and ψ is log-concave.
Next counterexample 2 explains the result in Theorem 4 (i) couldn't hold if α ≥ 1. As is seen in Figure 7, the curve of ∂ 2 F −1 Y 1:2 (F X 1:2 (x))/∂x 2 is not positive and this means that Y 1:2 ⩾̸ c X 1:2 . Hence the result in Theorem 4 couldn't hold if α ≥ 1. Theorem 4 indicates that the larger and the smaller order statistics, from dependent homogeneous lower-truncated Weibull samples, are more skewed to the right than those from dependent heterogeneous lower-truncated Weibull samples.

Stochastic comparisons based on the dispersive order
Aside from the convex transform order, which commonly measures skewness, the dispersive order also measures skewness from a different standpoint. Theorem 5 looks at the dispersive order of the smallest order statistic from dependent heterogeneous lower-truncated Weibull samples and dependent homogeneous lower-truncated Weibull samples with the identical Archimedean copulas.
From (4.5) and (4.6), we have f Y 1:n F −1 Y 1:n F X 1:n (x) ≤ f X 1:n (x). Hence the theorem follows. □ The following Corollary 4 derives from Theorem 5 and Theorem 6. Corollary 4 compares the smallest of two samples, one group from n dependent and heterogeneous lower-truncated Weibull samples and another group from n dependent and homogeneous lower-truncated Weibull samples with different Archimedean copulas.
Proof. Let Z i ∼ LT W(α, λ)(i = 1, 2, ..., n) and the associated Archimedean copula is with generator ψ 1 . Then from Theorem 5, we have Z 1:n ≥ disp X 1:n , and from Theorem 6, we have Y 1:n ≥ disp Z 1:n . □ In reliability, the smallest order statistics represents a series system. Corollary 4 indicates that the aging rate of dependent and homogeneous lower-truncated Weibull samples is usually faster than dependent and heterogeneous lower-truncated Weibull samples.

Applications
The study of ordering findings of order statistics is critical in many practical domains, including actuarial science, reliability theory, auction theory and multivariate statistics. In this part, we will look at how the established theoretical conclusions may be put to use. In industrial engineering areas. These are referred to as the parallel and series systems, respectively. Consider two series systems with n dependent components that are modeled using lower-truncated Weibull models.
More heterogeneous scale parameters in the weakly supermajority order yield a series system with a stochastically longer lifetime, according to Theorem 1. Similar findings may be drawn from Theorems 2 for more heterogeneous reciprocals of location parameters in terms of weak supermajority orders.

Conclusions
We studied the usual stochastic order, the hazard rate order, the dispersive order, and the convex transform order in this paper considering dependent and heterogeneous lower-truncated Weibull samples of both the biggest and smallest order statistics. It is highly significant to compare distinct extremes order statistics using the convex transform order, which may be used to examine the relative aging properties of systems from diverse perspectives. These findings generalize several previous findings in the literature. The fundamental restriction of the researched topic may be the absence of treatment of the kth order statistic, which has a significant practical impact. However, due to the complexities of modeling for statistically dependent order statistics, these intriguing problems remain unsolved and demand more investigation.
research was supported by the National Natural Science Foundation of China (No. 11861058 and No. 71471148).