New Improved Ranked Set Sampling Designs with an Application to Real Data

Amer Ibrahim Al-Omari and Ibrahim M.Almanjahie

1Department of Mathematics,Faculty of Science,Al al-Bayt University,Mafraq,Jordan

2Department of Mathematics,College of Science,King Khalid University,Abha,62529,Saudi Arabia

3Statistical Research and Studies Support Unit,King Khalid University,Abha,62529,Saudi Arabia

Abstract:This article proposes two new Ranked Set Sampling(RSS)designs for estimating the population parameters:Simple Z Ranked Set Sampling(SZRSS) and Generalized Z Ranked Set Sampling (GZRSS).These designs provide unbiased estimators for the mean of symmetric distributions.It is shown that for non-uniform symmetric distributions,the estimators of the mean under the suggested designs are more efficient than those obtained by RSS,Simple Random Sampling (SRS),extreme RSS and truncation based RSS designs.Also,the proposed RSS schemes outperform other RSS schemes and provide more efficient estimates than their competitors under imperfect rankings.The suggested mean estimators under perfect and imperfect rankings are more efficient than the linear regression estimator under SRS.Our proposed RSS designs are also extended to cover the estimation of the population median.Real data is used to examine wthe usefulness and efficiency of our estimators.

Keywords:Ranked set sampling;unbiased estimator;simple random sampling;mean squared error;efficiency;imperfect ranking

1 Introduction

The Ranked Set Sampling (RSS) is originally derived by McIntyre [1]as a new design to increase the efficiency of pasture and forage yields estimates for fixed sample units.The RSS is considered when the study variable can simply be ranked than quantified.Takahasi et al.[2],independently,introduced the background of the RSS design,mathematically.It is shown that mean estimator by the RSS is unbiased,and provides more efficient estimates than the simple random sampling mean estimator.Even when the measured observations are ranked with errors,the RSS still provides an unbiased estimator,but the imperfect ranking is generally better than ordering based on random [3].Stokes [4]considered the case of measuring the variable of interest and concluded that the study variable can be ranked by some concomitant variables.The competence of the estimator then depends on the relation between the study variables and the ancillary variables.With perfect ranking,the estimation based on RSS is more adequate as compared to the regression estimation based on SRS,especially,when the study variables and the ancillary variables are highly correlated (say |ρ|>0.85) [5].

In the last few decades,many applications and modifications of the RSS design have been proposed.Halls et al.[6]for considering an application of forage yields using RSS.Samawi et al.[7]introduced the Extreme RSS (ERSS) design.Al-Omari et al.[8]introduced ratio estimators of the population mean with missing values using RSS.Al-Omari [9]considered the median estimation based on double robust extreme RSS.Al-Saleh et al.[10]extended the work further and provided Multistage RSS (MSRSS) design.They proved that as the number of stages increases,the efficiency of the mean estimator under MSRSS increases and vice versa.A Robust L RSS procedure based on the idea of L estimators is suggested by Al-Naseer [11].Muttlak [12]introduced Median RSS (MRSS);he showed that it provides an unbiased estimator of the mean of symmetric distributions,and is more efficient than the SRS and RSS mean estimators.Jemain et al.[13]suggested multistage median RSS for estimating the population median and Jemain et al.[14]proposed some variations of RSS.Al-Omari [15]proposed ratio estimators of the population mean by considering ancillary information in SRS and median RSS and Al-Omari [16]considered the entropy estimation in RSS methods.Hossain et al.[17]suggested paired RSS for estimating the population mean.Shadid et al.[18]considered the BLUEs and BLIEs of the scale and location parameters together with the population mean using RSS.Al-Omari et al.[19]investigated the ratio estimation using a multi-stage median RSS approach.Al-Omari [20]proposed robust extreme RSS mean for mean estimation.Haq [21]proposed Shewhart control chart for monitoring process mean based on partially ordered judgment subset sampling.Haq et al.[22]suggested unbiased estimators for the basic linear regression model based on double RSS.Yu et al.[23]for investigating regression estimator in RSS.Al-Naseer et al.[24]proposed robust extreme RSS.Haq et al.[25]suggested some ratio estimators for the population mean in ERSS using two ancillary variables.Ozturk [26]studied sampling based on partially rankordered sets.Haq et al.[22]proposed the hybrid RSS method.Zamanzade et al.[27]introduced a new RSS estimation method for the population mean and variance.Haq [28]considered cluster sampling with hybrid RSS.Haq [29]studied the distribution function estimation under hybrid RSS.Al-Omari et al.[30,31]dealt with tests based on Laplace and logistic distributions.Al-Nasser et al.[32]studied information-theoretic weighted mean based on truncated RSS.Zamanzade et al.[33]used population proportion estimation in pair RSS.Haq [34]studied ordered partially subset sampling and consider the applications of this method to parametric inference.Al-Omari et al.[35]suggested a new RSS procedure called Truncation Based RSS (TBRSS),and showed that their estimator is unbaised of the population mean of symmetric distributions.Al-Nasser et al.[36]suggested minimax RSS method.Haq et al.[37]proposed the Hybrid ranked set sampling scheme.Wang et al.[38]investigated general ranked set sampling with cost consideration.Muttlak [39]introduced median ranked set sampling with concomitant variables and a comparison with ranked set sampling and regression estimators.For applications and new techniques based on RSS,we refer the readers to the references [40-44].

In this paper,we extended the work in this area and proposed two new improved RSS designs called the Simple Z Ranked Set Sampling (SZRSS) and the Generalized Z Ranked Set Sampling (GZRSS) methods.For some cases,SZRSS becomes a particular case of GZRSS design.The proposed sampling procedure estimator is unbiased of the population mean for symmetric distributions.It is shown,theoretically and numerically,that under perfect and imperfect rankings for symmetric non-uniform distributions,the proposed mean estimators under the GZRSS design are more efficient than those obtained by RSS and TBRSS.For asymmetric distributions,the proposed estimators based on GZRSS are more precise as compared to the estimators based on RSS and TBRSS.Also,we extended our sampling designs for estimating the population median.The efficiency of the suggested median estimators under GZRSS is better than that based on the RSS and TBRSS estimators,for symmetric non-uniform and asymmetric distributions.The GZRSS estimator of the population mean is investigated based on perfect and imperfect rankings,and is also compared to the SRS linear regression mean estimator.It is noteworthy that for small to moderate correlation between the auxiliary and study variables,the proposed estimators are more efficient than the SRS linear regression estimator of the population mean.

This paper is organized as follows.Some sampling methods are presented in Section 2.The proposed ZRSS designs are described in detail in Section 3.The problem of errors in ranking and a comparison with the SRS linear regression estimator is discussed in Section 4.The problem of estimating the population median is considered in Section 5.A detailed application to real data is given in Section 6,and finally,the paper is concluded in Section 7.

2 Sampling Methods

In this section,we explain some existing sampling schemes considered in this study.

2.1 Ranked Set Sampling

We describe the RSS design as follows:

Step 1:Given the value of sample size,saym,identifym2units from the corresponding population.

Step 2:These units are randomly allocated tomsets such that the size of each set ism.

Step 3:Now,rank the units within each set,this ranking can be done visually or by an inexpensive method with respect to the study variable.Then select the smallest ranked unit from the first set ofmunits.Similarly,select the second smallest ranked unit from the second set ofmunits.The procedure continues until the largest ranked unit is selected from the last set.This completes a cycle of a ranked set sample of sizem.

Step 4:For a large sample size,sayn,the above steps are repeatedrtimes until size of the sample becomesn=mr,forr≥1.

LetZbe the variable of interest with a distribution function (cdf)F(z)and a probability density function (pdf)f(z).Suppose thatZhas a meanµand a varianceσ2.LetZ1,Z2,...,Zmbe a SRS of sizemdrawn from the pdff(y),i.e.,Zi～f(z),fori=1,2,...,m.Then,the mean of SRS is denoted by(an unbiased estimator ofµ),and variance isSuppose thatZij,i,j=1,2,3,...,mbemindependent SRS each of sizem.LetZi(1:m),Zi(2:m),...,Zi(m:m)denotes the order statistics of theith sampleZi1,Zi2,...,Zim.Now,implement the RSS method tomselected samples.This gives a balanced RSS of sizem,Z1(1:m),Z2(2:m),...,Zm(m:m).The RSS mean estimator is denoted by=Assuming thatg(i:m)(z)be the pdf of theith order statisticZ(i:m),and noting that for eachi,Zi(i:m)Z(i:m),wheredstands for equality in distribution.The pdf of theZ(i:m)is given byg(i:m)(z)={1−F(z)}m−i f(z),−∞

2.2 Truncation Based RSS

As we mentioned in the introduction that Al-Omari et al.[35]derived the TBRSS design;its describtion is as follows:

Step 1:Choosemby SRS of sizemeach from the parent population.

Step 2:Within each choden sample,rank the units visually based on the variable of interest or by any inexpensive method.

Step 3:Define a coefficientδ=[αm],for 0≤α <0.5.Note that [t]denotes the integer part oft.

Step 4:Choose the minumum ranked unit from the firstδsamples and the maximum ranked unit from the lastδsamples.From the remainingm−2δ,choose theith ranked unit from theith sample fori=δ+1,...,m−δ.

Step 5:This finalizes a cycle of a TBRSS.Steps 1-4 are repeatedrtimes if needed to determine a sample of sizen=mr.

The corresponding estimator of population mean based on TBRSS is=The estimatorbecomes unbiased if the population is symmetric.For symmetric populations,the variance ofis=Note that forδ=0,1,For samples of odd sizes,whenδ=(m−1)/2,The TBRSS and ERSS become equivalent.For further details and application of this method,see Al-Omari et al.[35].

3 New Sampling Designs

This section introduces two new RSS methods;namely,simple Z ranked set sampling (SZRSS)and generalized Z ranked set sampling (GZRSS) designs.

3.1 Simple ZRSS Design

The SZRSS procedure for both even and odd samples is described as follows.To get an SZRSS ofmsize,selectmrandom samples each of sizem.Without yet knowing the values in the samples,rank the units within each sample based on any inexpensive or cost free method.

i) For evenm,choose the(i+1)th smallest ranked unit from the firstm/2 samples,fori=1,...,m/2.Similarly,choose(i−1)th the smallest ranked unit from the lastm/2 samples,fori=(m/2)+1,...,m.

ii) For odd sample sizem,choose the(i+1)th smallest ranked unit from the first(m−1)/2 samples,fori=1,...,(m−1)/2.Then select the median of the((m+1)/2)th sample.From the last(m−1)/2 samples,select the(i−1)th smallest ranked unit,fori=(m+3)/2,(m+5)/2,...,m.

This process provides a cycle of an SZRSS of sizemThe cycle are repeatedrtimes to determine the sizen=mr.The SZRSS estimator ofµfor an evenmis

Lemma 1:(i) For symmetric distributions,the estimator(J=E or O) of the population meanµis unbiased.(ii) Var()

Proof:

(i) For the estimator,given in Eq.(1),we haveFor any symmetric distribution,µ(i:m)−µ=µ−µ(m−i+1:m),fori=1,2,...,m.After some simplications,we can writeFollow the same process to prove that

(ii) The variance of Eq.(1) is defined as

For any symmetric distribution,i=1,2,...,m.With some algebraic operations,we can writeNote that the variance decreases asiincreases for symmetric (non-uniform) distributions,with minimum value occuring ati=[(m+1)/2],i.e.,≤fori,j=1,2,...,[(m+1)/2]withi≥j.Here,[t]represents the greatest integer value oft.Therefore,and hence,which completes the proof.Follow the same process to prove that

3.2 Generalized ZRSS Design

Now,we propose a generalized ZRSS (GZRSS) design.The steps of selecting a GZRSS are given below in which the Steps 1-3 are similar to the TBRSS method.

Step 1:Choosemsimple random samples,each with sizemselected from the corresponding population.

Step 2:Within each sample,rank the units visually with respect to the variable of interest or by inexpensive or cost free method.

Step 3:Define a coefficientδ=[αm],for 0≤α<0.5 and [t]symbolizes the integer value oft.

Step 4:From the firstδsamples,draw the(i+1)th smallest ranked unit.From the lastδsamples,draw the(i−1)th smallest ranked unit.But from the remainingm−2δsamples,draw theith ranked unit from theith sample fori=δ+1,...,m−δ.

Step 5:Previous steps finalize a cycle of a GZRSS of sizem.Steps 1-4 are donertimes if needed to determine a sample of sizen=mr.

LetZ11,Z12,...,Z1m,Z21,Z22,...,Z2m,...,Zi1,Zi2,...,Zim,...,Zm1,Zm2,...,Zmmbemindependent simple random samples each of sizem.The GZRSS estimator of population meanµbased on this sample is

and the corresponding variance isNote that forδ=0,we have

Lemma 2:For symmetric distributions about its population meanµwe have

Proof:

(i) Take the expectation of,given in Eq.(2),we haveAsµ(i:m)−µ=µ−µ(m−i+1:m),fori=1,2,...,m.Therefore,we can writewhich completes the proof.

(ii) Consider the estimator,given in Eq.(2),we have

For any symmetric distribution,i=1,2,...,m.After some simplication,we can writeAs explain above,for symmetric (non-uniform) distributions,σ2(i:m)≤σ2(j:m)fori,j=1,2,...,[(m+1)/2]withi≥j.

Therefore,σ2(1:m)>σ2(k+1:m)and hence,The equality is attained

whenk=0,which completes the proof.

As mentioned above,for symmetric (non-uniform) distributions,σ2(i:m)≤σ2(j:m)fori,j=1,2,...,[(m+1)/2]withi≥j.Therefore,all of the above differences are positive and henceThe equality is attained whenδ=0,which completes the proof.

In the case of symmetric of the parent distribution,theEffofwith respect tois defined byFor asymmetric populations,theEffwill beNow,to illustrate the method,some choices of the sample sizemand the coefficientδare considered for normal and Weibull distributions.

3.3 Examples

3.3.1 Normal Distribution

LetZ～N(0,1),where−∞

Case I:δ=0:The computed results for expectation and the variance of the estimator,given in Eq.(2),are,respectively,=0.0721.TheEffofwith respect tois given by

Case II:δ=1:The expectation ofis=0 and variance0.0613.TheEffofwith respect tois

Case III:δ=2:The expectation of the estimator given in Eq.(2) iswith varianceTheEffofwith respect tois3.3703.

3.3.2 Weibull Distribution

LetZ～Weibull(2,1),having pdff(z)=z >0.The cdf and pdf of theith ranked unit from a ranked set sample form=5,respectively,areG(i:5)(z)=BetaRegularizedandg(i:5)(z)=

Form=5,we haveBased on the order statistics,the means and variances areµ(1:5)=0.3963,µ(2:5)=0.6302,µ(3:5)=0.8479,µ(4:5)=1.0946,µ(5:5)=1.4619,andσ2(1:5)=0.0429,σ2(2:5)=0.0528,σ2(3:5)=0.0643σ2(4:5)=0.0850,σ2(5:5)=0.1459,respectively.Also,the cases below forδcan be treated as follows:

Case I:Consider δ=0:

Case II:Consider δ=1:Similarly,the mean (expectation) and variance ofare=0.8595 and=0.0135,respectively.As the estimator is not unbiased,therefore,the MSE ofis given byThe bias ofis=−0.0266.Therefore,=0.0143.TheEffofwith respect tois

Case III:Considerδ=2:The mean and variance ofare=0.8537 and=0.0132.Again the estimator is biased with=−0.0324 andTheEffofwith respect tois

Table 1:Exact Eff of mean estimators under symmetric distributions

Now,we consider the mean estimation for some symmetric distributions,and also for some asymmetric distributions.The exact relative efficiencies of our proposed estimators are presented in Tabs.1 and 2.

Tabs.1 and 2 show that,for symmetric distributions,the efficiency of the GZRSS increases as theδvalue increases except in the case of the uniform distribution.In the case of asymmetric distributions,generally,the efficiencies increase whenδincreases for 0 to 1,and they decrease function whenδ>1.For both asymmetric and symmetric distributions,the relative efficiency of mean estimators under GZRSS is an increasing function of the sample size.For all considered cases,GZRSS is more efficient than RSS and TBRSS except that TBRSS is more adequate than GZRSS when the considered distribution is standard uniform.

4 Errors in Ranking and Comparison with SRS Regression Estimator

We investigate the fulfillment of the suggested estimators for the mean under both GZRSS design and imperfect rankings.The suggested estimators under both rankings’ schemes are also compared with the SRS for the population mean based on the linear regression estimator.

4.1 Errors in Ranking

Accurate ranking increases the efficiency of the RSS.However,Dell et al.[3]show that even if the ranking has some errors,the estimator under RSS still remains unbiased and performs at least as well as the SRS estimator.Here,we study the performance of the estimators under the proposed RSS designs,when ranking has some errors.The mostly used RSS model to study the effect of errors in ranking is based on the ranking with respect to a concomitant variable that is correlated with the study variable.The efficiency of the estimator now depends on the correlation value between the study variableZand the concomitant or ancillary variableW.Stokes [4]suggested a model for imperfect ranking assuming that an ancillary variableWis available,can be simply measured and is correlated with the interest variableZ.For further details see Stokes [4],Patil et al.[5]and Muttlak [39].Stokes [4]imposed the following assumptions considered in developing the following model:

(i) The relationship betweenZand the regressorWis linear,

Table 2:Exact Eff comparison of mean estimators under asymmetric distributions

If(Z,W)follows the bivariate normal distribution,then both conditions are easily satisfied.Following Stokes [4]and (i),we can writeZ[i:m]=µZ+i=1,2,...,m,where whereρis the coefficient of correlation,σZandσWare the population standard deviations,µZandµWare the corresponding means.Note that the ranking of the auxiliary variableWis perfect whereas the ranking ofZis imperfect,i.e.,the ranking ofZhas some errors.Here,W(i:m)andZ[i:m]denote theith order statistic and theith judgment order statistic of a random sample of sizem.ξidenotes the error term with zero mean and a constant variance,i.e.,E(ξi)=0 and Var(ξi)=As SZRSS becomes a special case of GZRSS,therefore,we consider the estimator based on GZRSS.Now,the mean of the study variableZwith ranking based on the auxiliary variableWunder GZRSS can be written aswhereis unbiased estimaor ofµY,and it variance isNote that if we considerδ=0 in,then it becomes the simple RSS estimator of population mean.The efficiency ofwith respect tois

The exact efficiencies of our proposed estimators under GZRSS with respect to RSS and TBRSS are given in Tab.3.

It is clear from the results given in Tab.3 that,as the efficiencies under each design are a function of the correlation coefficientρ,i.e.,as the value ofρincreases,the relative efficiencies increase and vice versa.As expected,the increase in the sample size also increases the efficiency of the estimator under each of the RSS design.The proposed estimators are better than the existing counterparts.

4.2 Comparison with Regression Estimator Based on SRS

Patil et al.[5]compared the estimator of a population mean under RSS with the regression estimator based on SRS.It is shown that for a small correlation between the study variable and the ancillary variable,the RSS mean estimator is better than the regression estimator under SRS.In this section,we compare the performance of the proposed mean estimator under GZRSS with respect to the SRS regression estimator.It is assumed that the population mean of the ancillary variable is known.Following Muttlak [39],the linear regression ofZonWisZi=α+βWi+ξi,i=1,2,...,m,whereαandβare the intercept and slope of the regression line.Here,ξiis error term with zero mean.The linear regression estimator of the population meanµZwhenµWis known is

whereandare the corresponding sample mean ofZandW,based on an SRS of sizem.Note that,is the least square estimator of the slopeβof the regression line.Sukhatme and Sukhatme (1970) showed that the regression estimator given in Eq.(3) of the population meanµZis an unbiased estimator once the joint distribution ofZandWis a bivariate normal distribution.The variance ofis given byIn case of perfect ranking,theEffofrelative tois given bySimilarly,in case of imperfect ranking,theEffofrelative tois given by

In Tab.4,we provide exact relative efficiencies of the proposed estimators with respect to the classical linear regression estimator of mean.Note that the proposed mean estimator with perfect ranking under GZRSS outperforms other competitor estimators when the value ofρis less than 0.9.

Similarly,in Tab.5,we compared the performance of the suggested estimators under imperfect ranking with respect to the linear regression estimator.It is worth mentioning that even when there are errors in ranking,the proposed estimator is still more efficient than the linear regression estimator when the value ofρis less than 0.8.The efficiencies of the newly estimators are high based on perfect ranking as compared with the case of imperfect ranking.Note that here RE is a decreasing function of sample size because the performance of linear regression estimator is increasing with the increasing of the sample size.For all of the cases,GZRSS mean estimator always performs better than the TBRSS estimator.

Table 4:The Eff of the SRS linear regression estimator with respect to the GZRSS estimator based on perfect ranking

Table 5:The Eff of the SRS linear regression estimator with respect to the GZRSS estimator based on imperfect ranking

5 Estimation of Population Median

Estimation of the population median based on the sampling methods,studied in this paper,is presented in this section.LetQbe the population median andZ1,Z2,...,Zmbe an SRS of sizem.Then,the median estimator is given byfor oddm,andfor evenm.From RSS units of sizem,i.e.,Z1(1:m),Z2(2:m),...,Zm(m:m),the population median estimator based on the RSS isThe corresponding population median estimator based on the GZRSS is=medianThe efficiencies ofandwith respect to,are given bywhereh=GZRSS,RSS.The estimated MSE of any median estimator is defined ash=GZRSS,RSS,SRS.The median estimation of some symmetric and asymmetric distributions is considered here based on extensive Monte Carlo simulations.The obtained results are presented in Tabs.6-8.

Table 6:Eff comparison of median estimators under symmetric distributions

Table 7:Eff comparison of median estimators under asymmetric distributions

The results,given in Tabs.6 and 7,reveal that the attainment in efficiency determined by using the GZRSS method.For instance,whenm=7 andδ=2,the RE of the GZRSS is 3.5187 for estimating the median of the student’stdistribution.Also,GZRSS is more efficient than RSS and TBRSS based on the same sample size for a fixed value ofδ.To study the performance of the proposed median estimators under GZRSS for imperfect rankings,we have considered standard bivariate normal distribution.The relative efficiencies of the median estimators are obtained for different values of correlation coefficient using extensive Monte Carlo simulations and are displayed in Tab.8.

According to the results given in Tab.8,the median estimators under proposed designs are at least as efficient as compared with the SRS median estimator.Here,the relative efficiencies are also increasing function ofmandρ.The results under GZRSS are efficient as compared to RSS and TBRSS under perfect and imperfect rankings.

Table 8:Eff comparison of median estimators under bivariate normal distribution

6 An Application to Real Data

To illustrate the use of the GZRSS method in the field,a real data set is considered for both mean and median estimation.This real data set is considered by Platt et al.[45]and it is related to the height and diameter of 399 conifers (Pinus Palustris) trees.The data consists of 7 variables of which we have considered only 2 variables.Let the variable of interestZrepresents the height of the conifer tree measured in feet while the ancillary variableWis the diameter of the tree at breast height.In Tab.9,we provide the summary statistics of the data,and the corresponding plots of the data are displayed in Fig.1.

Table 9:Statistics summary of the trees data

Figure 1:List plot (left) and histogram (right) of the 399-tree data

Table 10:The Eff of estimating of the population mean and median of the study variable based on perfect and imperfect rankings

For the diameter and the height,the coefficients of skewness are 0.884 and 1.619 respectively,indicating that these data are non-symmetric.The MSEs for various estimators (under SRS,RSS,TBRSS and GZRSS methods) were calculated by one million iterations.The obtained results are summarized in Tab.10.The samples were drawn using SRS without replacement.The results given in Tab.10 are the mean and median estimation values of the trees’ heights under perfect and imperfect rankings.These results demonstrate that the GZRSS estimators are more efficient than their competitors.As we concluded in the above sections,the RE increases as sample size increases and vice versa.The perfect ranking provides efficient estimates than imperfect ranking.Also,the relative efficiencies under GZRSS in median estimation are greater than mean estimation because the data is asymmetrically distributed.The GZRSS is recommended for estimating the mean and median of the trees data.

7 Conclusions

We propose two new efficient RSS sampling methods for estimating the population mean and median.The proposed estimators based on the new designs are compared with their competitors using SRS,RSS and TBRSS techniques based on the same number of quantified units.It turns out that the GZRSS estimators of the population mean for symmetric populations are unbiased.It is worth mentioning that for non-uniform symmetric distributions,under perfect and imperfect rankings,the mean estimators under the proposed GZRSS are more efficient than those under SRS,RSS and TBRSS methods.We also compare the performance of the mean estimator under GZRSS with the SRS linear regression estimators.It is observed that for small and moderate correlation between the study and ancillary variables,the suggested estimators are more efficient than the SRS linear regression estimator for perfect and imperfect rankings.Therefore,we recommend the use of the proposed sampling methods over the existing RSS methods,considered here.The proposed methods,in this paper,can be considered in many real applications,such as mean estimation in case of missing data [46],quality control charts for monitoring the process mean [47],and in acceptance sampling plans [48,49].

Funding Statement:The authors extend their appreciation to Deanship of Scientific Research at King Khalid University for funding this work through Research Groups Program under Grant No.R.G.P.2/68/41.I.M.A.and A.I.A.received the grant.

Conflicts of Interest:The authors declare that they have no conflicts of interest to report regarding the present study.