Comparing two dependent groups via quantiles

This paper considers two general ways dependent groups might be compared based on quantiles. The first compares the quantiles of the marginal distributions. The second focuses on the lower and upper quantiles of the usual difference scores. Methods for comparing quantiles have been derived that typically assume that sampling is from a continuous distribution. There are exceptions, but generally, when sampling from a discrete distribution where tied values are likely, extant methods can perform poorly, even with a large sample size. One reason is that extant methods for estimating the standard error can perform poorly. Another is that quantile estimators based on a single-order statistic, or a weighted average of two-order statistics, are not necessarily asymptotically normal. Our main result is that when using the Harrell–Davis estimator, good control over the Type I error probability can be achieved in simulations via a standard percentile bootstrap method, even when there are tied values, provided the sample sizes are not too small. In addition, the two methods considered here can have substantially higher power than alternative procedures. Using real data, we illustrate how quantile comparisons can be used to gain a deeper understanding of how groups differ.     

Global comparisons of medians and other quantiles in a one-way design when there are tied values

For J ? 2 independent groups, the article deals with testing the global hypothesis that all J groups have a common population median or identical quantiles, with an emphasis on the quartiles. Classic rank-based methods are sometimes suggested for comparing medians, but it is well known that under general conditions they do not adequately address this goal. Extant methods based on the usual sample median are unsatisfactory when there are tied values except for the special case J = 2. A variation of the percentile bootstrap used in conjunction with the Harrell–Davis quantile estimator performs well in simulations. The method is illustrated with data from the Well Elderly 2 study.     

Efficient computation of the cdf of the maximal difference between a Brownian bridge and its concave majorant

In this paper, we describe two computational methods for calculating the cumulative distribution function and the upper quantiles of the maximal difference between a Brownian bridge and its concave majorant. The first method has two different variants that are both based on a Monte Carlo approach, whereas the second uses the Gaver–Stehfest (GS) algorithm for the numerical inversion of the Laplace transform. If the former method is straightforward to implement, it is very much outperformed by the GS algorithm, which provides a very accurate approximation of the cumulative distribution as well as its upper quantiles. Our numerical work has a direct application in statistics: the maximal difference between a Brownian bridge and its concave majorant arises in connection with a nonparametric test for monotonicity of a density or regression curve on [0,1]. Our results can be used to construct very accurate rejection region for this test at a given asymptotic level.     

Robust ANCOVA: Some Small-sample Results when there are Multiple Groups and Multiple Covariates

Numerous methods have been proposed for dealing with the serious practical problems associated with the conventional analysis of covariance method, with an emphasis on comparing two groups when there is a single covariate. Recently, Wilcox (2005a: section 11.8.2) outlined a method for handling multiple covariates that allows nonlinearity and heteroscedasticity. The method is readily extended to multiple groups, but nothing is known about its small-sample properties. This paper compares three variations of the method, each method based on one of three measures of location: means, medians and 20% trimmed means. The methods based on a 20% trimmed mean or median are found to avoid Type I error probabilities well above the nominal level, but the method based on medians can be too conservative in various situations; using a 20% trimmed mean gave the best results in terms of Type I errors. The methods are based in part on a running interval smoother approximation of the regression surface. Included are comments on required sample sizes that are relevant to the so-called curse of dimensionality.     

Assessing a hypothesis test for the difference between two quantiles from independent populations

An empirical distribution function estimator for the difference of order statistics from two independent populations can be used for inference between quantiles from these populations. The inferential properties of the approach are evaluated in a simulation study where different sample sizes, theoretical distributions, and quantiles are studied. Small to moderate sample sizes, tail quantiles, and quantiles which do not coincide with the expectation of an order statistic are identified as problematic for appropriate Type I error control.     

Exact nonparametric confidence,prediction and tolerance intervals based on multi-sample Type-II right censored data

Exact nonparametric inference based on ordinary Type-II right censored samples has been extended here to the situation when there are multiple samples with Type-II censoring from a common continuous distribution. It is shown that marginally, the order statistics from the pooled sample are mixtures of the usual order statistics with multivariate hypergeometric weights. Relevant formulas are then derived for the construction of nonparametric confidence intervals for population quantiles, prediction intervals, and tolerance intervals in terms of these pooled order statistics. It is also shown that this pooled-sample approach assists in achieving higher confidence levels when estimating large quantiles as compared to a single Type-II censored sample with same number of observations from a sample of comparable size. We also present some examples to illustrate all the methods of inference developed here.     

The Gamma–Kumaraswamy distribution: An alternative to Gamma distribution

In this article, a new generalization of the Kumaraswamy distribution, namely the Gamma–Kumaraswamy distribution, is defined and studied. Various properties of the Gamma–Kumaraswamy are obtained. The structural analysis of the distribution in this article includes limiting behavior, mode, quantiles, moments, skewness, kurtosis, Shannon’s entropy, and order statistics. The method of maximum likelihood estimation is proposed for estimating the model parameters. For illustrative purposes, two real datasets are analyzed as application of the Gamma–Kumaraswamy distribution.     

Estimation for Type II domain of attraction based on the W statistic

The paper presents an estimating equation approach to the estimation of high quantiles of a distribution in the Type II domain of attraction based on the k largest order statistics. The estimators are shown to be consistent. The method fits neatly into a general scheme for estimating high quantiles irrespective of the domain of attraction, which includes Wang's approach to optimally choosing k .     

Pareto Poisson–Lindley distribution with applications

A new lifetime distribution is introduced based on compounding Pareto and Poisson–Lindley distributions. Several statistical properties of the distribution are established, including behavior of the probability density function and the failure rate function, heavy- and long-right tailedness, moments, the Laplace transform, quantiles, order statistics, moments of residual lifetime, conditional moments, conditional moment generating function, stress–strength parameter, Rényi entropy and Song's measure. We get maximum-likelihood estimators of the distribution parameters and investigate the asymptotic distribution of the estimators via Fisher's information matrix. Applications of the distribution using three real data sets are presented and it is shown that the distribution fits better than other related distributions in practical uses.     

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].     

Distribution-free confidence intervals for quantiles in small samples

Estimators for quantiles based on linear combinations of order statistics have been proposed by Harrell and Davis(1982) and kaigh and Lachenbruch (1982). Both estimators have been demonstrated to be at least as efficient for small sample point estimation as an ordinary sample quantile estimator based on one or two order statistics: Distribution-free confidence intervals for quantiles can be constructed using either of the two approaches. By means of a simulation study, these confidence intervals have been compared with several other methods of constructing confidence intervals for quantiles in small samples. For the median, the Kaigh and Lachenbruch method performed fairly well. For other quantiles, no method performed better than the method which uses pairs of order statistics.     

A characterization of the Burr Type III and Type XII distributions through the method of percentiles and the Spearman correlation

A characterization of Burr Type III and Type XII distributions based on the method of percentiles (MOP) is introduced and contrasted with the method of (conventional) moments (MOM) in the context of estimation and fitting theoretical and empirical distributions. The methodology is based on simulating the Burr Type III and Type XII distributions with specified values of medians, inter-decile ranges, left-right tail-weight ratios, tail-weight factors, and Spearman correlations. Simulation results demonstrate that the MOP-based Burr Type III and Type XII distributions are substantially superior to their (conventional) MOM-based counterparts in terms of relative bias and relative efficiency.     

Some Results on Comparing the Quantiles of Dependent Groups

Recently, Lombard derived an extension of the Doksum–Sievers shift function to dependent groups. This article suggests using a particular numerical method for determining the critical value, reports on the ability of the method to control the probability of a Type I error when sample sizes are small, and it provides comparisons with methods aimed at comparing deciles. It is found that for continuous distributions, Lombard's method performs well and in particular has high power relative to the other two methods considered. But when tied values can occur, now it can have relatively poor power; a method based on the Harrell-Davis estimator is found to give more satisfactory results.     

QUANTILES OF SUMS AND EXPECTED VALUES OF ORDERED SUMS

Watson & Gordon (1986) investigated the relationship between the quantiles of a sum of independent continuous random variables and the sum of the individual quantiles. In this note some further results are obtained. Also corresponding relationships are developed for the expected values of the order statistics of a sum, and for the sum of the expected values of the individual order statistics.     

Parameter and quantile estimation for the three-parameter lognormal distribution based on statistics invariant to unknown location

Lognormal distribution is one of the popular distributions used for modelling positively skewed data, especially those encountered in economic and financial data. In this paper, we propose an efficient method for the estimation of parameters and quantiles of the three-parameter lognormal distribution, which avoids the problem of unbounded likelihood, by using statistics that are invariant to unknown location. Through a Monte Carlo simulation study, we then show that the proposed method performs well compared to other prominent methods in terms of both bias and mean-squared error. Finally, we present two illustrative examples.     

A New Goodness-of-Fit Test for Censored Data with an Application in Monitoring Processes

In this article, we propose a new goodness-of-fit test for Type I or Type II censored samples from a completely specified distribution. This test is a generalization of Michael's test for censored data, which is based on the empirical distribution and a variance stabilizing transformation. Using Monte Carlo methods, the distributions of the test statistics are analyzed under the null hypothesis. Tables of quantiles of these statistics are also provided. The power of the proposed test is studied and compared to that of other well-known tests also using simulation. The proposed test is more powerful in most of the considered cases. Acceptance regions for the PP, QQ, and Michael's stabilized probability plots are derived, which enable one to visualize which data contribute to the decision of rejecting the null hypothesis. Finally, an application in quality control is presented as illustration.     

Hybrid censoring schemes with exponential failure distribution

The mixture of Type I and Type I1 censoring schemes, called the hybrid censoring, is quite important in life–testing experiments. Epstein(1954, 1960) introduced this testing scheme and proposed a two–sided confidence interval to estimate the mean lifetime, θ, when the underlying lifetime distribution is assumed to be exponential. There are some two–sided confidence intervals and credible intervals proposed by Fairbanks et al. (1982) and Draper and Guttman (1987) respectively. In this paper we obtain the exact two–sided confidence interval of θ following the approach of Chen and Bhattacharya (1988). We also obtain the asymptotic confidence intervals in the Hybrid censoring case. It is important to observe that the results for Type I and Type II censoring schemes can be obtained as particular cases of the Hybrid censoring scheme. We analyze one data set and compare different methods by Monte Carlo simulations.     

Estimating pareto quantiles using two order statistics

Asymptotically best linear unbiased estimators of the population quantiles for the location-scale Pareto distribution with fixed shape parameter are obtained using two suitably chosen order statistics. Formulae for the appropriate order statistics, coefficients, variances, and asymptotic relative efficiencies (relative to the usual non-parametric estimator for quantiles) are given     

Nonparametric estimation of extreme conditional quantiles

The estimation of extreme conditional quantiles is an important issue in different scientific disciplines. Up to now, the extreme value literature focused mainly on estimation procedures based on independent and identically distributed samples. Our contribution is a two-step procedure for estimating extreme conditional quantiles. In a first step nonextreme conditional quantiles are estimated nonparametrically using a local version of [Koenker, R. and Bassett, G. (1978). Regression quantiles. Econometrica, 46, 33–50.] regression quantile methodology. Next, these nonparametric quantile estimates are used as analogues of univariate order statistics in procedures for extreme quantile estimation. The performance of the method is evaluated for both heavy tailed distributions and distributions with a finite right endpoint using a small sample simulation study. A bootstrap procedure is developed to guide in the selection of an optimal local bandwidth. Finally the procedure is illustrated in two case studies.     

NONPARAMETRIC TWO SAMPLE TESTS BASED ON AREA STATISTICS

ABSTRACT

Area statistics are sample versions of areas occurring in a probability plot of two distribution functions F and G. This paper presents a unified basis for five statistics of this type. They can be used for various testing problems in the framework of the two sample problem for independent observations, such as testing equality of distributions against inequality or testing stochastic dominance of distributions in one or either direction against nondominance. Though three of the statistics considered have already been suggested in literature, two of them are new and deserve our interest. The finite sample distributions of the statistics (under F=G) can be calculated via recursion formulae. Two tables with critical values of the new statistics are included. The asymptotic distribution of the properly normalized versions of the area statistics are functionals of the Brownian bridge. The distribution functions and quantiles thereof are obtained by Monte Carlo simulation. Finally, the power functions of the two new tests based on area statistics are compared to the power functions of the tests based on the corresponding supremum statistics, i.e., statistics of the Kolmogorov–Smirnov type.