# Some New Robust Estimators for Circular Logistic Regression Model with Applications on Meteorological and Ecological Data

Authors: Shokrya S. Alshqaq; Abdullah A. Ahmadini; Ali H. Abuzaid
Journal: Mathematical Problems in Engineering
 (2021)

---

## Abstract

Maximum likelihood estimation (MLE) is often used to estimate the parameters of the circular logistic regression model due to its efficiency under a parametric model. However, evidence has shown that the classical MLE extremely affects the parameter estimation in the presence of outliers. This article discusses the effect of outliers on circular logistic regression and extends four robust estimators, namely, Mallows, Schweppe, Bianco and Yohai estimator BY, and weighted BY estimators, to the circular logistic regression model. These estimators have been successfully used in linear logistic regression models for the same purpose. The four proposed robust estimators are compared with the classical MLE through simulation studies. They demonstrate satisfactory finite sample performance in the presence of misclassified errors and leverage points. Meteorological and ecological datasets are analyzed for illustration.


---

## Body

## 1. Introduction

Circular data arise whenever the values of a random variable can present the circumference of the unit circle. It is measured by angles with values between 0 and2π or 0° and 360°, for example, wind directions, animal navigation, or the values of any periodic phenomena such as a 24-hour clock or days of the year which can be converted to circular data [1]. The modeling of the relationship between circular variables is so-called “circular regression,” and it is classified into three main classes, namely, circular-circular, circular-linear, and linear-circular regression models [2].The applications on circular regression models are widely spread in many applied fields. Several regression models were proposed to predict a continuous circular variable from other circular or linear predictors [3–6].Logistic regression analysis is a useful statistical tool that analyses the relationship between a binary response and a predictor. The theory of logistic regression is well developed [7]. Daffaie and Khan [8] proposed the circular logistic model to predict a binary variable from a circular variable, such as modeling the rainfall (yes, no) and wind direction, fatal road accident (yes, no), and time of accident.The existence of outliers is a common problem in regression analysis. For the linear logistic model, Feser and Pia [9] showed that maximum likelihood estimation can be influenced by outliers. Croux et al. [10] found that the most dangerous outliers termed “bad leverage points” are misclassified observations that are outlying in the design space of predicted variables.Circular logistic regression is also subjected to the existence of outliers as shown in [11], where they have proposed an outlier detection procedure based on the penalized maximum likelihood and applied it on different real datasets.Several robust estimators that are less affected by outliers are proposed in the literature to improve the estimation performance in linear logistic regression models [12]. The authors in [13] introduced weights depending on the response and covariates and proposed the Mallows-type estimator. This estimator was analyzed deeply by Carroll and Pederson [14] using Mahalanobis distance to reduce the weights in terms of leverage, and Bianco and Yohai [15] proposed BY methods.Since the published work just considered the detection of outliers in the circular logistic regression model, this article attempts to overcome the problem of outliers in the circular logistic regression model by extending some robust estimators from the classical linear logistic to the circular logistic case.The rest of this paper is organized as follows. Section2 reviews the formulation of the circular logistic regression model and its parameter estimation via MLE. Section 3 presents the types of outliers in circular logistic regression and derives the proposed robust estimators in a logistic circular regression model. Section 4 discusses the effect of outliers on the circular logistic estimators by computing the influence functions. Section 5 investigates the performance of the considered robust estimators. Section 6 applies the proposed estimators on meteorological and ecological datasets. Section 7 provides the conclusion.


## 2. Model Formulation

A circular logistic regression describes the relationship between a binary response and circular predictors. It shows potential with various applications in the field of environmental sciences. The authors in [8] assumed that n binomial observations with a probability of success, pi=vi/ni, for i=1,…,k, depend on circular random variable ui, and the proposed model is given as follows:(1)logitpi=logpi1−pi=β0+β⋆cosui−u0,where β0 is the value of the logit (log odds) when ui−u0=90° and u0 is the angle where the logit reaches its highest value. Let β⋆=β12+β22 and u0=arctanβ2/β1, and equation (1) can be written as(2)logitpi=logpi1−pi=β0+β1cosui+β2sinui.Suppose that binomial data of the formvi successes out of ni trails are observations from a binomial distribution, the likelihood function is then given by(3)Lβ=∏i=1knivipivi1−pini−vi.Letηi=β0+β1cosui+β2sinui, and by using the exponential function, we obtain(4)pi1−pi=eηi⇒pi=eηi1+eηi.The maximum likelihood estimationMLE is classically used for parameter estimation and is defined by an objective function as(5)β^MLE=argmaxβ∑i=1klogLvi,ui,β,where(6)logLvi,ui,β=lognivi+vilogpi+ni−vilog1−pi.The maximum likelihood equations are given as follows:(7)∂logLβ∂β0=∑i=1kvi−∑i=1kni1+e−ηi,(8)∂logLβ∂β1=∑i=1kvicosui−∑i=1knicosui1+e−ηi,(9)∂logLβ∂β2=∑i=1kvisinui−∑i=1knisinui1+e−ηi.These equations are solved iteratively by using the Newton–Raphson method. Recently, Abuzaid and ElShekh Ahmed [11] used the penalized maximum likelihood estimator (PMLE) to identify outliers in the circular logistic regression model and investigated its performance via simulation. The following section discusses some possible robust estimators for the circular logistic regression model.


## 3. Robust Estimators for Circular Logistic Regression

### 3.1. Outliers in Circular Logistic Regression

This section distinguishes different cases of outlying observations in circular logistic regression, where outliers might occur in the dependent variable, independent variable, or both of them.For binary data, all thev’s values are either 0 or 1. Hence, an error in the vi direction can only occur as a transposition of 0 to 1 or vice versa. This type of outlier is known as residual outlier or misclassification-type error [16,17].A leverage outlier or leverage point occurs when the circular observation at positiond (e.g., ud) is contaminated as follows: ud⋆=ud+γπmod2π, where ud⋆ is the value after contamination and γ is the degree of contamination in the range 0≤γ≤1. A leverage point can be considered as a good leverage point when vi=1 with a large value of Pvi=1|ui while it is a bad leverage point when vi=0 with small value of Pvi=1|ui and vice versa. Abuzaid and ElShekh Ahmed [11] considered the misclassification-type error outlier in their simulation study without any investigation of the leverage point detection.Alternatively, robust estimators are the common methods used for handling the problem of outliers in logistic models. The following section presents four robust estimators for the parameters in the circular logistic model as follows.


### 3.2. Circular Mallows Class

The proposed estimator (MALLOWSC) is extended from the Mallows class in [18] to the circular logistic model for weighting the maximum likelihood estimator.AssumeF is a continuous and increasing distribution function and is given by(10)Fηi=e−ηi1+e−ηi.Then, the partial derivatives in (7)–(9) become(11)∂lβ∂β0=∑i=1kvi−Fηi,∂lβ∂β1=∑i=1kcosuivi−Fηi,∂lβ∂β2=∑i=1ksinuivi−Fηi,respectively. The robust estimates for the circular logistic regression model in equation (10) are given by the solution of β obtained by(12)∂lβ∂β0=∑i=1kwivi−Fηi−cui,β0=0,(13)∂lβ∂β1=∑i=1kwicosuivi−Fηi−cui,β1=0,(14)∂lβ∂β2=∑i=1kwisinuivi−Fηi+cui,β2=0,where wi are the weights that may depend on ui, vi, or both and c. is a correction function needed to ensure consistency. If wi=1 and c.=0, then equations (12)–(14) give the usual circular logistic regression estimate. If wi=wui,ηi and c.=0, then the weights depend only on ui.


### 3.3. Circular Schweppe Class

Stefanski [19] stated that wi is robust Mahalanobis distance for the ui vector that depends on the covariance matrix of the regression model, which is given by(15)Σ^=ui′Wui−1,where W is a diagonal matrix with wii=pi1−pi and 1−pi is the probability pi of vi=1; then, the Mahalanobis distance is given by(16)DM=ui−ui¯′Σ^−1ui−ui¯.Ifwi=wui,ηi,vi, then the estimator is the same as the linear logistic Schweppe class proposed in [13]. Here, wi depends on vi and circular ui. This estimator is called circular conditionally unbiased bounded influence function or CUBIFC estimator.


### 3.4. Circular BY Estimators

Bianco and Yohai [15] proposed BY methods for the linear logistic model, and in this section, we extend it to the circular case and referred to as BYC.Let theMLE be obtained by minimizing the deviance,(17)β^MLE=argminβ∑i=1kDiβ,where Diβ=−vilnFηi−1−viln1−Fηi. By replacing the deviance function in (17) with function ρ, the robust BYC estimator is defined by(18)β^BYC=argminβ∑i=1kρDFηi,vi+GFηi+G1−Fηi,where ρ is a bounded, differentiable, and nondecreasing function defined in [15] and given by(19)ρu=u−u22m,ifu≤m,m2,otherwise,where m is a positive number, ρ′u=ψu, and Gu=∫0uψ−lngdg.


### 3.5. Circular WBY Estimators

The extension ofBY estimator by including weights is reducing the influence of outliers in ui space. This weighted BY (WBY) estimator is extended to the circular logistic regression model and defined as(20)β^WBYC=argminβ∑i=1kwuiρDFηi,vi+GFηi+G1−Fηi,where the weights wui are distances which are computed using the minimum covariance determinant (MCD) estimator, to be a decreasing function of robust Mahalanobis distances (RDi), and given by (see [20])(21)wui=1,ifRDi2≤χp,0.9752,0,elsewhere.


## 3.1. Outliers in Circular Logistic Regression

This section distinguishes different cases of outlying observations in circular logistic regression, where outliers might occur in the dependent variable, independent variable, or both of them.For binary data, all thev’s values are either 0 or 1. Hence, an error in the vi direction can only occur as a transposition of 0 to 1 or vice versa. This type of outlier is known as residual outlier or misclassification-type error [16,17].A leverage outlier or leverage point occurs when the circular observation at positiond (e.g., ud) is contaminated as follows: ud⋆=ud+γπmod2π, where ud⋆ is the value after contamination and γ is the degree of contamination in the range 0≤γ≤1. A leverage point can be considered as a good leverage point when vi=1 with a large value of Pvi=1|ui while it is a bad leverage point when vi=0 with small value of Pvi=1|ui and vice versa. Abuzaid and ElShekh Ahmed [11] considered the misclassification-type error outlier in their simulation study without any investigation of the leverage point detection.Alternatively, robust estimators are the common methods used for handling the problem of outliers in logistic models. The following section presents four robust estimators for the parameters in the circular logistic model as follows.


## 3.2. Circular Mallows Class

The proposed estimator (MALLOWSC) is extended from the Mallows class in [18] to the circular logistic model for weighting the maximum likelihood estimator.AssumeF is a continuous and increasing distribution function and is given by(10)Fηi=e−ηi1+e−ηi.Then, the partial derivatives in (7)–(9) become(11)∂lβ∂β0=∑i=1kvi−Fηi,∂lβ∂β1=∑i=1kcosuivi−Fηi,∂lβ∂β2=∑i=1ksinuivi−Fηi,respectively. The robust estimates for the circular logistic regression model in equation (10) are given by the solution of β obtained by(12)∂lβ∂β0=∑i=1kwivi−Fηi−cui,β0=0,(13)∂lβ∂β1=∑i=1kwicosuivi−Fηi−cui,β1=0,(14)∂lβ∂β2=∑i=1kwisinuivi−Fηi+cui,β2=0,where wi are the weights that may depend on ui, vi, or both and c. is a correction function needed to ensure consistency. If wi=1 and c.=0, then equations (12)–(14) give the usual circular logistic regression estimate. If wi=wui,ηi and c.=0, then the weights depend only on ui.


## 3.3. Circular Schweppe Class

Stefanski [19] stated that wi is robust Mahalanobis distance for the ui vector that depends on the covariance matrix of the regression model, which is given by(15)Σ^=ui′Wui−1,where W is a diagonal matrix with wii=pi1−pi and 1−pi is the probability pi of vi=1; then, the Mahalanobis distance is given by(16)DM=ui−ui¯′Σ^−1ui−ui¯.Ifwi=wui,ηi,vi, then the estimator is the same as the linear logistic Schweppe class proposed in [13]. Here, wi depends on vi and circular ui. This estimator is called circular conditionally unbiased bounded influence function or CUBIFC estimator.


## 3.4. Circular BY Estimators

Bianco and Yohai [15] proposed BY methods for the linear logistic model, and in this section, we extend it to the circular case and referred to as BYC.Let theMLE be obtained by minimizing the deviance,(17)β^MLE=argminβ∑i=1kDiβ,where Diβ=−vilnFηi−1−viln1−Fηi. By replacing the deviance function in (17) with function ρ, the robust BYC estimator is defined by(18)β^BYC=argminβ∑i=1kρDFηi,vi+GFηi+G1−Fηi,where ρ is a bounded, differentiable, and nondecreasing function defined in [15] and given by(19)ρu=u−u22m,ifu≤m,m2,otherwise,where m is a positive number, ρ′u=ψu, and Gu=∫0uψ−lngdg.


## 3.5. Circular WBY Estimators

The extension ofBY estimator by including weights is reducing the influence of outliers in ui space. This weighted BY (WBY) estimator is extended to the circular logistic regression model and defined as(20)β^WBYC=argminβ∑i=1kwuiρDFηi,vi+GFηi+G1−Fηi,where the weights wui are distances which are computed using the minimum covariance determinant (MCD) estimator, to be a decreasing function of robust Mahalanobis distances (RDi), and given by (see [20])(21)wui=1,ifRDi2≤χp,0.9752,0,elsewhere.


## 4. Influence Function

Suppose two source populations withk-dimensional circular variables, which are both von Mises distribution with different means but the same concentration parameter, κ. Circular variable ui can arise from one of these populations(22)ui∼H1=viMμ1,κ,with probabilityp1,H0=viMμ0,κ,with probabilityp0,where i=1,…,k and p1+p0=1. Let the binary variable vi indicates the source population of the corresponding ui, then(23)vi∼1,with probabilityp1,0,with probabilityp0.Let the joint distribution of (ui, vi) be denoted by H and T be an estimator of the circular logistic parameters. Then, the influence function [21] is defined as(24)IFui,vi,T,H=limε⟶∞T1−εH+εΔui,vi−THε=∂∂εT1−εH+εΔui,vi|ε=0.IfT=MLE, as shown in equations (7)–(9), then the influence function is an unbounded function for spaces ui and vi. Specifically, a small amount of contamination in the training data due to the presence of possible outlier in ui or vi intensively affects the MLE, as shown in the simulation section.SupposeT=MALLOWSC, where the weights depend on the robust distance of observation ui, and robust distance RDi is equal to the Mahalanobis distance of ui to the center of the data cloud. This condition reduces the influence of outlying observations in the ui space. Thus, IFui,vi,MALLOWSC,H is bounded with respect to ui but unbounded with respect to vi. Similar conclusions as T=BYC can be derived.IfT=CUBIFC or T=WBYC, which adds a weight to BYC, then a fully bounded influence is obtained.


## 5. Simulation Study

### 5.1. Settings

This simulation aims to compare the robustness of the proposed robust estimators and the classicalMLE. The independent circular variable is generated from von Mises distribution with mean 60° and concentration parameter of κ  = 1, 2, 6, 10, and 15, with sample size N=200,250, and 300; a large sample size is chosen to avoid separation problems. The true parameter values are β=0,2,2.The simulation study is reported in a variety of situations. Initially, the data without contamination are simulated. The robust properties of all estimators with contaminated data are examined in three different ways. First, proportions are taken from the responses, andv is chosen randomly and changed from either 0 to 1 or 1 to 0. This process constitutes the misclassification-type error. For each contaminated case, 5%, 10%,20%, 30%, and 40% of the original data are contaminated. Second, the same proportions are taken to contaminate u with γ=0.5 for good leverage points. Finally, the same proportions are considered, and the generated data are contaminated with two types of outliers simultaneously. This process constitutes bad leverage points.Each simulation includes 1000 replications. The performance of the estimators is evaluated based on the bias and the median squared error (MSE) for each parameter, which are defined as follows:(25)bias=11000∑i=11000β^i−βi,MSE=median∑i=11000β^i−β¯i2.A good estimator has bias andMSE that are relatively small or close to zero.The simulation used the standard available “Robust” package in R to obtain the estimators and the” CircStats” package for generating the circular variableu.


### 5.2. Results

The bias andMSE of the five estimators are shown in Tables 1 to 7. Table 1 shows the results for uncontaminated data (i.e., clean data), where the biases and MSEs of all five estimators are fairly close to each other. However, the MALLOWSC’s MSE is larger than the others for large concentration parameter (κ  = 6, 10, and 15). Hence, MALLOWSC performs worse compared with the other estimators in this situation.Table 1
Bias andMSE of all methods for uncontaminated data.
κMethodsN=200N=250N=300BiasMSEBiasMSEBiasMSE1MLE0.00490.01660.00230.01270.00250.0114MALLOWSC0.00290.01940.00170.01400.00280.0123CBIFC0.00380.01660.00210.01280.00220.0114BYC0.00510.01790.00280.01350.00280.0117WBYC0.00500.01790.00300.01350.00300.01182MLE0.00640.02130.00120.01661.2e-030.0131MALLOWSC0.01130.14960.00550.06012.3e-030.0406CBIFC0.00620.02330.00100.01771.6e-030.0136BYC0.00830.02200.00240.01711.0e-030.0136WBYC0.00570.13960.00100.01824.2e-050.01436MLE0.01700.22210.01100.08490.00660.0642MALLOWSC0.10990.87340.07670.34990.04170.2006CBIFC0.01970.23460.01760.09630.01050.0749BYC0.00600.22480.00430.09080.00330.0641WBYC0.03070.24020.00070.03160.04710.063210MLE0.02160.13920.01620.11550.01680.0938MALLOWSC0.06340.36000.06240.27690.06090.2530CBIFC0.03410.15290.03190.13110.01600.1030BYC0.03260.14190.01730.11820.01770.0948WBYC0.07640.12640.07190.10900.07520.087115MLE0.08950.25010.07060.19920.03300.1789MALLOWSC0.15181.02600.11580.48440.14140.4132CBIFC0.07840.27490.07160.21370.05290.1945BYC0.08590.25600.06260.20350.03240.1805WBYC0.09340.97540.08350.19490.08160.1661Table 2
Bias andMSE of all methods for data with various percentages of misclassified errors (5%, 10%, and 20%).
κ% of misc errorMethodsN=200N=250N=300BiasMSEBiasMSEBiasMSE15MLE1.45702.19181.55702.20181.65812.3018MALLOWSC0.14000.30730.25100.41620.13670.3612CBIFC0.51480.44820.46480.34720.46360.3371BYC0.13440.37650.12330.27560.11440.2565WBYC0.17450.38420.08650.39450.08410.29512MLE1.55702.29181.65702.30181.75812.4018MALLOWSC0.24000.40730.35100.51620.23670.4612CBIFC0.61480.54820.56480.44720.56360.4371BYC0.23440.47650.22330.37560.21440.3565WBYC0.27450.48420.18650.49450.18410.39516MLE1.05702.18181.05702.20181.05812.3018MALLOWSC0.14000.31730.25100.62620.12670.5512CBIFC0.50480.42820.47480.33720.45360.2271BYC0.02440.32650.01330.26560.01440.2465WBYC0.06450.34420.07650.38450.07410.285110MLE1.55802.20181.66802.31281.76912.4128MALLOWSC0.03100.21830.14000.30520.02570.2502CBIFC0.40380.33720.35380.23620.35260.2261BYC0.02340.26550.01230.16460.00340.1455WBYC0.06450.27320.07550.28450.07310.184115MLE1.56812.20291.66812.31291.76922.4129MALLOWSC0.13110.31620.24010.40510.12560.3501CBIFC0.50370.43710.45370.33610.45250.3260BYC0.12330.36540.11220.26450.10330.2454WBYC0.16340.37310.07540.38340.07300.2840110MLE2.46486.10132.35376.00022.25266.1113MALLOWSC0.34840.20980.23930.10870.12920.0076CBIFC0.89331.45980.99220.94970.98110.9386BYC0.25420.17160.14310.06050.03200.0504WBYC0.05650.11200.04540.01100.03430.01212MLE2.36486.00132.25376.01022.15266.0113MALLOWSC0.24840.10980.13930.00870.02920.0066CBIFC0.69330.75980.89220.84970.88110.8386BYC0.15420.07160.04310.05050.02200.0404WBYC0.04650.01200.03540.00100.02430.00216MLE2.35486.00132.24376.01122.14266.0013MALLOWSC0.23840.11980.12930.01870.01920.0176CBIFC0.9833.64980.88220.83970.87110.8286BYC0.14420.06160.03310.15050.12200.1404WBYC0.14650.00200.13540.10100.12430.102110MLE2.57586.21232.46476.11122.14166.2223MALLOWSC0.23740.11880.12830.00770.01820.0166CBIFC0.88230.54880.88120.83870.87010.8276BYC0.14320.06060.03210.15150.12100.1414WBYC0.04550.00300.03440.00000.02330.001115MLE2.57596.21242.46486.11132.36376.2224MALLOWSC0.23730.10870.12820.00760.01810.0065CBIFC0.68220.74870.88110.83860.87000.8275BYC0.14310.06050.03210.05040.02100.0403WBYC0.04540.00100.03430.00000.02320.0010120MLE2.72887.47732.61777.36622.50667.2551MALLOWSC0.43090.64670.32080.53560.21070.4245CBIFC1.66031.80311.55022.70201.44012.6010BYC0.76881.32570.65771.21460.54661.1035WBYC0.07030.32170.06020.21060.05010.10052MLE3.72888.47733.61778.36623.50668.2551MALLOWSC0.53090.74670.42080.63560.31070.5245CBIFC2.66033.80312.55023.70202.44013.6010BYC0.86882.32570.75772.21460.64662.1035WBYC0.17030.42170.16020.31060.15010.20056MLE2.83887.58732.72777.47622.61667.3651MALLOWSC0.52090.75670.43080.64560.32070.5345CBIFC1.77032.91311.66022.81201.55012.7110BYC0.87881.43570.76771.32460.65661.2135WBYC0.18030.43170.17020.21060.16010.210510MLE3.93888.68733.82778.57623.51668.4651MALLOWSC0.41090.64670.32080.53560.21070.4245CBIFC0.67031.81310.56021.71200.45011.6110BYC0.76880.53570.65770.22460.54660.1135WBYC0.07030.32170.06020.10060.05010.100515MLE4.04889.79734.93779.68624.62669.5751MALLOWSC0.40080.63560.31070.52450.20060.4134CBIFC0.66021.80200.55011.70100.44001.6000BYC0.75770.52460.64660.21350.53550.1024WBYC0.06020.31060.05010.10050.04000.1004Table 3
Bias andMSE of all methods for data with various percentages of misclassified errors (30% and 40%).
κ% of misc errorMethodsN=200N=250N=300BiasMSEBiasMSEBiasMSE130MLE2.82478.00533.72477.00432.61476.0032MALLOWSC0.43540.56140.32540.65140.21540.5414CBIFC2.03184.16581.02083.15481.01083.0448BYC2.82588.01122.71588.00122.60588.0001WBYC0.37520.55600.26420.44500.15320.33402MLE3.82479.00534.72478.00433.61477.0032MALLOWSC0.53540.66140.42540.75140.31540.6414CBIFC3.03185.16582.02084.15482.01084.0448BYC3.82589.01123.71589.00123.60589.0001WBYC0.47520.65600.36420.54500.25320.43406MLE3.93479.11534.83478.11433.72477.1132MALLOWSC0.54640.67240.43640.76240.32640.6524CBIFC3.04285.17682.03184.16582.02184.0558BYC3.83689.02223.72689.01223.61689.0000WBYC0.48620.66700.37520.55600.26420.445010MLE3.93589.11644.83588.11543.72587.1143MALLOWSC0.32430.45030.21430.54030.10430.4303CBIFC2.02074.05471.01073.04371.00073.0337BYC2.71478.00012.60478.00012.50478.0000WBYC0.26410.66500.15310.33400.04210.223015MLE3.82589.00644.72588.00543.61587.0043MALLOWSC0.21430.34030.10430.43030.00430.3203CBIFC2.13074.16471.12073.15371.11073.1437BYC2.82478.11012.71478.11012.61478.1100WBYC0.37410.77500.26310.44400.15210.3330140MLE2.88388.53202.77278.42102.66168.3100MALLOWSC0.77160.88490.76050.87380.75040.8627CBIFC2.42875.93372.31775.02272.20665.0116BYC2.87718.31482.76608.20372.65508.1026WBYC0.05150.39610.04040.28700.03140.29806MLE2.99388.64202.88278.53102.77168.4200MALLOWSC0.88160.99490.87050.98380.86040.9727CBIFC2.53875.04372.42775.13272.31665.1216BYC2.98718.42482.87608.31372.76508.2126WBYC0.16150.40710.05150.29810.04250.209110MLE3.00499.54213.98289.63113.87179.5201MALLOWSC0.99261.00300.98171.09490.97151.0838CBIFC2.87535.37042.77425.27132.66315.1612BYC2.71898.48422.60878.37312.50768.2621WBYC0.51610.71400.15050.81290.25040.201315MLE3.11509.65323.09399.74223.98289.6312MALLOWSC1.09371.11411.98171.10501.06261.9747CBIFC2.98635.48172.88535.38242.77415.2723BYC2.82908.59532.71988.48422.61878.3732WBYC0.62720.82510.26160.92300.36150.3124Table 4
Bias andMSE of all methods for data with various percentages of good leverage points (5%, 10%, and 20%).
κ% of misc errorMethodsN=200N=250N=300BiasMSEBiasMSEBiasMSE15MLE0.80720.81220.81620.82320.82720.8342MALLOWSC0.81150.81880.82250.82980.81050.8188CBIFC0.81220.82190.82320.83290.83420.8439BYC0.81380.83560.82480.84660.83580.8576WBYC0.82090.86230.83190.87330.84290.88432MLE0.87200.82210.86210.83220.87220.8423MALLOWSC0.81510.88810.82520.89820.80510.8881CBIFC0.82210.81920.83220.82930.84230.8394BYC0.83810.85630.84820.86640.85830.8765WBYC0.80920.82360.81930.83370.82940.84386MLE0.81720.92220.92620.93320.93720.9442MALLOWSC0.92150.92880.93250.93980.92050.9288CBIFC0.92220.93190.93320.94290.94420.9539BYC0.92380.94560.93480.95660.94580.9676WBYC0.93090.97230.94190.98330.95290.994310MLE0.87120.82220.86220.83320.87320.8442MALLOWSC0.81520.88280.82530.89380.80250.8828CBIFC0.82220.81390.83320.82490.84420.8359BYC0.83280.85460.84380.86560.85480.8766WBYC0.80390.82730.81490.83830.82590.849315MLE0.82170.82230.82260.82330.82370.8244MALLOWSC0.82510.82880.83520.88390.85200.8228CBIFC0.82230.89310.82330.89420.82440.8953BYC0.88230.86450.88340.85660.88450.8667WBYC0.89300.83720.89410.88330.89520.8394110MLE0.80350.80960.71350.71960.70250.7086MALLOWSC0.80830.81830.71830.72820.72910.7393CBIFC0.80600.81430.71600.72430.72700.7353BYC0.81210.84160.72210.73160.73310.7426WBYC0.81180.85060.72180.76060.71280.70662MLE0.83550.89660.73550.79660.72550.7866MALLOWSC0.88330.88230.78330.78220.79110.7933CBIFC0.86000.84330.76000.74330.77000.7533BYC0.82110.81660.72110.71660.73110.7266WBYC0.81880.80660.71880.70660.72880.76666MLE0.91530.91960.81530.81690.81520.8168MALLOWSC0.91380.92380.82380.82280.82190.8339CBIFC0.91060.91340.81060.82340.82070.8335BYC0.91120.94610.82120.83610.83130.8462WBYC0.91810.95600.82810.86600.81820.806510MLE0.81450.81070.72450.72070.71350.7196MALLOWSC0.81930.82930.72930.73920.73020.7404CBIFC0.81700.82530.72700.73530.73800.7463BYC0.82310.85260.73310.74260.74410.7536WBYC0.82280.86160.73280.77160.72380.717615MLE0.82540.82070.73540.73700.72530.7269MALLOWSC0.82390.83390.73390.74290.74200.7540CBIFC0.82070.83350.73070.74350.74080.7536BYC0.83130.86620.74130.75620.75140.7663WBYC0.83820.87610.74820.78610.73830.7267120MLE0.79100.79030.60910.60940.61020.6105MALLOWSC0.79540.79790.60640.60890.61750.6190CBIFC0.79110.79220.68110.68220.67220.6733BYC0.80190.78670.71190.69670.72200.6078WBYC0.78670.81010.67670.70010.66560.72232MLE0.79110.79220.68110.68220.67220.6733MALLOWSC0.80190.78670.71190.69670.72200.6078CBIFC0.79100.79030.60910.60940.61020.6105BYC0.79540.79790.60640.60890.61750.6190WBYC0.77860.80110.66770.70100.65660.72326MLE0.79100.79030.60910.60940.61020.6105MALLOWSC0.79540.79790.60640.60890.61750.6190CBIFC0.79110.79220.68110.68220.67220.6733BYC0.80190.78670.71190.69670.72200.6078WBYC0.78670.81010.67670.70010.66560.722310MLE0.80200.80130.71010.71040.72120.7215MALLOWSC0.80640.80890.71740.71990.72850.7200CBIFC0.80210.80320.79210.79320.78320.7843BYC0.91290.89770.82290.70770.83300.7188WBYC0.89770.92110.78770.81110.77660.833315MLE0.81310.81240.72120.72150.73230.7326MALLOWSC0.81750.81900.72850.72000.73960.7311CBIFC0.81320.81430.70320.70430.79430.7954BYC0.92310.80880.83300.71880.82210.7299WBYC0.80880.93220.79880.82110.78770.8433Table 5
Bias andMSE of all methods for data with various percentages of good leverage points (30% and 40%).
κ% of misc errorMethodsN=200N=250N=300BiasMSEBiasMSEBiasMSE130MLE0.80890.83920.71890.74920.79890.7922MALLOWSC0.81240.84520.72240.75520.72440.7522CBIFC0.81110.84420.72110.75420.71120.7425BYC0.81500.86010.72500.77010.70520.7107WBYC0.81090.86320.72090.77320.70920.73272MLE0.81900.84030.72900.75030.70900.7033MALLOWSC0.82350.85620.73350.76630.73550.7633CBIFC0.82210.85530.73220.76530.72230.7536BYC0.82610.87120.73610.78120.71630.7218WBYC0.82100.87430.73100.78430.71030.74386MLE0.88900.89230.78910.79240.78990.7229MALLOWSC0.82410.85240.72420.75250.74420.7215CBIFC0.82310.84240.71120.74250.71210.7254BYC0.85010.80160.75020.70170.75200.7071WBYC0.80910.83260.70920.73270.79200.727310MLE0.89810.90120.79800.80130.79880.7337MALLOWSC0.83520.86350.73530.76360.75530.7326CBIFC0.83420.85350.72230.75360.72320.7365BYC0.86120.81270.76130.71280.76310.7182WBYC0.81020.84370.71030.74380.80310.738410MLE0.91890.91230.80890.90310.88890.8373MALLOWSC0.92350.93560.82350.83560.83550.8236CBIFC0.92340.93550.83220.83650.82230.8256BYC0.91620.92710.81360.82810.81360.8281WBYC0.90210.93740.80310.83840.91030.843815MLE0.99820.92320.89810.93110.88990.8734MALLOWSC0.93530.95640.83530.85640.85540.8363CBIFC0.93430.95540.82240.86540.82330.8563BYC0.96220.97130.83620.88130.83620.8813WBYC0.92110.97440.83110.88440.90320.8384140MLE0.81620.84210.72720.73110.73820.7421MALLOWSC0.82160.85190.73260.76290.74360.7739CBIFC0.81860.84580.72860.75580.73960.7658BYC0.84540.90450.75640.80350.76740.8145WBYC0.84630.90960.75730.91070.76830.92172MLE0.82730.85320.73830.74220.74930.7532MALLOWSC0.83270.86200.74370.77300.75460.7840CBIFC0.82970.85690.73970.76690.74070.7769BYC0.85650.91560.76740.81460.77850.8256WBYC0.85740.91070.76840.92180.77940.93286MLE0.83730.83260.78340.72250.79350.7326MALLOWSC0.82740.82070.73750.73080.74660.7409CBIFC0.89730.86960.79740.76970.70750.7698BYC0.86560.95620.77470.84620.78580.853WBYC0.87460.90720.78470.91830.79480.928410MLE0.91620.94210.84740.83110.83820.8421MALLOWSC0.92160.95110.83260.86210.84350.8731CBIFC0.91860.94580.82870.85580.83160.8650BYC0.94540.92670.85630.90350.86740.9145WBYC0.94630.92180.85740.91070.86840.921715MLE0.92730.95320.85850.84220.84930.8532MALLOWSC0.93270.96220.84370.87320.85460.8842CBIFC0.92970.95690.83980.86690.84270.8761BYC0.95660.93680.86740.91460.87850.9256WBYC0.93520.91070.84630.90060.85730.9106Table 6
Bias andMSE of all methods for data with various percentages of bad leverage points (5%, 10%, and 20%).
κ% of misc errorMethodsN=200N=250N=300BiasMSEBiasMSEBiasMSE15MLE1.45702.19181.34612.20281.23502.1017MALLOWSC0.14000.30730.25110.41840.24100.4293CBIFC0.51480.44820.62590.55930.51590.4472BYC0.13440.37650.14550.36540.13440.3543WBYC0.17450.38420.16340.37320.15230.36212MLE1.66812.38291.55722.41391.44612.3128MALLOWSC0.35110.51840.46220.62950.45210.6384CBIFC0.72590.65930.83680.76840.72680.6583BYC0.34550.58740.36660.57650.34550.5654WBYC0.38560.59530.37450.58430.36340.47326MLE1.55712.27191.44622.30291.33512.2018MALLOWSC0.24010.40740.35120.51850.34110.5274CBIFC0.61490.54830.72580.65740.61580.5473BYC0.23450.47640.25560.46550.23450.4544WBYC0.27460.48430.26350.47330.25240.362210MLE1.55602.27081.44512.30181.33402.2007MALLOWSC0.24000.40630.35010.51740.34000.5263CBIFC0.61380.54720.72470.65630.61470.5462BYC0.23340.47530.25450.46440.23340.4533WBYC0.27350.48320.26240.47220.25130.361115MLE1.44502.16071.33412.20071.22302.1106MALLOWSC0.13100.30520.24000.40630.23100.4152CBIFC0.50270.53610.61360.54520.50360.4351BYC0.12330.36430.14340.35330.12230.3422WBYC0.16240.37210.15130.36110.14020.2500110MLE2.46486.10132.35386.09032.24286.1813MALLOWSC0.34840.20980.23740.19880.23740.1988CBIFC1.19331.45981.08221.34871.07111.2376BYC0.25420.17160.14310.06050.03200.1514WBYC0.05650.11200.14540.00100.03430.01202MLE2.48466.13102.38356.03092.28246.1318MALLOWSC0.45840.20980.34740.28880.34740.2888CBIFC1.28331.56981.17221.45871.18111.3476BYC0.36420.28160.25310.17050.14200.2614WBYC0.06750.12300.15640.01200.04530.02206MLE2.49566.14202.39456.04192.29346.1428MALLOWSC0.46940.21880.35840.29980.35840.2998CBIFC1.29431.57881.18321.46971.19211.3586BYC0.37520.29260.26410.18150.15300.2724WBYC0.07850.13400.16740.02300.05630.032010MLE1.58665.03101.40355.25292.18246.2538MALLOWSC0.35050.10990.46950.18090.46730.3009CBIFC1.30531.68981.29421.57071.20311.3586BYC0.48620.30360.37510.29250.26400.3834WBYC0.08960.14510.17850.03410.06740.043115MLE1.69775.14211.51465.36302.29356.3549MALLOWSC0.46160.21880.57060.29180.57840.4108CBIFC1.41641.79991.39531.68181.31421.4697BYC0.59730.41470.48610.38360.37510.4945WBYC0.19850.16620.19960.05520.08740.0541120MLE2.72887.47732.93087.69832.09187.7093MALLOWSC0.43090.64670.54100.75780.55210.7689CBIFC1.66032.80311.55122.71200.76231.9231BYC0.76881.32570.88761.57320.69870.4621WBYC0.07030.32170.08140.33280.09250.34392MLE2.61787.36632.82187.58732.18087.6183MALLOWSC0.32190.53570.43000.64680.44110.6579CBIFC1.54132.71211.43222.60100.65131.8121BYC0.65781.21470.77661.46220.58770.3511WBYC0.06120.31060.07030.32170.08140.33286MLE2.72887.47732.93287.69832.29187.7293MALLOWSC0.43290.64670.54100.75780.55210.7689CBIFC1.65232.82311.54322.71200.76231.9231BYC0.76881.32570.88761.57320.69870.4621WBYC0.07230.32170.08120.31060.07030.321710MLE2.72887.47732.93287.69832.29187.7293MALLOWSC0.43290.64670.54100.75780.55210.7689CBIFC1.65232.82311.54322.71200.76231.9231BYC0.76881.32570.88761.57320.69870.4621WBYC0.07230.32170.08120.31060.07030.321715MLE3.61786.36631.84186.58731.18086.6183MALLOWSC0.54390.75670.65200.86780.66310.8789CBIFC2.54121.71200.75431.82311.85342.8120BYC1.65771.63681.79852.48431.58760.6510WBYC0.16120.41060.17010.41050.16120.4106Table 7
Bias andMSE of all methods for data with various percentages of bad leverage points (30% and 40%).
κ% of misc errorMethodsN=200N=250N=300BiasMSEBiasMSEBiasMSE130MLE2.82478.00532.71378.00421.82489.1153MALLOWSC0.43540.56140.32430.45030.21320.3402CBIFC2.03184.16582.02074.05481.13185.1659BYC2.82588.01122.71478.00011.60367.9880WBYC0.37520.55600.36410.44510.25300.33402MLE2.71368.11422.72268.11311.93599.2242MALLOWSC0.54430.67030.43320.56020.32210.4501CBIFC2.02074.05472.13064.16371.24075.2748BYC2.93478.12012.82358.11001.71258.0970WBYC0.48410.66510.47300.55400.36200.44306MLE2.60258.00312.61158.00201.82489.1131MALLOWSC0.43320.56120.32210.45110.21100.3410CBIFC2.11164.14362.02154.05261.13165.1637BYC2.82368.01102.71248.00111.60148.0861WBYC0.37300.55400.36210.44300.25110.332110MLE1.82479.22533.83379.22422.05609.9352MALLOWSC0.65540.78340.54430.67330.43320.5632CBIFC2.33384.37582.24374.27481.35385.2859BYC2.04588.03322.73468.33221.82368.2983WBYC0.59520.77620.58430.66520.47330.554315MLE1.71379.11433.72279.11322.14509.8242MALLOWSC0.54440.67240.43330.56230.32220.4522CBIFC2.22284.26482.13274.16381.24285.1749BYC2.03478.02212.62368.22121.71268.1873WBYC0.48420.66520.47330.55420.36230.4433140MLE2.88388.53202.77288.42101.66177.3101MALLOWSC0.77160.88490.66050.77380.65140.7627CBIFC2.42875.93372.31765.82261.52876.9337BYC2.87718.31482.78808.40372.68708.3126WBYC0.05150.39610.04040.28700.13130.17612MLE3.99489.65303.88389.53202.77278.4211MALLOWSC0.66060.77390.55151.66280.54040.6517CBIFC3.31776.82273.42866.93362.41777.8227BYC2.87718.31482.78808.40372.79708.4216WBYC0.06260.30720.05150.29810.14240.18726MLE3.00399.76413.99499.64312.88388.5322MALLOWSC0.77070.88490.66260.77390.45150.7628CBIFC3.42886.93383.53976.94472.52877.9337BYC2.98818.42682.89918.51483.68617.3105WBYC0.17370.41830.16260.30920.25350.298310MLE2.11288.95313.10389.75202.99278.6411MALLOWSC0.88060.99380.55150.88280.56040.8717CBIFC3.53786.09273.64086.85362.63767.0426BYC2.08708.31792.78028.62593.79507.4204WBYC0.28260.52720.27150.41810.36240.317215MLE2.22398.96423.21499.76312.00178.6522MALLOWSC0.99160.90490.66250.89390.67140.8828CBIFC3.64886.00383.75186.86472.74867.0315BYC2.17608.32882.89128.63693.80607.4315WBYC0.39360.53830.38150.42920.47340.3283Tables2 and 3 show the results of data with misclassified errors. The bias and MSE of the MLE estimates are immediately affected by 5% misclassified type error. The results suggest that the MLE becomes biased with 5% contamination, CUBIFC with 20% contamination, and BYC with 30% contamination. MALLOWSC and WBYC are good robust estimators and considered as the best methods.As shown in Tables4 and 5, small differences are found between the classical MLE and the robust methods when contaminated data with good leverage points are used, where the biases and MSEs for all methods are relatively small. Hence, good leverage points do not affect the data.The results in the case of bad leverage points change the performance of estimators intensively. As shown in Tables6 and 7, the CUBIFC estimator can only withstand up to 5% contamination. The BYC estimator can tolerate up to 30% contamination. However, the performance of WBYC and MALLOWSC is equally good: they can withstand up to 40% contamination.


## 5.1. Settings

This simulation aims to compare the robustness of the proposed robust estimators and the classicalMLE. The independent circular variable is generated from von Mises distribution with mean 60° and concentration parameter of κ  = 1, 2, 6, 10, and 15, with sample size N=200,250, and 300; a large sample size is chosen to avoid separation problems. The true parameter values are β=0,2,2.The simulation study is reported in a variety of situations. Initially, the data without contamination are simulated. The robust properties of all estimators with contaminated data are examined in three different ways. First, proportions are taken from the responses, andv is chosen randomly and changed from either 0 to 1 or 1 to 0. This process constitutes the misclassification-type error. For each contaminated case, 5%, 10%,20%, 30%, and 40% of the original data are contaminated. Second, the same proportions are taken to contaminate u with γ=0.5 for good leverage points. Finally, the same proportions are considered, and the generated data are contaminated with two types of outliers simultaneously. This process constitutes bad leverage points.Each simulation includes 1000 replications. The performance of the estimators is evaluated based on the bias and the median squared error (MSE) for each parameter, which are defined as follows:(25)bias=11000∑i=11000β^i−βi,MSE=median∑i=11000β^i−β¯i2.A good estimator has bias andMSE that are relatively small or close to zero.The simulation used the standard available “Robust” package in R to obtain the estimators and the” CircStats” package for generating the circular variableu.


## 5.2. Results

The bias andMSE of the five estimators are shown in Tables 1 to 7. Table 1 shows the results for uncontaminated data (i.e., clean data), where the biases and MSEs of all five estimators are fairly close to each other. However, the MALLOWSC’s MSE is larger than the others for large concentration parameter (κ  = 6, 10, and 15). Hence, MALLOWSC performs worse compared with the other estimators in this situation.Table 1
Bias andMSE of all methods for uncontaminated data.
κMethodsN=200N=250N=300BiasMSEBiasMSEBiasMSE1MLE0.00490.01660.00230.01270.00250.0114MALLOWSC0.00290.01940.00170.01400.00280.0123CBIFC0.00380.01660.00210.01280.00220.0114BYC0.00510.01790.00280.01350.00280.0117WBYC0.00500.01790.00300.01350.00300.01182MLE0.00640.02130.00120.01661.2e-030.0131MALLOWSC0.01130.14960.00550.06012.3e-030.0406CBIFC0.00620.02330.00100.01771.6e-030.0136BYC0.00830.02200.00240.01711.0e-030.0136WBYC0.00570.13960.00100.01824.2e-050.01436MLE0.01700.22210.01100.08490.00660.0642MALLOWSC0.10990.87340.07670.34990.04170.2006CBIFC0.01970.23460.01760.09630.01050.0749BYC0.00600.22480.00430.09080.00330.0641WBYC0.03070.24020.00070.03160.04710.063210MLE0.02160.13920.01620.11550.01680.0938MALLOWSC0.06340.36000.06240.27690.06090.2530CBIFC0.03410.15290.03190.13110.01600.1030BYC0.03260.14190.01730.11820.01770.0948WBYC0.07640.12640.07190.10900.07520.087115MLE0.08950.25010.07060.19920.03300.1789MALLOWSC0.15181.02600.11580.48440.14140.4132CBIFC0.07840.27490.07160.21370.05290.1945BYC0.08590.25600.06260.20350.03240.1805WBYC0.09340.97540.08350.19490.08160.1661Table 2
Bias andMSE of all methods for data with various percentages of misclassified errors (5%, 10%, and 20%).
κ% of misc errorMethodsN=200N=250N=300BiasMSEBiasMSEBiasMSE15MLE1.45702.19181.55702.20181.65812.3018MALLOWSC0.14000.30730.25100.41620.13670.3612CBIFC0.51480.44820.46480.34720.46360.3371BYC0.13440.37650.12330.27560.11440.2565WBYC0.17450.38420.08650.39450.08410.29512MLE1.55702.29181.65702.30181.75812.4018MALLOWSC0.24000.40730.35100.51620.23670.4612CBIFC0.61480.54820.56480.44720.56360.4371BYC0.23440.47650.22330.37560.21440.3565WBYC0.27450.48420.18650.49450.18410.39516MLE1.05702.18181.05702.20181.05812.3018MALLOWSC0.14000.31730.25100.62620.12670.5512CBIFC0.50480.42820.47480.33720.45360.2271BYC0.02440.32650.01330.26560.01440.2465WBYC0.06450.34420.07650.38450.07410.285110MLE1.55802.20181.66802.31281.76912.4128MALLOWSC0.03100.21830.14000.30520.02570.2502CBIFC0.40380.33720.35380.23620.35260.2261BYC0.02340.26550.01230.16460.00340.1455WBYC0.06450.27320.07550.28450.07310.184115MLE1.56812.20291.66812.31291.76922.4129MALLOWSC0.13110.31620.24010.40510.12560.3501CBIFC0.50370.43710.45370.33610.45250.3260BYC0.12330.36540.11220.26450.10330.2454WBYC0.16340.37310.07540.38340.07300.2840110MLE2.46486.10132.35376.00022.25266.1113MALLOWSC0.34840.20980.23930.10870.12920.0076CBIFC0.89331.45980.99220.94970.98110.9386BYC0.25420.17160.14310.06050.03200.0504WBYC0.05650.11200.04540.01100.03430.01212MLE2.36486.00132.25376.01022.15266.0113MALLOWSC0.24840.10980.13930.00870.02920.0066CBIFC0.69330.75980.89220.84970.88110.8386BYC0.15420.07160.04310.05050.02200.0404WBYC0.04650.01200.03540.00100.02430.00216MLE2.35486.00132.24376.01122.14266.0013MALLOWSC0.23840.11980.12930.01870.01920.0176CBIFC0.9833.64980.88220.83970.87110.8286BYC0.14420.06160.03310.15050.12200.1404WBYC0.14650.00200.13540.10100.12430.102110MLE2.57586.21232.46476.11122.14166.2223MALLOWSC0.23740.11880.12830.00770.01820.0166CBIFC0.88230.54880.88120.83870.87010.8276BYC0.14320.06060.03210.15150.12100.1414WBYC0.04550.00300.03440.00000.02330.001115MLE2.57596.21242.46486.11132.36376.2224MALLOWSC0.23730.10870.12820.00760.01810.0065CBIFC0.68220.74870.88110.83860.87000.8275BYC0.14310.06050.03210.05040.02100.0403WBYC0.04540.00100.03430.00000.02320.0010120MLE2.72887.47732.61777.36622.50667.2551MALLOWSC0.43090.64670.32080.53560.21070.4245CBIFC1.66031.80311.55022.70201.44012.6010BYC0.76881.32570.65771.21460.54661.1035WBYC0.07030.32170.06020.21060.05010.10052MLE3.72888.47733.61778.36623.50668.2551MALLOWSC0.53090.74670.42080.63560.31070.5245CBIFC2.66033.80312.55023.70202.44013.6010BYC0.86882.32570.75772.21460.64662.1035WBYC0.17030.42170.16020.31060.15010.20056MLE2.83887.58732.72777.47622.61667.3651MALLOWSC0.52090.75670.43080.64560.32070.5345CBIFC1.77032.91311.66022.81201.55012.7110BYC0.87881.43570.76771.32460.65661.2135WBYC0.18030.43170.17020.21060.16010.210510MLE3.93888.68733.82778.57623.51668.4651MALLOWSC0.41090.64670.32080.53560.21070.4245CBIFC0.67031.81310.56021.71200.45011.6110BYC0.76880.53570.65770.22460.54660.1135WBYC0.07030.32170.06020.10060.05010.100515MLE4.04889.79734.93779.68624.62669.5751MALLOWSC0.40080.63560.31070.52450.20060.4134CBIFC0.66021.80200.55011.70100.44001.6000BYC0.75770.52460.64660.21350.53550.1024WBYC0.06020.31060.05010.10050.04000.1004Table 3
Bias andMSE of all methods for data with various percentages of misclassified errors (30% and 40%).
κ% of misc errorMethodsN=200N=250N=300BiasMSEBiasMSEBiasMSE130MLE2.82478.00533.72477.00432.61476.0032MALLOWSC0.43540.56140.32540.65140.21540.5414CBIFC2.03184.16581.02083.15481.01083.0448BYC2.82588.01122.71588.00122.60588.0001WBYC0.37520.55600.26420.44500.15320.33402MLE3.82479.00534.72478.00433.61477.0032MALLOWSC0.53540.66140.42540.75140.31540.6414CBIFC3.03185.16582.02084.15482.01084.0448BYC3.82589.01123.71589.00123.60589.0001WBYC0.47520.65600.36420.54500.25320.43406MLE3.93479.11534.83478.11433.72477.1132MALLOWSC0.54640.67240.43640.76240.32640.6524CBIFC3.04285.17682.03184.16582.02184.0558BYC3.83689.02223.72689.01223.61689.0000WBYC0.48620.66700.37520.55600.26420.445010MLE3.93589.11644.83588.11543.72587.1143MALLOWSC0.32430.45030.21430.54030.10430.4303CBIFC2.02074.05471.01073.04371.00073.0337BYC2.71478.00012.60478.00012.50478.0000WBYC0.26410.66500.15310.33400.04210.223015MLE3.82589.00644.72588.00543.61587.0043MALLOWSC0.21430.34030.10430.43030.00430.3203CBIFC2.13074.16471.12073.15371.11073.1437BYC2.82478.11012.71478.11012.61478.1100WBYC0.37410.77500.26310.44400.15210.3330140MLE2.88388.53202.77278.42102.66168.3100MALLOWSC0.77160.88490.76050.87380.75040.8627CBIFC2.42875.93372.31775.02272.20665.0116BYC2.87718.31482.76608.20372.65508.1026WBYC0.05150.39610.04040.28700.03140.29806MLE2.99388.64202.88278.53102.77168.4200MALLOWSC0.88160.99490.87050.98380.86040.9727CBIFC2.53875.04372.42775.13272.31665.1216BYC2.98718.42482.87608.31372.76508.2126WBYC0.16150.40710.05150.29810.04250.209110MLE3.00499.54213.98289.63113.87179.5201MALLOWSC0.99261.00300.98171.09490.97151.0838CBIFC2.87535.37042.77425.27132.66315.1612BYC2.71898.48422.60878.37312.50768.2621WBYC0.51610.71400.15050.81290.25040.201315MLE3.11509.65323.09399.74223.98289.6312MALLOWSC1.09371.11411.98171.10501.06261.9747CBIFC2.98635.48172.88535.38242.77415.2723BYC2.82908.59532.71988.48422.61878.3732WBYC0.62720.82510.26160.92300.36150.3124Table 4
Bias andMSE of all methods for data with various percentages of good leverage points (5%, 10%, and 20%).
κ% of misc errorMethodsN=200N=250N=300BiasMSEBiasMSEBiasMSE15MLE0.80720.81220.81620.82320.82720.8342MALLOWSC0.81150.81880.82250.82980.81050.8188CBIFC0.81220.82190.82320.83290.83420.8439BYC0.81380.83560.82480.84660.83580.8576WBYC0.82090.86230.83190.87330.84290.88432MLE0.87200.82210.86210.83220.87220.8423MALLOWSC0.81510.88810.82520.89820.80510.8881CBIFC0.82210.81920.83220.82930.84230.8394BYC0.83810.85630.84820.86640.85830.8765WBYC0.80920.82360.81930.83370.82940.84386MLE0.81720.92220.92620.93320.93720.9442MALLOWSC0.92150.92880.93250.93980.92050.9288CBIFC0.92220.93190.93320.94290.94420.9539BYC0.92380.94560.93480.95660.94580.9676WBYC0.93090.97230.94190.98330.95290.994310MLE0.87120.82220.86220.83320.87320.8442MALLOWSC0.81520.88280.82530.89380.80250.8828CBIFC0.82220.81390.83320.82490.84420.8359BYC0.83280.85460.84380.86560.85480.8766WBYC0.80390.82730.81490.83830.82590.849315MLE0.82170.82230.82260.82330.82370.8244MALLOWSC0.82510.82880.83520.88390.85200.8228CBIFC0.82230.89310.82330.89420.82440.8953BYC0.88230.86450.88340.85660.88450.8667WBYC0.89300.83720.89410.88330.89520.8394110MLE0.80350.80960.71350.71960.70250.7086MALLOWSC0.80830.81830.71830.72820.72910.7393CBIFC0.80600.81430.71600.72430.72700.7353BYC0.81210.84160.72210.73160.73310.7426WBYC0.81180.85060.72180.76060.71280.70662MLE0.83550.89660.73550.79660.72550.7866MALLOWSC0.88330.88230.78330.78220.79110.7933CBIFC0.86000.84330.76000.74330.77000.7533BYC0.82110.81660.72110.71660.73110.7266WBYC0.81880.80660.71880.70660.72880.76666MLE0.91530.91960.81530.81690.81520.8168MALLOWSC0.91380.92380.82380.82280.82190.8339CBIFC0.91060.91340.81060.82340.82070.8335BYC0.91120.94610.82120.83610.83130.8462WBYC0.91810.95600.82810.86600.81820.806510MLE0.81450.81070.72450.72070.71350.7196MALLOWSC0.81930.82930.72930.73920.73020.7404CBIFC0.81700.82530.72700.73530.73800.7463BYC0.82310.85260.73310.74260.74410.7536WBYC0.82280.86160.73280.77160.72380.717615MLE0.82540.82070.73540.73700.72530.7269MALLOWSC0.82390.83390.73390.74290.74200.7540CBIFC0.82070.83350.73070.74350.74080.7536BYC0.83130.86620.74130.75620.75140.7663WBYC0.83820.87610.74820.78610.73830.7267120MLE0.79100.79030.60910.60940.61020.6105MALLOWSC0.79540.79790.60640.60890.61750.6190CBIFC0.79110.79220.68110.68220.67220.6733BYC0.80190.78670.71190.69670.72200.6078WBYC0.78670.81010.67670.70010.66560.72232MLE0.79110.79220.68110.68220.67220.6733MALLOWSC0.80190.78670.71190.69670.72200.6078CBIFC0.79100.79030.60910.60940.61020.6105BYC0.79540.79790.60640.60890.61750.6190WBYC0.77860.80110.66770.70100.65660.72326MLE0.79100.79030.60910.60940.61020.6105MALLOWSC0.79540.79790.60640.60890.61750.6190CBIFC0.79110.79220.68110.68220.67220.6733BYC0.80190.78670.71190.69670.72200.6078WBYC0.78670.81010.67670.70010.66560.722310MLE0.80200.80130.71010.71040.72120.7215MALLOWSC0.80640.80890.71740.71990.72850.7200CBIFC0.80210.80320.79210.79320.78320.7843BYC0.91290.89770.82290.70770.83300.7188WBYC0.89770.92110.78770.81110.77660.833315MLE0.81310.81240.72120.72150.73230.7326MALLOWSC0.81750.81900.72850.72000.73960.7311CBIFC0.81320.81430.70320.70430.79430.7954BYC0.92310.80880.83300.71880.82210.7299WBYC0.80880.93220.79880.82110.78770.8433Table 5
Bias andMSE of all methods for data with various percentages of good leverage points (30% and 40%).
κ% of misc errorMethodsN=200N=250N=300BiasMSEBiasMSEBiasMSE130MLE0.80890.83920.71890.74920.79890.7922MALLOWSC0.81240.84520.72240.75520.72440.7522CBIFC0.81110.84420.72110.75420.71120.7425BYC0.81500.86010.72500.77010.70520.7107WBYC0.81090.86320.72090.77320.70920.73272MLE0.81900.84030.72900.75030.70900.7033MALLOWSC0.82350.85620.73350.76630.73550.7633CBIFC0.82210.85530.73220.76530.72230.7536BYC0.82610.87120.73610.78120.71630.7218WBYC0.82100.87430.73100.78430.71030.74386MLE0.88900.89230.78910.79240.78990.7229MALLOWSC0.82410.85240.72420.75250.74420.7215CBIFC0.82310.84240.71120.74250.71210.7254BYC0.85010.80160.75020.70170.75200.7071WBYC0.80910.83260.70920.73270.79200.727310MLE0.89810.90120.79800.80130.79880.7337MALLOWSC0.83520.86350.73530.76360.75530.7326CBIFC0.83420.85350.72230.75360.72320.7365BYC0.86120.81270.76130.71280.76310.7182WBYC0.81020.84370.71030.74380.80310.738410MLE0.91890.91230.80890.90310.88890.8373MALLOWSC0.92350.93560.82350.83560.83550.8236CBIFC0.92340.93550.83220.83650.82230.8256BYC0.91620.92710.81360.82810.81360.8281WBYC0.90210.93740.80310.83840.91030.843815MLE0.99820.92320.89810.93110.88990.8734MALLOWSC0.93530.95640.83530.85640.85540.8363CBIFC0.93430.95540.82240.86540.82330.8563BYC0.96220.97130.83620.88130.83620.8813WBYC0.92110.97440.83110.88440.90320.8384140MLE0.81620.84210.72720.73110.73820.7421MALLOWSC0.82160.85190.73260.76290.74360.7739CBIFC0.81860.84580.72860.75580.73960.7658BYC0.84540.90450.75640.80350.76740.8145WBYC0.84630.90960.75730.91070.76830.92172MLE0.82730.85320.73830.74220.74930.7532MALLOWSC0.83270.86200.74370.77300.75460.7840CBIFC0.82970.85690.73970.76690.74070.7769BYC0.85650.91560.76740.81460.77850.8256WBYC0.85740.91070.76840.92180.77940.93286MLE0.83730.83260.78340.72250.79350.7326MALLOWSC0.82740.82070.73750.73080.74660.7409CBIFC0.89730.86960.79740.76970.70750.7698BYC0.86560.95620.77470.84620.78580.853WBYC0.87460.90720.78470.91830.79480.928410MLE0.91620.94210.84740.83110.83820.8421MALLOWSC0.92160.95110.83260.86210.84350.8731CBIFC0.91860.94580.82870.85580.83160.8650BYC0.94540.92670.85630.90350.86740.9145WBYC0.94630.92180.85740.91070.86840.921715MLE0.92730.95320.85850.84220.84930.8532MALLOWSC0.93270.96220.84370.87320.85460.8842CBIFC0.92970.95690.83980.86690.84270.8761BYC0.95660.93680.86740.91460.87850.9256WBYC0.93520.91070.84630.90060.85730.9106Table 6
Bias andMSE of all methods for data with various percentages of bad leverage points (5%, 10%, and 20%).
κ% of misc errorMethodsN=200N=250N=300BiasMSEBiasMSEBiasMSE15MLE1.45702.19181.34612.20281.23502.1017MALLOWSC0.14000.30730.25110.41840.24100.4293CBIFC0.51480.44820.62590.55930.51590.4472BYC0.13440.37650.14550.36540.13440.3543WBYC0.17450.38420.16340.37320.15230.36212MLE1.66812.38291.55722.41391.44612.3128MALLOWSC0.35110.51840.46220.62950.45210.6384CBIFC0.72590.65930.83680.76840.72680.6583BYC0.34550.58740.36660.57650.34550.5654WBYC0.38560.59530.37450.58430.36340.47326MLE1.55712.27191.44622.30291.33512.2018MALLOWSC0.24010.40740.35120.51850.34110.5274CBIFC0.61490.54830.72580.65740.61580.5473BYC0.23450.47640.25560.46550.23450.4544WBYC0.27460.48430.26350.47330.25240.362210MLE1.55602.27081.44512.30181.33402.2007MALLOWSC0.24000.40630.35010.51740.34000.5263CBIFC0.61380.54720.72470.65630.61470.5462BYC0.23340.47530.25450.46440.23340.4533WBYC0.27350.48320.26240.47220.25130.361115MLE1.44502.16071.33412.20071.22302.1106MALLOWSC0.13100.30520.24000.40630.23100.4152CBIFC0.50270.53610.61360.54520.50360.4351BYC0.12330.36430.14340.35330.12230.3422WBYC0.16240.37210.15130.36110.14020.2500110MLE2.46486.10132.35386.09032.24286.1813MALLOWSC0.34840.20980.23740.19880.23740.1988CBIFC1.19331.45981.08221.34871.07111.2376BYC0.25420.17160.14310.06050.03200.1514WBYC0.05650.11200.14540.00100.03430.01202MLE2.48466.13102.38356.03092.28246.1318MALLOWSC0.45840.20980.34740.28880.34740.2888CBIFC1.28331.56981.17221.45871.18111.3476BYC0.36420.28160.25310.17050.14200.2614WBYC0.06750.12300.15640.01200.04530.02206MLE2.49566.14202.39456.04192.29346.1428MALLOWSC0.46940.21880.35840.29980.35840.2998CBIFC1.29431.57881.18321.46971.19211.3586BYC0.37520.29260.26410.18150.15300.2724WBYC0.07850.13400.16740.02300.05630.032010MLE1.58665.03101.40355.25292.18246.2538MALLOWSC0.35050.10990.46950.18090.46730.3009CBIFC1.30531.68981.29421.57071.20311.3586BYC0.48620.30360.37510.29250.26400.3834WBYC0.08960.14510.17850.03410.06740.043115MLE1.69775.14211.51465.36302.29356.3549MALLOWSC0.46160.21880.57060.29180.57840.4108CBIFC1.41641.79991.39531.68181.31421.4697BYC0.59730.41470.48610.38360.37510.4945WBYC0.19850.16620.19960.05520.08740.0541120MLE2.72887.47732.93087.69832.09187.7093MALLOWSC0.43090.64670.54100.75780.55210.7689CBIFC1.66032.80311.55122.71200.76231.9231BYC0.76881.32570.88761.57320.69870.4621WBYC0.07030.32170.08140.33280.09250.34392MLE2.61787.36632.82187.58732.18087.6183MALLOWSC0.32190.53570.43000.64680.44110.6579CBIFC1.54132.71211.43222.60100.65131.8121BYC0.65781.21470.77661.46220.58770.3511WBYC0.06120.31060.07030.32170.08140.33286MLE2.72887.47732.93287.69832.29187.7293MALLOWSC0.43290.64670.54100.75780.55210.7689CBIFC1.65232.82311.54322.71200.76231.9231BYC0.76881.32570.88761.57320.69870.4621WBYC0.07230.32170.08120.31060.07030.321710MLE2.72887.47732.93287.69832.29187.7293MALLOWSC0.43290.64670.54100.75780.55210.7689CBIFC1.65232.82311.54322.71200.76231.9231BYC0.76881.32570.88761.57320.69870.4621WBYC0.07230.32170.08120.31060.07030.321715MLE3.61786.36631.84186.58731.18086.6183MALLOWSC0.54390.75670.65200.86780.66310.8789CBIFC2.54121.71200.75431.82311.85342.8120BYC1.65771.63681.79852.48431.58760.6510WBYC0.16120.41060.17010.41050.16120.4106Table 7
Bias andMSE of all methods for data with various percentages of bad leverage points (30% and 40%).
κ% of misc errorMethodsN=200N=250N=300BiasMSEBiasMSEBiasMSE130MLE2.82478.00532.71378.00421.82489.1153MALLOWSC0.43540.56140.32430.45030.21320.3402CBIFC2.03184.16582.02074.05481.13185.1659BYC2.82588.01122.71478.00011.60367.9880WBYC0.37520.55600.36410.44510.25300.33402MLE2.71368.11422.72268.11311.93599.2242MALLOWSC0.54430.67030.43320.56020.32210.4501CBIFC2.02074.05472.13064.16371.24075.2748BYC2.93478.12012.82358.11001.71258.0970WBYC0.48410.66510.47300.55400.36200.44306MLE2.60258.00312.61158.00201.82489.1131MALLOWSC0.43320.56120.32210.45110.21100.3410CBIFC2.11164.14362.02154.05261.13165.1637BYC2.82368.01102.71248.00111.60148.0861WBYC0.37300.55400.36210.44300.25110.332110MLE1.82479.22533.83379.22422.05609.9352MALLOWSC0.65540.78340.54430.67330.43320.5632CBIFC2.33384.37582.24374.27481.35385.2859BYC2.04588.03322.73468.33221.82368.2983WBYC0.59520.77620.58430.66520.47330.554315MLE1.71379.11433.72279.11322.14509.8242MALLOWSC0.54440.67240.43330.56230.32220.4522CBIFC2.22284.26482.13274.16381.24285.1749BYC2.03478.02212.62368.22121.71268.1873WBYC0.48420.66520.47330.55420.36230.4433140MLE2.88388.53202.77288.42101.66177.3101MALLOWSC0.77160.88490.66050.77380.65140.7627CBIFC2.42875.93372.31765.82261.52876.9337BYC2.87718.31482.78808.40372.68708.3126WBYC0.05150.39610.04040.28700.13130.17612MLE3.99489.65303.88389.53202.77278.4211MALLOWSC0.66060.77390.55151.66280.54040.6517CBIFC3.31776.82273.42866.93362.41777.8227BYC2.87718.31482.78808.40372.79708.4216WBYC0.06260.30720.05150.29810.14240.18726MLE3.00399.76413.99499.64312.88388.5322MALLOWSC0.77070.88490.66260.77390.45150.7628CBIFC3.42886.93383.53976.94472.52877.9337BYC2.98818.42682.89918.51483.68617.3105WBYC0.17370.41830.16260.30920.25350.298310MLE2.11288.95313.10389.75202.99278.6411MALLOWSC0.88060.99380.55150.88280.56040.8717CBIFC3.53786.09273.64086.85362.63767.0426BYC2.08708.31792.78028.62593.79507.4204WBYC0.28260.52720.27150.41810.36240.317215MLE2.22398.96423.21499.76312.00178.6522MALLOWSC0.99160.90490.66250.89390.67140.8828CBIFC3.64886.00383.75186.86472.74867.0315BYC2.17608.32882.89128.63693.80607.4315WBYC0.39360.53830.38150.42920.47340.3283Tables2 and 3 show the results of data with misclassified errors. The bias and MSE of the MLE estimates are immediately affected by 5% misclassified type error. The results suggest that the MLE becomes biased with 5% contamination, CUBIFC with 20% contamination, and BYC with 30% contamination. MALLOWSC and WBYC are good robust estimators and considered as the best methods.As shown in Tables4 and 5, small differences are found between the classical MLE and the robust methods when contaminated data with good leverage points are used, where the biases and MSEs for all methods are relatively small. Hence, good leverage points do not affect the data.The results in the case of bad leverage points change the performance of estimators intensively. As shown in Tables6 and 7, the CUBIFC estimator can only withstand up to 5% contamination. The BYC estimator can tolerate up to 30% contamination. However, the performance of WBYC and MALLOWSC is equally good: they can withstand up to 40% contamination.


## 6. Application to Real Data

Two real datasets are considered to illustrate the performance of the proposed robust estimators. The results of the estimated coefficients with their standard errors (SEs) and p values are presented. A good estimator has small SE and the lower p values. The analyzed datasets are available from the corresponding author upon request.


### 6.1. Meteorological Data

A total of 365 daily observations of rainfall (yes, no) and wind directions are measured in degrees at 0.00 a.m. from January 1 to December 31, 2019. The data are obtained from the Palestinian Meteorological Department for the Hebron meteorology station in the West Bank, Palestine, with geographical coordinates (latitude =31°32 N and longitude =  35°06 W) and 1005 m above sea level. Forty-eight rainy days are reported in 2019 in Hebron governorate. Figure 1 shows the circular plot of wind direction data in Hebron.Figure 1
Circular plot of daily wind directions in Hebron, 2019.Figure2 suggests no outlier is found in the data (clean dataset). Abuzaid and ElShekh Ahmed [11] have not identified any outliers by their PMLE method and refereed that to the low concentration parameter of the wind direction (κ^=0.883).Figure 2
Scatter plot of meteorological data.Table8 shows the estimated parameters and estimated SE for the different estimators. As observed in Table 8, all the five estimators are fairly close to each other. Thus, the model of the relation between the rainfall as a binary response and the wind directions as a circular predictor of the considered dataset is given by(26)logitpi=logpi1−pi=−0.7894+0.2508cosui+0.1789sinui.Table 8
Meteorological data: estimated parameters,SE, and p values.
EstimationInterceptcoswindsinwindMethodEstimationSEp valueEstimationSEp valueEstimationSEp valueMLE−0.78940.1149≤0.0010.25080.16270.01230.17890.16200.0269MALLOWSC−0.78860.1149≤0.0010.25230.16270.01210.18020.16190.0266CBIFC−0.78940.1149≤0.0010.25080.16270.01230.17890.16200.0269BYC−0.78870.1148≤0.0010.25070.16850.01370.16840.15710.0284WBYC−0.78750.1149≤0.0010.25130.16300.01230.16660.16230.0305


### 6.2. Ecological Data

Leaf inclination angles are extremely important to model the light transmission through the vegetation canopies and to indirectly quantify the canopy attributes such as leaf area index and G-function. Two hundred measurements of leaf inclination angles ofBetula pendula trees were measured on May 9, 2013, and their canopies either (top or bottom) were recorded [22].Figure3 shows the circular plot of leaf inclination angles, and Figure 4 shows the scatter plot of the data. Six observations far from the mass of directions were identified as outliers by Abuzaid and ElShekh Ahmed [11]. Table 9 presents the parameter estimates, SE, and p values for the various robust estimators where the CUBIFC and MALLOWSC estimates are reasonably close to the MLE. However, BYC and WBYC have the smallest SE values. Therefore, BYC and WBYC estimators give the best results for this dataset. Thus, the model of the relation between the vegetation canopy (v) as a binary response and the leaf inclination angles (u) as a circular predictor of the considered dataset is given by(27)logitpi=logpi1−pi=−6.8431+6.5681cosui+4.5951sinui.Figure 3
Circular plot for leaf inclination angles (N=200).Figure 4
Scatter plot of ecological data.Table 9
Ecological data: estimated parameters, standard errors, andp values.
EstimationInterceptcosanglesinangleMethodEstimationSEp valueEstimationSEp valueEstimationSEp valueMLE−0.01000.25230.0944−0.01060.19990.0958−0.14720.20260.0468MALLOWSC−0.00630.25240.0964−0.01570.20000.0937−0.14280.20280.0481CBIFC−0.01000.25230.0943−0.01060.19990.0957−0.14720.20260.0467BYC−6.84310.1524≤0.0016.56810.1487≤0.0014.59510.13260.0200WBYC−6.81410.1423≤0.0016.53810.1484≤0.0014.56510.10260.0100


## 6.1. Meteorological Data

A total of 365 daily observations of rainfall (yes, no) and wind directions are measured in degrees at 0.00 a.m. from January 1 to December 31, 2019. The data are obtained from the Palestinian Meteorological Department for the Hebron meteorology station in the West Bank, Palestine, with geographical coordinates (latitude =31°32 N and longitude =  35°06 W) and 1005 m above sea level. Forty-eight rainy days are reported in 2019 in Hebron governorate. Figure 1 shows the circular plot of wind direction data in Hebron.Figure 1
Circular plot of daily wind directions in Hebron, 2019.Figure2 suggests no outlier is found in the data (clean dataset). Abuzaid and ElShekh Ahmed [11] have not identified any outliers by their PMLE method and refereed that to the low concentration parameter of the wind direction (κ^=0.883).Figure 2
Scatter plot of meteorological data.Table8 shows the estimated parameters and estimated SE for the different estimators. As observed in Table 8, all the five estimators are fairly close to each other. Thus, the model of the relation between the rainfall as a binary response and the wind directions as a circular predictor of the considered dataset is given by(26)logitpi=logpi1−pi=−0.7894+0.2508cosui+0.1789sinui.Table 8
Meteorological data: estimated parameters,SE, and p values.
EstimationInterceptcoswindsinwindMethodEstimationSEp valueEstimationSEp valueEstimationSEp valueMLE−0.78940.1149≤0.0010.25080.16270.01230.17890.16200.0269MALLOWSC−0.78860.1149≤0.0010.25230.16270.01210.18020.16190.0266CBIFC−0.78940.1149≤0.0010.25080.16270.01230.17890.16200.0269BYC−0.78870.1148≤0.0010.25070.16850.01370.16840.15710.0284WBYC−0.78750.1149≤0.0010.25130.16300.01230.16660.16230.0305


## 6.2. Ecological Data

Leaf inclination angles are extremely important to model the light transmission through the vegetation canopies and to indirectly quantify the canopy attributes such as leaf area index and G-function. Two hundred measurements of leaf inclination angles ofBetula pendula trees were measured on May 9, 2013, and their canopies either (top or bottom) were recorded [22].Figure3 shows the circular plot of leaf inclination angles, and Figure 4 shows the scatter plot of the data. Six observations far from the mass of directions were identified as outliers by Abuzaid and ElShekh Ahmed [11]. Table 9 presents the parameter estimates, SE, and p values for the various robust estimators where the CUBIFC and MALLOWSC estimates are reasonably close to the MLE. However, BYC and WBYC have the smallest SE values. Therefore, BYC and WBYC estimators give the best results for this dataset. Thus, the model of the relation between the vegetation canopy (v) as a binary response and the leaf inclination angles (u) as a circular predictor of the considered dataset is given by(27)logitpi=logpi1−pi=−6.8431+6.5681cosui+4.5951sinui.Figure 3
Circular plot for leaf inclination angles (N=200).Figure 4
Scatter plot of ecological data.Table 9
Ecological data: estimated parameters, standard errors, andp values.
EstimationInterceptcosanglesinangleMethodEstimationSEp valueEstimationSEp valueEstimationSEp valueMLE−0.01000.25230.0944−0.01060.19990.0958−0.14720.20260.0468MALLOWSC−0.00630.25240.0964−0.01570.20000.0937−0.14280.20280.0481CBIFC−0.01000.25230.0943−0.01060.19990.0957−0.14720.20260.0467BYC−6.84310.1524≤0.0016.56810.1487≤0.0014.59510.13260.0200WBYC−6.81410.1423≤0.0016.53810.1484≤0.0014.56510.10260.0100


## 7. Conclusion

This paper aimed to propose robust estimators for the circular logistic regression model. We have compared the performance of theMLE and the proposed robust estimators for the circular logistic model under clean and contaminated datasets. The findings indicate that the MLE tends to bias in the presence of misclassified error and bad leverage points. A good performance is obtained for the proposed robust estimator in estimating the parameters in such models. Some robust estimators, such as BYC and WBYC, show superiority over others depending on the type of contamination.



---

*Source: 9944363-2021-05-26.xml*