Combining multi‐observer information in partially rank‐ordered judgment post‐stratified and ranked set samples

This paper develops two sampling designs to create artificially stratified samples. These designs use a small set of experimental units to determine their relative ranks without measurement. In each set, the units are ranked by all available observers (rankers), with ties whenever the units cannot be ranked with high confidence. The rankings from all the observers are then combined in a meaningful way to create a single weight measure. This weight measure is used to create judgment strata in both designs. The first design constructs the strata through judgment post‐stratification after the data has been collected. The second design creates the strata before any measurements are made on the experimental units. The paper constructs estimators and confidence intervals, and develops testing procedures for the mean and median of the underlying distribution based on these sampling designs. We show that the proposed sampling designs provide a substantial improvement over their competitor designs in the literature. The Canadian Journal of Statistics 41: 304–324; 2013 © 2013 Statistical Society of Canada     

Inference in the presence of ranking error in ranked set sampling

The author proposes inference techniques for ranked set sample data in the presence of judgment ranking errors. He bases his analysis on the models of Bohn & Wolfe (1994) and Frey (2007a, b), of which parameters are estimated by minimizing a distance measure. He then uses the fitted models to calibrate confidence intervals and tests. He shows the validity of his approach through simulation and illustrates its application through the construction of distribution‐free confidence intervals for the median area of apple tree leaves covered by a spray.     

Pseudo‐empirical likelihood ratio confidence intervals for complex surveys

The authors show how an adjusted pseudo‐empirical likelihood ratio statistic that is asymptotically distributed as a chi‐square random variable can be used to construct confidence intervals for a finite population mean or a finite population distribution function from complex survey samples. They consider both non‐stratified and stratified sampling designs, with or without auxiliary information. They examine the behaviour of estimates of the mean and the distribution function at specific points using simulations calling on the Rao‐Sampford method of unequal probability sampling without replacement. They conclude that the pseudo‐empirical likelihood ratio confidence intervals are superior to those based on the normal approximation, whether in terms of coverage probability, tail error rates or average length of the intervals.     

Coverage‐adjusted Confidence Intervals for a Binomial Proportion

We consider the classic problem of interval estimation of a proportion p based on binomial sampling. The ‘exact’ Clopper–Pearson confidence interval for p is known to be unnecessarily conservative. We propose coverage adjustments of the Clopper–Pearson interval that incorporate prior or posterior beliefs into the interval. Using heatmap‐type plots for comparing confidence intervals, we show that the coverage‐adjusted intervals have satisfying coverage and shorter expected lengths than competing intervals found in the literature.     

Prediction of order statistics and record values based on ordered ranked set sampling

In this paper, we consider two-sample prediction problems. First, based on ordered ranked set sampling (ORSS) introduced by Balakrishnan and Li [Ordered ranked set samples and applications to inference. Ann Inst Statist Math. 2006;58:757–777], we obtain prediction intervals for order statistics from a future sample and compare the results with the one based on the usual-order statistics. Next, we construct prediction intervals for record values from a future sequence based on ORSS and compare the results with the one based on an another independent record sequence developed recently by Ahmadi and Balakrishnan [Prediction of order statistics and record values from two independent sequences. Statistics. 2010;44:417–430].     

Empirical likelihood‐based inferences for a low income proportion

Low income proportion is an important index in comparisons of poverty in countries around the world. The stability of a society depends heavily on this index. An accurate and reliable estimation of this index plays an important role for government's economic policies. In this paper, the authors study empirical likelihood‐based inferences for a low income proportion under the simple random sampling and stratified random sampling designs. It is shown that the limiting distributions of the empirical likelihood ratios for the low income proportion are the scaled chi‐square distributions. The authors propose various empirical likelihood‐based confidence intervals for the low income proportion. Extensive simulation studies are conducted to evaluate the relative performance of the normal approximation‐based interval, bootstrap‐based intervals, and the empirical likelihood‐based intervals. The proposed methods are also applied to analyzing a real economic survey income dataset. The Canadian Journal of Statistics 39: 1–16; 2011 ©2011 Statistical Society of Canada     

Empirical Likelihood Inference for the Rao‐Hartley‐Cochran Sampling Design

The Hartley‐Rao‐Cochran sampling design is an unequal probability sampling design which can be used to select samples from finite populations. We propose to adjust the empirical likelihood approach for the Hartley‐Rao‐Cochran sampling design. The approach proposed intrinsically incorporates sampling weights, auxiliary information and allows for large sampling fractions. It can be used to construct confidence intervals. In a simulation study, we show that the coverage may be better for the empirical likelihood confidence interval than for standard confidence intervals based on variance estimates. The approach proposed is simple to implement and less computer intensive than bootstrap. The confidence interval proposed does not rely on re‐sampling, linearization, variance estimation, design‐effects or joint inclusion probabilities.     

Nonparametric Prediction Intervals Based on Ranked Set Samples

In this article, we develop nonparametric prediction intervals based on generalized ranked set samples using conditional as well as unconditional approaches. The predictions are developed for order statistics from a future sample as well as for order statistics from a future balanced ranked set sample. The effects of ranking errors on the coverage probabilities of these prediction intervals are also examined.     

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.     

Optimal sign test for quantiles in ranked set samples

This paper considers the one-sample sign test for population quantiles in general ranked set sampling, and proposes a weighted sign test because observations with different ranks are not identically distributed. It is shown analytically that optimal weight always improves the Pitman efficiency for all distributions. For each quantile, the sampling allocation that maximizes the sign test efficacy is identified and shown to not depend on the population distribution. Moreover, distribution-free confidence intervals for quantiles based on ordered values of optimal ranked set samples are discussed.     

An improved ranked set two‐sample mann‐whitney‐wilcoxon test

The authors present an improved ranked set two‐sample Mann‐Whitney‐Wilcoxon test for a location shift between samples from two distributions F and G. They define a function that measures the amount of information provided by each observation from the two samples, given the actual joint ranking of all the units in a set. This information function is used as a guide for improving the Pitman efficacy of the Mann‐Whitney‐Wilcoxon test. When the underlying distributions are symmetric, observations at their mode(s) must be quantified in order to gain efficiency. Analogous results are provided for asymmetric distributions.     

Confidence intervals for the mean of a population containing many zero values under unequal‐probability sampling

In many applications, a finite population contains a large proportion of zero values that make the population distribution severely skewed. An unequal‐probability sampling plan compounds the problem, and as a result the normal approximation to the distribution of various estimators has poor precision. The central‐limit‐theorem‐based confidence intervals for the population mean are hence unsatisfactory. Complex designs also make it hard to pin down useful likelihood functions, hence a direct likelihood approach is not an option. In this paper, we propose a pseudo‐likelihood approach. The proposed pseudo‐log‐likelihood function is an unbiased estimator of the log‐likelihood function when the entire population is sampled. Simulations have been carried out. When the inclusion probabilities are related to the unit values, the pseudo‐likelihood intervals are superior to existing methods in terms of the coverage probability, the balance of non‐coverage rates on the lower and upper sides, and the interval length. An application with a data set from the Canadian Labour Force Survey‐2000 also shows that the pseudo‐likelihood method performs more appropriately than other methods. The Canadian Journal of Statistics 38: 582–597; 2010 © 2010 Statistical Society of Canada     

New imperfect rankings models for ranked set sampling

Ranked set sampling is a sampling approach that leads to improved statistical inference in situations where the units to be sampled can be ranked relative to each other prior to formal measurement. This ranking may be done either by subjective judgment or according to an auxiliary variable, and it need not be completely accurate. In fact, results in the literature have shown that no matter how poor the quality of the ranking, procedures based on ranked set sampling tend to be at least as efficient as procedures based on simple random sampling. However, efforts to quantify the gains in efficiency for ranked set sampling procedures have been hampered by a shortage of available models for imperfect rankings. In this paper, we introduce a new class of models for imperfect rankings, and we provide a rigorous proof that essentially any reasonable model for imperfect rankings is a limit of models in this class. We then describe a specific, easily applied method for selecting an appropriate imperfect rankings model from the class.     

Inference in Targeted Group‐Sequential Covariate‐Adjusted Randomized Clinical Trials

This article is devoted to the construction and asymptotic study of adaptive, group‐sequential, covariate‐adjusted randomized clinical trials analysed through the prism of the semiparametric methodology of targeted maximum likelihood estimation. We show how to build, as the data accrue group‐sequentially, a sampling design that targets a user‐supplied optimal covariate‐adjusted design. We also show how to carry out sound statistical inference based on such an adaptive sampling scheme (therefore extending some results known in the independent and identically distributed setting only so far), and how group‐sequential testing applies on top of it. The procedure is robust (i.e. consistent even if the working model is mis‐specified). A simulation study confirms the theoretical results and validates the conjecture that the procedure may also be efficient.     

Likelihood‐based Inference with Missing Data Under Missing‐at‐Random

Likelihood‐based inference with missing data is challenging because the observed log likelihood is often an (intractable) integration over the missing data distribution, which also depends on the unknown parameter. Approximating the integral by Monte Carlo sampling does not necessarily lead to a valid likelihood over the entire parameter space because the Monte Carlo samples are generated from a distribution with a fixed parameter value. We consider approximating the observed log likelihood based on importance sampling. In the proposed method, the dependency of the integral on the parameter is properly reflected through fractional weights. We discuss constructing a confidence interval using the profile likelihood ratio test. A Newton–Raphson algorithm is employed to find the interval end points. Two limited simulation studies show the advantage of the Wilks inference over the Wald inference in terms of power, parameter space conformity and computational efficiency. A real data example on salamander mating shows that our method also works well with high‐dimensional missing data.     

NONPARAMETRIC META‐ANALYSIS FOR MINIMAL‐REPAIR SYSTEMS

We consider the situation that repair times of several identically structured technical systems are observed. As an example of such data we discuss the Boeing air conditioner data, consisting of successive failures of the air conditioning system of each member of a fleet of Boeing jet airplanes. The repairing process is assumed to be performed according to a minimal‐repair strategy. This reflects the idea that only those operations are accomplished that are absolutely necessary to restart the system after a failure. The ‘after‐repair‐state’ of the system is the same as it was shortly before the failure. Clearly, the observed repair times contain valuable information about the repair times of an identically structured system put into operation in the future. Thus, for statistical analysis and prediction, it is certainly favourable to take into account all repair times from each system. The resulting pooled sample is used to construct nonparametric prediction intervals for repair times of a future minimal‐repair system. To illustrate our results we apply them to the above‐mentioned data set. As expected, the maximum coverage probabilities of prediction intervals based on two samples exceed those based on one sample. We show that the relative gain for a two‐sample prediction over a one‐sample prediction can be substantial. One of the advantages of the present approach is that it allows nonparametric prediction intervals to be constructed directly. This provides a beneficial alternative to existing nonparametric methods for minimal‐repair systems that construct prediction intervals via the asymptotic distribution of quantile estimators. Moreover, the prediction intervals presented here are exact regardless of the sample size.     

On Entropy-Based Test of Exponentiality in Ranked Set Sampling

When the sampling units can be easily ranked than quantified, ranked set sampling (RSS) is a viable alternative to the traditional simple random sampling (SRS). Much effort has been made for modifying basic RSS protocol with the aim of deriving more efficient estimators of the population attributes. Entropy has been seminal in developing measures of distributional disparities as a tool for statistical inference. This article is concerned with testing exponentiality based on sample entropy under some RSS-based designs. A simulation study shows that the proposed tests possess good power properties against several alternatives as compared with the ordinary test based on SRS.     

An affine‐invariant multivariate sign test for cluster correlated data

The author presents a multivariate location model for cluster correlated observations. He proposes an affine‐invariant multivariate sign statistic for testing the value of the location parameter. His statistic is an adaptation of that proposed by Randles (2000). The author shows, under very mild conditions, that his test statistic is asymptotically distributed as a chi‐squared random variable under the null hypothesis. In particular, the test can be used for skewed populations. In the context of a general multivariate normal model, the author obtains values of his test's Pitman asymptotic efficiency relative to another test based on the overall average. He shows that there is an improvement in the relative performance of the new test as soon as intra‐cluster correlation is present Even in the univariate case, the new test can be very competitive for Gaussian data. Furthermore, the statistic is easy to compute, even for large dimensional data. The author shows through simulations that his test performs well compared to the average‐based test. He illustrates its use with real data.     

Estimation of the means of the bivariate normal using moving extreme ranked set sampling with concomitant variable

The estimation of the means of the bivariate normal distribution, based on a sample obtained using a modification of the moving extreme ranked set sampling technique (MERSS) is considered. The modification involves using a concomitant random variable. Nonparametric-type methods as well as the maximum likelihood estimation are considered. The estimators obtained are compared to their counterparts based on simple random sampling (SRS). It appears that the suggested estimators are more efficient. Also, MERSS with concomitant variable is easier to use in practice than the usual ranked set sampling (RSS) with concomitant variable. The issue of robustness of the procedure is addressed. Real trees data set is used for illustration.     

A Further Study on Reliable Life Estimation under Ranked Set Sampling

This article considers nonparametric estimation of reliable life based on ranked set sampling and its properties. It is proven analytically that the large sample efficiency of the reliable life estimator under the balanced ranked set sampling is higher than that under the simple random sampling of the same size, but the relative efficiency damps away as the reliable life moves away from the median on both directions. To improve the efficiency for the estimation of extreme reliable life, we then propose a reliable life estimator under a modified ranked set sampling protocol, its strong consistency and asymptotic normality are established. The proposed sampling is shown to be superior to the balanced ranked set sampling, and the relative advantage improves as the reliable life moves away from median. Finally, results of simulation studies for small sample as well as an application to a real data set are presented to illustrate some of the theoretical findings.