An Affine-Invariant Generalization of the Wilcoxon Signed-Rank Test for the Bivariate Location Problem

This paper proposes an affine‐invariant test extending the univariate Wilcoxon signed‐rank test to the bivariate location problem. It gives two versions of the null distribution of the test statistic. The first version leads to a conditionally distribution‐free test which can be used with any sample size. The second version can be used for larger sample sizes and has a limiting χ₂² distribution under the null hypothesis. The paper investigates the relationship with a test proposed by Jan & Randles (1994). It shows that the Pitman efficiency of this test relative to the new test is equal to 1 for elliptical distributions but that the two tests are not necessarily equivalent for non‐elliptical distributions. These facts are also demonstrated empirically in a simulation study. The new test has the advantage of not requiring the assumption of elliptical symmetry which is needed to perform the asymptotic version of the Jan and Randles test.     

Quantile inference based on partially rank-ordered set samples

This paper develops statistical inference for population quantiles based on a partially rank-ordered set (PROS) sample design. A PROS sample design is similar to a ranked set sample with some clear differences. This design first creates partially rank-ordered subsets by allowing ties whenever the units in a set cannot be ranked with high confidence. It then selects a unit for full measurement at random from one of these partially rank-ordered subsets. The paper develops a point estimator, confidence interval and hypothesis testing procedure for the population quantile of order p. Exact, as well as asymptotic, distribution of the test statistic is derived. It is shown that the null distribution of the test statistic is distribution-free, and statistical inference is reasonably robust against possible ranking errors in ranking process.     

Partitioning Anderson's statistic for tied data

Anderson (Biometrics 15 (1959) 582) proposed a χ²-type statistic for the nonparametric analysis of a randomized blocks design with no ties in the data. In this paper, we propose an Anderson statistic that allows for ties in the data. We show that the asymptotic distribution of the statistic under the null hypothesis of no treatment effect is a χ² distribution. Under weak assumptions on the tie structure it is shown that the degrees of freedom for the asymptotic distribution is unchanged compared to the untied case. An extended analysis based on a partition of the statistic into independent components is suggested. The first component is shown to equal the Friedman rank statistic corrected for ties. The subsequent components allow for the detection of dispersion effects, higher order effects and differences in distribution. A simulation study is given and the new analysis is applied to a sensory evaluation data set.     

Fast computation of the exact null distribution of Spearman's ρ and Page's L statistic for samples with and without ties

We present a new algorithm for computing the exact null distribution of the Spearman rank correlation statistic ρ, which also works in the case of ties. The algorithm is based on symmetries in the representation of the probability generating function as a permanent with monomial entries. We present new critical values for sample sizes 19⩽n⩽22. Finally, we show how to derive the exact null distribution of Page's L statistic from the null distribution of ρ.     

More powerful sign test using median ranked set sample: Finite sample power comparison

Sign test using median ranked set samples (MRSS) is introduced and investigated. We show that, this test is more powerful than the sign tests based on simple random sample (SRS) and ranked set sample (RSS) for finite sample size. It is found that, when the set size of MRSS is odd, the null distribution of the MRSS sign test is the same as the sign test obtained by using SRS. The exact null distributions and the power functions, in case of finite sample sizes, of these tests are derived. Also, the asymptotic distribution of the MRSS sign tests are derived. Numerical comparison of the MRSS sign test power with the power of the SRS sign test and the RSS sign test is given. Illustration of the procedure, using real data set of bilirubin level in Jaundice babies who stay in neonatal intensive care is introduced.     

Unbiased Estimation of the Distribution Function of an Exponential Population Using Order Statistics with Application in Ranked Set Sampling

In this paper we consider the problem of unbiased estimation of the distribution function of an exponential population using order statistics based on a random sample. We present a (unique) unbiased estimator based on a single, say ith, order statistic and study some properties of the estimator for i = 2. We also indicate how this estimator can be utilized to obtain unbiased estimators when a few selected order statistics are available as well as when the sample is selected following an alternative sampling procedure known as ranked set sampling. It is further proved that for a ranked set sample of size two, the proposed estimator is uniformly better than the conventional nonparametric unbiased estimator, further, for a general sample size, a modified ranked set sampling procedure provides an unbiased estimator uniformly better than the conventional nonparametric unbiased estimator based on the usual ranked set sampling procedure.     

Multistage ranked set sampling

The superiority of using ranked set sampling, for estimating the mean of a population, over simple random sampling, is well established. This technique is useful when visual ordering of a small set of size (m) can be done easily and fairly accurately, but exact measurement of an observation is difficult and expensive. It is noted that for many distributions, an increase in the efficiency of ranked set sampling can be achieved by increasing the set size m. However, in practice, m should be kept very small so that visual ranking errors will not destroy the gain in efficiency. In this paper, multistage ranked set sampling is considered as a generalization of ranked set sampling, that results in an increase of the efficiency for fixed value of m. Steady state efficiency, the limiting efficiency as the number of stages approaches infinity, varies from one distribution to another. It is shown that this efficiency is always larger than 1, close to m² for symmetric distributions and equal to m² for the uniform distribution. Some real applications of the technique are discussed. Data on olive yield of olive trees is collected to illustrate the technique.     

A Bayesian analysis for the Wilcoxon signed-rank statistic

A Bayesian analysis is provided for the Wilcoxon signed-rank statistic (T⁺). The Bayesian analysis is based on a sign-bias parameter φ on the (0, 1) interval. For the case of a uniform prior probability distribution for φ and for small sample sizes (i.e., 6 ? n ? 25), values for the statistic T⁺ are computed that enable probabilistic statements about φ. For larger sample sizes, approximations are provided for the asymptotic likelihood function P(T⁺|φ) as well as for the posterior distribution P(φ|T⁺). Power analyses are examined both for properly specified Gaussian sampling and for misspecified non Gaussian models. The new Bayesian metric has high power efficiency in the range of 0.9–1 relative to a standard t test when there is Gaussian sampling. But if the sampling is from an unknown and misspecified distribution, then the new statistic still has high power; in some cases, the power can be higher than the t test (especially for probability mixtures and heavy-tailed distributions). The new Bayesian analysis is thus a useful and robust method for applications where the usual parametric assumptions are questionable. These properties further enable a way to do a generic Bayesian analysis for many non Gaussian distributions that currently lack a formal Bayesian model.     

Computations of Bayesian t-values for multiple comparisons

Results from a simulation study of the power of eight statistics for testing that a sample is form a uniform distribution on the unit interval are reported. Power is given for each statistic against four classes if alternatives. The statistics studied include the discrete Pearson chi-square with ten and twenty cells, X² ₁₀ and X² ₂₀; Kolmogorov-smirov, D; Cramer-Von Mises, W²; Watson, U²; Anderson-Darling, A; Greenwood. G;and a new statistic called O A modified form of each of these statistic is also studied by first transforming the sample using a transformation given by Durbin. On the basis of the results observed in this study, the Watson U² statistic is recommended as a general test for uniformity.     

Wilcoxon's signed‐rank statistic: what null hypothesis and why it matters

In statistical literature, the term ‘signed‐rank test’ (or ‘Wilcoxon signed‐rank test’) has been used to refer to two distinct tests: a test for symmetry of distribution and a test for the median of a symmetric distribution, sharing a common test statistic. To avoid potential ambiguity, we propose to refer to those two tests by different names, as ‘test for symmetry based on signed‐rank statistic’ and ‘test for median based on signed‐rank statistic’, respectively. The utility of such terminological differentiation should become evident through our discussion of how those tests connect and contrast with sign test and one‐sample t‐test. Published 2014. This article is a U.S. Government work and is in the public domain in the USA.     

Estimation of the shape and scale parameters of Pareto distribution using ranked set sampling

In the current paper, the estimation of the shape and location parameters α and c, respectively, of the Pareto distribution will be considered in cases when c is known and when both are unknown. Simple random sampling (SRS) and ranked set sampling (RSS) will be used, and several traditional and ad hoc estimators will be considered. In addition, the estimators of α, when c is known using an RSS version based on the order statistic that maximizes the Fisher information for a fixed set size, will be considered. These estimators will be compared in terms of their biases and mean square errors. The estimators based on RSS can be real competitors against those based on SRS.     

Some Pitman closeness properties pertinent to symmetric populations

In this paper, we focus on Pitman closeness probabilities when the estimators are symmetrically distributed about the unknown parameter θ. We first consider two symmetric estimators θ?₁ and θ?₂ and obtain necessary and sufficient conditions for θ?₁ to be Pitman closer to the common median θ than θ?₂. We then establish some properties in the context of estimation under the Pitman closeness criterion. We define Pitman closeness probability which measures the frequency with which an individual order statistic is Pitman closer to θ than some symmetric estimator. We show that, for symmetric populations, the sample median is Pitman closer to the population median than any other independent and symmetrically distributed estimator of θ. Finally, we discuss the use of Pitman closeness probabilities in the determination of an optimal ranked set sampling scheme (denoted by RSS) for the estimation of the population median when the underlying distribution is symmetric. We show that the best RSS scheme from symmetric populations in the sense of Pitman closeness is the median and randomized median RSS for the cases of odd and even sample sizes, respectively.     

Nonparametric confidence intervals for location in time series data

We extend the confidence interval construction procedure for location for symmetric iid data using the one-sample Wilcoxon signed rank statistic (T⁺) to stationary time series data. We propose a normal approximation procedure when explicit knowledge of the underlying dependence structure/distribution is unknown. By conducting extensive simulations from linear and nonlinear time series models, we show that the extended procedure is a strong contender for use in the construction of confidence intervals in time series analysis. Finally we demonstrate real application implementations in two case studies.     

A note on the nonparametric test based on theL 1-version of the Cramér-von Mises statistic

We consider the test based on theL ₁-version of the Cramér-von Mises statistic for the nonparametric two-sample problem. Some quantiles of the exact distribution under H₀ of the test statistic are computed for small sample sizes. We compare the test in terms of power against general alternatives to other two-sample tests, namely the Wilcoxon rank sum test, the Smirnov test and the Cramér-von Mises test in the case of unbalanced small sample sizes. The computation of the power is rather complicated when the sample sizes are unequal. Using Monte Carlo power estimates it turns out that the Smirnov test is more sensitive to non stochastically ordered alternatives than the new test. And under location-contamination alternatives the power estimates of the new test and of the competing tests are equal.     

One-and two-sample sign tests for ranked set sample selective designs

Statistical inference based on a ranked set sample depends very much on the location of the quantified observations. A selective design which determines the location of the quantified observations in a ranked set sample is introduced. The paper investigates the effects of selective designs on one and two sample sign test statistics. The Pitman efficiencies of one- and two sample sign tests are calculated for selective designs and compared with ranked set samples of the same size. If the design quantifies observations at the center points, then the proposed procedure is superior to a ranked set sample of the same size in the sense of Pitman efficiency. Some practical problems are addressed for the two-sample sign test.     

Sign test based on k-tuple general ranked set samples

ABSTRACT

The sign test based on the k-tuple ranked set samples is discussed here. We first derive the distribution of the k-tuple ranked set sample sign test statistic, and then the asymptotic distribution is also obtained. We then compare its performance with its counterparts based on simple random sample and classical ranked set sample. The asymptotic relative efficiency and the power are then derived. Finally, the effect of imperfect ranking on the procedure is assessed.     

Asymptotic Expansion of the Non Null Distribution of the Two-Sample t-Statistic Under Non Normality with Application in Power Comparison

This article studies the non null distribution of the two-sample t-statistic, or Welch statistic, under non normality. The asymptotic expansion of the non null distribution is derived up to n ^?1, where n is the pooled sample size, under general conditions. It is used to compare the power with that obtained by normal theory method. A simple technique is recommended to achieve more power through a monotone transformation in practice.     

Some simple nonparametric methods to test for perfect ranking in ranked set sampling

A lot of research on ranked set sampling (RSS) is based on the assumption that the ranking is perfect. Hence, it is necessary to develop some tests that could be used to validate this assumption of perfect ranking. In this paper, we introduce some simple nonparametric methods for this purpose. We specifically define three test statistics, _Nk,_Sk$N_k,S_k$ and _Ak$A_k$, based on one-cycle RSS, which are all associated with the ordered ranked set sample (ORSS). We then derive the exact null distributions and exact power functions of all these tests. Next, by using the sum or the maximum of each statistic over all cycles, we propose six test statistics for the case of multi-cycle RSS. We compare the performance of all these tests with that of the Kolmogorov–Smirnov test statistic proposed earlier by Stokes and Sager [1988. Characterization of a ranked-set sample with application to estimating distribution functions. J. Amer. Statist. Assoc. 83, 35–42] and display that all proposed test statistics are more powerful. Finally, we present an example to illustrate the test procedures discussed here.     

A useful property of order statistics and its application to UMVU estimation

Consider a random sample of sizen drawn from a continuous parent distributionF. A basic and useful known property associated with such sample is the following: the conditional distribution of thej ^th order statistic given a valuet of thei ^th order statistics, (j>i), coincides with the distribution of the(j?i) ^th order statistic in a sample of size (n?i) drawn from the parent distributionF truncated at the left att. In this article we mention some applications of this property, and provide a new application to the construction of an Uniformly Minimum Variance Unbiased (UMVU) estimator in the case of two-truncation parameters family of distributions.     

Ordered Alternatives and the 1 1/2-Tail Test

Under proper conditions, two independent tests of the null hypothesis of homogeneity of means are provided by a set of sample averages. One test, with tail probability P ₁, relates to the variation between the sample averages, while the other, with tail probability P ₂, relates to the concordance of the rankings of the sample averages with the anticipated rankings under an alternative hypothesis. The quantity G = P ₁ P ₂ is considered as the combined test statistic and, except for the discreteness in the null distribution of P ₂, would correspond to the Fisher statistic for combining probabilities. Illustration is made, for the case of four means, on how to get critical values of G or critical values of P ₁ for each possible value of P ₂, taking discreteness into account. Alternative measures of concordance considered are Spearman's ρ and Kendall's τ. The concept results, in the case of two averages, in assigning two-thirds of the test size to the concordant tail, one-third to the discordant tail.