Mathematical Problems in Engineering

Mathematical Problems in Engineering / 2021 / Article
Special Issue

Robust Estimation Methods in the Presence of Extreme Observations

View this Special Issue

Research Article | Open Access

Volume 2021 |Article ID 1953546 | https://doi.org/10.1155/2021/1953546

Sehar Saleem, Rehan Ahmad Khan Sherwani, Muhammad Amin, Maryam Khalid, Nouman Ali, "Development of New Robust Optimal Score Function for the Weibull Distributed Error Term in Multilevel Models", Mathematical Problems in Engineering, vol. 2021, Article ID 1953546, 15 pages, 2021. https://doi.org/10.1155/2021/1953546

Development of New Robust Optimal Score Function for the Weibull Distributed Error Term in Multilevel Models

Academic Editor: Ishfaq Ahmad
Received07 Apr 2021
Revised08 May 2021
Accepted15 May 2021
Published25 May 2021

Abstract

A popular robust estimation technique for linear models is the rank-based method as an alternative to the ordinary least square (OLS) and restricted maximum likelihood (REML) in the presence of extreme observations. This method is applied in machine reliability analysis and quantum engineering, especially in artificial intelligence and optimization problems where outliers are commonly observed. This technique is also extended for the multilevel model, where the shape of error distribution contributes a significant role in more efficient estimation. In this study, we proposed the Weibull score function for the Weibull distributed error terms in the multilevel model. The efficiency of the proposed score function is compared with the existing Wilcoxon score function and the traditional method REML via Monte Carlo simulations after adding simulated extreme observations. For small values of shape parameter in Weibull distribution of error term showing the presence of outliers, the Weibull score function was found to be efficient as compared to the Wilcoxon and REML methods. However, for a large value of shape parameter, Wilcoxon score appeared either equally efficient than the Weibull score function. REML is observed least precise in all situations. These findings are verified through a real application on test scores data, with a small value of shape parameter, and the Weibull score function turned out the most efficient.

1. Introduction

The rank-based method is an alternative robust estimation method over conventional estimation techniques of the linear models such as restricted maximum likelihood (REML) and ordinary least square (OLS) when errors are asymmetric due to the presence of outliers. Rank-based fitting works by assigning ranks to the residuals based on a pseudonorm similar to the Euclidean squared norm used in the OLS providing for a robust fit [1]. The rank estimation method is also a popular parameter estimation method in reliability analysis, machine learning, and artificial intelligence, where the structure of error term is positively skewed. The reason for skewness may be due to the presence of one or more extreme observations. For example, if one wants to study the life (reliability) batteries, then some batteries maybe failed at an initial time and some are failed after a long time. Both these times are extreme observations. These observations are sometimes also called outliers. So, in this situation, we need a ranked estimation method as a robust estimation method alternative to REML and OLS.

Garcia et al. [2] developed a new rank-based constraint handling technique for the solution of engineering optimization problems. The comparison with other techniques revealed the accuracy and robustness of the rank-based technique. Neelima et al. [3] presented a rank-based skeleton extraction algorithm that shows improved performance based on speed, time, rate, and connectivity than various methods in the literature. Li et al. [4] presented a new powerful rank-based procedure for the EWMA control scheme for online sequential monitoring. This procedure can be used to monitor the location and scale parameters of a continuous distribution. From the simulation study, the rank-based method appeared quite robust to nonnormally distributed data and efficient in detecting various process shifts. Chen et al. [5] advocated that rank regression is a more objective approach for dealing with nonnormal data due to the presence of outliers. This finding is evidenced by comparing the rank-based method with traditional parametric and semiparametric regression models through Monte Carlo simulation and real data examples. Liu et al. [6] proposed estimators of nonparametric coefficients in robust estimation for the varying coefficient partially functional linear regression model (VCPFLM) based on B-spline approximation and rank regression. The proposed estimators are found robust for nonnormal error distribution compared with OLS through a simulation study and real-life dataset. Henry et al. [7] compared rank-based estimation and OLS for the Granger and Lee asymmetric price transmission model. Both estimation methods appeared equivalent for normal data irrespective of the sample size. On the other side, for the data affected by outliers, OLS estimates turned inaccurate. However, the rank-based estimation is robust to outliers in large sample size and produced unbiased estimates. The results of the simulation study confirmed the outperformance of rank-based estimation for asymmetric data with outliers and almost equivalent performance for normal datasets in the context of asymmetric price transmission modeling. The proper choice of a score function leads to asymptotically efficient estimators. The efficiency of robust rank-based analyses is compared with the traditional REML method for normal and nonnormal error distributions. For nonnormal error term in linear models, rank-based estimators attain more efficiency than REML but produce quite a similar efficiency for normal error term [1]. Fixed effect estimates, standard errors, and estimates of variance components are compared through a large simulation study. Rank-based analyses appeared much more powerful and efficient than REML. McKean and Hettmansperger presented a review of the development of rank-based theory for linear models and various extensions [8]. Kloke et al. [9] extended this rank theory for the linear models with cluster-correlated errors for example mixed models and nonlinear models. A complete unified methodology of rank-based analyses is presented in [10]. Usually, rank-based analyses are a robust and efficient estimation method. Nevertheless, this efficiency can be optimized based on the information available regarding the distribution of the error term. The rank-based analysis turns out to be more accurate and efficient when the selected score function is close to the form of ‘f’. If the probability density function (pdf) of u is f (u), then the expression for optimal score function can be obtained through some mathematical derivation [8]. Thus, the prior knowledge of the distribution of error term would assist in proper choice of score function leading towards efficient regression estimators. The underlying shapes of error distributions can be symmetric, right-tailed, left-tailed, light-, moderate-, and heavy-tailed. One common probability distribution for right-skewed data is the Weibull distribution, considered the case of the right-tailed error distribution. The shape of the distribution is right-skewed and controlled by two parameters such as scale () and shape . The level of skewness of the Weibull distribution decreases as the value of the shape parameter increases. Weibull distribution mostly occurs in the reliability of products based on design, manufacturing, and development in prelaunch and postlaunch stages [11]. The rank-based estimation method is highly efficient than OLS and REML in the case of the nonnormal error term. The dependent error structure in the error term spoils the normality. The efficiency of the rank-based method can be optimized by selecting the appropriate score function for Weibull distributed error term. In peridynamics, the variability in the brittle fracture can be modeled well by Weibull distribution [12]. Weibull score function would provide the most efficient estimates of standard error. Some research work related to the rank-based estimation method is available in the literature. Mckean and Kloke [13] developed a family of optimal score functions for a skew-normal family of distributions. Additionally, the validity and efficiency of rank-based estimators are compared with MLE and Wilcoxon score through simulation study for skew-normal and contaminated normal distribution. Mckean and Sievers [14] proposed rank score functions suitable for skewed error distributions, specifically for generalized F family of error distributions. These score functions are bounded and provide a robust and powerful rank analysis for a linear model with application to model the lifetime data. Zhen et al. [15] proposed a new score function of intuitionistic fuzzy values (IFVs). The new score function is capable to overcome the problems of other score functions for ranking IFVs. Lalancette et al. [16] provided a review of rank-based estimation methods for asymptotic dependence and independence in linear models. Rank-based M-estimators are proposed, and asymptotic normality is established under weak regularity conditions. Watcharotone et al. [17] developed a robust rank-based picked-point analysis optimizable for heavy-tailed or skewed distributions for the analysis of covariance with heteroscedastic slopes. This technique is compared with least-square based picked-point analyses for normal and nonnormal models by Monte Carlo simulation. Rank-based analyses appeared valid and more powerful than the least-square method over all the situations, although losing little efficiency for normal models. Kloke and McKean [18] developed R package for computing rank-based estimators and their associated inferences. This package includes a library of several score functions suitable for different shapes of error distributions and provides detailed rank analyses. Shao et al. [19] provided ordered treatment methods based on maximum likelihood and robust estimation for one-factor randomized group design including a vector of the covariate for the cluster-correlated model. The proposed estimation method produced higher power and showed robustness against outliers through theoretical and simulation results. Bindele et al. [20] proposed the estimator for rank estimation of regression coefficients in a single-index regression model. The conditions are established for consistency and asymptotic normality of the estimators. The efficiency and robustness of the rank estimator are compared with semiparametric least-square by using Monte Carlo simulation. A real-life example is used as an illustration that rank regression fixes model nonlinearity in the presence of outliers. Abebe and McKean [21] discussed the estimation of parameters in the linear regression model. They considered the estimator which minimizes weighted Wilcoxon dispersion function and established its asymptotic properties. The robustness, efficiency, and validity of these estimates are verified over several normal and nonnormal error distributions. Terpstra and McKean [22] demonstrated robust analyses based on the calculation of weighted Wilcoxon (WW) for linear models. An R suite of function is developed for WW estimation and testing of hypotheses, diagnostic measures, and residual analysis is available. Bindele [23] considered rank estimator of the general linear regression model. The author established strong consistency of the rank estimator is under mild conditions. Abebe et al. [24] presented a rank-based fitting procedure for repeated measurement design that is based on replacing a norm based on a score function for the Euclidean norm. By using a simulation study, proposed estimators are proved efficient, asymptotically normal, and valid, also a real dataset with inherent hierarchy is used for illustration. Černý et al. [25] studied algorithms for minimization of Jaeckel’s dispersion function to reach the robust rank estimator that is insensitive to outliers. A new two-stage algorithm is developed merging the benefits of two algorithms already used in the literature, approximate and exact algorithms. The behavior of this two-stage algorithm is illustrated using a computational experiment, and convergence and exact result is achieved. Cerny et al. [26] focused on the optimization algorithms of rank estimators by dividing these in two major classes: continuous and convex objective function (CCC), and another class (GEN). For CCC unconditionally polynomial algorithm and for GEN, the enumerative algorithm is proposed that is efficient and superior to other algorithms in the literature. Dutta and Datta [27] developed rank-based weighted estimating equations appropriate for intracluster group size is informative. Additionally, an aligned rank-sum test based on covariate adjusted outcomes is constructed. The authors provided asymptotic distributions, and test-statistics of rank estimators are presented. If the informativeness is present, the significance of selecting proper weights is shown through simulation studies. The superiority and robustness of proposed method is shown in comparison to traditional mixed models in clustered data by using real-life datasets. Dutta and Datta [28] presented a novel extension of the rank-sum test for the scenario where group-specific marginal distributions are based on intracluster group size. The performance of proposed test is compared with the Wilcoxon rank-sum test and classical signed-rank test by using the simulation study of informative intracluster group size. It is observed that the proposed test maintained correct size and attained more power [27, 28]. Xie et al. [29] developed a rank test based on the rank score function using functional principal component analysis. Also, asymptotic properties of the test are established under local and alternative hypotheses. The proposed test attained good size and power in the results of simulation study.

The literature showed that no one still worked on the multilevel modeling when the error term follows a Weibull distribution. So, in this article, we proposed a score function for the multilevel model when the error term follows a Weibull distribution to optimize rank-based analyses. Additionally, we compared the performance of the proposed score function with the Wilcoxon score and the traditional REML method.

2. Description of Multilevel Model and Rank Theory

In this section, a brief overview of the random intercept multilevel model, rank theory, and optimal score function is described.

2.1. Random Intercept Model

The multilevel models generate cluster-correlated errors due to the hierarchical structure in data and are a special case of mixed effect models. A general two-level random intercept multilevel model is defined as [30]where is an vector of responses measured at level-1 for the individual I in cluster j (i subscript refers to individual-level variation and j refers to group-level variation), is the design matrix of fixed effects, and is the design matrix of random effects. is the vector of fixed effect, is the vector of individual-level error terms for group j, and is the vector of random effects for cluster j and following normal distribution, i.e., and , respectively. It is assumed that the error terms follow to be independently normally distributed across individuals with distribution function F and density function f.

2.2. Brief Overview of Rank Theory

In the least-square estimation, the vector is obtained by minimizing the Euclidean norm; similarly, rank-based estimation is accomplished by minimizing a pseudonorm. Thus, the robust estimate of can be defined as [9]

It is efficient and robust in Y-space (Kloke and McKean, 2012). Al-Shomrani (2003) defined this pseudonorm bywhere R denotes the rank of among and these are invariant to constant shifts, is the dependent variable for the individual i in a cluster j, and is the corresponding matrix of covariates.

A suitable selection of score function according to the shape of error distribution is required for optimal analysis. The score function is a nondecreasing square-integrable function bounded in the interval (0, 1). The score function can be standardized such that and . Scores are calculated as , where t = 1,…, N, and N is the total sample size, and these scores sum equal to zero, i.e., .

The asymptotic distribution of the rank-based estimator is defined as [13]

Hence, minimizing the variance of the is equivalent to minimizing . The parameter is given by [13]

The comparison based on asymptotic relative efficiency can be obtained as , where is the variance of the REML estimator that is a traditional method to estimate multilevel models based on the maximization of the likelihood function.

2.3. Optimal Scores

The rank-based analyses are based on the selection of a score function. If the underlying distribution of the error term is known, then [31] showed that the following expression can be used to derive the optimal score function as

This expression calculates scores that lead to fully efficient rank-based estimates where and are probability density function (pdf) and cumulative density function (cdf), respectively, of the respective error distributions. Wilcoxon score function is the optimal score for the logistic distribution of the error term. The expression for the Wilcoxon score is .

According to Mckean and Kloke [13], from equations (5) and (6), can be explained as

In the last expression, is the correlation coefficient and is the Fisher Information matrix. Hence, minimization of scale estimator can be achieved by maximizing by taking  = 1, i.e., . Thus, (6) is a compact expression of the optimal score function. Only the knowledge about the form of f (x) is required to make the rank-based estimator fully efficient [13].

3. Weibull Score Function

In this section, after presenting a brief description of the Weibull distribution, a path diagram to reach the optimal score function for the Weibull distribution is shown. This segment includes all the important steps in completing the derivation of the score function.

3.1. Proposed Score Function

The path to propose optimal score function is briefly described in a flow diagram in Figure 1.

3.2. Weibull Error Distribution

The Weibull distribution is a continuous right-skewed distribution. It was first identified by Fréchet in 1928 [32] and later on discussed in detail by Swedish mathematician Waloddi Weibull [33]. This distribution is extensively used in survival analysis, reliability analysis, and extreme value theory. The pdf of Weibull distribution is of the form:where and skewness is controlled by shape parameter typically ranges between 0.5 and 8.0. In this study, we have fixed the value of the scale parameter and varied the value of the shape parameter from 1.0 to 4.0, capable to change the distribution shape from skewed to symmetric. Hence, the pdf simplifies to the following form:

In this article, we are concerned with the random intercept multilevel model (1), in which the level-1 and level-2 random errors follow the Weibull distribution. To develop an optimal score function for some specific value of shape parameter , it requires only knowledge about the form of pdf.

3.3. Derivation of Score Function

The pdf of Weibull distribution for the error term is defined bywhere is the shape parameter and is the scale parameter. Let us assume , then the pdf becomes .

The quantile function of Weibull distribution for is given by

The optimal score function can be derived by the following expression:

Let us assume that for convenience in derivation, where ,

After substituting the value of x,

The score function is

The derivative of the score function is required for implementation in the R package:

The plot of the score function is suitable for Weibull distribution for the shape parameter and is shown in Figure 2. The left panel comprises pdf’s and the right panel displays its corresponding optimal score functions. It is noticeable that the pdf is extremely right-skewed for the smallest value . The skewness in distribution tends to become symmetric with increasing the value of . Hence, when the value of , then its pdf is closer to symmetric distribution.

4. Simulation

The simulated model is as follows:where is the response variable and are the explanatory variables generated from uniform distribution, i.e., measured at level-1. and are random errors at the group level and individual level, respectively, simulated with . The number of groups is denoted by j, and the vector has observations, where is the number of observations. The score function is evaluated for multilevel models at different sample sizes. Level-1 sample size (N1) was taken 5, 10, 20, and 40 within each level-2 sample size (N2) 5, 10, 20, and 40 leading to the range of total sample size N = 25-1600. The simulation is repeated 1000 times on R software. The number of covariates included in the model is considered to be and . The true regression parametric values are selected according to [13]. The efficiency of the Weibull score function is compared with the Wilcoxon score function and the traditional method REML. To assess the performance of the Weibull score function, we use bias, MSE, and precision as the performance evaluation criteria. Absolute estimation of bias for the fixed effects of covariates and the random effect is calculated as , where is the vector of estimates. The estimates of MSE is calculated as . Three estimates of precision are calculated by taking a ratio of scale estimators obtained by a rank-based method with Weibull score, Wilcoxon score, and REML as .

Two kinds of error distributions are reflected in this simulation study: Weibull distribution and contaminated Weibull distribution. The contaminated Weibull error model is generated aswhere is a random variable with a binomial distribution (1,  = 0.20). The is the Weibull random error with the shape parameter . is a random variable from Weibull distribution . A similar setting is used for group-level error too. For robust analysis through rank theory, the score is selected suitable for in Weibull distribution. Five values of close to are considered to show a comparative analysis of optimal score function, i.e., .

5. Results and Discussion

Table 1 and Table 2 summarize the results of the simulation study with estimated bias and MSE of . The estimates of bias and MSE of the multilevel model including =3 and =6 covariates, when , are computed through the rank-based method, Wilcoxon score function, and REML. After obtaining three estimates of precisions, it is worthy to notice that the rank-based method is most efficient than REML, either applied through Weibull or Wilcoxon score function. Auda et al. [1] debated the same results that the rank-based method is generally highly efficient even if few outliers are present in the data. The efficiency of regression estimates is many times better than MLE for skewed data. Therefore, low MSE is observed for the rank-based method either applied through Wilcoxon or Weibull score function than MLE.


REMLWeibullWilcoxonPrecision
N2N1BiasMSEBiasMSEBiasMSEWeib/REMLWeib/WilWil/REML

552.4164.8672.0933.8662.2074.6220.7940.8360.95
102.4065.8572.0594.3552.1565.0640.7440.860.865
202.4229.2412.036.8842.1297.780.7450.8850.842
402.3715.4181.97811.3182.07112.6170.7340.8970.818

1052.5413.3162.1792.5192.2382.8580.760.8820.862
102.5173.9532.1633.0182.2213.3180.7630.910.839
202.5525.8172.1994.4872.2524.750.7710.9450.817
402.579.882.2097.6722.2588.0380.7770.9540.814

2052.6232.6812.2862.0752.3212.2360.7740.9280.834
102.5652.8992.2112.2742.2532.3660.7850.9610.816
202.5973.9752.2583.132.2833.2330.7870.9680.813
402.6015.8472.2614.5782.2834.6990.7830.9740.804

4052.6072.2632.2671.7692.2911.8390.7820.9620.813
102.6872.652.3582.0962.3722.150.7910.9750.811
202.5942.9062.2652.2872.282.3270.7870.9830.801
402.6444.0342.3133.2112.333.2410.7960.9910.804


REMLWeibullWilcoxonPrecision
N2N1BiasMSEBiasMSEBiasMSEWeib/REMLWeib/WilWil/REML

552.7419.382.4783.4252.3524.1840.8870.8571.035
102.8549.6082.6324.7532.4174.6950.8120.8940.908
202.72510.4952.5836.3152.3616.1730.8040.9460.85
402.81416.9522.64512.2142.43211.9020.8210.9810.837

1053.0124.8412.8012.7162.5532.6430.8520.9190.927
102.9374.9772.7823.0562.5232.9160.8290.9590.864
202.9587.6562.8416.012.5765.8130.8311.0110.822
402.9899.7882.8057.642.5447.3190.8431.0320.817

2052.9463.232.8511.9572.5741.8420.8330.9610.867
10332.7992.2032.5242.0870.8341.0080.828
203.0634.1072.8913.1482.6082.9940.8391.0310.813
402.9955.1782.8424.1822.5513.9340.8571.060.808

4053.0522.1992.8691.5182.5941.4270.8311.0010.83
103.1173.0222.9752.442.6782.3480.8451.0310.82
203.0332.8712.872.242.572.1290.8581.0530.815
403.0483.8642.8743.1422.5742.9870.8621.0720.805

Furthermore, the smallest bias and MSE are achieved through the Weibull score function as compared to the Wilcoxon, for with =3 case. Although for  = 6 covariates, for group size 5, Weibull score is better than Wilcoxon but for group size larger than 10, Weibull and Wilcoxon are equally precise.

Table 3 and Table 4 comprise the estimates of bias and MSE of for =3 and =6 covariates in the multilevel model when . It is observed that in all the simulation scenarios, Weibull score function is more efficient than the Wilcoxon score. Additionally, it is seen that REML is least efficient than rank-based estimation with Weibull and Wilcoxon score function. Table 5 and Table 6 display the estimates of bias and MSE of for =3 and =6 covariates in the multilevel model . It is observed that for =3 against all the sample sizes, the Weibull score function appeared more efficient than the Wilcoxon score. Additionally, it is seen that REML is least efficient than rank-based estimation with Weibull and Wilcoxon score function. However, for =6, both the score functions produce quite a similar fit and all three techniques are equally precise. Table 7 and Table 8 show the estimated bias and MSE of for =3 and =6 covariates in the multilevel model . It is observed that for =3 against all the sample sizes, the Weibull score function performed efficiently than the Wilcoxon score function. Additionally, it is seen that REML is equally efficient as rank-based estimation with Weibull and Wilcoxon score function. However, for =6, both the score functions produce quite a similar fit and all three techniques are almost equally precise. As the value of the shape parameter increases the distribution shape becomes symmetric, hence the Wilcoxon score function starts becoming the optimal score.


REMLWeibullWilcoxonPrecision
N2N1BiasMSEBiasMSEBiasMSEWeib/REMLWeib/WilWil/REML

552.1482.9531.9952.852.0593.1110.9651.0540.916
102.1314.0541.9773.7262.0434.0140.9190.990.928
202.1566.6561.9985.9372.0636.3790.8920.9580.931
402.13711.591.97810.0362.04110.850.8660.9360.925

1052.232.1622.0972.0032.1362.120.9270.9810.945
102.2542.852.1132.5792.1472.690.9050.9440.959
202.2644.1812.1243.7582.1613.890.8990.930.966
402.277.0272.1276.2682.1636.4830.8920.9220.967

2052.3131.8242.181.6732.2071.7250.9170.9460.97
102.3152.2112.1832.0142.2032.0540.9110.9290.98
202.3092.8922.1812.6122.1992.6540.9030.9180.984
402.3044.3022.1763.8732.1933.9390.90.9160.983

4052.3161.6422.2021.4882.2091.5230.9070.9270.978
102.3251.8352.2011.6782.2131.6930.9140.9220.991
202.2892.12.1751.8962.1811.9160.9030.9130.989
402.3252.9072.2062.6322.2142.6570.9050.9140.991


REMLWeibullWilcoxonPrecision
N2N1BiasMSEBiasMSEBiasMSEWeib/REMLWeib/WilWil/REML

552.2072.2432.3422.5352.1742.5640.9890.9891.143
102.1522.6482.2652.7492.0822.7550.9980.9981.04
202.1633.9752.2844.0552.0633.941.0291.0290.991
402.2247.2582.3687.1872.1336.8551.0481.0480.944

1052.3391.4572.5051.5132.2651.4871.0171.0171.021
102.3251.7922.491.8292.221.7511.0451.0450.977
202.3142.552.4632.5892.2162.4021.0781.0780.942
402.3074.1452.4564.182.2023.8491.0861.0860.929

2052.321.1232.4981.1442.2161.0911.0491.0490.972
102.3471.2862.5181.3012.2411.211.0751.0750.941
202.3321.6932.5121.6982.2271.5691.0821.0820.927
402.3762.6012.5432.6222.2662.3971.0941.0940.922

4052.3871.0092.56712.2770.9481.0551.0550.94
102.3961.1262.5871.1192.2921.041.0761.0760.924
202.3591.2732.5291.282.2491.171.0941.0940.919
402.3721.7282.5511.7342.2621.5841.0951.0950.917


REMLWeibullWilcoxonPrecision
N2N1BiasMSEBiasMSEBiasMSEWeib/REMLWeib/WilWil/REML

552.0792.4570.2250.5272.0832.7171.0720.971.106
102.093.5370.2240.6832.1073.6910.9940.9531.044
202.0985.7960.241.1532.0966.0050.9810.9471.036
402.09610.2960.2371.9512.10410.530.9640.9431.023

1052.1991.9080.0780.2842.2051.9671.020.9891.031
102.1932.5510.0940.3692.1852.61.0010.9821.019
202.2133.8240.0680.5122.2113.8390.9850.9811.004
402.2146.3490.1120.8572.2086.3760.9840.981.004

2052.2761.6740.0330.2312.2631.6760.9970.9961.001
102.2862.0150.0480.262.2742.0130.9920.9930.999
202.2562.6280.0520.3442.2452.6050.9870.9950.991
402.2653.9640.0460.4632.2523.9340.9840.9920.992

4052.2881.5110.0340.2032.2741.4970.9870.9960.991
102.2811.6850.0380.1932.2681.6580.980.9960.984
202.2942.0280.0180.2122.2742.0010.9820.9950.987
402.3022.7370.0110.3052.2842.7030.9840.9960.988


REMLWeibullWilcoxonPrecision
N2N1BiasMSEBiasMSEBiasMSEWeib/REMLWeib/WilWil/REML

552.1451.6160.2230.562.1471.9551.0921.0681.21
102.1812.2190.2140.7212.1762.4291.0891.0861.095
202.153.480.2361.1472.1533.6531.0391.0851.05
402.1396.2080.2362.0842.1436.4131.0121.0761.033

1052.2641.2120.0730.2962.2661.2941.0651.0911.068
102.2781.5610.1030.3422.2681.6141.0231.0861.034
202.252.2730.1440.5472.2452.3051.0191.0141.014
402.2593.790.0870.852.2613.8111.0171.091.006

2052.3371.0240.0780.2542.3271.0361.0071.0951.012
102.311.2020.0340.2352.3021.2011.0931.0940.999
202.3191.5790.0530.3272.311.5721.0991.1040.996
402.3282.3910.0630.4962.3172.3751.0971.1040.993

4052.3380.9080.0220.22.3210.9061.0811.0840.998
102.351.0170.0230.2042.3351.0061.0821.0930.989
202.3341.2130.030.2252.3171.1981.0791.0930.988
402.3451.613-0.0010.2982.3281.5921.0861.10.987


REMLWeibullWilcoxonPrecision
N2N1BiasMSEBiasMSEBiasMSEWeib/REMLWeib/WilWil/REML

552.1522.390.2420.3072.2192.7151.1210.9871.136
102.153.5630.2230.3822.223.9131.0580.9641.098
202.1265.7670.2360.712.2016.231.0280.9521.08
402.15310.5970.2281.2232.22511.4521.010.9351.081

1052.2741.9950.1030.1452.3242.11.0440.9911.053
102.272.6260.0950.1662.3142.7611.0410.991.051
202.2593.9010.1170.2612.3074.0681.0340.9921.043
402.2746.5480.1050.4022.3186.8521.0290.9831.046

2052.3171.7170.040.1072.3441.7651.0270.9991.028
102.3232.0480.0390.1152.3472.1091.0230.9931.03
202.3232.7310.0360.1462.3482.8041.0220.9951.027
402.3254.1130.0420.2312.3544.2111.0140.9911.024

4052.3531.5740.020.0932.3671.6011.0160.9991.017
102.3491.7430.0350.0992.3641.7671.0130.9991.014
202.3522.0970.0170.1112.3692.1251.0110.9981.013
402.3462.7870.0340.1362.362.8281.0120.9981.015


REMLWeibullWilcoxonPrecision
N2N1 BiasMSEBiasMSEBiasMSEWeib/REMLWeib/WilWil/REML

552.1671.5330.2140.2862.382.0041.0871.0961.092
102.1822.1390.2350.4172.482.5721.0121.0811.013
202.1823.5140.220.6272.4584.1171.0721.0651.1
402.1846.240.2271.2142.4647.2271.0581.0681.084

1052.3231.2270.0960.1552.6391.4321.0671.0911.069
102.3241.5910.1110.1792.6251.8531.0651.0961.062
202.3312.3540.1070.2642.6512.7011.0471.0951.048
402.3043.8570.1010.442.6224.4271.0481.0981.046

2052.3691.0210.0440.1042.6671.1651.0411.0981.039
102.3751.230.0380.1092.6891.3871.0281.0951.03
202.381.6380.0580.1452.691.8421.0251.0951.027
402.3682.4370.0520.2112.6762.7481.0281.1011.025

4052.3950.9360.0420.0932.7071.0381.0101.0881.019
102.4171.0580.0440.0942.7291.1691.0151.0871.016
202.41.2510.0180.1022.7181.3831.0781.0891.015
402.3821.6430.020.1472.6931.8271.0971.0961.015

5.1. Contaminated Weibull Distributed Error Term

A contaminated Weibull distribution of the error terms is also considered. The contamination model is described in (11) that is designed to include outliers and make the distribution shape more skewed. Table 9 consists of bias and MSE of when with =3, it is obvious from estimates of precision that the Weibull score function is efficient than the Wilcoxon score function, and the rank-based method by using Wilcoxon and Weibull is far more efficient than REML. McKean and Kloke [13] compared the rank-based fit by using the skew-normal score function for skew-normal data and compared the empirical relative efficiencies with MLE and Wilcoxon score. It is observed that the most efficient estimates are obtained for the value of parameters for which the score function is derived. Although empirical efficiencies are not significantly different from empirical efficiencies of a nearby range of parameter values. In the same manner, it is observed that estimates with the Weibull score function are most efficient than for neighborhood parametric values. Also, the Weibull score function appeared most efficient than Wilcoxon and MLE.


REMLWeibullWilcoxonPrecision
N2N1BiasMSEBiasMSEBiasMSEWeib/REMLWeib/WilWil/REML

552.73811.7822.3655.6572.4536.0570.480.9340.514
102.71611.1932.2036.2582.2836.5270.5590.9590.583
202.68416.9962.23411.0332.31911.4980.6490.960.676
402.66926.4632.18417.7852.27218.6390.6720.9540.704

1052.9055.8692.4133.3972.4513.4390.5790.9880.586
102.9949.8382.5147.6592.5547.7030.7790.9940.783
202.7747.6372.3395.0972.3735.1890.6670.9820.68
402.93317.4612.47313.6422.51213.830.7810.9860.792

2053.1426.4412.6784.8292.6994.8570.750.9940.754
102.9074.5782.4313.1452.4513.1510.6870.9980.688
202.9345.542.4153.892.4313.9210.7020.9920.708
403.0099.4542.527.1552.5377.2040.7570.9930.762

4053.1426.4412.6784.8292.6994.8570.750.9940.754
102.9074.5782.4313.1452.4513.1510.6870.9980.688
202.9345.542.4153.892.4313.9210.7020.9920.708
403.0099.4542.527.1552.5377.2040.7570.9930.762

From Table 10 for =6, it is clear that the Weibull score is more precise than REML but Wilcoxon scores in almost equal precisely as REML. Moreover, the Weibull score is efficient than the Wilcoxon score function as it produces low estimates of bias and MSE. Table 11 comprises the estimates for =3 when , and it is observed that Weibull and Wilcoxon score functions are almost equally precise. However, the rank-based method appeared efficient in some instances than REML. In Table 12, for =6, at three occasions of sample size, Weibull appeared efficient. REML is found precise in some combinations of level-1 and level-2 sample sizes.


REMLWeibullWilcoxonPrecision
N2N1BiasMSEBiasMSEBiasMSEWeib/REMLWeib/WilWil/REML

559.382.4783.4252.3524.1840.2490.3650.8190.819
109.6082.6324.7532.4174.6950.0370.4951.0121.012
2010.4952.5836.3152.3616.1730.4790.6021.0231.023
4016.9522.64512.2142.43211.9020.40.7211.0261.026

1054.8412.8012.7162.5532.6431.0090.5611.0281.028
104.9772.7823.0562.5232.9160.6270.6141.0481.048
207.6562.8416.012.5765.8131.5660.7851.0341.034
409.7882.8057.642.5447.3190.3940.7811.0441.044

2053.232.8511.9572.5741.8420.2030.6061.0621.062
1032.7992.2032.5242.0870.0570.7341.0551.055
204.1072.8913.1482.6082.9940.5220.7671.0511.051
405.1782.8424.1822.5513.9341.6040.8081.0631.063

4052.1992.8691.5182.5941.4270.8940.6911.0641.064
103.0222.9752.442.6782.3481.1370.8071.0391.039
202.8712.872.242.572.1290.360.781.0521.052
403.8642.8743.1422.5742.9870.7790.8131.0521.052


REMLWeibullWilcoxonPrecision
N2N1BiasMSEBiasMSEBiasMSEWeib/REMLWeib/WilWil/REML

552.2072.2430.2121.2372.3422.5351.130.9891.143
102.1522.6480.1831.6182.2652.7491.0380.9981.04
202.1633.9750.2452.7372.2844.0551.021.0290.991
402.2247.2580.2284.3382.3687.1870.991.0480.944

1052.3391.4570.1080.7352.5051.5131.0381.0171.021
102.3251.7920.0760.8122.491.8291.0211.0450.977
202.3142.550.0861.2642.4632.5891.0151.0780.942
402.3074.1450.0931.962.4564.181.0081.0860.929

2052.321.1230.0330.5612.4981.1441.0191.0490.972
102.3471.2860.0730.6872.5181.3011.0121.0750.941
202.3321.6930.0590.8462.5121.6981.0031.0820.927
402.3762.6010.0411.2662.5432.6221.0081.0940.922

4052.3871.0090.0360.5622.56710.9911.0550.94
102.3961.1260.0180.4922.5871.1190.9941.0760.924
202.3591.2730.0030.5592.5291.281.0051.0940.919
402.3721.728-0.0110.7642.5511.7341.0031.0950.917


REMLWeibullWilcoxonPrecision
N2N1BiasMSEBiasMSEBiasMSEWeib/REMLWeib/WilWil/REML

552.640.4851.0922.3862.981.1252.321.8411.26
102.4092.2910.9962.9762.9442.7551.2031.5880.757
201.6631.8570.090.0411.8982.2721.2231.1060.581
402.79611.0790.3831.3192.71511.8221.0671.1020.968

1051.2261.6950.2640.1051.5741.6440.970.7561.283
102.372.0850.2180.0952.951.9970.9581.3070.733
202.7365.702-0.2220.1482.7655.7821.0141.011.004
402.4244.440.5281.3932.7424.4170.9951.010.985

2052.0840.287-1.0041.2592.4170.3141.0941.6880.648
105.0865.119-0.6820.6984.9145.1421.0040.9931.011
202.6361.7120.2980.1773.0221.7311.0111.230.822
402.1222.0220.8662.2512.3752.11.0391.0560.984

4052.3351.7041.3291.9883.111.270.7450.9380.795
102.6771.154-0.7070.6252.3411.1160.9671.2860.752
201.4260.521-0.5110.3921.540.5691.0921.2080.904
401.690.8580.320.2051.8920.9271.081.2630.855

Table 13 displays the estimated bias and MSE of when with =3. It is obvious that the Weibull score function computes the smallest estimates, leads to the most efficient analysis. Rank-based analysis by using the Wilcoxon score function and REML produces almost equal fit. Table 14, for =6, shows that Weibull score is better than Wilcoxon but in most instances, both score functions produce quite equal fit. Similar behavior is seen between the efficiency of rank-based methods and REML. Table 15, presents the estimates of bias and MSE of when with =3. In most of the sample sizes, the Weibull score is better than the Wilcoxon score and sometimes both techniques performed equally well. Rank-based estimation and REML are found equally precise. Table 16, for =6, indicates that the Weibull score is as equally efficient as the Wilcoxon score for all sample sizes. Similar behavior is seen between the efficiency of rank-based methods and REML. Both techniques produce quite a similar fit. For =3, the Weibull score is found efficient than Wilcoxon, whereas for =6, both score functions showed almost the same precision.


REMLWeibullWilcoxonPrecision
N2N1BiasMSEBiasMSEBiasMSEWeib/REMLWeib/WilWil/REML

552.0852.4590.1970.5212.1152.6921.0570.9651.095
102.0663.4410.230.7542.0663.6030.9980.9531.047
202.0995.7710.2371.2572.1065.9550.990.961.032
402.10210.3790.2081.92.11510.5540.9670.9511.017

1052.2051.9270.1190.3442.2031.9841.0130.9841.03
102.2082.5690.1030.3812.2012.6120.9950.9791.017
202.223.8420.0890.5072.2153.8740.9870.9791.008
402.2166.3960.10.8952.2096.4290.9850.981.005

2052.2751.7010.0290.212.2651.7070.9940.9911.004
102.2762.0070.040.2652.2662.0040.9890.990.999
202.2592.6340.0380.2942.2482.6120.9860.9940.992
402.2683.990.0390.5062.2563.9540.9830.9920.991

4052.2791.5090.0370.2052.2561.5020.9930.9970.995
102.2841.6820.0140.2012.2661.6660.9890.9980.99
202.3072.0540.0040.2482.2932.0240.9840.9990.985
402.2912.6970.0240.2812.2732.6590.9830.9970.986


REMLWeibullWilcoxonPrecision
N2N1 BiasMSEBiasMSEBiasMSEWeib/REMLWeib/WilWil/REML

552.40.7720.6020.7262.2390.7550.9780.8391.166
102.5531.2630.2140.1383.1291.5431.0221.0891.028
201.7522.340.2450.32.0193.2491.0881.0031.084
402.5459.3570.7985.7243.119.1850.9820.8671.033

1051.4730.3991.3032.5461.3480.7011.0570.8871.98
102.1511.254-0.3620.2622.2051.1550.9211.0940.771
201.7911.3170.0180.0012.8471.8331.0921.7420.799
402.1443.688-0.060.0182.4613.5920.9740.9241.054

2051.7650.9030.020.0012.0340.6850.7590.7940.956
101.5840.9280.0920.0131.5770.8750.9431.0960.86
202.1411.046-0.3530.252.3341.1821.0301.0630.971
402.3352.362-0.060.0112.6012.9251.0381.0911.04

4051.9070.947-0.5010.2822.4550.8540.9021.0250.801
102.110.8410.2960.1092.4760.9241.0991.0540.952
202.2221.2260.8811.1652.5241.3791.0251.0231.002
402.3211.6150.0860.0152.6151.7581.0891.0791.009


REMLWeibullWilcoxonPrecision
N2N1BiasMSEBiasMSEBiasMSEWeib/REMLWeib/WilWil/REML

552.1572.3760.230.2872.1672.6131.10.9791.124
102.1433.5790.2420.4322.1383.7531.0490.9551.098
202.1485.870.220.6752.1515.9951.0210.9441.082
402.1410.5540.2381.2312.13510.6411.0080.9361.077

1052.26720.0950.1512.282.1131.0571.0011.056
102.2682.6210.1040.192.2912.7231.0390.9911.048
202.253.8390.1110.2592.2673.9731.0350.9881.048
402.2646.5120.0960.4222.2846.6861.0270.9851.042

2052.3221.7070.0640.1062.3281.7671.0351.0061.029
102.3262.0630.0490.1132.3412.1041.020.9941.026
202.3252.7470.0450.1462.342.8011.020.9941.025
402.3334.140.0470.2132.3514.2061.0160.9931.023

4052.3421.5660.0070.0882.3591.5841.0110.9961.015
102.3391.7280.0310.12.351.7521.0140.9991.015
202.3332.070.0220.0992.3412.1021.01511.015
402.3532.8030.0340.1392.3632.8341.0110.9971.014


REMLWeibullWilcoxonPrecision
N2N1BiasMSEBiasMSEBiasMSEWeib/REMLWeib/WilWil/REML

552.1861.5240.2290.2962.1861.5241.0231.0951.208
102.1922.1510.2240.3962.1922.1511.0111.0761.125
202.2033.5360.2320.692.2033.5361.0741.0651.103
402.216.4280.2321.1732.216.4281.0511.061.086

1052.3071.2070.1090.1542.3071.2071.0811.0941.08
102.3221.5640.1050.1822.3221.5641.0681.0961.059
202.3162.340.1060.252.3162.341.0521.0981.049
402.3183.8920.1050.4412.3183.8921.0461.0951.047

2052.3781.0480.0320.1042.3781.0481.0241.0941.028
102.3631.2290.0650.112.3631.2291.0271.0961.028
202.3621.6260.0360.1452.3621.6261.0391.0991.028
402.3772.4620.0480.2042.3772.4621.0281.1011.024

4052.3870.9320.0260.0922.3870.9321.0111.0921.017
102.3941.0320.0090.0911.05-0.0021.0081.0891.017
202.391.2380.030.1081.2550.0191.0071.0921.014
402.4041.6690.020.1421.6930.011.0111.0961.014

5.2. Real Data Application

The dataset “CASchools” taken from [34] encloses the information on test performance, school characteristics, and student demographic backgrounds for all 420 K-6 and K-8 districts in California with data available for 1998 and 1999. The reading score is a response variable with three explanatory variables included in the model. These include the percentage of students who are English learners (student’s demographic variable), the expenditure per student spent by the school (the US $ 1,000), and the student-teacher ratio is computed by working on school characteristics. The data are analyzed to verify the simulation results because the response variable is Weibull distributed with shape ( = 2.60825) and scale parameter ( = 6.67620). A multilevel model has been applied to these data through REML, rank-based method with Wilcoxon, and Weibull score functions. The estimates of fixed effects, SE, and value are summarized in Table 17.


Rank-based (Weibull)Rank-based (Wilcoxon)REML

EffectsEstimatesSE-valueEstimatesSE-valueEstimatesSE-value
Constant644.06317.344 ≤ 0.001641.58116.193 ≤ 0.001677.58918.507 ≤ 0.001
Std teach ratio-0.2070.6220.740−0.5020.5940.403−0.8960.6070.028
English-0.7400.036 ≤ 0.001−0.7760.040 ≤ 0.001−0.7310.044 ≤ 0.001
Expend0.0050.0020.0190.0060.0020.0010.0010.0020.328
 = 14.923 = 16.176 = 17.701

As the number of students who are English learners (English is not their first language) increases, the average score on reading declines significantly. The improvement in the student-teacher ratio remains unable to change reading scores considerably. The amount of school spent on each student is significant in rank-based estimation but nonsignificant in REML. By using the Weibull score function, minimum SE is observed for all explanatory variables. The estimate of scale parameter is the smallest for Weibull ( = 14.923) among all three methods.

6. Conclusion

We proposed a score function of the rank-based estimation for the multilevel model when the error terms follow the Weibull distribution. To reach the optimal analysis, we compared the performance of Weibull and Wilcoxon score functions with the traditional method REML through a simulation study and a real application. A random intercept multilevel model is analyzed by a rank-based method including three and six covariates for level-1 and level-2 sample sizes under Weibull and contaminated Weibull distributed error terms. The estimates of bias, MSE, and precision of three methods, i.e., Weibull, Wilcoxon, and REML, are computed for each simulation scenario. For estimates of with =3, the Weibull score function is found to be more precise as compared to the Wilcoxon score for , whereas for =6, both score functions showed the same precision. The REML is appeared least efficient among all three methods for all simulation scenarios. Similarly, for contaminated Weibull distribution with , quite similar behavior is observed. For =3, the Weibull score function is mostly efficient than the Wilcoxon score function whereas for =6, both score functions show equal efficiency. For a small value of , the Weibull score function is performed as the most efficient estimation method because the distribution shows the skewed shape. On the other hand, for large values of , the Wilcoxon score function outperformed due to the closed symmetric shape. Moreover, real application results also show that the Weibull score function produces a smaller variation as compared to the Wilcoxon score and REML methods.

7. Limitations and Future Recommendations

Limitations: Weibull score function is derived only for the Weibull distribution of the error term. It is derived for the cluster-correlated error term in linear models. In this study, the Weibull score function is applied to analyze the random intercept multilevel models only.

Future recommendations: the side versions of the Weibull score function could be derived for special cases of Weibull distributions. It could be derived for other dependent error structures like in the time-series AR (1) model of error terms. Its application could be made valid for random intercept and random slope multilevel models.

Data Availability

No data were used to support the findings of this study.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

  1. H. A. Auda, J. W. McKean, J. D. Kloke, and M. Sadek, “A Monte Carlo study of REML and robust rank-based analyses for the random intercept mixed model,” Communications in Statistics - Simulation and Computation, vol. 48, no. 3, pp. 837–860, 2017. View at: Publisher Site | Google Scholar
  2. R. d. P. Garcia, B. S. L. P. de Lima, A. C. d. C. Lemonge, and B. P. Jacob, “A rank-based constraint handling technique for engineering design optimization problems solved by genetic algorithms,” Computers & Structures, vol. 187, pp. 77–87, 2017. View at: Publisher Site | Google Scholar
  3. N. Neelima, D. A. Srikrishna, and D. K. Gangadhara Rao, “Rank-based approach for extracting Unit Pixel Width skeleton,” Materials Science and Engineering, vol. 981, no. 3, Article ID 032005, 2020. View at: Publisher Site | Google Scholar
  4. Z. Li, M. Xie, and M. Zhou, “Rank-based EWMA procedure for sequentially detecting changes of process location and variability,” Quality Technology & Quantitative Management, vol. 15, no. 3, pp. 354–373, 2018. View at: Publisher Site | Google Scholar
  5. T. Chen, W. Tang, Y. Lu, and X. Tu, “Rank regression: an alternative regression approach for data with outliers,” Shanghai Archives of Psychiatry, vol. 26, no. 5, pp. 310–315, 2014. View at: Publisher Site | Google Scholar
  6. W. Liu, J. Sun, and J. Yang, “Rank-based estimation in varying coefficient partially functional linear regression models,” Communications in Statistics, 2020. View at: Google Scholar
  7. A. D.-G. Henry, “RANK-BASED estimation for COBB-DOUGLAS Modelling in the presence OF outliers,” Russian Journal of Agricultural and Socio-Economic Sciences, vol. 78, no. 6, 2018. View at: Google Scholar
  8. J. W. Mckean and T. P. Hettmansperger, “Rank-based analysis of linear models and beyond: a review,” Robust Rank-Based and Nonparametric Methods, pp. 1–24, 2016. View at: Publisher Site | Google Scholar
  9. J. D. Kloke, J. W. McKean, and M. M. Rashid, “Rank-based estimation and associated Inferences for linear models with cluster correlated errors,” Journal of the American Statistical Association, vol. 104, no. 485, pp. 384–390, 2009. View at: Google Scholar
  10. T. P. Hettmansperger and J. W. McKean, Robust Nonparametric Statistical Methods, CRC Press, Boca Rotan, FL, USA, Second Edi edition, 2010.
  11. C.-D. Lai, D. N. Murthy, and M. Xie, “Weibull distributions and their applications,” Weibull Distributions and Their Applications Springer Handbook of Engineering Statistics, Springer, Berlin, Germany, 2006. View at: Publisher Site | Google Scholar
  12. L. D. Jones, L. J. Vandeperre, T. A. Haynes, and M. R. Wenman, “Theory and application of Weibull distributions to 1D peridynamics for brittle solids,” Computer Methods in Applied Mechanics and Engineering, vol. 363, p. 112903, 2020. View at: Publisher Site | Google Scholar
  13. J. W. McKean and J. D. Kloke, “Efficient and adaptive rank-based fits for linear models with skew-normal errors,” Journal of Statistical Distributions and Applications, vol. 1, no. 1, pp. 1–18, 2014. View at: Publisher Site | Google Scholar
  14. J. W. McKean and G. L. Sievers, “Rank scores suitable for analyses of linear models under asymmetric error distributions,” Technometrics, vol. 31, no. 2, pp. 207–218, 1989. View at: Publisher Site | Google Scholar
  15. S. Zeng, S.-M. Chen, and L.-W. Kuo, “Multiattribute decision making based on novel score function of intuitionistic fuzzy values and modified VIKOR method,” Information Sciences, vol. 488, pp. 76–92, 2019. View at: Publisher Site | Google Scholar
  16. M. Lalancette, S. Engelke, and S. Volgushev, Rank-based Estimation under Asymptotic Dependence and Independence, with Applications to Spatial Extremes, Cornell University, Ithaca, NY, USA, 2020.
  17. K. Watcharotone, J. W. Mckean, and B. E. Huitema, Communications in Statistics-Simulation and Computation A Monte Carlo Study of Traditional and Rank-Based Picked-Points Analyses A Monte Carlo Study of Traditional and Rank-Based Picked-Points Analyses, Taylor & Francis, Abingdon, UK, 2017.
  18. J. D. Kloke and J. W. McKean, “Rfit: rank-based estimation for linear models,” The R Journal, vol. 4, no. 2, pp. 57–64, 2012. View at: Publisher Site | Google Scholar
  19. Y. Shao, J. W. Mckean, and B. E. Huitema, “Traditional and rank-based tests for ordered alternatives in a cluster corrrelated model,” Psychometrika, vol. 85, no. 3, pp. 531–554, 2020. View at: Publisher Site | Google Scholar
  20. H. F. Bindele, A. Abebe, and K. N. Meyer, “General rank-based estimation for regression single index models,” Annals of the Institute of Statistical Mathematics, vol. 70, no. 5, pp. 1115–1146, 2018. View at: Publisher Site | Google Scholar
  21. A. Abebe and J. W. McKean, “Weighted Wilcoxon estimators in nonlinear regression,” Australian & New Zealand Journal of Statistics, vol. 55, no. 4, pp. 401–420, 2013. View at: Publisher Site | Google Scholar
  22. J. T. Terpstra and J. W. Mckean, “Rank-based analyses of linear models using R,” Journal of Statistical Software, vol. 14, 2005. View at: Google Scholar
  23. H. F. Bindele, “Strong consistency of the general rank estimator,” Communications in Statistics - Theory and Methods, vol. 46, no. 2, pp. 532–539, 2017. View at: Publisher Site | Google Scholar
  24. A. Abebe, J. W. McKean, J. D. Kloke, and Y. K. Bilgic, “Iterated reweighted rank-based estimates for GEE models,” Robust Rank-Based and Nonparametric Methods, vol. 168, pp. 61–79, 2016. View at: Publisher Site | Google Scholar
  25. M. Černý, M. Rada, and O. Sokol, “Rank estimators for robust regression: approximate algorithms, exact algorithms, and two-stage methods,” in Lecture Notes in Computer Science (including subseries Lecture Notes in Artificial Intelligence and Lecture Notes in Bioinformatics), pp. 162–173, Springer, Berlin, Germany, 2020. View at: Google Scholar
  26. M. Cerny, M. Rada, J. Antoch, and M. Hladik, A Class of Optimization Problems Motivated by Rank Estimators in Robust Regression, Taylor & Francis, Abingdon, UK, 2019.
  27. S. Dutta and S. Datta, “Rank‐based inference for covariate and group effects in clustered data in presence of informative intra‐cluster group size,” Statistics in Medicine, vol. 37, no. 30, pp. 4807–4822, 2018. View at: Publisher Site | Google Scholar
  28. S. Dutta and S. Datta, “A rank-sum test for clustered data when the number of subjects in a group within a cluster is informative,” Biometrics, vol. 72, no. 2, pp. 432–440, 2016. View at: Publisher Site | Google Scholar
  29. T. Xie, R. Cao, and P. Yu, “Rank-based test for partial functional linear regression models,” Journal of Systems Science and Complexity, vol. 33, no. 5, pp. 1571–1584, 2020. View at: Publisher Site | Google Scholar
  30. M. A. Scott, J. S. Simonoff, and B. D. Marx, The SAGE Handbook of Multilevel Modeling, Sage, Thousands Oaks, CA, USA, 2013.
  31. J. Hájek and Z. Sidák, Theory of Rank Tests, Academic Press, Cambridge, MA, USA, 1967.
  32. M. Fréchet, “Sur la loi de probabilité de l’écart maximum,” in Annales de la societe Polonaise de Mathematique, Polish Mathematical Society, Warsaw, Poland, 1928. View at: Google Scholar
  33. W. Weibull, “A statistical distribution function of wide applicability,” Journal of Applied Mechanics, vol. 18, pp. 290–293, 1951. View at: Google Scholar
  34. J. H. Stock and M. W. Watson, Introduction to Econometrics, Pearson Education, London, UK, 2015.

Copyright © 2021 Sehar Saleem et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Related articles

No related content is available yet for this article.
 PDF Download Citation Citation
 Download other formatsMore
 Order printed copiesOrder
Views343
Downloads327
Citations

Related articles

No related content is available yet for this article.

Article of the Year Award: Outstanding research contributions of 2021, as selected by our Chief Editors. Read the winning articles.