Submit manuscript...
eISSN: 2378-315X

Biometrics & Biostatistics International Journal

Research Article Volume 1 Issue 3

A more efficient nonparametric test of symmetry based on overlapping coefficient

Hani M Samawi, Robert Vogel

Jiann-Ping Hsu College of Public Health, Georgia Southern University, USA

Correspondence: Hani M Samawi, Jiann-Ping Hsu College of Public Health, Georgia Southern University, USA, Tel 912 478 1345

Received: November 14, 2014 | Published: December 16, 2014

Citation: Samawi HM, Vogel R. A more efficient nonparametric test of symmetry based on overlapping coefficient. Biom Biostat Int J. 2014;1(3):75-83. DOI: 10.15406/bbij.2014.01.00015

Download PDF

Abstract

In this paper we provide a more efficient nonparametric test of symmetry based on the empirical overlap coefficient using kernel density estimation applied to an extreme order statistics, namely extreme ranked set sampling. Our simulation investigation reveals that our proposed test of symmetry is at least as powerful as currently available tests of symmetry. Intensive simulation is conducted to examine the power of the proposed test. An illustration is provided using cardiac output and body weight of neonates in a neonatal intensive care unit.

Keywords: test of symmetry, power of the test, bootstrap method, overlap coefficients, weitzman’s measure, extreme ranked set sample, kernel density estimation, AMS, 62G10

Abbreviations

ERSS, extreme ranked set sample; SRS, simple random sample; RSS, ranked set sampling

Introduction

Parametric and some nonparametric statistical inferences and modeling are valid only under certain assumptions. One of most common assumptions in the literature is that of symmetry of the underlying distribution. If the underlying distribution is not symmetric the question becomes how to define the appropriate location and scale measures. Thus to choose the appropriate statistical analysis, we need to check for underlying assumptions, including symmetry. Most tests of symmetry available in the literature typically have low statistical power and fail to detect a small but meaningful asymmetry in the population. Examples of those tests have been suggested by Butler,1 Rothman and Woodroofe,2 Hill and Rao,3 McWilliams4 and Ozturk.5 McWilliams’s4 runs test of symmetry is more powerful than those provided by Butler,1 Rothman and Woodroofe,2 Hill and Rao3 against various asymmetric alternatives. Tajuddin,6 proposed a test for symmetry based on the Wilcoxon two-sample test and found his test to be more powerful than the runs test.

Baklizi7 suggested a runs test of symmetry based on the conditional distribution and demonstrated that it performed slightly better than the unconditional test by McWilliams.4 Baklizi‘s test is also very robust for misspecification of the median. Modarres and Gastwirth8-9 provided a modification to McWilliams4 runs test based on Wilcoxon scores to weigh the runs. Their procedure improved the power for testing symmetry when the center of the distribution is known. However, their test did not perform well when asymmetry is focused on regions close to the median. Samawi10 investigated the use of extreme ranked set sample (ERSS). Samawi et al.11 used (ERSS) to provide a more powerful runs test of symmetry. Finally, Samawi et al.12 used the overlap coefficient to test for symmetry and showed that their test procedure is competitive with the other available tests of symmetry. This paper uses ERSS to provide a more powerful overlap coefficient test of symmetry.

The overlap measure (OVL) is defined as the area of intersection of the graphs of two probability density functions. It measures the similarity, which is the agreement or the closeness of the two probability distributions. The OVL measure was originally introduced by Weitzman.13 Recently, several authors including Bradley and Piantadosi 14, Inman and Bradley 15, Clemons,16 Reiser and Faraggi,17 Clemons and Bradley,18 Mulekar and Mishra,19 Al-Saidy, et al.,20 Schmid and Schmidt,21 Al-Saleh and Samawi,22 and Samawi and Al-Saleh23 considered this measure. The sampling behavior of a nonparametric estimator of using naive kernel density estimation was examined by Clemons and Bradley,18 using Monte Carlo and bootstrap techniques.

Let f 1 ( x )  and  f 2 ( x ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzamaaBaaaleaacaaIXaaabeaakiaacIcacaWG4bGaaiykaiaabccacaqGHbGaaeOBaiaabsgacaqGGaGaamOzamaaBaaaleaacaaIYaaabeaakiaacIcacaWG4bGaaiykaaaa@446C@ be two probability density functions. Assume samples of observations are drawn from continuous distributions. The overlap measure used in the literature is defined by Weitzman’s Measure (1970) as Δ = min { f 1 ( x ) , f 2 ( x ) d x . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaeyypa0Zaa8qaaeaaciGGTbGaaiyAaiaac6gacaGG7bGaamOzamaaBaaaleaacaaIXaaabeaakiaacIcacaWG4bGaaiykaiaacYcacaWGMbWaaSbaaSqaaiaaikdaaeqaaOGaaiikaiaadIhacaGGPaaaleqabeqdcqGHRiI8aOGaaeyFaiaabccacaWGKbGaamiEaiaac6caaaa@4D97@

The overlap measure of two densities assumes values between 0 and 1. An overlap value close to 0 indicates extreme inequality of the two density functions, and an overlap value of 1 indicates exact equality. In most statistical applications the data used is assumed to consist of a simple random sample (SRS). Cost savings of quantifying sampling units can be achieved by using ranked set sampling (RSS) methods as described by McIntyre24 to estimate the population mean. The procedure introduced by McIntyre was later called RSS. As a variation of RSS, an extreme ranked set sample (ERSS) is introduced and investigated by Samawi et al.10

First we describe the RSS procedure as follows: First, randomly sample a group of sampling units, of size r 2 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOCamaaCaaaleqabaGaaGOmaaaaaaa@39E5@ from the target population. Randomly partition the group into disjoint subsets each having a pre-assigned size r. In most practical situations, the size r will be two, three or four. Rank the elements in each subset by a suitable method of ranking such as prior information, visual inspection or by the subject-matter experimenter himself. The ith order statistic from the i-th subset, X i ( i ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiwaSWaaSbaaeaacaWGPbGaamikaiaadMgacaWGPaaabeaaaaa@3C45@ , i = 1, …, r, will be quantified (actual measurement). Finally, X 1 ( 1 ) , X 2 ( 2 ) , ... , X r ( r ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiwaSWaaSbaaeaacaaIXaGaamikaiaaigdacaWGPaaabeaakiaadYcacaWGybWcdaWgaaqaaiaaikdacaWGOaGaaGOmaiaadMcaaeqaaOGaamilaiaad6cacaWGUaGaamOlaiaadYcacaWGybWcdaWgaaqaaiaadkhacaWGOaGaamOCaiaadMcaaeqaaaaa@484D@ constitutes the RSS. This will represents one complete cycle. The procedure can be repeated m-times as needed, to get a RSS of size n=Mr. A detailed explanation of uni variate RSS and its variations may be found in Kaur et al.,25 Patil et al.,26 Kaur et al.27 and Sinha 28. However, an extreme ranked set sample (ERSS) of size n=2m, is described by Samawi et al.,10 is similar to RSS procedure except that we quantify only the minima or the maxima from each subject to get an ERSS as X 1 ( 1 ) , X 2 ( 1 ) , ... , X m ( 1 ) ; X 1 ( r ) , X 2 ( r ) , ... , X m ( r ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiwaSWaaSbaaeaacaaIXaGaamikaiaaigdacaWGPaaabeaakiaadYcacaWGybWcdaWgaaqaaiaaikdacaWGOaGaaGymaiaadMcaaeqaaOGaamilaiaad6cacaWGUaGaamOlaiaadYcacaWGybWcdaWgaaqaaiaad2gacaWGOaGaaGymaiaadMcaaeqaaiaacUdakiaadIfalmaaBaaabaGaaGymaiaadIcacaWGYbGaamykaaqabaGccaWGSaGaamiwaSWaaSbaaeaacaaIYaGaamikaiaadkhacaWGPaaabeaakiaadYcacaWGUaGaamOlaiaad6cacaWGSaGaamiwaSWaaSbaaeaacaWGTbGaamikaiaadkhacaWGPaaabeaaaaa@598E@ .Now consider testing the null hypothesis of symmetry for an underlying absolutely continuous distribution F ( . ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOraiaacIcacaGGUaGaaiykaaaa@3ADB@ with density denoted by f ( . ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzaiaacIcacaGGUaGaaiykaaaa@3AFB@ : H 0 : f ( x ) = f ( x )   MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamisamaaBaaaleaacaaIWaaabeaakiaacQdacaWGMbGaaiikaiaadIhacaGGPaGaeyypa0JaamOzaiaacIcacqGHsislcaWG4bGaaiykaiaabccaaaa@4398@ versus  H a : f ( x ) f ( x ) ; MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeODaiaabwgacaqGYbGaae4CaiaabwhacaqGZbGaaeiiaiaadIeadaWgaaWcbaGaamyyaaqabaGccaGG6aGaamOzaiaacIcacaWG4bGaaiykaiabgcMi5kaadAgacaGGOaGaeyOeI0IaamiEaiaacMcacaGG7aaaaa@4AFE@ for some x. It is clear that under the null hypothesis of symmetry, if we let f 1 ( x ) = f ( x )  and   f 2 ( x ) = f ( x ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzamaaBaaaleaacaaIXaaabeaakiaacIcacaWG4bGaaiykaiabg2da9iaadAgacaGGOaGaamiEaiaacMcacaqGGaGaaeyyaiaab6gacaqGKbGaaeiiaiaabccacaWGMbWaaSbaaSqaaiaaikdaaeqaaOGaaiikaiaadIhacaGGPaGaeyypa0JaamOzaiaacIcacqGHsislcaWG4bGaaiykaaaa@4E8A@ then the overlap measure is equal to one ( Δ = 1 ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaiikaiabfs5aejabg2da9iaaigdacaGGPaaaaa@3C85@ , which will be our focus in this paper. Samawi et al. 12 used the overlap measure Δ to develop a new test of symmetry based on kernel density estimation of Δ. The availability of kernel density estimation in some of statistical software also, motivated us to use Δ when the sample in hand is ERSS. This paper will introduce a powerful test of symmetry based on ERSS overlap measure. The overlap test of symmetry using ERSS and its asymptotic properties are introduced is Section 2. A simulation study is given in Section 3. Illustrations using cardiac output data from neonates along with final comments are given in Section 4.

Test of Symmetry Based on ERSS Overlap Measure Let be X 1 ( 1 ) , X 2 ( 1 ) , ... , X m ( 1 ) ; X 1 ( r ) , X 2 ( r ) , ... , X m ( r ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiwaSWaaSbaaeaacaaIXaGaamikaiaaigdacaWGPaaabeaakiaadYcacaWGybWcdaWgaaqaaiaaikdacaWGOaGaaGymaiaadMcaaeqaaOGaamilaiaad6cacaWGUaGaamOlaiaadYcacaWGybWcdaWgaaqaaiaad2gacaWGOaGaaGymaiaadMcaaeqaaiaacUdakiaadIfalmaaBaaabaGaaGymaiaadIcacaWGYbGaamykaaqabaGccaWGSaGaamiwaSWaaSbaaeaacaaIYaGaamikaiaadkhacaWGPaaabeaakiaadYcacaWGUaGaamOlaiaad6cacaWGSaGaamiwaSWaaSbaaeaacaWGTbGaamikaiaadkhacaWGPaaabeaaaaa@598E@ an ERSS random sample from an absolutely continuous and differentiable distribution F ( . ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaajugGbiaadAeacaGGOaGaaiOlaiaacMcaaaa@3C32@ having known median. Without loss of generality, we will assume the median to be zero. When the median or the center of the distribution is unknown, the data can be centered by a consistent estimate of the median. The implications on the asymptotic properties resulting from centering the data on a consistent estimator of the median are not intuitively clear. Therefore, further investigations are needed to study the robustness of the proposed test of symmetry and compare it with other available tests of symmetry in case of an unknown median. In this paper we will discuss only the case when the median of the underlying distribution is assumed known. Consider the test for symmetry H 0 : f ( x ) = f ( x )   MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamisamaaBaaaleaacaaIWaaabeaakiaacQdacaWGMbGaaiikaiaadIhacaGGPaGaeyypa0JaamOzaiaacIcacqGHsislcaWG4bGaaiykaiaabccaaaa@4398@ versus  H a : f ( x ) f ( x ) ; for some  x . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeODaiaabwgacaqGYbGaae4CaiaabwhacaqGZbGaaeiiaiaadIeadaWgaaWcbaGaamyyaaqabaGccaGG6aGaamOzaiaacIcacaWG4bGaaiykaiabgcMi5kaadAgacaGGOaGaeyOeI0IaamiEaiaacMcacaGG7aGaaeOzaiaab+gacaqGYbGaaeiiaiaabohacaqGVbGaaeyBaiaabwgacaqGGaGaamiEaiaac6caaaa@5483@ Under the assumption of symmetry, F ( x ) = 1 F ( x )   MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOraiaacIcacaWG4bGaaiykaiabg2da9iaaigdacqGHsislcaWGgbGaaiikaiabgkHiTiaadIhacaGGPaGaaeiiaaaa@4285@ . Let f ( 1 ) ( x ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzamaaBaaaleaacaGGOaGaaGymaiaacMcaaeqaaOGaaiikaiaadIhacaGGPaaaaa@3D90@ be the density function of the first order statistics X ( 1 ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiwamaaBaaaleaacaGGOaGaaGymaiaacMcaaeqaaaaa@3B22@ and f ( r ) ( x ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzamaaBaaaleaacaGGOaGaamOCaiaacMcaaeqaaOGaaiikaiaadIhacaGGPaaaaa@3DCC@ be the density function of the rth order statistics X ( r ) 9 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiwamaaBaaaleaacaGGOaGaamOCaiaacMcaaeqaaOGaaGyoaaaa@3A18@ from random samples of size r respectively. Under the assumption of symmetry, it can be shown that f ( 1 ) ( x ) = f ( r ) ( x ) . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzamaaBaaaleaacaGGOaGaaGymaiaacMcaaeqaaOGaaiikaiabgkHiTiaadIhacaGGPaGaeyypa0JaamOzamaaBaaaleaacaGGOaGaamOCaiaacMcaaeqaaOGaaiikaiaadIhacaGGPaGaaiOlaaaa@45FC@ If we let f 1 ( x ) = f ( 1 ) ( x )  and   f 2 ( x ) = f ( r ) ( x ) , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzamaaBaaaleaacaaIXaaabeaakiaacIcacaWG4bGaaiykaiabg2da9iaadAgadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG4bGaaiykaiaabccacaqGHbGaaeOBaiaabsgacaqGGaGaaeiiaiaadAgadaWgaaWcbaGaaGOmaaqabaGccaGGOaGaamiEaiaacMcacqGH9aqpcaWGMbWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaamiEaiaacMcacaGGSaaaaa@540A@ the null hypothesis of symmetry is equivalent to H o : f ( 1 ) ( x ) = f ( r ) ( x ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamisamaaBaaaleaacaWGVbaabeaakiaacQdacaWGMbWaaSbaaSqaaiaacIcacaaIXaGaaiykaaqabaGccaGGOaGaeyOeI0IaamiEaiaacMcacqGH9aqpcaWGMbWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaamiEaiaacMcaaaa@47FF@ and under the null hypothesis Δ = 1. MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaeyypa0JaaGymaiaac6caaaa@3BDE@ Therefore, an equivalent hypothesis for testing the symmetry is H 0 : Δ = 1 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamisamaaBaaaleaacaaIWaaabeaakiaacQdacqqHuoarcqGH9aqpcaaIXaaaaa@3DA7@ versus  H a : Δ < 1. MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaeODaiaabwgacaqGYbGaae4CaiaabwhacaqGZbGaaeiiaiaadIeadaWgaaWcbaGaamyyaaqabaGccaGG6aGaeuiLdqKaeyipaWJaaGymaiaac6caaaa@44E0@ we propose using Δ ^ E R S S MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafuiLdqKbaKaadaWgaaWcbaGaamyraiaadkfacaWGtbGaam4uaaqabaaaaa@3CF8@ as it is a consistent nonparametric estimator of Δ using ERSS. Under the null hypothesis of symmetry and some mild regularity assumptions, which will be discussed later in this paper, we will derive the asymptotic distribution Δ ^ E R S S MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafuiLdqKbaKaadaWgaaWcbaGaamyraiaadkfacaWGtbGaam4uaaqabaaaaa@3CF8@ to use it as a test of symmetry, say: z 0 N ( 0 , 1 ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOEamaaBaaaleaacaaIWaaabeaakiablYJi6iaad6eacaGGOaGaaGimaiaacYcacaaIXaGaaiykaaaa@3F6E@ , (1)

for large n=2m. An asymptotic significant test procedure at level α MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaaccaGae8xSdegaaa@3796@ is to reject H 0 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamisamaaBaaaleaacaaIWaaabeaaaaa@39B8@ if z 0 < z α MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOEamaaBaaaleaacaaIWaaabeaakiabgYda8iabgkHiTiaadQhadaWgaaWcbaGaeqySdegabeaaaaa@3EAF@ , z α MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOEamaaBaaaleaacqaHXoqyaeqaaaaa@3ACF@ where is the upper α MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqySdegaaa@39A4@ percentile of the standard normal distribution.

Kernel estimation of ∆using ERSS

Based on the results of Schmid and Schmidt21 and Anderson et al.29 we will study the asymptotic properties using ERSS. Using one of the several available nonparametric density estimation procedures, see for example Wegman,30-31 Van Kerm32 and Chen and Kelton,33 one can use the overlap coefficient estimators for inferential purposes.

Let X 1 ( 1 ) , X 2 ( 1 ) , ... , X m ( 1 ) ; X 1 ( r ) , X 2 ( r ) , ... , X m ( r ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiwaSWaaSbaaeaacaaIXaGaamikaiaaigdacaWGPaaabeaakiaadYcacaWGybWcdaWgaaqaaiaaikdacaWGOaGaaGymaiaadMcaaeqaaOGaamilaiaad6cacaWGUaGaamOlaiaadYcacaWGybWcdaWgaaqaaiaad2gacaWGOaGaaGymaiaadMcaaeqaaiaacUdakiaadIfalmaaBaaabaGaaGymaiaadIcacaWGYbGaamykaaqabaGccaWGSaGaamiwaSWaaSbaaeaacaaIYaGaamikaiaadkhacaWGPaaabeaakiaadYcacaWGUaGaamOlaiaad6cacaWGSaGaamiwaSWaaSbaaeaacaWGTbGaamikaiaadkhacaWGPaaabeaaaaa@598E@ be ERSS random sample from a differentiable distribution F ( . ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaajugGbiaadAeacaGGOaGaaiOlaiaacMcaaaa@3C32@ having known median. Without a loss of generality assume the median to be zero. Let X 1 ( 1 ) , X 2 ( 1 ) , ... , X m ( 1 ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiwaSWaaSbaaeaacaaIXaGaamikaiaaigdacaWGPaaabeaakiaadYcacaWGybWcdaWgaaqaaiaaikdacaWGOaGaaGymaiaadMcaaeqaaOGaamilaiaad6cacaWGUaGaamOlaiaadYcacaWGybWcdaWgaaqaaiaad2gacaWGOaGaaGymaiaadMcaaeqaaaaa@480B@ denote a random sample of minimums in ERSS and X 1 ( r ) , X 2 ( r ) , ... , X m ( r ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiwaSWaaSbaaeaacaaIXaGaamikaiaadkhacaWGPaaabeaakiaadYcacaWGybWcdaWgaaqaaiaaikdacaWGOaGaamOCaiaadMcaaeqaaOGaamilaiaad6cacaWGUaGaamOlaiaadYcacaWGybWcdaWgaaqaaiaad2gacaWGOaGaamOCaiaadMcaaeqaaaaa@48BF@ denote a random sample of maximums in ERSS, where n = 2 m MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGUbGaeyypa0JaaGOmaiaad2gaaaa@3C5F@ . We will use a kernel function K that satisfies the condition
K ( x ) d x = 1 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWdXaqaaiaadUeacaGGOaGaamiEaiaacMcacaWGKbGaamiEaiabg2da9iaaigdaaSqaaiabgkHiTiabg6HiLcqaaiabg6HiLcqdcqGHRiI8aaaa@457D@ (2)

The kernel K is normally considered as a symmetric density function with mean 0 and finite variance; an example is the standard normal density. The kernel estimators of f ( 1 ) ( w i )  and  f ( r ) ( w i )    , i = 1 , 2 , ... , C MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGMbWaaSbaaSqaaiaacIcacaaIXaGaaiykaaqabaGccaGGOaGaeyOeI0Iaam4DamaaBaaaleaacaWGPbaabeaakiaacMcacaqGGaGaaeyyaiaab6gacaqGKbGaaeiiaiaadAgadaWgaaWcbaGaaiikaiaadkhacaGGPaaabeaakiaacIcacaWG3bWaaSbaaSqaaiaadMgaaeqaaOGaaiykaiaabccacaqGGaGaaiilaiaadMgacqGH9aqpcaaIXaGaaiilaiaaikdacaGGSaGaaiOlaiaac6cacaGGUaGaaiilaiaadoeaaaa@558E@ ,respectively are
f ^ ( 1 ) ( w i ) = 1 m h   j = 1 m K ( w i X j ( 1 ) h ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaaceWGMbGbaKaadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bWaaSbaaSqaaiaadMgaaeqaaOGaaiykaiabg2da9maalaaabaGaaGymaaqaaiaad2gacaWGObWaaSbaaSqaaiabgkHiTaqabaaaaOGaaeiiamaaqahabaGaam4samaabmaabaWaaSaaaeaacqGHsislcaWG3bWaaSbaaSqaaiaadMgaaeqaaOGaeyOeI0IaamiwamaaBaaaleaacaWGQbGaaiikaiaaigdacaGGPaaabeaaaOqaaiaadIgadaWgaaWcbaGaeyOeI0cabeaaaaaakiaawIcacaGLPaaaaSqaaiaadQgacqGH9aqpcaaIXaaabaGaamyBaaqdcqGHris5aaaa@5846@ (3)
and f ^ ( r ) ( w i ) = 1 m h +   j = 1 m K ( w i X j ( r ) h + ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaaceWGMbGbaKaadaWgaaWcbaGaaiikaiaadkhacaGGPaaabeaakiaacIcacaWG3bWaaSbaaSqaaiaadMgaaeqaaOGaaiykaiabg2da9maalaaabaGaaGymaaqaaiaad2gacaWGObWaaSbaaSqaaiabgUcaRaqabaaaaOGaaeiiamaaqahabaGaam4samaabmaabaWaaSaaaeaacaWG3bWaaSbaaSqaaiaadMgaaeqaaOGaeyOeI0IaamiwamaaBaaaleaacaWGQbGaaiikaiaadkhacaGGPaaabeaaaOqaaiaadIgadaWgaaWcbaGaey4kaScabeaaaaaakiaawIcacaGLPaaaaSqaaiaadQgacqGH9aqpcaaIXaaabaGaamyBaaqdcqGHris5aaaa@56CE@ (4) where C MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGdbaaaa@3980@ is the number of bins that depends on the sample size. In practice, we suggest to take C = integer of  m MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGdbGaeyypa0JaaeyAaiaab6gacaqG0bGaaeyzaiaabEgacaqGLbGaaeOCaiaabccacaqGVbGaaeOzaiaabccadaGcaaqaaiaad2gaaSqabaaaaa@4537@ (or the default setting of the software). Also, h  and  h + MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGObWaaSbaaSqaaiabgkHiTaqabaGccaqGGaGaaeyyaiaab6gacaqGKbGaaeiiaiaadIgadaWgaaWcbaGaey4kaScabeaaaaa@40C5@ are the bandwidths of the kernel estimators satisfying the conditions that h , h + > 0 , h , h + 0  and ( m h , m h + ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGObWaaSbaaSqaaiabgkHiTaqabaGccaGGSaGaamiAamaaBaaaleaacqGHRaWkaeqaaOGaeyOpa4JaaGimaiaacYcacaWGObWaaSbaaSqaaiabgkHiTaqabaGccaGGSaGaamiAamaaBaaaleaacqGHRaWkaeqaaOGaeyOKH4QaaGimaiaabccacaqGHbGaaeOBaiaabsgacaqGGaGaaeikaiaad2gacaWGObWaaSbaaSqaaiabgkHiTaqabaGccqGHsgIRcqGHEisPcaGGSaGaamyBaiaadIgadaWgaaWcbaGaey4kaScabeaakiabgkziUkabg6HiLkaacMcaaaa@5A1A@ as m . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGTbGaeyOKH4QaeyOhIuQaaiOlaaaa@3DBA@ There are many choices of the bandwidths ( h , h + ) , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaGGOaGaamiAamaaBaaaleaacqGHsislaeqaaOGaaiilaiaadIgadaWgaaWcbaGaey4kaScabeaakiaacMcacaGGSaaaaa@3F86@ however, in our procedure we use [34] Silverman’s [33] suggestion as follows: Using the normal distribution as the parametric family, the bandwidths of the kernel estimators are h = 0.9 A ( m ) 1 / 5  and  h + = 0.9 A + ( m ) 1 / 5 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGObWaaSbaaSqaaiabgkHiTaqabaGccqGH9aqpcaaIWaGaaiOlaiaaiMdacaWGbbWaaSbaaSqaaiabgkHiTaqabaGccaGGOaGaamyBaiaacMcadaahaaWcbeqaaiabgkHiTiaaigdacaGGVaGaaGynaaaakiaabccacaqGHbGaaeOBaiaabsgacaqGGaGaamiAamaaBaaaleaacqGHRaWkaeqaaOGaeyypa0JaaGimaiaac6cacaaI5aGaamyqamaaBaaaleaacqGHRaWkaeqaaOGaaiikaiaad2gacaGGPaWaaWbaaSqabeaacqGHsislcaaIXaGaai4laiaaiwdaaaaaaa@562E@ , (5)

where A MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGbbWaaSbaaSqaaiabgkHiTaqabaaaaa@3A97@ =min{standard deviation of ( X 1 ( 1 ) , X 2 ( 1 ) , ... , X m ( 1 ) ) , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaiikaiaadIfalmaaBaaabaGaaGymaiaadIcacaaIXaGaamykaaqabaGccaWGSaGaamiwaSWaaSbaaeaacaaIYaGaamikaiaaigdacaWGPaaabeaakiaadYcacaWGUaGaamOlaiaad6cacaWGSaGaamiwaSWaaSbaaeaacaWGTbGaamikaiaaigdacaWGPaaabeaakiaacMcacaGGSaaaaa@4A1E@ interquantile range of ( X 1 ( 1 ) , X 2 ( 1 ) , ... , X m ( 1 ) ) / 1.349 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaiikaiaadIfalmaaBaaabaGaaGymaiaadIcacaaIXaGaamykaaqabaGccaWGSaGaamiwaSWaaSbaaeaacaaIYaGaamikaiaaigdacaWGPaaabeaakiaadYcacaWGUaGaamOlaiaad6cacaWGSaGaamiwaSWaaSbaaeaacaWGTbGaamikaiaaigdacaWGPaaabeaakiaacMcacaGGVaGaaGymaiaac6cacaaIZaGaaGinaiaaiMdaaaa@4DCC@ } and A + MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGbbWaaSbaaSqaaiabgUcaRaqabaaaaa@3A8C@ =min{standard deviation of ( X 1 ( r ) , X 2 ( r ) , ... , X m ( r ) ) , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaiikaiaadIfalmaaBaaabaGaaGymaiaadIcacaWGYbGaamykaaqabaGccaWGSaGaamiwaSWaaSbaaeaacaaIYaGaamikaiaadkhacaWGPaaabeaakiaadYcacaWGUaGaamOlaiaad6cacaWGSaGaamiwaSWaaSbaaeaacaWGTbGaamikaiaadkhacaWGPaaabeaakiaacMcacaGGSaaaaa@4AD2@ interquantile range of ( ( X 1 ( r ) , X 2 ( r ) , ... , X m ( r ) ) / 1.349 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaiikaiaadIfalmaaBaaabaGaaGymaiaadIcacaWGYbGaamykaaqabaGccaWGSaGaamiwaSWaaSbaaeaacaaIYaGaamikaiaadkhacaWGPaaabeaakiaadYcacaWGUaGaamOlaiaad6cacaWGSaGaamiwaSWaaSbaaeaacaWGTbGaamikaiaadkhacaWGPaaabeaakiaacMcacaGGVaGaaGymaiaac6cacaaIZaGaaGinaiaaiMdaaaa@4E80@ }. These were found to be adequate choices of the bandwidth for minimizing the integrated mean squared error (IMSE),
I M S E = E [ f ^ ( i ) ( x ) f ( i ) ( x ) ] 2 d x ,  where  i = 1 , r . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGjbGaamytaiaadofacaWGfbGaeyypa0Zaa8qaaeaacaWGfbGaai4waiqadAgagaqcamaaBaaaleaacaGGOaGaamyAaiaacMcaaeqaaOGaaiikaiaadIhacaGGPaGaeyOeI0IaamOzamaaBaaaleaacaGGOaGaamyAaiaacMcaaeqaaOGaaiikaiaadIhacaGGPaGaaiyxamaaCaaaleqabaGaaGOmaaaakiaadsgacaWG4bGaaiilaiaabccacaqG3bGaaeiAaiaabwgacaqGYbGaaeyzaiaabccacaWGPbGaeyypa0JaaGymaiaacYcacaWGYbGaaiOlaaWcbeqab0Gaey4kIipaaaa@5C7F@ (6)

The bins used are as follows: Let R 1 = M a x ( R ( 1 ) , R ( r ) ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGsbWaaSbaaSqaaiaaigdaaeqaaOGaeyypa0JaamytaiaadggacaWG4bGaaiikaiaadkfadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacYcacaWGsbWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGPaaaaa@46C2@ , where R ( 1 ) = r a n g e ( X 1 ( 1 ) , X 2 ( 1 ) , ... , X m ( 1 ) ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGsbWaaSbaaSqaaiaacIcacaaIXaGaaiykaaqabaGccqGH9aqpcaWGYbGaamyyaiaad6gacaWGNbGaamyzaiaacIcacaWGybWcdaWgaaqaaiaaigdacaWGOaGaaGymaiaadMcaaeqaaOGaamilaiaadIfalmaaBaaabaGaaGOmaiaadIcacaaIXaGaamykaaqabaGccaWGSaGaamOlaiaad6cacaWGUaGaamilaiaadIfalmaaBaaabaGaamyBaiaadIcacaaIXaGaamykaaqabaGccaGGPaaaaa@52EE@ and R ( r ) = r a n g e ( X 1 ( r ) , X 2 ( r ) , ... , X m ( r ) ) . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGsbWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccqGH9aqpcaWGYbGaamyyaiaad6gacaWGNbGaamyzaiaacIcacaWGybWcdaWgaaqaaiaaigdacaWGOaGaamOCaiaadMcaaeqaaOGaamilaiaadIfalmaaBaaabaGaaGOmaiaadIcacaWGYbGaamykaaqabaGccaWGSaGaamOlaiaad6cacaWGUaGaamilaiaadIfalmaaBaaabaGaamyBaiaadIcacaWGYbGaamykaaqabaGccaGGPaGaaiOlaaaa@5490@ The bins will be selected as w i = w i 1 + δ x , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWG3bWaaSbaaSqaaiaadMgaaeqaaOGaeyypa0Jaam4DamaaBaaaleaacaWGPbGaeyOeI0IaaGymaaqabaGccqGHRaWkcqaH0oazdaWgaaWcbaGaamiEaaqabaGccaGGSaaaaa@4410@ where i = 2 , ... , C MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGPbGaeyypa0JaaGOmaiaacYcacaGGUaGaaiOlaiaac6cacaGGSaGaam4qaaaa@3FA6@ , w 1 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWG3bWaaSbaaSqaaiaaigdaaeqaaaaa@3A9B@ is an initial value chosen based on the minimum value used in R’s calculation and δ x = R 1 C . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH0oazdaWgaaWcbaGaamiEaaqabaGccqGH9aqpdaWcaaqaaiaadkfadaWgaaWcbaGaaGymaaqabaaakeaacaWGdbaaaiaac6caaaa@3FE8@

Using the aforementioned kernel estimator the nonparametric kernel estimator of Δ MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaaccaGae8hLdqeaaa@396F@ is given by
Δ ^ E R S S = min ( f ^ ( 1 ) ( w i ) , f ^ ( r ) ( w i ) ) d w , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafuiLdqKbaKaadaWgaaWcbaGaamyraiaadkfacaWGtbGaam4uaaqabaGccqGH9aqpdaWdraqaaiGac2gacaGGPbGaaiOBaiaacIcaceWGMbGbaKaadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bWaaSbaaSqaaiaadMgaaeqaaOGaaiykaiaacYcaceWGMbGbaKaadaWgaaWcbaGaaiikaiaadkhacaGGPaaabeaakiaacIcacaWG3bWaaSbaaSqaaiaadMgaaeqaaOGaaiykaiaacMcacaWGKbGaam4DaaWcbaGaeSyhHekabeqdcqGHRiI8aOGaaiilaaaa@57A1@ (7)

which can be approximated by a trapezoidal rule, resulting in
Δ ^ E R S S i = 1 C δ x 2 [ min ( f ^ ( 1 ) ( w i ) , f ^ ( r ) ( w i ) ) + min ( f ^ ( 1 ) ( w i 1 ) , f ^ ( r ) ( w i 1 ) ) ] . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafuiLdqKbaKaadaWgaaWcbaGaamyraiaadkfacaWGtbGaam4uaaqabaGccqGHijYUdaaeWbqaamaalaaabaGaeqiTdq2aaSbaaSqaaiaadIhaaeqaaaGcbaGaaGOmaaaadaWadaqaaiGac2gacaGGPbGaaiOBaiaacIcaceWGMbGbaKaadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bWaaSbaaSqaaiaadMgaaeqaaOGaaiykaiaacYcaceWGMbGbaKaadaWgaaWcbaGaaiikaiaadkhacaGGPaaabeaakiaacIcacaWG3bWaaSbaaSqaaiaadMgaaeqaaOGaaiykaiaacMcacqGHRaWkciGGTbGaaiyAaiaac6gacaGGOaGabmOzayaajaWaaSbaaSqaaiaacIcacaaIXaGaaiykaaqabaGccaGGOaGaeyOeI0Iaam4DamaaBaaaleaacaWGPbGaeyOeI0IaaGymaaqabaGccaGGPaGaaiilaiqadAgagaqcamaaBaaaleaacaGGOaGaamOCaiaacMcaaeqaaOGaaiikaiaadEhadaWgaaWcbaGaamyAaiabgkHiTiaaigdaaeqaaOGaaiykaiaacMcaaiaawUfacaGLDbaacaGGUaaaleaacaWGPbGaeyypa0JaaGymaaqaaiaadoeaa0GaeyyeIuoaaaa@75E6@ Asymptotic properties of Δ ^ E R S S MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafuiLdqKbaKaadaWgaaWcbaGaamyraiaadkfacaWGtbGaam4uaaqabaaaaa@3CF8@

The nonparametric kernel estimator of Δ MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqeaaa@396B@ ( Δ ^ E R S S ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaaiikaiqbfs5aezaajaWaaSbaaSqaaiaadweacaWGsbGaam4uaiaadofaaeqaaOGaaiykaaaa@3E5B@ is based on the univariate kernel for density estimation, K : MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGlbGaaiOoaiabl2riHkabgkziUkabl2riHcaa@3F13@ . Some of the necessary regularity conditions imposed on the univariate kernel for density estimation sees for example Silverman,34 Wand and Jones35 and Schmid and Schmidt21 are stated below:
1. K ( z ) d z = 1. 9 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWdraqaaiaadUeacaGGOaGaamOEaiaacMcacaWGKbGaamOEaiabg2da9iaaigdacaGGUaaaleaacqWIDesOaeqaniabgUIiYdGccaaI5aaaaa@4482@
2. z β K ( z ) d z = 0  for any  β = 1 , ... ,  d 1 , and  | z | d K ( z ) d z < , 1 < d . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaadaWdraqaaiaadQhadaahaaWcbeqaaiabek7aIbaakiaadUeacaGGOaGaamOEaiaacMcacaWGKbGaamOEaiabg2da9iaaicdacaqGGaGaaeOzaiaab+gacaqGYbGaaeiiaiaabggacaqGUbGaaeyEaiaabccacqaHYoGycqGH9aqpcaaIXaGaaiilaiaac6cacaGGUaGaaiOlaiaacYcacaqGGaGaaeizaaWcbaGaeSyhHekabeqdcqGHRiI8aOGaeyOeI0IaaGymaiaabYcacaqGGaGaaeyyaiaab6gacaqGKbGaaeiiamaapebabaGaaiiFaiaadQhacaGG8bWaaWbaaSqabeaacaWGKbaaaOGaam4saiaacIcacaWG6bGaaiykaiaadsgacaWG6bGaeyipaWJaeyOhIuQaaiilaiaaigdacqGH8aapcaWGKbGaeyicI4SaeSyfHukaleaacqWIDesOaeqaniabgUIiYdGccaGGUaaaaa@724A@
3. R = K 2 ( z ) d z < . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGsbGaeyypa0Zaa8qeaeaacaWGlbWaaWbaaSqabeaacaaIYaaaaOGaaiikaiaadQhacaGGPaGaamizaiaadQhacqGH8aapcqGHEisPcaGGUaaaleaacqWIDesOaeqaniabgUIiYdaaaa@4739@
4. h , h + > 0 , h , h + 0  , ( m h , m h + )  and ( m h log m , m h + log m ) . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGObWaaSbaaSqaaiabgkHiTaqabaGccaGGSaGaamiAamaaBaaaleaacqGHRaWkaeqaaOGaeyOpa4JaaGimaiaacYcacaWGObWaaSbaaSqaaiabgkHiTaqabaGccaGGSaGaamiAamaaBaaaleaacqGHRaWkaeqaaOGaeyOKH4QaaGimaiaabccacaqGSaGaaeiiaiaabIcacaWGTbGaamiAamaaBaaaleaacqGHsislaeqaaOGaeyOKH4QaeyOhIuQaaiilaiaad2gacaWGObWaaSbaaSqaaiabgUcaRaqabaGccqGHsgIRcqGHEisPcaGGPaGaaeiiaiaabggacaqGUbGaaeizaiaabccacaqGOaWaaSaaaeaacaWGTbGaamiAamaaBaaaleaacqGHsislaeqaaaGcbaGaciiBaiaac+gacaGGNbGaamyBaaaacqGHsgIRcqGHEisPcaGGSaWaaSaaaeaacaWGTbGaamiAamaaBaaaleaacqGHRaWkaeqaaaGcbaGaciiBaiaac+gacaGGNbGaamyBaaaacqGHsgIRcqGHEisPcaGGPaGaaiOlaaaa@7322@

To show consistency of Δ ^ E R S S MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafuiLdqKbaKaadaWgaaWcbaGaamyraiaadkfacaWGtbGaam4uaaqabaaaaa@3CF8@ , we use some of the kernel density asymptotic properties from Silverman33 and Wand and Jones34 under simple random sample (SRS). Under the assumptions 1-4 and assuming that the density f : MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGMbGaaiOoaiabl2riHkabgkziUkabl2riHcaa@3F2E@ is continuous at each w i MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWG3bWaaSbaaSqaaiaadMgaaeqaaaaa@3ACE@ , i=1, 2, …, C , and F(x) is absolutely continuous. Note that f ( 1 ) ( x ) = r [ 1 F ( x ) ] r 1 f ( x )  and  f ( r ) ( x ) = r [ F ( x ) ] r 1 f ( x ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzamaaBaaaleaacaGGOaGaaGymaiaacMcaaeqaaOGaaiikaiaadIhacaGGPaGaeyypa0JaamOCaiaacUfacaaIXaGaeyOeI0IaamOraiaacIcacaWG4bGaaiykaiaac2fadaahaaWcbeqaaiaadkhacqGHsislcaaIXaaaaOGaamOzaiaacIcacaWG4bGaaiykaiaabccacaqGHbGaaeOBaiaabsgacaqGGaGaamOzamaaBaaaleaacaGGOaGaamOCaiaacMcaaeqaaOGaaiikaiaadIhacaGGPaGaeyypa0JaamOCaiaacUfacaWGgbGaaiikaiaadIhacaGGPaGaaiyxamaaCaaaleqabaGaamOCaiabgkHiTiaaigdaaaGccaWGMbGaaiikaiaadIhacaGGPaaaaa@62EB@ are also absolutely continuous, therefore, the following apply:
B i a s ( f ^ ( 1 ) ( w i ) ) = o ( 1 )  and   B i a s ( f ^ ( r ) ( w i ) ) = o ( 1 ) + MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGcbGaamyAaiaadggacaWGZbGaaiikaiqadAgagaqcamaaBaaaleaacaGGOaGaaGymaiaacMcaaeqaaOGaaiikaiabgkHiTiaadEhadaWgaaWcbaGaamyAaaqabaGccaGGPaGaaiykaiabg2da9iaad+gacaGGOaGaaGymaiaacMcadaWgaaWcbaGaeyOeI0cabeaakiaabccacaqGHbGaaeOBaiaabsgacaqGGaGaaeiiaiaadkeacaWGPbGaamyyaiaadohacaGGOaGabmOzayaajaWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaam4DamaaBaaaleaacaWGPbaabeaakiaacMcacaGGPaGaeyypa0Jaam4BaiaacIcacaaIXaGaaiykamaaBaaaleaacqGHRaWkaeqaaaaa@6027@ (8)
V a r ( f ^ ( 1 ) ( w i ) ) = f ( 1 ) ( w i ) m h K 2 ( z ) d z + o ( 1 m h ) and   V a r ( f ^ ( r ) ( w i ) ) = f ( r ) ( w i ) m h + K 2 ( z ) d z + o ( 1 m h + ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakqaabeqaaiaadAfacaWGHbGaamOCaiaacIcaceWGMbGbaKaadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bWaaSbaaSqaaiaadMgaaeqaaOGaaiykaiaacMcacqGH9aqpdaWcaaqaaiaadAgadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bWaaSbaaSqaaiaadMgaaeqaaOGaaiykaaqaaiaad2gacaWGObWaaSbaaSqaaiabgkHiTaqabaaaaOWaa8qeaeaacaWGlbWaaWbaaSqabeaacaaIYaaaaOGaaiikaiaadQhacaGGPaGaamizaiaadQhacqGHRaWkcaWGVbGaaiikamaalaaabaGaaGymaaqaaiaad2gacaWGObWaaSbaaSqaaiabgkHiTaqabaaaaaqaaiabl2riHcqab0Gaey4kIipakiaabMcacaqGGaGaaeyyaiaab6gacaqGKbGaaeiiaiaabccaaeaacaWGwbGaamyyaiaadkhacaGGOaGabmOzayaajaWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaam4DamaaBaaaleaacaWGPbaabeaakiaacMcacaGGPaGaeyypa0ZaaSaaaeaacaWGMbWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaam4DamaaBaaaleaacaWGPbaabeaakiaacMcaaeaacaWGTbGaamiAamaaBaaaleaacqGHRaWkaeqaaaaakmaapebabaGaam4samaaCaaaleqabaGaaGOmaaaakiaacIcacaWG6bGaaiykaiaadsgacaWG6bGaey4kaSIaam4BaiaacIcadaWcaaqaaiaaigdaaeaacaWGTbGaamiAamaaBaaaleaacqGHRaWkaeqaaaaaaeaacqWIDesOaeqaniabgUIiYdGccaqGPaGaaeOlaiaabccaaaaa@8CF8@ (9) Since both variances converge to zero and for h , h + > 0 , h , h + 0  and ( m h , m h + ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGObWaaSbaaSqaaiabgkHiTaqabaGccaGGSaGaamiAamaaBaaaleaacqGHRaWkaeqaaOGaeyOpa4JaaGimaiaacYcacaWGObWaaSbaaSqaaiabgkHiTaqabaGccaGGSaGaamiAamaaBaaaleaacqGHRaWkaeqaaOGaeyOKH4QaaGimaiaabccacaqGHbGaaeOBaiaabsgacaqGGaGaaeikaiaad2gacaWGObWaaSbaaSqaaiabgkHiTaqabaGccqGHsgIRcqGHEisPcaGGSaGaamyBaiaadIgadaWgaaWcbaGaey4kaScabeaakiabgkziUkabg6HiLkaacMcaaaa@5A1A@ as m MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGTbGaeyOKH4QaeyOhIukaaa@3D08@ , then f ^ ( 1 ) ( w i ) P f ( 1 ) ( w i )  and  f ^ ( r ) ( w i ) P f ( r ) ( w i ) , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaaceWGMbGbaKaadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bWaaSbaaSqaaiaadMgaaeqaaOGaaiykaiabgkziUoaaCaaaleqabaGaamiuaaaakiaadAgadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bWaaSbaaSqaaiaadMgaaeqaaOGaaiykaiaabccacaqGHbGaaeOBaiaabsgacaqGGaGabmOzayaajaWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaam4DamaaBaaaleaacaWGPbaabeaakiaacMcacqGHsgIRdaahaaWcbeqaaiaadcfaaaGccaWGMbWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaam4DamaaBaaaleaacaWGPbaabeaakiaacMcacaGGSaaaaa@6086@ for all continuous points, w i MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWG3bWaaSbaaSqaaiaadMgaaeqaaaaa@3ACE@ of f ( 1 )  and  f ( r ) . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzamaaBaaaleaacaGGOaGaaGymaiaacMcaaeqaaOGaaeiiaiaabggacaqGUbGaaeizaiaabccacaWGMbWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGUaaaaa@435F@ Also, if f ( i ) ( . )   , i = 1 , r MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVG0dg9vqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGMbWaaSbaaSqaaiaacIcacaWGPbGaaiykaaqabaGccaGGOaGaaiOlaiaacMcacaqGGaGaaiilaiaadMgacqGH9aqpcaaIXaGaaiilaiaadkhaaaa@412B@ are uniformly continuous, then the kernel densities estimate are strongly consistent. As in Samawi et al.12 we can redefine Δ MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqeaaa@3758@ as follows: for any two numbers a and b min ( a , b ) = a + b 2 | a b | 2 , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaciyBaiaacMgacaGGUbGaaiikaiaadggacaGGSaGaamOyaiaacMcacqGH9aqpdaWcaaqaaiaadggacqGHRaWkcaWGIbaabaGaaGOmaaaacqGHsisldaWcaaqaaiaacYhacaWGHbGaeyOeI0IaamOyaiaacYhaaeaacaaIYaaaaiaacYcaaaa@4A51@ then Δ = 1 1 2 | f ( 1 ) ( x ) f ( r ) ( x ) d x . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqKaeyypa0JaaGymaiabgkHiTmaalaaabaGaaGymaaqaaiaaikdaaaWaa8qaaeaacaGG8bGaamOzamaaBaaaleaacaGGOaGaaGymaiaacMcaaeqaaOGaaiikaiabgkHiTiaadIhacaGGPaGaeyOeI0IaamOzamaaBaaaleaacaGGOaGaamOCaiaacMcaaeqaaOGaaiikaiaadIhacaGGPaaaleqabeqdcqGHRiI8aOGaaeiFaiaabccacaWGKbGaamiEaiaac6caaaa@520B@ Thus Δ ^ E R S S i = 1 C δ x 2 [ min ( f ^ ( 1 ) ( w i ) , f ^ ( r ) ( w i ) ) + min ( f ^ ( 1 ) ( w i 1 ) , f ^ ( r ) ( w i 1 ) ) ] , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafyiLdqKbaKaadaWgaaWcbaGaamyraiaadkfacaWGtbGaam4uaaqabaGccqGHijYUdaaeWbqaamaalaaabaGaeqiTdq2aaSbaaSqaaiaadIhaaeqaaaGcbaGaaGOmaaaadaWadaqaaiGac2gacaGGPbGaaiOBaiaacIcaceWGMbGbaKaadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bWaaSbaaSqaaiaadMgaaeqaaOGaaiykaiaacYcaceWGMbGbaKaadaWgaaWcbaGaaiikaiaadkhacaGGPaaabeaakiaacIcacaWG3bWaaSbaaSqaaiaadMgaaeqaaOGaaiykaiaacMcacqGHRaWkciGGTbGaaiyAaiaac6gacaGGOaGabmOzayaajaWaaSbaaSqaaiaacIcacaaIXaGaaiykaaqabaGccaGGOaGaeyOeI0Iaam4DamaaBaaaleaacaWGPbGaeyOeI0IaaGymaaqabaGccaGGPaGaaiilaiqadAgagaqcamaaBaaaleaacaGGOaGaamOCaiaacMcaaeqaaOGaaiikaiaadEhadaWgaaWcbaGaamyAaiabgkHiTiaaigdaaeqaaOGaaiykaiaacMcaaiaawUfacaGLDbaacaGGSaaaleaacaWGPbGaeyypa0JaaGymaaqaaiaadoeaa0GaeyyeIuoaaaa@75E5@ can be written as
Δ ^ E R S S = 1 1 2 | f ^ K ( 1 ) ( w ) f ^ K ( r ) ( w ) | d w      1 1 2 i = 1 C δ x 2 [ | f ^ ( 1 ) ( w i ) f ^ ( r ) ( w i ) | + | f ^ ( 1 ) ( w i 1 ) f ^ ( r ) ( w i 1 ) |   ] . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGceaqabeaacuqHuoargaqcamaaBaaaleaacaWGfbGaamOuaiaadofacaWGtbaabeaakiabg2da9iaaigdacqGHsisldaWcaaqaaiaaigdaaeaacaaIYaaaamaapeaabaGaaiiFaiqadAgagaqcamaaBaaaleaacaWGlbGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bGaaiykaiabgkHiTiqadAgagaqcamaaBaaaleaacaWGlbGaaiikaiaadkhacaGGPaaabeaakiaacIcacaWG3bGaaiykaiaacYhacaWGKbGaam4DaaWcbeqab0Gaey4kIipaaOqaaiaabccacaqGGaGaeyisISRaaeiiaiaaigdacqGHsisldaWcaaqaaiaaigdaaeaacaaIYaaaamaaqahabaWaaSaaaeaacqaH0oazdaWgaaWcbaGaamiEaaqabaaakeaacaaIYaaaamaadmaabaGaaiiFaiqadAgagaqcamaaBaaaleaacaGGOaGaaGymaiaacMcaaeqaaOGaaiikaiabgkHiTiaadEhadaWgaaWcbaGaamyAaaqabaGccaGGPaGaeyOeI0IabmOzayaajaWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaam4DamaaBaaaleaacaWGPbaabeaakiaacMcacaGG8bGaey4kaSIaaiiFaiqadAgagaqcamaaBaaaleaacaGGOaGaaGymaiaacMcaaeqaaOGaaiikaiabgkHiTiaadEhadaWgaaWcbaGaamyAaiabgkHiTiaaigdaaeqaaOGaaiykaiabgkHiTiqadAgagaqcamaaBaaaleaacaGGOaGaamOCaiaacMcaaeqaaOGaaiikaiaadEhadaWgaaWcbaGaamyAaiabgkHiTiaaigdaaeqaaOGaaiykaiaacYhacaqGGaaacaGLBbGaayzxaaaaleaacaWGPbGaeyypa0JaaGymaaqaaiaadoeaa0GaeyyeIuoakiaac6caaaaa@90DF@

Using the above results and under the null hypothesis of symmetry, i.e. f ( 1 ) ( w i ) = f ( r ) ( w i ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGMbWaaSbaaSqaaiaacIcacaaIXaGaaiykaaqabaGccaGGOaGaeyOeI0Iaam4DamaaBaaaleaacaWGPbaabeaakiaacMcacqGH9aqpcaWGMbWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaam4DamaaBaaaleaacaWGPbaabeaakiaacMcaaaa@4843@ , then,
| f ^ ( 1 ) ( w ) f ^ ( r ) ( w ) | = | ( f ^ ( 1 ) ( w ) f ( 1 ) ( w ) ) + ( f ( r ) ( w ) f ^ ( r ) ( w ) ) |                              < | f ^ ( 1 ) ( w ) f ( 1 ) ( w ) | + | f ( r ) ( w ) f ^ ( r ) ( w ) | P 0 , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGceaqabeaacaGG8bGabmOzayaajaWaaSbaaSqaaiaacIcacaaIXaGaaiykaaqabaGccaGGOaGaeyOeI0Iaam4DaiaacMcacqGHsislceWGMbGbaKaadaWgaaWcbaGaaiikaiaadkhacaGGPaaabeaakiaacIcacaWG3bGaaiykaiaacYhacqGH9aqpcaGG8bGaaiikaiqadAgagaqcamaaBaaaleaacaGGOaGaaGymaiaacMcaaeqaaOGaaiikaiabgkHiTiaadEhacaGGPaGaeyOeI0IaamOzamaaBaaaleaacaGGOaGaaGymaiaacMcaaeqaaOGaaiikaiabgkHiTiaadEhacaGGPaGaaiykaiabgUcaRiaacIcacaWGMbWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaam4DaiaacMcacqGHsislceWGMbGbaKaadaWgaaWcbaGaaiikaiaadkhacaGGPaaabeaakiaacIcacaWG3bGaaiykaiaacMcacaGG8baabaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiabgYda8iaacYhaceWGMbGbaKaadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bGaaiykaiabgkHiTiaadAgadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bGaaiykaiaacYhacqGHRaWkcaGG8bGaamOzamaaBaaaleaacaGGOaGaamOCaiaacMcaaeqaaOGaaiikaiaadEhacaGGPaGaeyOeI0IabmOzayaajaWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaam4DaiaacMcacaGG8bGaeyOKH46aaWbaaSqabeaacaWGqbaaaOGaaGimaiaacYcaaaaa@9EF5@ (10)

To prove that Δ ^ E R S S P Δ MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafuiLdqKbaKaadaWgaaWcbaGaamyraiaadkfacaWGtbGaam4uaaqabaGcdaGdKaWcbaGaamiuaaqabOGaayPKHaGaeuiLdqeaaa@40D5@ , where Δ MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeuiLdqeaaa@396B@ =1 under the null hypothesis, we need to show that
| Δ ^ E R S S Δ | p 0 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaqWaaeaacuqHuoargaqcamaaBaaaleaacaWGfbGaamOuaiaadofacaWGtbaabeaakiabgkHiTiabfs5aebGaay5bSlaawIa7amaaoqcaleaacaWGWbaabeGccaGLsgcacaaIWaaaaa@45BE@ .Now
| Δ ^ E R S S Δ | = 1 2 | | f ^ ( 1 ) ( w ) f ^ ( r ) ( w ) | | d w | | f ( 1 ) ( w ) f ( r ) ( w ) | d w                  1 2 [ | | f ^ ( 1 ) ( w ) f ^ ( r ) ( w ) | | f ( 1 ) ( w ) f ( r ) ( w ) | | ] d w                 1 2 | [ f ^ ( 1 ) ( w ) f ^ ( r ) ( w ) ] [ f ( 1 ) ( w ) f ( r ) ( w ) ] | d w                 1 2 | f ^ ( 1 ) ( w ) f ( 1 ) ( w ) | d w + 1 2 | f ( r ) ( w ) f ^ ( r ) ( w ) | d w . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGceaqabeaadaabdaqaaiqbfs5aezaajaWaaSbaaSqaaiaadweacaWGsbGaam4uaiaadofaaeqaaOGaeyOeI0IaeuiLdqeacaGLhWUaayjcSdGaeyypa0ZaaSaaaeaacaaIXaaabaGaaGOmaaaadaabdaqaamaapeaabaGaaiiFaiqadAgagaqcamaaBaaaleaacaGGOaGaaGymaiaacMcaaeqaaOGaaiikaiabgkHiTiaadEhacaGGPaGaeyOeI0IabmOzayaajaWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaam4DaiaacMcacaGG8baaleqabeqdcqGHRiI8aaGccaGLhWUaayjcSdGaamizaiaadEhacqGHsisldaabdaqaamaapeaabaGaaiiFaiaadAgadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bGaaiykaiabgkHiTiaadAgadaWgaaWcbaGaaiikaiaadkhacaGGPaaabeaakiaacIcacaWG3bGaaiykaaWcbeqab0Gaey4kIipaaOGaay5bSlaawIa7aiaadsgacaWG3baabaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiabgsMiJoaalaaabaGaaGymaaqaaiaaikdaaaWaa8qaaeaadaWadaqaamaaemaabaGaaiiFaiqadAgagaqcamaaBaaaleaacaGGOaGaaGymaiaacMcaaeqaaOGaaiikaiabgkHiTiaadEhacaGGPaGaeyOeI0IabmOzayaajaWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaam4DaiaacMcacaGG8bGaeyOeI0IaaiiFaiaadAgadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bGaaiykaiabgkHiTiaadAgadaWgaaWcbaGaaiikaiaadkhacaGGPaaabeaakiaacIcacaWG3bGaaiykaiaacYhaaiaawEa7caGLiWoaaiaawUfacaGLDbaaaSqabeqaniabgUIiYdGccaWGKbGaam4DaaqaaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiabgsMiJoaalaaabaGaaGymaaqaaiaaikdaaaWaa8qaaeaadaabdaqaaiaacUfaceWGMbGbaKaadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bGaaiykaiabgkHiTiqadAgagaqcamaaBaaaleaacaGGOaGaamOCaiaacMcaaeqaaOGaaiikaiaadEhacaGGPaGaaiyxaiabgkHiTiaacUfacaWGMbWaaSbaaSqaaiaacIcacaaIXaGaaiykaaqabaGccaGGOaGaeyOeI0Iaam4DaiaacMcacqGHsislcaWGMbWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaam4DaiaacMcacaGGDbaacaGLhWUaayjcSdaaleqabeqdcqGHRiI8aOGaamizaiaadEhaaeaacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacqGHKjYOdaWcaaqaaiaaigdaaeaacaaIYaaaamaapeaabaGaaiiFaiqadAgagaqcamaaBaaaleaacaGGOaGaaGymaiaacMcaaeqaaOGaaiikaiabgkHiTiaadEhacaGGPaGaeyOeI0IaamOzamaaBaaaleaacaGGOaGaaGymaiaacMcaaeqaaOGaaiikaiabgkHiTiaadEhacaGGPaGaaiiFaaWcbeqab0Gaey4kIipakiaadsgacaWG3bGaey4kaSYaaSaaaeaacaaIXaaabaGaaGOmaaaadaWdbaqaaiaacYhacaWGMbWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaam4DaiaacMcacqGHsislceWGMbGbaKaadaWgaaWcbaGaaiikaiaadkhacaGGPaaabeaakiaacIcacaWG3bGaaiykaiaacYhaaSqabeqaniabgUIiYdGccaWGKbGaam4Daiaac6caaaaa@1033@ Since, f ^ ( 1 ) ( w ) p f ( 1 ) ( w )  and  f ^ ( r ) ( w ) p f ( r ) ( w ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGabmOzayaajaWaaSbaaSqaaiaacIcacaaIXaGaaiykaaqabaGccaGGOaGaeyOeI0Iaam4DaiaacMcadaGdKaWcbaGaamiCaaqabOGaayPKHaGaamOzamaaBaaaleaacaGGOaGaaGymaiaacMcaaeqaaOGaaiikaiabgkHiTiaadEhacaGGPaGaaeiiaiaabggacaqGUbGaaeizaiaabccaceWGMbGbaKaadaWgaaWcbaGaaiikaiaadkhacaGGPaaabeaakiaacIcacaWG3bGaaiykamaaoqcaleaacaWGWbaabeGccaGLsgcacaWGMbWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaam4DaiaacMcaaaa@59BB@ for all continuous points w of f ( 1 )  and  f ( r ) , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqiVu0Je9sqqrpepC0xbbL8F4rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzamaaBaaaleaacaGGOaGaaGymaiaacMcaaeqaaOGaaeiiaiaabggacaqGUbGaaeizaiaabccacaWGMbWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGSaaaaa@414A@ we have
lim m | Δ ^ E R S S Δ | 1 2 lim m | f ^ ( 1 ) ( w ) f ( 1 ) ( w ) | d w + 1 2 lim m | f ( r ) ( w ) f ^ ( r ) ( w ) | d w . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaCbeaeaaciGGSbGaaiyAaiaac2gaaSqaaiaad2gacqGHsgIRcqGHEisPaeqaaOWaaqWaaeaacuqHuoargaqcamaaBaaaleaacaWGfbGaamOuaiaadofacaWGtbaabeaakiabgkHiTiabfs5aebGaay5bSlaawIa7aiabgsMiJoaalaaabaGaaGymaaqaaiaaikdaaaWaaCbeaeaaciGGSbGaaiyAaiaac2gaaSqaaiaad2gacqGHsgIRcqGHEisPaeqaaOWaa8qaaeaacaGG8bGabmOzayaajaWaaSbaaSqaaiaacIcacaaIXaGaaiykaaqabaGccaGGOaGaeyOeI0Iaam4DaiaacMcacqGHsislcaWGMbWaaSbaaSqaaiaacIcacaaIXaGaaiykaaqabaGccaGGOaGaeyOeI0Iaam4DaiaacMcacaGG8baaleqabeqdcqGHRiI8aOGaamizaiaadEhacqGHRaWkdaWcaaqaaiaaigdaaeaacaaIYaaaamaaxababaGaciiBaiaacMgacaGGTbaaleaacaWGTbGaeyOKH4QaeyOhIukabeaakmaapeaabaGaaiiFaiaadAgadaWgaaWcbaGaaiikaiaadkhacaGGPaaabeaakiaacIcacaWG3bGaaiykaiabgkHiTiqadAgagaqcamaaBaaaleaacaGGOaGaamOCaiaacMcaaeqaaOGaaiikaiaadEhacaGGPaGaaiiFaaWcbeqab0Gaey4kIipakiaadsgacaWG3bGaaiOlaaaa@853F@

Now for a given ε > 0 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyTduMaeyOpa4JaaGimaaaa@3B6E@ , let A = { w : | f ^ ( 1 ) ( w ) f ( 1 ) ( w ) | > ε / 2 }  and  A 1 = { w : | f ( r ) ( w ) f ^ ( r ) ( w ) | > ε / 2 } MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyqaiabg2da9iaacUhacaWG3bGaaiOoaiaacYhaceWGMbGbaKaadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bGaaiykaiabgkHiTiaadAgadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bGaaiykaiaacYhacqGH+aGpcqaH1oqzcaGGVaGaaGOmaiaac2hacaqGGaGaaeyyaiaab6gacaqGKbGaaeiiaiaadgeadaWgaaWcbaGaaGymaaqabaGccqGH9aqpcaGG7bGaam4DaiaacQdacaGG8bGaamOzamaaBaaaleaacaGGOaGaamOCaiaacMcaaeqaaOGaaiikaiaadEhacaGGPaGaeyOeI0IabmOzayaajaWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaam4DaiaacMcacaGG8bGaeyOpa4JaeqyTduMaai4laiaaikdacaGG9baaaa@6EB4@ , we have
lim m | Δ ^ E R S S Δ | 1 2 lim m | f ^ ( 1 ) ( w ) f ( 1 ) ( w ) | d w + 1 2 lim m | f ( r ) ( w ) f ^ ( r ) ( w ) | d w                       = 1 2 lim m { A | f ^ ( 1 ) ( w ) f ( 1 ) ( w ) | d w + A c | f ^ ( 1 ) ( w ) f ( 1 ) ( w ) | d w                             + A 1 | f ( r ) ( w ) f ^ ( r ) ( w ) | d w + A 1 c | f ( r ) ( w ) f ^ ( r ) ( w ) | d w } , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGceaqabeaadaWfqaqaaiGacYgacaGGPbGaaiyBaaWcbaGaamyBaiabgkziUkabg6HiLcqabaGcdaabdaqaaiqbfs5aezaajaWaaSbaaSqaaiaadweacaWGsbGaam4uaiaadofaaeqaaOGaeyOeI0IaeuiLdqeacaGLhWUaayjcSdGaeyizIm6aaSaaaeaacaaIXaaabaGaaGOmaaaadaWfqaqaaiGacYgacaGGPbGaaiyBaaWcbaGaamyBaiabgkziUkabg6HiLcqabaGcdaWdbaqaaiaacYhaceWGMbGbaKaadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bGaaiykaiabgkHiTiaadAgadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bGaaiykaiaacYhaaSqabeqaniabgUIiYdGccaWGKbGaam4DaiabgUcaRmaalaaabaGaaGymaaqaaiaaikdaaaWaaCbeaeaaciGGSbGaaiyAaiaac2gaaSqaaiaad2gacqGHsgIRcqGHEisPaeqaaOWaa8qaaeaacaGG8bGaamOzamaaBaaaleaacaGGOaGaamOCaiaacMcaaeqaaOGaaiikaiaadEhacaGGPaGaeyOeI0IabmOzayaajaWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaam4DaiaacMcacaGG8baaleqabeqdcqGHRiI8aOGaamizaiaadEhaaeaacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacqGH9aqpdaWcaaqaaiaaigdaaeaacaaIYaaaamaaxababaGaciiBaiaacMgacaGGTbaaleaacaWGTbGaeyOKH4QaeyOhIukabeaakiaacUhadaWdrbqaaiaacYhaceWGMbGbaKaadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bGaaiykaiabgkHiTiaadAgadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG3bGaaiykaiaacYhaaSqaaiaadgeaaeqaniabgUIiYdGccaWGKbGaam4DaiabgUcaRmaapefabaGaaiiFaiqadAgagaqcamaaBaaaleaacaGGOaGaaGymaiaacMcaaeqaaOGaaiikaiabgkHiTiaadEhacaGGPaGaeyOeI0IaamOzamaaBaaaleaacaGGOaGaaGymaiaacMcaaeqaaOGaaiikaiabgkHiTiaadEhacaGGPaaaleaacaWGbbWaaWbaaWqabeaacaWGJbaaaaWcbeqdcqGHRiI8aOGaaiiFaiaadsgacaWG3baabaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaey4kaSYaa8quaeaacaGG8bGaamOzamaaBaaaleaacaGGOaGaamOCaiaacMcaaeqaaOGaaiikaiaadEhacaGGPaGaeyOeI0IabmOzayaajaWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaam4DaiaacMcacaGG8baaleaacaWGbbWaaSbaaWqaaiaaigdaaeqaaaWcbeqdcqGHRiI8aOGaamizaiaadEhacqGHRaWkdaWdrbqaaiaacYhacaWGMbWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaam4DaiaacMcacqGHsislceWGMbGbaKaadaWgaaWcbaGaaiikaiaadkhacaGGPaaabeaakiaacIcacaWG3bGaaiykaaWcbaGaamyqamaaDaaameaacaaIXaaabaGaam4yaaaaaSqab0Gaey4kIipakiaacYhacaWGKbGaam4Daiaac2hacaGGSaaaaaa@073F@
then
lim m | Δ ^ E R S S Δ | 1 2 lim m { P ( { w : | f ^ ( 1 ) ( w ) f ( 1 ) ( w ) | > ε / 2 } ) + ε / 2                             + P ( { w : | f ( r ) ( w ) f ^ ( r ) ( w ) | > ε / 2 } ) + ε / 2 } . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGceaqabeaadaWfqaqaaiGacYgacaGGPbGaaiyBaaWcbaGaamyBaiabgkziUkabg6HiLcqabaGcdaabdaqaaiqbfs5aezaajaWaaSbaaSqaaiaadweacaWGsbGaam4uaiaadofaaeqaaOGaeyOeI0IaeuiLdqeacaGLhWUaayjcSdGaeyizIm6aaSaaaeaacaaIXaaabaGaaGOmaaaadaWfqaqaaiGacYgacaGGPbGaaiyBaaWcbaGaamyBaiabgkziUkabg6HiLcqabaGccaGG7bGaamiuaiaacIcacaGG7bGaam4DaiaacQdacaGG8bGabmOzayaajaWaaSbaaSqaaiaacIcacaaIXaGaaiykaaqabaGccaGGOaGaeyOeI0Iaam4DaiaacMcacqGHsislcaWGMbWaaSbaaSqaaiaacIcacaaIXaGaaiykaaqabaGccaGGOaGaeyOeI0Iaam4DaiaacMcacaGG8bGaeyOpa4JaeqyTduMaai4laiaaikdacaGG9bGaaiykaiabgUcaRiabew7aLjaac+cacaaIYaaabaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaey4kaSIaamiuaiaacIcacaGG7bGaam4DaiaacQdacaGG8bGaamOzamaaBaaaleaacaGGOaGaamOCaiaacMcaaeqaaOGaaiikaiaadEhacaGGPaGaeyOeI0IabmOzayaajaWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaam4DaiaacMcacaGG8bGaeyOpa4JaeqyTduMaai4laiaaikdacaGG9bGaaiykaiabgUcaRiabew7aLjaac+cacaaIYaGaaiyFaiaac6caaaaa@A3B5@

Clearly lim m | Δ ^ E R S S Δ | 0 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaCbeaeaaciGGSbGaaiyAaiaac2gaaSqaaiaad2gacqGHsgIRcqGHEisPaeqaaOWaaqWaaeaacuqHuoargaqcamaaBaaaleaacaWGfbGaamOuaiaadofacaWGtbaabeaakiabgkHiTiabfs5aebGaay5bSlaawIa7aiabgkziUkaaicdaaaa@4C81@ in probability, i.e. | Δ ^ E R S S Δ | p 0 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaqWaaeaacuqHuoargaqcamaaBaaaleaacaWGfbGaamOuaiaadofacaWGtbaabeaakiabgkHiTiabfs5aebGaay5bSlaawIa7amaaoqcaleaacaWGWbaabeGccaGLsgcacaaIWaaaaa@45BE@ . Hence

Δ^ERSSpΔ.

The asymptotic distribution of

Δ ^ E R S S MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafuiLdqKbaKaadaWgaaWcbaGaamyraiaadkfacaWGtbGaam4uaaqabaaaaa@3CF8@ under the null hypothesis, using the results derived by Anderson et al.29 is as follows: Let f 1 ( x ) = f ( 1 ) ( x )  and  f 2 ( x ) = f ( r ) ( x ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzamaaBaaaleaacaaIXaaabeaakiaacIcacaWG4bGaaiykaiabg2da9iaadAgadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG4bGaaiykaiaabccacaqGHbGaaeOBaiaabsgacaqGGaGaamOzamaaBaaaleaacaaIYaaabeaakiaacIcacaWG4bGaaiykaiabg2da9iaadAgadaWgaaWcbaGaaiikaiaadkhacaGGPaaabeaakiaacIcacaWG4bGaaiykaaaa@52B7@ ,
C f 1 , f 2 = { x : f 1 ( x ) = f 2 ( x ) > 0 } , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBaaaleaacaWGMbWaaSbaaWqaaiaaigdaaeqaaSGaaiilaiaadAgadaWgaaadbaGaaGOmaaqabaaaleqaaOGaeyypa0Jaai4EaiaadIhacqGHiiIZcqWIDesOcaGG6aGaamOzamaaBaaaleaacaaIXaaabeaakiaacIcacaWG4bGaaiykaiabg2da9iaadAgadaWgaaWcbaGaaGOmaaqabaGccaGGOaGaamiEaiaacMcacqGH+aGpcaaIWaGaaiyFaiaacYcaaaa@5102@ , C f 1 = { x : f 1 ( x ) < f 2 ( x ) > 0 } MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBaaaleaacaWGMbWaaSbaaWqaaiaaigdaaeqaaaWcbeaakiabg2da9iaacUhacaWG4bGaeyicI4SaeSyhHeQaaiOoaiaadAgadaWgaaWcbaGaaGymaaqabaGccaGGOaGaamiEaiaacMcacqGH8aapcaWGMbWaaSbaaSqaaiaaikdaaeqaaOGaaiikaiaadIhacaGGPaGaeyOpa4JaaGimaiaac2haaaa@4DC1@ and C f 2 = { x : f 1 ( x ) > f 2 ( x ) > 0 } . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaam4qamaaBaaaleaacaWGMbWaaSbaaWqaaiaaikdaaeqaaaWcbeaakiabg2da9iaacUhacaWG4bGaeyicI4SaeSyhHeQaaiOoaiaadAgadaWgaaWcbaGaaGymaaqabaGccaGGOaGaamiEaiaacMcacqGH+aGpcaWGMbWaaSbaaSqaaiaaikdaaeqaaOGaaiikaiaadIhacaGGPaGaeyOpa4JaaGimaiaac2hacaGGUaaaaa@4E78@ Let n 1 = n 2 = m MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOBamaaBaaaleaacaaIXaaabeaakiabg2da9iaad6gadaWgaaWcbaGaaGOmaaqabaGccqGH9aqpcaWGTbaaaa@3ECC@ , h = h + = h MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiAamaaBaaaleaacqGHsislaeqaaOGaeyypa0JaamiAamaaBaaaleaacqGHRaWkaeqaaOGaeyypa0JaamiAaaaa@3F13@ p 0 = P ( X C f 1 , f 2 ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiCamaaBaaaleaacaaIWaaabeaakiabg2da9iaadcfacaGGOaGaamiwaiabgIGiolaadoeadaWgaaWcbaGaamOzamaaBaaameaacaaIXaaabeaaliaacYcacaWGMbWaaSbaaWqaaiaaikdaaeqaaaWcbeaakiaacMcaaaa@44EA@ , p 1 = P ( X C f 1 ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiCamaaBaaaleaacaaIXaaabeaakiabg2da9iaadcfacaGGOaGaamiwaiabgIGiolaadoeadaWgaaWcbaGaamOzamaaBaaameaacaaIXaaabeaaaSqabaGccaGGPaaaaa@425C@ and p 2 = P ( X C f 2 ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamiCamaaBaaaleaacaaIYaaabeaakiabg2da9iaadcfacaGGOaGaamiwaiabgIGiolaadoeadaWgaaWcbaGaamOzamaaBaaameaacaaIYaaabeaaaSqabaGccaGGPaaaaa@425E@ .Under the above assumption, we have the following asymptotic result:
m ( Δ ^ Δ ) a m N ( 0 , v ) , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaOaaaeaacaWGTbaaleqaaOGaaiikaiqbfs5aezaajaGaeyOeI0IaeyiLdqKaaiykaiabgkHiTiaadggadaWgaaWcbaGaamyBaaqabaGccqGHshI3caWGobGaaiikaiaaicdacaGGSaGaamODaiaacMcacaGGSaaaaa@48D8@ where, v = p 0 σ 0 2 + σ 1 2 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamODaiabg2da9iaadchadaWgaaWcbaGaaGimaaqabaGccqaHdpWCdaqhaaWcbaGaaGimaaqaaiaaikdaaaGccqGHRaWkcqaHdpWCdaqhaaWcbaGaaGymaaqaaiaaikdaaaaaaa@43A4@ , σ 1 2 = p 1 ( 1 p 1 ) + p 2 ( 1 p 2 ) , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4Wdm3aa0baaSqaaiaaigdaaeaacaaIYaaaaOGaeyypa0JaamiCamaaBaaaleaacaaIXaaabeaakiaacIcacaaIXaGaeyOeI0IaamiCamaaBaaaleaacaaIXaaabeaakiaacMcacqGHRaWkcaWGWbWaaSbaaSqaaiaaikdaaeqaaOGaaiikaiaaigdacqGHsislcaWGWbWaaSbaaSqaaiaaikdaaeqaaOGaaiykaiaacYcaaaa@4BAA@ a m = R h C f 1 , f 2 f ( r ) 1 / 2 ( x ) d x . E ( min { Z 1 , Z 2 } ) , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamyyamaaBaaaleaacaWGTbaabeaakiabg2da9maakaaabaWaaSaaaeaacaWGsbaabaGaamiAaaaaaSqabaGcdaWdraqaaiaadAgadaqhaaWcbaGaaiikaiaadkhacaGGPaaabaGaaGymaiaac+cacaaIYaaaaOGaaiikaiaadIhacaGGPaGaamizaiaadIhaaSqaaiaadoeadaWgaaadbaGaamOzamaaBaaabaGaaGymaiaacYcaaeqaaiaadAgadaWgaaqaaiaaikdaaeqaaaqabaaaleqaniabgUIiYdGccaGGUaGaamyraiaacIcaciGGTbGaaiyAaiaac6gacaGG7bGaamOwamaaBaaaleaacaaIXaaabeaakiaacYcacaWGAbWaaSbaaSqaaiaaikdaaeqaaOGaaiyFaiaacMcacaGGSaaaaa@5AE5@
σ 0 2 = R T 0 ( cov ( min { Z 1 , Z 2 } , min { ρ ( t ) Z 1 + 1 ρ ( t ) 2 Z 3 , ρ ( t ) Z 2 + 1 ρ ( t ) 2 Z 4 } ) d t , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeq4Wdm3aa0baaSqaaiaaicdaaeaacaaIYaaaaOGaeyypa0JaamOuamaapebabaWaaeWaaeaaciGGJbGaai4BaiaacAhacaGGOaGaciyBaiaacMgacaGGUbGaai4EaiaadQfadaWgaaWcbaGaaGymaaqabaGccaGGSaGaamOwamaaBaaaleaacaaIYaaabeaakiaac2hacaGGSaGaciyBaiaacMgacaGGUbGaai4Eaiabeg8aYjaacIcacaWG0bGaaiykaiaadQfadaWgaaWcbaGaaGymaaqabaGccqGHRaWkdaGcaaqaaiaaigdacqGHsislcqaHbpGCcaGGOaGaamiDaiaacMcadaahaaWcbeqaaiaaikdaaaaabeaakiaadQfadaWgaaWcbaGaaG4maaqabaGccaGGSaGaeqyWdiNaaiikaiaadshacaGGPaGaamOwamaaBaaaleaacaaIYaaabeaakiabgUcaRmaakaaabaGaaGymaiabgkHiTiabeg8aYjaacIcacaWG0bGaaiykamaaCaaaleqabaGaaGOmaaaaaeqaaOGaamOwamaaBaaaleaacaaI0aaabeaakiaac2haaiaawIcacaGLPaaaaSqaaiaadsfadaWgaaadbaGaaGimaaqabaaaleqaniabgUIiYdGccaWGKbGaamiDaiaacYcaaaa@76CB@
ρ ( t ) = 1 R K ( u ) K ( u + t ) d u ,   R = K 2 ( u ) d u , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeqyWdiNaaiikaiaadshacaGGPaGaeyypa0ZaaSaaaeaacaaIXaaabaGaamOuaaaadaWdraqaaiaadUeacaGGOaGaamyDaiaacMcacaWGlbGaaiikaiaadwhacqGHRaWkcaWG0bGaaiykaiaadsgacaWG1bGaaiilaaWcbaGaeSyhHekabeqdcqGHRiI8aOGaaeiiaiaadkfacqGH9aqpdaWdraqaaiaadUeadaahaaWcbeqaaiaaikdaaaGccaGGOaGaamyDaiaacMcacaWGKbGaamyDaiaacYcaaSqaaiabl2riHcqab0Gaey4kIipaaaa@599A@ Z 1 , Z 2 , Z 3  and  Z 4 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOwamaaBaaaleaacaaIXaaabeaakiaacYcacaWGAbWaaSbaaSqaaiaaikdaaeqaaOGaaiilaiaadQfadaWgaaWcbaGaaG4maaqabaGccaqGGaGaaeyyaiaab6gacaqGKbGaaeiiaiaadQfadaWgaaWcbaGaaGinaaqabaaaaa@44A3@ are independent standard normal variables and T 0 = { t : | t | < 1 } . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamivamaaBaaaleaacaaIWaaabeaakiabg2da9iaacUhacaWG0bGaeyicI4SaeSyhHeQaaiOoaiaacYhacaWG0bGaaiiFaiabgYda8iaaigdacaGG9bGaaiOlaaaa@46E9@ However, under the null hypothesis that f 1 ( x ) = f ( 1 ) ( x ) = f 2 ( x ) = f ( r ) ( x ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaamOzamaaBaaaleaacaaIXaaabeaakiaacIcacaWG4bGaaiykaiabg2da9iaadAgadaWgaaWcbaGaaiikaiaaigdacaGGPaaabeaakiaacIcacqGHsislcaWG4bGaaiykaiabg2da9iaadAgadaWgaaWcbaGaaGOmaaqabaGccaGGOaGaamiEaiaacMcacqGH9aqpcaWGMbWaaSbaaSqaaiaacIcacaWGYbGaaiykaaqabaGccaGGOaGaamiEaiaacMcaaaa@4FBB@ then the above result is reduced to
m ( Δ ^ E R S S Δ ) H 0 a m N ( 0 , σ 0 2 ) . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaWaaOaaaeaacaWGTbaaleqaaOGaaiikaiqbfs5aezaajaWaaSbaaSqaaiaadweacaWGsbGaam4uaiaadofaaeqaaOGaeyOeI0IaeyiLdqKaaiykamaaBaaaleaacaWGibWaaSbaaWqaaiaaicdaaeqaaaWcbeaakiabgkHiTiaadggadaWgaaWcbaGaamyBaaqabaGccqGHshI3caWGobGaaiikaiaaicdacaGGSaGaeq4Wdm3aa0baaSqaaiaaicdaaeaacaaIYaaaaOGaaiykaiaac6caaaa@50CB@

Simulation study

To get some insights about the performance of our new test of symmetry based on Δ ^ E R S S MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafuiLdqKbaKaadaWgaaWcbaGaamyraiaadkfacaWGtbGaam4uaaqabaaaaa@3CF8@ we conducted the following simulation. We compared our proposed test of symmetry with its counterpart used by McWilliams,4 Modarres and Gastwirth 8 and Samawi et al.12 (using overlap measure) tests of symmetry.

McWilliams4 runs test is described as follows: For any random sample of size n, let Y ( 1 ) , Y ( 2 ) , ...   , Y ( n ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGzbWaaSbaaSqaaiaacIcacaaIXaGaaiykaaqabaGccaGGSaGaamywamaaBaaaleaacaGGOaGaaGOmaiaacMcacaGGSaGaaiOlaiaac6cacaGGUaGaaeiiaiaacYcaaeqaaOGaamywamaaBaaaleaacaGGOaGaamOBaiaacMcaaeqaaaaa@4728@ denote the sample values ordered from the smallest to largest according to their absolute value (signs are retained), and S 1 , S 2 , ... , S n MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGtbWaaSbaaSqaaiaaigdaaeqaaOGaaiilaiaadofadaWgaaWcbaGaaGOmaiaacYcacaGGUaGaaiOlaiaac6cacaGGSaaabeaakiaadofadaWgaaWcbaGaamOBaaqabaaaaa@4268@ denote indicator variables designating the sign of the Y ( j ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGzbWaaSbaaSqaaiaacIcacaWGQbGaaiykaaqabaaaaa@3C0A@ values [ S j = 1  if  Y ( j )  is nonnegative, 0 otherwise MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGtbWaaSbaaSqaaiaadQgaaeqaaOGaeyypa0JaaGymaiaabccacaqGPbGaaeOzaiaabccacaWGzbWaaSbaaSqaaiaacIcacaWGQbGaaiykaaqabaGccaqGGaGaaeyAaiaabohacaqGGaGaaeOBaiaab+gacaqGUbGaaeOBaiaabwgacaqGNbGaaeyyaiaabshacaqGPbGaaeODaiaabwgacaqGSaGaaeiiaiaabcdacaqGGaGaae4BaiaabshacaqGObGaaeyzaiaabkhacaqG3bGaaeyAaiaabohacaqGLbaaaa@5B71@ ]. Thus, the test statistic used for testing symmetry is R * MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGsbWaaWbaaSqabeaacaGGQaaaaaaa@3A6A@ = the number of runs in S 1 , S 2 , ... , S n MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGtbWaaSbaaSqaaiaaigdaaeqaaOGaaiilaiaadofadaWgaaWcbaGaaGOmaiaacYcacaGGUaGaaiOlaiaac6cacaGGSaaabeaakiaadofadaWgaaWcbaGaamOBaaqabaaaaa@4268@ sequence= 1 + j = 2 n I j MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaaliaaigdacqGHRaWkkmaaqahabaGaamysamaaBaaaleaacaWGQbaabeaaaeaacaWGQbGaeyypa0JaaGOmaaqaaiaad6gaa0GaeyyeIuoaaaa@422E@ , where .
I j = { 0    if  S j = S j 1 1    if  S j S j 1   MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGjbWaaSbaaSqaaiaadQgaaeqaaOGaeyypa0ZaaiqaaeaafaqabeGabaaabaGaaGimaiaabccacaqGGaGaaeiiaiaabMgacaqGMbGaaeiiaiaadofadaWgaaWcbaGaamOAaaqabaGccqGH9aqpcaWGtbWaaSbaaSqaaiaadQgacqGHsislcaaIXaaabeaaaOqaaiaaigdacaqGGaGaaeiiaiaabccacaqGPbGaaeOzaiaabccacaWGtbWaaSbaaSqaaiaadQgaaeqaaOGaeyiyIKRaam4uamaaBaaaleaacaWGQbGaeyOeI0IaaGymaaqabaGccaqGGaaaaaGaay5Eaaaaaa@55C3@

The test is to reject the null hypothesis if R * MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWGsbWaaWbaaSqabeaacaGGQaaaaaaa@3A6A@ is smaller than a critical value ( c α ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaGGOaGaam4yamaaBaaaleaacqaHXoqyaeqaaOGaaiykaaaa@3CCE@ at level of significant α MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaHXoqyaaa@3A57@ . However, the Modarres and Gastwirth8 test is M p = 1 + j = 2 + n p n ϕ ( j ) I j MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaaliaad2eadaWgaaadbaGaamiCaaqabaWccqGH9aqpcaaIXaGaey4kaSIcdaaeWbqaaiabew9aMjaacIcacaWGQbGaaiykaiaadMeadaWgaaWcbaGaamOAaaqabaaabaGaamOAaiabg2da9iaaikdacqGHRaWkcaWGUbWaaSbaaWqaaiaadchaaeqaaaWcbaGaamOBaaqdcqGHris5aaaa@4C45@ , where,
ϕ ( j ) = { j n p        if j> n p 0           otherwise,   MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaHvpGzcaGGOaGaamOAaiaacMcacqGH9aqpdaGabaqaauaabeqaceaaaeaacaWGQbGaeyOeI0IaamOBamaaBaaaleaacaWGWbaabeaakiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGPbGaaeOzaiaabccacaqGQbGaaeOpaiaad6gadaWgaaWcbaGaamiCaaqabaaakeaacaaIWaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaabccacaqGGaGaaeiiaiaab+gacaqG0bGaaeiAaiaabwgacaqGYbGaae4DaiaabMgacaqGZbGaaeyzaiaabYcaaaaacaGL7baacaqGGaaaaa@5E87@
and  n p  is an interger . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaqGHbGaaeOBaiaabsgacaqGGaGaamOBamaaBaaaleaacaWGWbaabeaakiaabccacaqGPbGaae4CaiaabccacaqGHbGaaeOBaiaabccacaqGPbGaaeOBaiaabshacaqGLbGaaeOCaiaabEgacaqGLbGaaeOCaiaab6caaaa@4BFE@

If p=0, terms are Wilcoxon scores. Otherwise, they are percentile-modified scores. The Modarres and Gastwirth9 test is the hybrid test of sign test in the first stage and a percentile-modified two-sample Wilcoxon test in the second stage ( W 0.80 ) . MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaaliaacIcacaWGxbWaaSbaaWqaaiaaicdacaGGUaGaaGioaiaaicdaaeqaaSGaaiykaiaac6caaaa@3ECA@

In this simulation, SAS version 9.2 {proc kde; method=srot} is used. The generalized lambda distribution see, Ramberg and Schmeiser36 is used in our simulation with the following set of parameters:
1- λ 1 = 0 , λ 2 = 0.197454 , λ 3 = 0.134915 , λ 4 = 0.134915 ,  (Symmetric) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH7oaBdaWgaaWcbaGaaGymaaqabaGccqGH9aqpcaaIWaGaaiilaiabeU7aSnaaBaaaleaacaaIYaaabeaakiabg2da9iaaicdacaGGUaGaaGymaiaaiMdacaaI3aGaaGinaiaaiwdacaaI0aGaaiilaiabeU7aSnaaBaaaleaacaaIZaaabeaakiabg2da9iaaicdacaGGUaGaaGymaiaaiodacaaI0aGaaGyoaiaaigdacaaI1aGaaiilaiabeU7aSnaaBaaaleaacaaI0aaabeaakiabg2da9iaaicdacaGGUaGaaGymaiaaiodacaaI0aGaaGyoaiaaigdacaaI1aGaaiilaiaabccacaqGOaGaae4uaiaabMhacaqGTbGaaeyBaiaabwgacaqG0bGaaeOCaiaabMgacaqGJbGaaeykaaaa@66DA@
2- λ 1 = 0 , λ 2 = 1 , λ 3 = 1.4 , λ 4 = 0.25 ,   MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH7oaBdaWgaaWcbaGaaGymaaqabaGccqGH9aqpcaaIWaGaaiilaiabeU7aSnaaBaaaleaacaaIYaaabeaakiabg2da9iaaigdacaGGSaGaeq4UdW2aaSbaaSqaaiaaiodaaeqaaOGaeyypa0JaaGymaiaac6cacaaI0aGaaiilaiabeU7aSnaaBaaaleaacaaI0aaabeaakiabg2da9iaaicdacaGGUaGaaGOmaiaaiwdacaGGSaGaaeiiaaaa@5154@
3- λ 1 = 0 , λ 2 = 1 , λ 3 = 0.00007 , λ 4 = 0.1 , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH7oaBdaWgaaWcbaGaaGymaaqabaGccqGH9aqpcaaIWaGaaiilaiabeU7aSnaaBaaaleaacaaIYaaabeaakiabg2da9iaaigdacaGGSaGaeq4UdW2aaSbaaSqaaiaaiodaaeqaaOGaeyypa0JaaGimaiaac6cacaaIWaGaaGimaiaaicdacaaIWaGaaG4naiaacYcacqaH7oaBdaWgaaWcbaGaaGinaaqabaGccqGH9aqpcaaIWaGaaiOlaiaaigdacaGGSaaaaa@52DB@
4- λ 1 = 3.586508 , λ 2 = 0.04306 , λ 3 = 0.025213 , λ 4 = 0.094029 , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH7oaBdaWgaaWcbaGaaGymaaqabaGccqGH9aqpcaaIZaGaaiOlaiaaiwdacaaI4aGaaGOnaiaaiwdacaaIWaGaaGioaiaacYcacqaH7oaBdaWgaaWcbaGaaGOmaaqabaGccqGH9aqpcaaIWaGaaiOlaiaaicdacaaI0aGaaG4maiaaicdacaaI2aGaaiilaiabeU7aSnaaBaaaleaacaaIZaaabeaakiabg2da9iaaicdacaGGUaGaaGimaiaaikdacaaI1aGaaGOmaiaaigdacaaIZaGaaiilaiabeU7aSnaaBaaaleaacaaI0aaabeaakiabg2da9iaaicdacaGGUaGaaGimaiaaiMdacaaI0aGaaGimaiaaikdacaaI5aGaaiilaaaa@60E5@
5- λ 1 = 0 , λ 2 = 1 , λ 3 = 0.0075 , λ 4 = 0.03 , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH7oaBdaWgaaWcbaGaaGymaaqabaGccqGH9aqpcaaIWaGaaiilaiabeU7aSnaaBaaaleaacaaIYaaabeaakiabg2da9iabgkHiTiaaigdacaGGSaGaeq4UdW2aaSbaaSqaaiaaiodaaeqaaOGaeyypa0JaeyOeI0IaaGimaiaac6cacaaIWaGaaGimaiaaiEdacaaI1aGaaiilaiabeU7aSnaaBaaaleaacaaI0aaabeaakiabg2da9iabgkHiTiaaicdacaGGUaGaaGimaiaaiodacaGGSaaaaa@55A9@
6- λ 1 = 0.116734 , λ 2 = 0.351663 , λ 3 = 0.13 , λ 4 = 0.16 , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH7oaBdaWgaaWcbaGaaGymaaqabaGccqGH9aqpcqGHsislcaaIWaGaaiOlaiaaigdacaaIXaGaaGOnaiaaiEdacaaIZaGaaGinaiaacYcacqaH7oaBdaWgaaWcbaGaaGOmaaqabaGccqGH9aqpcqGHsislcaaIWaGaaiOlaiaaiodacaaI1aGaaGymaiaaiAdacaaI2aGaaG4maiaacYcacqaH7oaBdaWgaaWcbaGaaG4maaqabaGccqGH9aqpcqGHsislcaaIWaGaaiOlaiaaigdacaaIZaGaaiilaiabeU7aSnaaBaaaleaacaaI0aaabeaakiabg2da9iabgkHiTiaaicdacaGGUaGaaGymaiaaiAdacaGGSaaaaa@5F67@
7- λ 1 = 0 , λ 2 = 1 , λ 3 = 0.1 , λ 4 = 0.18 , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH7oaBdaWgaaWcbaGaaGymaaqabaGccqGH9aqpcaaIWaGaaiilaiabeU7aSnaaBaaaleaacaaIYaaabeaakiabg2da9iabgkHiTiaaigdacaGGSaGaeq4UdW2aaSbaaSqaaiaaiodaaeqaaOGaeyypa0JaeyOeI0IaaGimaiaac6cacaaIXaGaaiilaiabeU7aSnaaBaaaleaacaaI0aaabeaakiabg2da9iabgkHiTiaaicdacaGGUaGaaGymaiaaiIdacaGGSaaaaa@5376@
8- λ 1 = 0 , λ 2 = 1 , λ 3 = 0.001 , λ 4 = 0.13 , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH7oaBdaWgaaWcbaGaaGymaaqabaGccqGH9aqpcaaIWaGaaiilaiabeU7aSnaaBaaaleaacaaIYaaabeaakiabg2da9iabgkHiTiaaigdacaGGSaGaeq4UdW2aaSbaaSqaaiaaiodaaeqaaOGaeyypa0JaeyOeI0IaaGimaiaac6cacaaIWaGaaGimaiaaigdacaGGSaGaeq4UdW2aaSbaaSqaaiaaisdaaeqaaOGaeyypa0JaeyOeI0IaaGimaiaac6cacaaIXaGaaG4maiaacYcaaaa@54E5@
9- λ 1 = 0 , λ 2 = 1 , λ 3 = 0.0001 , λ 4 = 0.17. MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaH7oaBdaWgaaWcbaGaaGymaaqabaGccqGH9aqpcaaIWaGaaiilaiabeU7aSnaaBaaaleaacaaIYaaabeaakiabg2da9iabgkHiTiaaigdacaGGSaGaeq4UdW2aaSbaaSqaaiaaiodaaeqaaOGaeyypa0JaeyOeI0IaaGimaiaac6cacaaIWaGaaGimaiaaicdacaaIXaGaaiilaiabeU7aSnaaBaaaleaacaaI0aaabeaakiabg2da9iabgkHiTiaaicdacaGGUaGaaGymaiaaiEdacaGGUaaaaa@55A5@ To generate the observations we used x i = λ 1 + 1 λ 2 ( u i λ 3 ( 1 u i ) λ 4 ,   i = 1 , ... , m , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWG4bWaaSbaaSqaaiaadMgaaeqaaOGaeyypa0Jaeq4UdW2aaSbaaSqaaiaaigdaaeqaaOGaey4kaSYaaSaaaeaacaaIXaaabaGaeq4UdW2aaSbaaSqaaiaaikdaaeqaaaaakiaacIcacaWG1bWaa0baaSqaaiaadMgaaeaacqaH7oaBdaWgaaadbaGaaG4maaqabaaaaOGaeyOeI0IaaiikaiaaigdacqGHsislcaWG1bWaaSbaaSqaaiaadMgaaeqaaOGaaiykamaaCaaaleqabaGaeq4UdW2aaSbaaWqaaiaaisdaaeqaaaaakiaacYcacaqGGaGaamyAaiabg2da9iaaigdacaGGSaGaaiOlaiaac6cacaGGUaGaaiilaiaad2gacaGGSaaaaa@5A3C@ where u i MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaWG1bWaaSbaaSqaaiaadMgaaeqaaaaa@3ACC@ a uniform random number is. The significance level considered is α = 0.05 , MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaHXoqycqGH9aqpcaaIWaGaaiOlaiaaicdacaaI1aGaaiilaaaa@3EF2@ with sample sizes n=30, 50, and 100. Our simulation is based on 1,000 simulated samples. It is clear that 95% and 99% confidence intervals of the true probability of type I error under the null hypothesis with α = 0.05 MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacqaHXoqycqGH9aqpcaaIWaGaaiOlaiaaicdacaaI1aaaaa@3E42@ are (0.0457, 0.0543) and (0.0435, 0.0575) respectively. Note that in the below tables, the values of skewnees ( α 3 ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaaliaacIcacqaHXoqydaWgaaadbaGaae4maaqabaWccaqGPaaaaa@3CA8@ and kurtosis ( α 4 ) MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0di9LrFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbbG8FasPYRqj0=yi0dXdbba9pGe9xq=JbbG8A8frFve9Fve9Ff0dmeaabaqaciGacaGaaeqabaWaaeaaeaaakeaacaGGOaWccqaHXoqydaWgaaadbaGaaGinaaqabaWccaGGPaaaaa@3CB1@ are from McWilliams.4

Table 1a first showcases the estimated probability of type I error. Our test is an asymptotic test with a slight bias in Δ MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGaeyiLdqeaaa@396C@ estimation and in the variance estimation for a small sample size. For sample sizes more than 30, the test seems to have an estimated probability of type I error close to the nominal value 0.05. Table 1a and Table 1b show that using Δ ^ E R S S MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafuiLdqKbaKaadaWgaaWcbaGaamyraiaadkfacaWGtbGaam4uaaqabaaaaa@3CF8@ based test is more powerful than McWilliams4 and Baklizi.7 In all cases, our proposed procedure is even more efficient than the tests of symmetry proposed by Modarres and Gastwirth,8 Modarres and Gastwirth9 and Samawi et al.12 In all of the tests within the comparison, the power of all tests of symmetry increases as the sample size increases. Finally the power of Δ ^ E R S S MathType@MTEF@5@5@+=feaaguart1ev2aaatCvAUfeBSjuyZL2yd9gzLbvyNv2CaerbuLwBLnhiov2DGi1BTfMBaeXatLxBI9gBaerbd9wDYLwzYbItLDharqqtubsr4rNCHbGeaGqk0Jf9crFfpeea0xh9v8qiW7rqqrFfpeea0xe9Lq=Jc9vqaqpepm0xbba9pwe9Q8fs0=yqaqpepae9pg0FirpepeKkFr0xfr=xfr=xb9adbaqaaeGaciGaaiaabeqaamaabaabaaGcbaGafuiLdqKbaKaadaWgaaWcbaGaamyraiaadkfacaWGtbGaam4uaaqabaaaaa@3CF8@ based tests increases as the set size r increases.

Illustration using noninvasive measurement of cardiac output by electrical velocimetry in neonate data

The samples selected in our illustration is from a study designed to evaluate the effectiveness of a new technology, Electrical Velocimetry (E.V.) for a non-invasive cardiac output (CO) and stroke volume (SV) in neonates.37 One of the research questions is whether the CO measure is the same for low birth-weight infants and non-low birth-weight infants. Low birth weight is defined as less than 1.5 kg. Thus we compared CO for neonates with birth-weight less than 1.5kg to neonates with birth-weight greater or equal to 1.5kg.

As it is frequently the case for this type of study, the underlying distribution is assumed “normal”, or at least symmetric. In either case, a test of symmetry is almost never considered in determining how to proceed in the analysis. Based on the conclusions of a test of symmetry, the analyst can choose the most powerful test for location. However, before deciding on the test procedure, we need to check the assumption of symmetry of underlying distribution of CO for the premature and term infants, with birth weight less than 1.5 gk and with birth weight greater or equal to 1.5 kg.

Distribution

n

R*

Δ^SRS

Samawi et al.12

M0.25 W0.80

Δ^ERSS

r=2

Δ^ERSS

r=3

Δ^ERSS

r=4

Δ^ERSS

r=5

(1)

λ1=0,λ2=0.197454,λ3=0.134915,λ4=0.134915, α3=0,α4=3.0

30

0.047

0.069

0.053

0.055

0.046

0.056

0.056

0.044

50

0.050

0.054

0.048

0.048

0.054

0.056

0.049

0.049

100

0.064

0.053

0.050

0.052

0.048

0.053

0.046

0.051

(2)

λ1=0,λ2=1,λ3=1.4,λ4=0.25 α3=0.5,α4=2.2

30

0.297

0.495

0.583

0.656

0.751

0.906

0.973

0.999

50

0.476

0.836

0.846

0.949

0.997

1.000

1.000

1.000

100

0.776

0.999

0.990

0.999

1.000

1.000

1.000

1.000

(3)

λ1=0,λ2=1,λ3=0.00007,λ4=0.1,α3=1.5,α4=5.8

30

0.438

0.852

0.761

0.762

0.960

0.975

0.999

0.999

50

0.683

0.966

0.950

0.992

1.000

1.000

1.000

1.000

100

0.927

1.000

0.999

1.000

1.000

1.000

1.000

1.000

(4)

λ1=3.586508,λ2=0.04306,λ3=0.025213,λ4=0.094029α3=0.9,α4=4.2

30

0.117

0.375

0.172

0.280

0.384

0.398

0.413

0.482

50

0.131

0.512

0.251

0.544

0.689

0.706

0.766

0.808

100

0.223

0.767

0.414

0.883

0.929

0.940

0.958

0.985

(5)

λ1=0,λ2=1,λ3=0.0075,λ4=0.03,α3=1.5,α4=7.5

30

0.145

0.459

0.234

0.407

0.484

0.569

0.616

0.716

50

0.192

0.580

0.356

0.736

0.832

0.846

0.889

0.921

100

0.338

0.846

0.588

0.972

0.983

0.985

0.991

0.997

Table 1A Probability of Type I Error under the Null Hypothesis. (α=0.05 )

Table 2a consists of two selected SRS samples of CO for neonates with birth weight less than 1.5 gk and neonates with birth weight more than 1.5 kg. Also, Table 2a consists of two selected ERSS samples of CO for both neonates with birth weight less than 1.5 kg and neonates with birth weight more than or equal to1.5 kg. Since CO measure and birth weight are positively correlated, the ranking was performed based on birth weights. ERSS samples in Table 2a consist of first half as the first order statistics and the second half as the third order statistics (the maximum).

Table 2b has the results of the runs and overlap tests of symmetry for the underlying distribution for CO patients. From all samples, we reject the assumption that the underlying distribution is symmetric. Table 2b shows the results of the Mann-Whitney test for two-independent samples. Table 2c shows that there is a significant difference on average in the CO measures between the low birth weight neonates (less than 1.5kg) and the non-low-birth-weight neonates (greater than or equal to 1.5kg).

Based on our simulation and real data example, the proposed test of symmetry based on ERSS sample overlap measure, appears to outperform the other tests of symmetry in the literature in terms of power. Our test is more sensitive to detect a slight asymmetry in the underlying distribution than other tests proposed in the literature. Drawing an ERSS is easier than the ordinary RSS and other RSS variations. Also, the kernel density estimation literature is very rich and many of the proposed and the improved methods are available on statistical software, such as SAS™, S-plus, Stata and R. Since overlap measures can be used in multivariate cases as well as in univariate cases, our proposed test of symmetry can be extended to multivariate cases for diagonal symmetry, conditional symmetry and other types of symmetry. In addition, our test procedure and kernel density estimation are valid under large sample size (n>30) and some regular conditions such as light tail underlying distribution functions. However, our simulation indicates that our procedure is still perform better than other test even for a sample size n=30 or larger and different underlying distributions.

Case #

n

R* Δ^SRS M0.25 W0.80

Δ^ERSS

r=2

Δ^ERSS

r=3

Δ^ERSS

r=4

Δ^ERSS

r=5

(6)

λ1=0.116734,λ2=0.351663,λ3=0.13,λ4=0.16,α3=0.8,α4=11.4

30

0.050

0.155

0.055

0.068

0.253

0.304

0.327

0.431

50

0.056

0.166

0.060

0.077

0.500

0.717

0.739

0.747

100

0.051

0.207

0.068

0.130

0.793

0.875

0.911

0.942

(7)

λ1=0,λ2=1,λ3=0.1,λ4=0.18,α3=2.0,α4=21.2

30

0.090

0.196

0.096

0.166

0.325

0.410

0.458

0.520

50

0.097

0.236

0.125

0.284

0.656

0.667

0.758

0.835

100

0.124

0.354

0.176

0.589

0.891

0.946

0.963

0.974

(8)

λ1=0,λ2=1,λ3=0.001,λ4=0.13,α3=3.16,α4=23.8

30

0.534

1.000

0.830

0.806

1.000

1.000

1.000

1.000

50

0.744

1.000

0.972

0.995

1.000

1.000

1.000

1.000

100

0.972

1.000

1.000

1.000

1.000

1.000

1.000

1.000

(9)

λ1=0,λ2=1,λ3=0.0001,λ4=0.17α3=3.88,α4=40.7

30

0.560

1.000

0.865

0.808

1.000

1.000

1.000

1.000

50

0.816

1.000

0.985

0.997

1.000

1.000

1.000

1.000

100

0.976

1.000

1.000

1.000

1.000

1.000

1.000

1.000

Table 1B (continue) Power of Overlap based test and Run Tests under Alternative Hypotheses. (α=0.05)
*
Results are taken from Modarres and Gastwirth (1996) and Modarres and Gastwirth (1998) respectively

SRS

ERSS

Birth-weight <1.5

Birth-weight 1.5

Birth-weight <1.5

Birth-weight 1.5

0.05

0.21

0.050

0.100

0.06

0.16

0.050

0.120

0.06

0.12

0.050

0.132

0.05

0.20

0.060

0.220

0.05

0.13

0.060

0.150

0.06

0.22

0.060

0.228

0.06

0.22

0.090

0.178

0.08

0.15

0.080

0.182

0.10

0.22

0.080

0.215

0.08

0.20

0.080

0.158

0.08

0.23

0.070

0.155

0.08

0.25

0.070

0.218

0.07

0.23

0.080

0.208

0.11

0.16

0.128

0.220

0.08

0.22

0.128

0.350

0.13

0.16

0.208

0.350

0.16

0.22

0.165

0.360

0.12

0.33

0.162

0.420

0.18

0.26

0.188

0.270

0.16

0.31

0.188

0.440

0.16

0.29

0.262

0.510

0.19

0.27

0.372

0.510

0.37

0.51

0.358

0.520

0.16

0.51

0.222

0.540

0.18

0.52

0.202

0.520

0.11

0.52

0.182

0.520

 

0.53

 

0.530

 

0.73

 

0.730

Table 2A Selected samples of CO data

 

CO Measure

N

 

 

Runs test

 

 

P-value

OVL test

P-value

SRS

Birth-weight <1.5

26

6

 

0.003

-6.414

 

<0.00001

Birth-weight 1.5

28

4

 

<0.00001

-4.356

 

<0.00001

ERSS

Birth-weight <1.5

26

2

 

<0.00001

-14.540

 

<0.00001

Birth-weight 1.5

28

2

 

<0.00001

-4.729

 

<0.00001

Table 2B Runs test of symmetry with summary statistics

Table 2B has the results of the runs and overlap tests of symmetry for the underlying distribution for CO patients. From all samples, we reject the assumption that the underlying distribution is symmetric

 

Mann-Whitney Utest (difference of medians of CO between<1.5 kg and1.5kg weight)

Sample Type

Test

P-value

SRS

59.5

<0.00001

ERSS

118

<0.00001

Table 2C Mann-Whitney test for two-Independent Samples

Shows that there is a significant difference on average in the CO measures between the low birth weight neonates (less than 1.5kg) and the non-low-birth-weight neonates (greater than or equal to 1.5kg)

Acknowledgment

The authors would like to thank the reviewers for their valuable comments which improve the manuscript.

Conflicts of interest

Author declares that there are no conflicts of interest.

References

  1. Butler CC. A test for symmetry using sample distribution function. The Annals of Mathematical Statistics. 1969;40:2211‒2214.
  2. Rothman ED, Woodroofe M. A Cramer‒Von Mises type statistic for testing symmetry. The Annals of Mathematical Statistics. 1972;43(6):2035‒2038.
  3. Hill DL, Rao PV. Test of Symmmetry based on Cramer‒Von Mises statistics. Biometrika. 1977;64(3):489‒494.
  4. McWilliams TP. A distribution–free test of symmetry based on a runs statistic. Journal of the American Statistical Association. 1990;85(412):1130‒1133.
  5. Ozturk O. A nonparametric test of symmetry versus asymmetry for ranked‒set samples. Communication in Statistics‒ Theory and Methods. 2001;30(10):2117‒2133.
  6. Tajuddin IH. Distribution‒Free test for symmetry based on Wilcoxon two‒sample test. Journal of Applied Statistics. 1994;21(5):409‒415.
  7. Baklizi A. A conditional distribution free runs test for symmetry. Nonparametric Statistics. 2003;15(6):713‒718.
  8. Modarres R, Gastwirth JL. A modified runs test of symmetry. Statistics & Probability Letters. 1996;31:107‒112.
  9. Modarres R, Gastwirth JL. Hybrid test for the hypothesis of symmetry. Journal of Applied Statistics. 1998;25(6):777‒783.
  10. Samawi HM, Ahmed MS, Abu‒Dayyeh W. Estimating the population mean using extreme ranked set sampling. Biometrical Journal. 1996;38(5):577‒586.
  11. Samawi HM, Vogel R, Senkowski CK. On Distribution‒Free Runs Test for Symmetry Using Extreme Ranked Set Sampling With an Application Involving Base Deficit Score. Journal of Statistical Theory and Practice. 2010;4(2):289‒301.
  12. Samawi HM, Helu A, Vogel R. A Nonparametric Test of Symmetry Based on the Overlapping Coefficient. Journal of Applied Statistics. 2011;38(5):885‒898.
  13. Weitzman MS. Measures of overlap of income distributions of white and Negro families in the United States. Technical paper No. 22, Department of Commerce, Bureau of Census, Washington. USA, 1970.
  14. Bradley EL, Piantadosi S. The overlapping coefficient as a measure of Agreement between distributions. Technical Report, Department of Biostatisticsand Biomathematics, University of Alabama at Birmingham, Birmingham, USA, 1982.
  15. Inman HF, Bradley EL. The Overlapping coefficient as a measure of agreement between probability distributions and point estimation of the overlap of two normal densities. Communications in Statistics ‒ Theory and Methods. 1989;18(10):3851‒3874.
  16. Clemons TE. The overlapping coefficient for two normal probability functions with unequal variances. Unpublished Thesis, Department of Biostatistics, University of Alabama at Birmingham, Birmingham, USA, 1996.
  17. Reiser B, Faraggi D. Confidence intervals for the overlapping coefficient: the normal equal variance case. The statistician. 1999;48(3):413‒418.
  18. Clemons TE, Bradley EL. A nonparametric measure of the overlapping coefficient. Computational Statistics & Data Analysis. 2000;34(1):51‒61.
  19. Mulekar MS, Mishra SN. Confidence interval estimation of overlap: Equal means case. Computational Statistics & Data Analysis. 2000;34(2):121‒137.
  20. Al‒Saidy O, Samawi HM, Al‒Saleh MF. Inference on Overlapping Coefficients under the Wei bull Distribution: Equal Shape Parameter. ESAIM: Probability and Statistics (P&S). 2005;9:p206‒219.
  21. Schmid F, Schmidt A. Nonparametric estimation of the coefficient of overlapping‒theory and empirical application. Computational Statistics & Data Analysis. 2006;50(6):1583‒1596.
  22. Al‒Saleh MF, Samawi HM. Inference on Overlapping Coefficients in Two Exponential Populations. Journal of Modern Statistical Science. 2007;6(2):503‒516.
  23. Samawi HM, Al‒Saleh MF. Inference of Overlapping Coefficients in two Exponential Populations using Ranked set Sample. Communication of Korean of Statistical Society. 2008,15(2):147‒159.
  24. McIntyre GA. A method for unbiased selective sampling, using ranked sets. Australian Journal of Agriculture Research. 1952;39(4):385‒390.
  25. Kaur A, Patil GP, Sinha AK, et al. Ranked set sampling: Anannotated bibliography. Environmental Ecological Statistics. 1995;2: 25‒54.
  26. Patil GP, Sinha Ak, Tillie C. Ranked set sampling: a bibliography. Environmental Ecological Statistics. 1999;6:91‒98.
  27. Kaur A, Patil GP, Taillie C. Optimal allocation for symmetric distributions in ranked set sampling. Ann Inst Statist Math. 2000;52(2):239‒254.
  28. Sinha, Arun K. On some recent developments in ranked set sampling. Bulletin of Informatics and Cybernetics. 2005;37:137‒160.
  29. Anderson G, Linton O, Whang Y. Nonparametric estimation of a polarization measure. University of Toronto, The London School of Economics and Seoul National University. Discussion Paper. 2010.
  30. Wegman EJ. Nonparametric probability density estimation: a summary of available methods. Technometrics. 1972;14(3):533‒546.
  31. Wegman EJ. Density estimation. Encyclopedia of Statistical Sciences. (Vol. 2), Wiley, New York, 1982;Pp.309‒315.
  32. Van Kerm P. Adaptive kernel density estimation. The Stata Journal. 2003;3(2):148‒156.
  33. Chen EJ, Kelton WD. Empirical Evaluation of Data‒Based Density Estimation. Proceedings of the Winter Simulation Conference. 2006;333‒341.
  34. Silverman BW. Density estimation for statistics and data analysis. Chapman & Hall, London, 1986.
  35. Wand MP, Jones MC. Kernel Smoothing, Chapman & Hall, London, 1995.
  36. Ramberg JS, Schmeiser BW. An approximate method for generating Asymmetric random variables. Communications of the ACM. 1974;17:78‒82.
  37. Norozi K, C Beck C, Osthaus WA, et al. Electrical velocimetry for measuring cardiac output in children with congenital heart disease. Oxford Journals Medicine & Health BJA. 2008;100(1):88‒94.
Creative Commons Attribution License

©2014 Samawi, et al. This is an open access article distributed under the terms of the, which permits unrestricted use, distribution, and build upon your work non-commercially.