In this article, we study the complete convergence and the complete moment convergence for negatively dependent (ND) random variables under sub-linear expectations. Under proper conditions of the moment of random variables, we establish the complete convergence and the complete moment convergence. As applications, we obtain the Marcinkiewcz-Zygmund type strong law of large numbers of ND random variables under sub-linear expectations. The results here generalize the corresponding ones in classic probability space to those under sub-linear expectations.
Citation: Mingzhou Xu, Xuhang Kong. Note on complete convergence and complete moment convergence for negatively dependent random variables under sub-linear expectations[J]. AIMS Mathematics, 2023, 8(4): 8504-8521. doi: 10.3934/math.2023428
Related Papers:
[1]
M. G. M. Ghazal, Yusra A. Tashkandy, Oluwafemi Samson Balogun, M. E. Bakr .
Exponentiated extended extreme value distribution: Properties, estimation, and applications in applied fields. AIMS Mathematics, 2024, 9(7): 17634-17656.
doi: 10.3934/math.2024857
[2]
Tabassum Naz Sindhu, Zawar Hussain, Naif Alotaibi, Taseer Muhammad .
Estimation method of mixture distribution and modeling of COVID-19 pandemic. AIMS Mathematics, 2022, 7(6): 9926-9956.
doi: 10.3934/math.2022554
[3]
M. G. M. Ghazal .
Modified Chen distribution: Properties, estimation, and applications in reliability analysis. AIMS Mathematics, 2024, 9(12): 34906-34946.
doi: 10.3934/math.20241662
[4]
Refah Alotaibi, Hassan Okasha, Hoda Rezk, Abdullah M. Almarashi, Mazen Nassar .
On a new flexible Lomax distribution: statistical properties and estimation procedures with applications to engineering and medical data. AIMS Mathematics, 2021, 6(12): 13976-13999.
doi: 10.3934/math.2021808
[5]
Haidy A. Newer, Bader S Alanazi .
Bayesian estimation and prediction for linear exponential models using ordered moving extremes ranked set sampling in medical data. AIMS Mathematics, 2025, 10(1): 1162-1182.
doi: 10.3934/math.2025055
[6]
Abdulhakim A. Al-Babtain, Amal S. Hassan, Ahmed N. Zaky, Ibrahim Elbatal, Mohammed Elgarhy .
Dynamic cumulative residual Rényi entropy for Lomax distribution: Bayesian and non-Bayesian methods. AIMS Mathematics, 2021, 6(4): 3889-3914.
doi: 10.3934/math.2021231
[7]
Bing Long, Zaifu Jiang .
Estimation and prediction for two-parameter Pareto distribution based on progressively double Type-II hybrid censored data. AIMS Mathematics, 2023, 8(7): 15332-15351.
doi: 10.3934/math.2023784
[8]
M. R. Irshad, S. Aswathy, R. Maya, Amer I. Al-Omari, Ghadah Alomani .
A flexible model for bounded data with bathtub shaped hazard rate function and applications. AIMS Mathematics, 2024, 9(9): 24810-24831.
doi: 10.3934/math.20241208
[9]
Abdullah Ali H. Ahmadini, Amal S. Hassan, Ahmed N. Zaky, Shokrya S. Alshqaq .
Bayesian inference of dynamic cumulative residual entropy from Pareto Ⅱ distribution with application to COVID-19. AIMS Mathematics, 2021, 6(3): 2196-2216.
doi: 10.3934/math.2021133
[10]
Emrah Altun, Mustafa Ç. Korkmaz, M. El-Morshedy, M. S. Eliwa .
The extended gamma distribution with regression model and applications. AIMS Mathematics, 2021, 6(3): 2418-2439.
doi: 10.3934/math.2021147
Abstract
In this article, we study the complete convergence and the complete moment convergence for negatively dependent (ND) random variables under sub-linear expectations. Under proper conditions of the moment of random variables, we establish the complete convergence and the complete moment convergence. As applications, we obtain the Marcinkiewcz-Zygmund type strong law of large numbers of ND random variables under sub-linear expectations. The results here generalize the corresponding ones in classic probability space to those under sub-linear expectations.
1.
Introduction
The main objective of statistical inference is to draw conclusions about a certain population based on a random sample of that population. Conclusions may be about the functional form of the distribution for the population of interest, estimates of its characteristics, testing statistical hypotheses regarding parameters or making predictions about future events based on current knowledge.
One of the essential elements of reliability analysis is testing the product failure times under normal operating conditions. Products are frequently highly dependable, and extensive applied testing is prohibitively costly. The use of complete samples is not a good choice in statistical analysis when the product is highly reliable and applied testing is prohibitively expensive. In such situations, statisticians have recommended many types of censoring to save time and money. Type Ⅱ-right censorship is one of these censoring types and has been more popular over the past few decades. It is expected for this type that a few lower-order statistics are observed and that subsequent inferences are then required, such as determining the distribution from which the data is collected, estimating the distribution parameters and making predictions about unobserved events. Here, the emphasis is on estimate for censored samples as well as prediction problems.
There are many contributions to the estimation of complete and censored data. Excellent references for this subject, among others, include Balakrishnan and Cramer [1], Casella and Berger [2], Lehmann and Casella [3] and the references therein.
Prediction of future observations is a major concern in many real-world problems. The prediction for ordered random variables (RVs) is frequently used in industrial applications and survival studies to estimate the future prevalence of defective products. In lifetimes testing experiments, the interval and point predictions are also helpful in determining the best censoring strategy. To put it more precisely, we can place identical items in a given life-test experiment and wait until a manageable number of failed items are obtained in relation to cost and time. We can therefore predict the failure times of the survivor items using these observed failure times. This prediction enables the experimenter to choose an appropriate censoring scheme. Such censoring schemes may be Type Ⅰ or Type Ⅱ censoring, progressive Type Ⅱ censoring, or hybrid censoring. The prediction has additional benefits that help us decide whether or not the life test needs to be accelerated and when a future failed item should cause the life test to end. Many authors have considered the prediction issue in statistical literature from both a theoretical and an applied standpoint. Ahsanullah [4] studied the linear prediction of record values of the two-parameter exponential distribution. Bayesian predictions of both ordinary order statistics (OOSs) and generalized order statistics (GOSs) were studied in AL-Hussaini [5] and AL-Hussaini and Ahmed [6], respectively. In Aly [7], two-point predictors of the fractional kth upper record values from the exponential distribution are given. Some predictive and reconstructive results of dual GOSs from the inverse Weibull distribution are obtained in Aly [8] via the pivotal quantities approach. The problem of predicting future lifetimes from the Weibull probability model for a simple step-stress plan under the Khamis-Higgins model is studied by Amleh and Raqab [9]. Barakat et al. [10] suggested a new method for predicting future order statistics based on some of their distributional properties. Moreover, by applying the cumulative hazard transformation, Barakat et al. [11] constructed nonparametric prediction intervals for GOSs, and their exact coverage probabilities were determined. For all distributions that are attracted to the Weibull distribution, El-Adll et al. [12] have constructed different asymptotic prediction intervals for future lower extreme order statistics. The maximum likelihood predictor was first introduced by Kaminsky and Rhodin [13]. In life-testing experiments with fixed sample sizes, Lawless [14] suggested a pivotal quantity for constructing a prediction interval for an order statistic. Lingappaiah [15] suggested a different pivotal quantity for the same purpose. Raqab [16] obtained an optimal prediction interval for m-GOSs based on a one-parameter exponential distribution. Concise reviews that are related to this topic and include several prominent prediction results may be found in the books by Aitcheson and Dunsmore [17], David and Nagaraja [18], Geisser [19] and Kaminsky and Nelson [20].
The main contribution of the paper is presented in Section 3. More specifically, based on Rényi representation and its extensions, a new least squares method for estimation and prediction is introduced, relying on the cumulative hazard function. The proposed method is formulated in a general manner so that it can be applied to any model of ordered random variables. Furthermore, the method is used for complete or type Ⅱ right-censored data and simplifies the computations significantly for those distributions with an explicit form for the cumulative hazard function.
The rest of this paper is organized as follows: In Section 2, some necessary preliminary results that are related to the present work are presented. A new least-squares method for estimation and prediction is presented in Section 3. The methodology presented in Section 3 is validated in Section 4 through comprehensive simulation studies. In Section 5, two real data sets are analyzed for illustrative and comparison purposes.
2.
Preliminary results
As a unified model of various ascending models of ordered RVs, Kamps [22] established the GOSs, based on any continuous distribution function (CDF) F, via their joint probability density function (PDF)
on the cone F−1(0+)<x1≤⋯≤xn<F−1(1), where Cr=r∏i=1γi,r=1,2,...,n,F−1 denotes the quantile function of F. The model parameters are defined by the vector γ_=(γ1,...,γn), where γr=k+n−1∑i=r(mi+1)>0,i=1,2,...,n,withγn=k>0and(m1,...,mn−1)∈Rn−1. Particular models, such as OOSs, sequential order statistics (SOSs), progressive type Ⅱ censored order statistics, standard kth record values and Pfeifer's record, result from specific choices of the model parameters, γ1,...,γn. When mi=m for all i=1,2,...n, extensive asymptotic results for bootstrapping m-GOSs are obtained by Barakat et al. [21].
According to the restriction γi≠γj for i≠j, a wide subclass of GOSs-model, excluding record values, have discussed by Kamps and Cramer [23], where some important distributional properties have been considered. In particular, the PDF of the rth GOS X∗r:n is given by
f∗r:n(x)=Crr∑i=1ai(r)[1−F(x)]γi−1f(x),−∞<x<∞,
(2.2)
and the CDF of X∗r:n is
F∗r:n(x)=1−Crr∑i=1ai(r)γi[1−F(x)]γi,−∞<x<∞,
(2.3)
where
ai(r)=r∏j=1,j≠i(γj−γi)−1,1≤i≤r≤n.
In view of the probability integral transform, the RVs, U∗r:n=F−1(X∗r:n),r=1,2,...,n, are uniform GOSs. Consequently, it can be shown that
The ordinary least squares method of estimation for complete samples was originally proposed by Swain et al. [24]. The method is based on minimizing the function
LF(Θ|˜xn):=n∑i=1(F(xi;Θ)−E[Ui:n])2,
(2.8)
with respect to the unknown parameter vector Θ=(θ1,θ2,...,θℓ), where ˜xn=(x1,x2,...,xn) is an observed ordered sample. The weighted least squares estimates can be accomplished by minimizing the function
with respect to Θ. Several authors used the least squares method and the weighted least squares method. For estimating the parameters from different distributions; among them are Gupta and Kundu [25] and Kundu and Raqab [26].
El-Adll and Aly [27] have extended the above results to type Ⅱ censored samples. Namely, based on the first r observed OOSs, x1,x2,...,xr, approximate least square estimates of Θ can be obtained by minimizing the least squares function
˜xr=(x1,x2,...,xr) are observed values of the GOSs X∗1:n,X∗2:n,...,X∗r:n,E[U∗i:n] and Var(U∗i:n) are given by (2.4) and (2.5), respectively.
Kamps [22] extended Rényi's representation [29], to GOSs model (see also Barakat et al. [28]). This representation enables us to express the rth GOS based on the exponential distribution as a linear combination of independent and identically distributed (iid) RVs from the EXP(1). Namely,
X∗r:nd=r∑i=1Ziγi,r=1,2,...,n,
(2.13)
where Z1,...,Zn are iid RVs from the EXP(1) and "Ud=V" means that the RVs U and V have the same CDF.
3.
A new least squares method for estimation and prediction
Minimizing the least squares and weighted least squares functions in (2.11) and (2.12) is not always a simple problem. In this section, a novel and efficient method for solving such problems is proposed.
3.1. Least squares estimation via the cumulative hazard function
The cumulative hazard function of any RV X with CDF F is defined by
H(x)=−log[1−F(x)],−∞<x<∞.
Evidently, the function H(x) is a nonnegative and nondecreasing function on x. Suppose now that X∗i:n:=X(i,n,˜m,k) denotes the ith GOS from a continuous CDF F. Therefore, the RV's
Z∗i:n:=H(X∗i:n),i=1,2,...,n,
(3.1)
are GOSs from the standard exponential distribution (denoted by EXP(1). Clearly, if the CDF F depends on an unknown vector of parameters Θ=(θ1,θ2,...,θℓ), with ℓ≥1, then H depends on the same vector of parameters. Hence, we can develop a least square method based on H. In view of Kamps [22] and Barakat et al. [28], we have
Consequently, the parameters can be estimated by minimizing the least squares function
L∗H(Θ|˜xn):=n∑i=1(H(xi;Θ)−μ∗i:n)2,
(3.3)
where −∞<x1<x2<...<xn<∞ are observed values of the GOSs X∗1:n,X∗2:n,...,X∗n:n. Moreover, the weighted least squares estimates based on the cumulative hazard function can be obtained by minimizing the weighted least squares function.
WL∗H(Θ|˜xn):=n∑i=1w∗i(H(xi;Θ)−μ∗i:n)2.
(3.4)
According to the extended Rényi's representation (2.13), we have
where z1,...,zn is an observed random sample from the EXP(1). Clearly, E[∑ij=r+1Zj−1γj]=0. Therefore, we can approximate the sum of the last term by its mean. Consequently, we can suggest the modified least squares function.
An approximate least squares estimate of Θ based on the first r observed GOSs, ˜xr:=(x1,x2,...,xr) for r≤n, can be obtained by minimizing L∗H,r(Θ|˜xr) with respect to Θ. Similarly, minimizing the function
with respect to Θ and xs. Similarly, a weighted least squares predictor can be derived by choosing the weights wi defined by (2.12). As we proceeded in (3.5) and (3.6), an approximate point predictor of the unobserved sth GOS can be obtained via the minimization of the proposed predictive least squares function
By the same manner, approximate weighted least squares estimates of Θ and xs based on x1,x2,...,xr are derived via minimizing the predictive weighed least squares function
One advantage of choosing the cumulative hazard function transformation is that it always follows the standard exponential distribution. In addition, its mean and variance can simply be computed for any model of ordered random variables, and they do not depend on unknown parameters. Moreover, as a quick comparison between the proposed method and one of the most well-known methods of estimation, we find that the Maximum likelihood estimation (MLE) for the distribution parameters may be difficult to obtain in certain cases. Particularly when the support of the distribution is unknown. Also, the MLEs may not be robust enough to depart from the assumed distribution. These considerations motivated the least-squares approach to be used.
4.
Numerical simulation experiments
In this section, simulation experiments are carried out to compare the proposed method with different estimation and prediction techniques. We are mainly interested in some important probability distributions that can be widely applied in survival analysis, reliability theory and life-testing experiments. For brevity, we compare only three methods, namely the maximum likelihood, the ordinary weighted least squares and the weighted least squares via the cumulative hazard function. The following assumes that the first r (r≤n) elements of the GOSs based on a continuous distribution F are observed and used to estimate the unknown distribution parameters and predict some future observations. In this simulation, 10,000 independent random samples are generated from F and then used in the proposed estimation and prediction methods.
4.1. Estimation
In view of Kamps [22], the likelihood function based on the first r elements of the GOSs, which is the joint PDF of X∗1:n,X∗2:n,...,X∗r:n, is given by
The maximum likelihood estimates (MLEs) of the unknown parameters Θ=(θ1,θ2,...,θℓ) can be accomplished via maximizing L∗(Θ|˜xr) by solving the nonlinear equations ∂log(L∗(Θ|˜xr))∂θj=0,j=1,2,...,ℓ, numerically. The ordinary weighted least squares estimates (WLSEs) are obtained by minimizing the function WL∗F,r(Θ|˜xr) given in (2.12), through solving the nonlinear equations ∂WL∗F,r(Θ|˜xr)∂θj=0,j=1,2,...,ℓ, numerically. Finally, the modified weighted least squares estimates (MWLSEs) can be obtained by minimizing the function WL∗H,r(Θ|˜xr) in (3.7), via solving the nonlinear equations ∂WL∗H,r(Θ|˜xr)∂θj=0,j=1,2,...,ℓ, numerically. The OOSs model, where γi=n−i+1, is primarily the focus of the next sections of the study.
We use the root mean square error to compare different estimators or predictors. As a result of its sensitivity to extreme values, the root-mean-square error may not even exist. Pitman's measure of closeness (Pitman [30]) is also used for comparing estimators or predictors. According to Keating et al. [31], who provided several inspiring instances and examples, Pitman's measure of closeness is an efficient criterion for selecting among estimators. Pitman's measure of closeness is widely applied to assess estimators and predictors by several authors. Pitman's measure is used by Balakrishnan et al. [32] for comparing different point predictors for type Ⅱ censored data that follows an exponential distribution in one sample and two sample cases. Nagaraja [33] used the mean square error and Pitman's measure of closeness to compare the best linear predictor with the best linear invariant predictor for the record value and order statistic. Raqab et al. [34] compared different point predictors of progressively censored units using Pitman's measure of closeness.
Let ˆθ1 and ˆθ2 be two estimators of the same parameter θ. Then, ˆθ1 is Pitman closer than ˆθ2 if
PMC(ˆθ1,ˆθ2):=P(|ˆθ1−θ|<|ˆθ2−θ|)⩾0.5,
for all values of θ, with strict inequality for at least one value of θ.
Pitman's measure of closeness is applied to explore which of two estimators is the closest (in probability) to the true value of a parameter.
4.1.1. Two parameter-exponential distribution
If the CDF of the RV X is given by
F(x;θ,σ)=1−e−(x−θσ),x≥θ,σ>0,
it is said to have a two-parameter exponential distribution, with a location parameter θ and a scale parameter σ. The MLEs of θ and σ based on the first r elements of GOSs are given in El-Adll [35] by
ˆθ=X∗1:nandˆσ=1rr∑i=2γi(X∗i:n−X∗i−1:n).
(4.2)
The WLSEs, ˜θ and ˜σ, are obtained by minimizing the function
subject to σ>0 and θ−a>0, where a is a suitable real number chosen, provided that it is less than the minimum of the data. Minimization can be accomplished by equating the first partial derivatives with respect to θ and σ with zero and then solving the resulting equations, numerically. The results are presented in Tables 1 and 2.
Table 1.
Different estimates of the parameters for the exponential distribution, EXP (2.5,10), along with their associated RMSEs.
The main reason for choosing the two-parameter exponential distribution in the simulation study is not only its theoretical importance but also because it has explicit forms for the maximum likelihood estimator of its parameters, and consequently, it can be considered a reference in comparison with the proposed method.
4.1.2. Three-parameter Weibull distribution
We consider the three-parameter Weibull distribution with CDF
F(x;α,σ,θ)=1−exp[−(x−θσ)α],x≥θ,α,σ>0,
where α,σ and θ are the shape, scale and location parameters, respectively. As in the exponential distribution, the MLEs, WLSEs and MWLSEs of the parameters are obtained numerically and shown in Tables 3 and 4.
Table 3.
Different estimates of the Weibull distribution parameters, Weibull(5,10,0.65), along with their associated RMSEs.
4.1.3. The two-parameter modified Kies exponential distribution
Al-Babtain et al. [36] introduced a new family of continuous probability distributions, which they called the "new modified Kies family." They discussed the two-parameter MKi-exponential (MKiExp) distribution as a special case in detail. In the same paper, they demonstrated the practical importance and flexibility of fitting several types of real data. The CDF of MKiExp distribution is given by
F(x;α,σ)=1−exp[−(exσ−1)α],x≥0,α,σ>0,
where α and σ are the shape and scale parameters, respectively. Abd El-Raheem et al. [37] consider the estimation problem of multiple constant-stress tests for progressive type-Ⅱ censored MKiExp data with binomial removals.
The problem of determining numerical estimates in this work is reduced to a suitable minimization problem subject to the following constraints: α>0 and σ>0. The results are given in Tables 5 and 6.
Table 5.
Different estimates of parameters for the MKiEx (0.5,10) distribution along with their associated RMSEs.
For OOSs we have μ∗i:n=∑ij=1(n−j+1)−1 and w∗i=[∑ij=1(n−j+1)−2]−1,i=1,2,...,n.
Remark 4.2.
By setting α=1, the above results are specialized to the two-parameter exponential distribution. The prediction problem of the two-parameter exponential distribution is discussed in details in Barakat et al. [38]. Moreover, θ=0 yields the two-parameter Weibull distribution.
4.2.2. Predicting future observations from the modified Kies exponential distribution
As in the Weibull distribution, the prediction of a future order statistic from the two-parameter modified Kies exponential distribution can be accomplished by minimizing the following three functions:
The minimization problem is subject to the constraints σ>0,α>0 and xs>xr.
Remark 4.3.
The RMSE and PMC are obtained numerically via a simulation, and all computations are performed through Mathematica 13.1.
In view of the simulation studies given above, the following comments are extracted:
1. In all cases, for the scale and shape parameters as well as the point predictors, the RMSEs decreased as r increased.
2. In most cases, the RMSE of the MLEs is smaller than the RMSEs of both the WLSEs and MWLSEs of the parameters.
3. For the Exp(θ,σ), according to PMC shown in Table 2, the MLE is the best followed by the MWLSE, which is followed by the WLSE for estimating the location parameter θ, while for estimating the scale parameter σ, the MWLSE is the best followed by the WLSE, which is followed by the MLE.
4. For the Weibull distribution, from Table 3, it is noted that
5. RMSE(ˆθ)<RMSE(˜θ)<RMSE(θ∗), RMSE(σ∗)≤RMSE(˜σ)<RMSE(ˆσ) and RMSE(α∗)<RMSE(ˆα) and RMSE(˜α)<RMSE(ˆα).
6. In view of PMC, Table 4 reveals that the WLSE is the best for the location parameter, while the MWLSEs of the scale and shape parameters are the best whenever n−r is small, but the WLSEs are the best whenever n−r is large. Moreover, for complete samples (i.e., r=n), the MWLSEs of the location, scale and shape parameters are the best.
7. According to the results presented in Tables 5 and 6, there is no preference for one method over the others when estimating the MKiExp distribution parameters.
8. In most cases, it is noted that the MWLSP and WLSP are better than the MLP, according to both the RMSE and PMC (see Tables 7–10).
Table 7.
Three point predictors of future unobserved order statistics from the exponential distribution with parameters θ=2 and σ=10, as well as their associated estimates and RMSEs.
Table 8.
Comparing different point predictors and their corresponding estimates of parameters for the two-parameter exponential with location parameter θ=2 and scale parameter σ=10, via PMC.
Table 9.
Three point predictors of unobserved future order statistics from the Weibull distribution with parameters σ=20 and α=1.5, as well as their associated estimates and RMSEs.
Table 10.
Comparing different point predictors and their corresponding estimates of parameters for the Weibull distribution with σ=20 and α=1.5, via PMC.
The windscreen on a large aircraft is a complex piece of equipment comprised of several layers. Failures of these items typically involve damage to or delamination of the heating system's nonstructural outer ply. These failures do not result in damage to the aircraft but do require replacement. Data of this type is incomplete in that all failure times have not yet been observed and may include failures to date of a particular model or combination of models. Murthy et al. [39] reported failure and service times for a specific model windscreen in Table 16.11 on page 297. The data represents 84 observed failure times for a specific windscreen device. Al-Babtain et al. [36] have shown that the MKiExp is an appropriate model for fitting this data.
We assume that the first r failure times have been observed, and we apply our prediction methods to predict the remaining failure times in two different scenarios. Three-point predictors, along with their related estimates and their relative errors, are obtained in the first scenario. Alternatively, to avoid complications in computations, in the second scenario, we first estimate the parameters based on the first r observed failure times and then compute three different point predictors for the future sth failure. In order to assess the prediction results, we compute the relative error (RE) for each point predictor. Recall that the RE is defined by RE=100×|˜Xs:n−Xs:n|Xs:n, where ˜Xs:n denotes the point predictor and Xs:n is the exact value of the quantity to be predicted. The results are presented in Table 11.
Table 11.
Different predictors, their associated estimations of parameters and the relative errors for xs:n,s=r+1,...,n(2) based on the first r order statistics with r=47 for the above failure times.
The second data set is reported and analyzed by Hoel [40] and Abu El Azm [41] among others. According to these data, male mice received a radiation dose of 300 roentgen at an age of 5–6 weeks. Each mouse's cause of death was identified by autopsy as either thymic lymphoma, reticulum cell sarcoma or other causes. Reticulum cell sarcoma is designated as cause 1 in this instance, and the other two causes of death are merged to make cause 2. There were 77 observations in this set of data, 39 of which are attributable to the second cause of death, while 38 are attributed to the first. For analysis purposes, we consider the following observations that are due to the first cause of death: 317,318,399,495,525,536,549,552,554,557,558,571,586,594,596,605,612,621,628,631,636,643,647,648,649,661,663,666,670,695,697,700,705,712,713,738,748,753.
It has been shown that the Makeham-Gompertz distribution is quite adequate for fitting reticulum cell sarcomas (e.g., Hoel [40]). Arguing as in the first data set, the prediction results are shown in Table 12.
Table 12.
Three different point predictors, the corresponding estimates and the relative errors of future OOSs xs:n,s=r+1,...,n based on the first r∈{22,26,30} OOS for Example 2.
The real data analysis given the the two preceding applications revealed that the maximum RE appears for the WLS method, followed by the MWLS method, followed by the ML method.
6.
Conclusions
In this article, by using the cumulative hazard function, a new least squares method for estimation and prediction has been proposed. The method is presented in a general setup so that it can be applied to any model of GOSs. A simulation study and numerical comparisons based on the RMSEs and PMCs have been performed through three important probability distributions. For applicability, two examples of real data sets are provided to illustrate the prescribed method. The comparisons reveal that the method is comparable with the ML and WLS methods, in the sense that there is no obvious preference for one method over the others for all estimation and prediction situations. Moreover, analyzing the real data revealed that the second scenario, in which we first estimate the unknown distribution parameters using Type Ⅱ right censoring and then predict future unobserved failures, performs better than the first scenario for the three methods.
Use of AI tools declaration
The authors declare they have not used Artificial Intelligence (AI) tools in the creation of this article.
Acknowledgments
The authors would like to express their sincere gratitude to the anonymous reviewers for their insightful comments and recommendations that raised the caliber of this paper.
This study is supported via funding from Prince Sattam Bin Abdulaziz University project number (PSAU/2023/R/1444).
Conflict of interest
The authors declare no conflict of interest.
References
[1]
S. G. Peng, G-expectation, G-Brownian motion and related stochastic calculus of Itô type, In: Stochastic analysis and applications, Abel Symposia, Vol. 2, Springer, Berlin, Heidelberg, 2007. https://doi.org/10.1007/978-3-540-70847-6_25
L. X. Zhang, Exponential inequalities under the sub-linear expectations with applications to laws of the iterated logarithm, Sci. China Math., 59 (2016), 2503–2526. https://doi.org/10.1007/s11425-016-0079-1 doi: 10.1007/s11425-016-0079-1
[4]
L. X. Zhang, Rosenthal's inequalities for independent and negatively dependent random variables under sub-linear expectations with applications, Sci. China Math., 59 (2016), 751–768. https://doi.org/10.1007/s11425-015-5105-2 doi: 10.1007/s11425-015-5105-2
[5]
L. X. Zhang, Donsker's invariance principle under the sub-linear expectation with an application to Chung's law of the iterated logarithm, Commun. Math. Stat., 3 (2015), 187–214. https://doi.org/10.1007/s40304-015-0055-0 doi: 10.1007/s40304-015-0055-0
[6]
J. P. Xu, L. X. Zhang, Three series theorem for independent random variables under sub-linear expectations with applications, Acta Math. Sin., English Ser., 35 (2019), 172–184. https://doi.org/10.1007/s10114-018-7508-9 doi: 10.1007/s10114-018-7508-9
[7]
J. P. Xu, L. X. Zhang, The law of logarithm for arrays of random variables under sub-linear expectations, Acta Math. Appl. Sin. English Ser., 36 (2020), 670–688. https://doi.org/10.1007/s10255-020-0958-8 doi: 10.1007/s10255-020-0958-8
[8]
Q. Y. Wu, Y. Y. Jiang, Strong law of large numbers and Chover's law of the iterated logarithm under sub-linear expectations, J. Math. Anal. Appl., 460 (2018), 252–270. https://doi.org/10.1016/j.jmaa.2017.11.053 doi: 10.1016/j.jmaa.2017.11.053
[9]
L. X. Zhang, J. H. Lin, Marcinkiewicz's strong law of large numbers for nonlinear expectations, Stat. Probab. Lett., 137 (2018), 269–276. https://doi.org/10.1016/j.spl.2018.01.022 doi: 10.1016/j.spl.2018.01.022
[10]
H. Y. Zhong, Q. Y. Wu, Complete convergence and complete moment convergence for weighted sums of extended negatively dependent random variables under sub-linear expectation, J. Inequal. Appl., 2017 (2017), 261. https://doi.org/10.1186/s13660-017-1538-1 doi: 10.1186/s13660-017-1538-1
[11]
Z. J. Chen, Strong laws of large numbers for sub-linear expectations, Sci. China Math., 59 (2016), 945–954. https://doi.org/10.1007/s11425-015-5095-0 doi: 10.1007/s11425-015-5095-0
[12]
X. C. Chen, Q. Y. Wu, Complete convergence and complete integral convergence of partial sums for moving average process under sub-linear expectations, AIMS Math., 7 (2022), 9694–9715. https://doi.org/10.3934/math.2022540 doi: 10.3934/math.2022540
[13]
L. X. Zhang, Strong limit theorems for extended independent and extended negatively dependent random variables under sub-linear expectations, Acta Math. Sci. English Ser., 42 (2022), 467–490. https://doi.org/10.1007/s10473-022-0203-z doi: 10.1007/s10473-022-0203-z
[14]
F. Hu, Z. J. Chen, D. F. Zhang, How big are the increments of G-Brownian motion, Sci. China Math., 57 (2014), 1687–1700. https://doi.org/10.1007/s11425-014-4816-0 doi: 10.1007/s11425-014-4816-0
[15]
F. Q. Gao, M. Z. Xu, Large deviations and moderate deviations for independent random variables under sublinear expectations, Sci. China Math., 41 (2011), 337–352. https://doi.org/10.1360/012009-879 doi: 10.1360/012009-879
[16]
A. Kuczmaszewska, Complete convergence for widely acceptable random variables under the sublinear expectations, J. Math. Anal. Appl., 484 (2020), 123662. https://doi.org/10.1016/j.jmaa.2019.123662 doi: 10.1016/j.jmaa.2019.123662
[17]
M. Z. Xu, K. Cheng, Convergence for sums of iid random variables under sublinear expectations, J. Inequal. Appl., 2021 (2021), 157. https://doi.org/10.1186/s13660-021-02692-x doi: 10.1186/s13660-021-02692-x
[18]
M. Z. Xu, K. Cheng, How small are the increments of G-Brownian motion, Stat. Probab. Lett., 186 (2022), 109464. https://doi.org/10.1016/j.spl.2022.109464 doi: 10.1016/j.spl.2022.109464
[19]
M. Z. Xu, K. Cheng, Note on precise asymptotics in the law of the iterated logarithm under sublinear expectations, Math. Probl. Eng., 2022 (2022), 6058563. https://doi.org/10.1155/2022/6058563 doi: 10.1155/2022/6058563
[20]
M. Z. Xu, K. Cheng, W. K. Yu, Complete convergence for weighted sums of negatively dependent random variables under sub-linear expectations, AIMS Math., 7 (2022), 19998–20019. https://doi.org/10.3934/math.20221094 doi: 10.3934/math.20221094
[21]
M. Yao, B. Q. Xiao, Equivalent conditions of complete convergence and complete moment convergence for END random variables, Chin. Ann. Math. Ser. B, 39 (2018), 83–96. https://doi.org/10.1007/s11401-018-1053-9 doi: 10.1007/s11401-018-1053-9
[22]
P. L. Hsu, H. Robbins, Complete convergence and the law of large numbers, Proc. Natl. Acad. Sci., 33 (1947), 25–31. https://doi.org/10.1073/pnas.33.2.25 doi: 10.1073/pnas.33.2.25
[23]
Y. S. Chow, On the rate of moment convergence of sample sums and extremes, Bull. Inst. Math. Acad. Sin., 16 (1988), 177–201.
[24]
M. H. Ko, Complete moment convergence of moving average process generated by a class of random variables, J. Inequal. Appl., 2015 (2015), 225. https://doi.org/10.1186/s13660-015-0745-x doi: 10.1186/s13660-015-0745-x
[25]
B. Meng, D. C. Wang, Q. Y. Wu, Convergence of asymptotically almost negatively associated random variables with random coefficients, Commun. Stat.-Theor. M., 2021 (2021), 1963457. https://doi.org/10.1080/03610926.2021.1963457 doi: 10.1080/03610926.2021.1963457
[26]
S. M. Hosseini, A. Nezakati, Complete moment convergence for the dependent linear processes with random coefficients, Acta Math. Sin., English Ser., 35 (2019), 1321–1333. https://doi.org/10.1007/s10114-019-8205-z doi: 10.1007/s10114-019-8205-z
[27]
B. Meng, D. C. Wang, Q. Y. Wu, Complete convergence and complete moment convergence for weighted sums of extended negatively dependent random variables, Commun. Stat.-Theor. M., 51 (2022), 3847–3863. https://doi.org/10.1080/03610926.2020.1804587 doi: 10.1080/03610926.2020.1804587
[28]
H. M. Srivastava, B. B. Jena, S. K. Paikray, A certain class of statistical probability convergence and its applications to approximation theorems, Appl. Anal. Discrete Math., 14 (2020), 579–598. https://doi.org/10.2298/AADM190220039S doi: 10.2298/AADM190220039S
[29]
H. M. Srivastava, B. B. Jena, S. K. Paikray, Statistical probability convergence via the deferred Nörlund mean and its applications to approximation theorems, RACSAM, 114 (2020), 1–14. https://doi.org/10.1007/s13398-020-00875-7 doi: 10.1007/s13398-020-00875-7
[30]
H. M. Srivastava, B. B. Jena, S. K. Paikray, U. K. Misra, Statistically and relatively modular deferred-weighted summability and Korovkin-type approximation theorems, Symmetry, 11 (2019), 1–20. https://doi.org/10.3390/sym11040448 doi: 10.3390/sym11040448
[31]
H. M. Srivastava, B. B. Jena, S. K. Paikray, Deferred Cesàro statistical probability convergence and its applications to approximation theorems, J. Nonlinear Convex Anal., 20 (2019), 1777–1792.
[32]
H. M. Srivastava, B. B. Jena, S. K. Paikray, Statistical deferred Nörlund summability and Korovkin-type approximation theorem, Mathematics, 8 (2020), 1–11. https://doi.org/10.3390/math8040636 doi: 10.3390/math8040636
[33]
S. K. Paikray, P. Parida, S. A. Mohiuddine, A certain class of relatively equi-statistical fuzzy approximation theorems, Eur. J. Pure Appl. Math., 13 (2020), 1212–1230. https://doi.org/10.29020/nybg.ejpam.v13i5.3711 doi: 10.29020/nybg.ejpam.v13i5.3711
[34]
L. X. Zhang, On the laws of the iterated logarithm under sub-linear expectations, Probab. Uncertain. Qua., 6 (2021), 409–460. https://doi.org/10.3934/puqr.2021020 doi: 10.3934/puqr.2021020
[35]
S. H. Sung, Moment inequalities and complete monent convergence, J. Inequal. Appl., 2009 (2009), 271265. https://doi.org/10.1155/2009/271265 doi: 10.1155/2009/271265
Mingzhou Xu, Xuhang Kong. Note on complete convergence and complete moment convergence for negatively dependent random variables under sub-linear expectations[J]. AIMS Mathematics, 2023, 8(4): 8504-8521. doi: 10.3934/math.2023428
Mingzhou Xu, Xuhang Kong. Note on complete convergence and complete moment convergence for negatively dependent random variables under sub-linear expectations[J]. AIMS Mathematics, 2023, 8(4): 8504-8521. doi: 10.3934/math.2023428
Table 7.
Three point predictors of future unobserved order statistics from the exponential distribution with parameters θ=2 and σ=10, as well as their associated estimates and RMSEs.
Table 8.
Comparing different point predictors and their corresponding estimates of parameters for the two-parameter exponential with location parameter θ=2 and scale parameter σ=10, via PMC.
Table 9.
Three point predictors of unobserved future order statistics from the Weibull distribution with parameters σ=20 and α=1.5, as well as their associated estimates and RMSEs.
Table 10.
Comparing different point predictors and their corresponding estimates of parameters for the Weibull distribution with σ=20 and α=1.5, via PMC.
Table 11.
Different predictors, their associated estimations of parameters and the relative errors for xs:n,s=r+1,...,n(2) based on the first r order statistics with r=47 for the above failure times.
Table 12.
Three different point predictors, the corresponding estimates and the relative errors of future OOSs xs:n,s=r+1,...,n based on the first r∈{22,26,30} OOS for Example 2.