Comparison of nonparametric point estimators for interquantile differences in moderate sized samples

Non-parametric estimators for interquantile differences in moderate sized samples are identified, including a simple method based on sample quantiles, an extension of the Harrell and Davis (1982) method, an extension of the Kaigh and Lachenbruch (1982) approach, and methods based on Yang (1985). For those estimators compared via a simulation study (Harrell-Davis and sample quantile), using relative bias and relative efficiency as judgment criteria, the Harrell-Davis method is judged to be preferable for interquantile estimation with data sets of moderate size.     

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.     

Confidence intervals based on L-moments for quantiles of the GP and GEV distributions with application to market-opening asset prices data

In a ground-breaking paper published in 1990 by the Journal of the Royal Statistical Society, J.R.M. Hosking defined the L-moment of a random variable as an expectation of certain linear combinations of order statistics. L-moments are an alternative to conventional moments and recently they have been used often in inferential statistics. L-moments have several advantages over the conventional moments, including robustness to the the presence of outliers, which may lead to more accurate estimates in some cases as the characteristics of distributions. In this contribution, asymptotic theory and L-moments are used to derive confidence intervals of the population parameters and quantiles of the three-parametric generalized Pareto and extreme-value distributions. Computer simulations are performed to determine the performance of confidence intervals for the population quantiles based on L-moments and to compare them to those obtained by traditional estimation techniques. The results obtained show that they perform well in comparison to the moments and maximum likelihood methods when the interest is in higher quantiles, or even best. L-moments are especially recommended when the tail of the distribution is rather heavier and the sample size is small. The derived intervals are applied to real economic data, and specifically to market-opening asset prices.     

QUANTILE ESTIMATION AND PROBABILITY COVERAGES

When estimating population quantiles via a random sample from an unknown continuous distribution function it is well known that a pair of order statistics may be used to set a confidence interval for any single desired, population quantile. In this paper the technique is generalized so that more than one pair of order statistics may be used to obtain simultaneous confidence intervals for the various quantiles that might be required. The generalization immediately extends to the problem of obtaining interval estimates for quantile intervals. Distributions of the ordered and unordered probability coverages of these confidence intervals are discussed as are the associated distributions of linear combinations of the coverages.     

Saddlepoint Approximations for the Difference of Order Statistics and Studentized Sample Quantiles

For a sample from a given distribution the difference of two order statistics and the Studentized quantile are statistics whose distribution is needed to obtain tests and confidence intervals for quantiles and quantile differences. This paper gives saddlepoint approximations for densities and saddlepoint approximations of the Lugannani–Rice form for tail probabilities of these statistics. The relative errors of the approximations are n ⁻¹ uniformly in a neighbourhood of the parameters and this uniformity is global if the densities are log-concave.     

Nonparametric asymptotic confidence intervals for extreme quantiles

In this paper, we propose new asymptotic confidence intervals for extreme quantiles, that is, for quantiles located outside the range of the available data. We restrict ourselves to the situation where the underlying distribution is heavy-tailed. While asymptotic confidence intervals are mostly constructed around a pivotal quantity, we consider here an alternative approach based on the distribution of order statistics sampled from a uniform distribution. The convergence of the coverage probability to the nominal one is established under a classical second-order condition. The finite sample behavior is also examined and our methodology is applied to a real dataset.     

Monte Carlo Approximations of the Quantiles of a Sample Statistic

We consider in this article the problem of numerically approximating the quantiles of a sample statistic for a given population, a problem of interest in many applications, such as bootstrap confidence intervals. The proposed Monte Carlo method can be routinely applied to handle complex problems that lack analytical results. Furthermore, the method yields estimates of the quantiles of a sample statistic of any sample size though Monte Carlo simulations for only two optimally selected sample sizes are needed. An analysis of the Monte Carlo design is performed to obtain the optimal choices of these two sample sizes and the number of simulated samples required for each sample size. Theoretical results are presented for the bias and variance of the numerical method proposed. The results developed are illustrated via simulation studies for the classical problem of estimating a bivariate linear structural relationship. It is seen that the size of the simulated samples used in the Monte Carlo method does not have to be very large and the method provides a better approximation to quantiles than those based on an asymptotic normal theory for skewed sampling distributions.     

Comparing two independent groups via the lower and upper quantiles

The most common strategy for comparing two independent groups is in terms of some measure of location intended to reflect the typical observation. However, it can be informative and important to compare the lower and upper quantiles as well, but when there are tied values, extant techniques suffer from practical concerns reviewed in the paper. For the special case where the goal is to compare the medians, a slight generalization of the percentile bootstrap method performs well in terms of controlling Type I errors when there are tied values [Wilcox RR. Comparing medians. Comput. Statist. Data Anal. 2006;51:1934–1943]. But our results indicate that when the goal is to compare the quartiles, or quantiles close to zero or one, this approach is highly unsatisfactory when the quantiles are estimated using a single order statistic or a weighted average of two order statistics. The main result in this paper is that when using the Harrell–Davis estimator, which uses all of the order statistics to estimate a quantile, control over the Type I error probability can be achieved in simulations, even when there are tied values, provided the sample sizes are not too small. It is demonstrated that this method can also have substantially higher power than the distribution free method derived by Doksum and Sievers [Plotting with confidence: graphical comparisons of two populations. Biometrika 1976;63:421–434]. Data from two studies are used to illustrate the practical advantages of the method studied here.     

Empirical Likelihood Inference for Population Quantiles with Unbalanced Ranked Set Samples

We consider nonparametric interval estimation for the population quantiles based on unbalanced ranked set samples. We derived the large sample distribution of the empirical log likelihood ratio statistic for the quantiles. Approximate intervals for quantiles are obtained by inverting the likelihood ratio statistic. The performance of the empirical likelihood interval is investigated and compared with the performance of the intervals based on the ranked set sample order statistics.     

Direct estimation of the percentile 'p-value' for the one-sample median test

In this paper we outline and illustrate an easy-to-use inference procedure for directly calculating the approximate bootstrap percentile-type p-value for the one-sample median test, i.e. we calculate the bootstrap p -value without resampling, by using a fractional order statistics based approach. The method parallels earlier work on fractionalorder-statistics-based non-parametric bootstrap percentile-type confidence intervals for quantiles. Monte Carlo simulation studies are performed, which illustrate that the fractional-order-statistics-based approach to the one-sample median test has accurate type I error control for small samples over a wide range of distributions; is easy to calculate; and is preferable to the sign test in terms of type I error control and power. Furthermore, the fractional-order-statistics-based median test is easily generalized to testing that any quantile has some hypothesized value; for example, tests for the upper or lower quartile may be performed using the same framework.     

Confidence intervals for a two-parameter exponential distribution: One- and two-sample problems

The problems of interval estimating the mean, quantiles, and survival probability in a two-parameter exponential distribution are addressed. Distribution function of a pivotal quantity whose percentiles can be used to construct confidence limits for the mean and quantiles is derived. A simple approximate method of finding confidence intervals for the difference between two means and for the difference between two location parameters is also proposed. Monte Carlo evaluation studies indicate that the approximate confidence intervals are accurate even for small samples. The methods are illustrated using two examples.     

COMPARISON OF FEATURE SIGNIFICANCE QUANTILE APPROXIMATIONS

Curve estimates and surface estimates often contain features such as inclines, bumps or ridges which may signify an underlying structural mechanism. However, spurious features are also a common occurrence and it is important to identify those features that are statistically significant. A method has been developed recently for recognising feature significance based on the derivatives of the function estimate. It requires simultaneous confidence intervals and tests, which in turn require quantiles for the maximal deviation statistics. This paper reviews and compares various approximations to these quantiles. Applying upcrossing‐probability theory to this problem yields better quantile approximations than the use of an independent blocks method.     

Simultaneous Confidence Intervals for Ratios of Fixed Effect Parameters in Linear Mixed Models

In multiple comparisons of fixed effect parameters in linear mixed models, treatment effects can be reported as relative changes or ratios. Simultaneous confidence intervals for such ratios had been previously proposed based on Bonferroni adjustments or multivariate normal quantiles accounting for the correlation among the multiple contrasts. We propose Fieller-type intervals using multivariate t quantiles and the application of Markov chain Monte Carlo techniques to sample from the joint posterior distribution and construct percentile-based simultaneous intervals. The methods are compared in a simulation study including bioassay problems with random intercepts and slopes, repeated measurements designs, and multicenter clinical trials.     

Small sample asymptotic inference for the coefficient of variation: normal and nonnormal models

In applied statistics, the coefficient of variation is widely calculated and interpreted even when the sample size of the data set is very small. However, confidence intervals for the coefficient of variation are rarely reported. One of the reasons is the exact confidence interval for the coefficient of variation, which is given in Lehmann (Testing Statistical Hypotheses, 2nd Edition, Wiley, New York, 1996), is very difficult to calculate. Various asymptotic methods have been proposed in literature. These methods, in general, require the sample size to be large. In this article, we will apply a recently developed small sample asymptotic method to obtain approximate confidence intervals for the coefficient of variation for both normal and nonnormal models. These small sample asymptotic methods are very accurate even for very small sample size. Numerical examples are given to illustrate the accuracy of the proposed method.     

Nonparametric confidence intervals for quantiles and quantile intervals from inversely sampled record breaking data in a multi-sample plan

There are a number of situations in which an observation is retained only if it is a record value, which include studies in industrial quality control experiments, destructive stress testing, meteorology, hydrology, seismology, athletic events and mining. When the number of records is fixed in advance, the data are referred to as inversely sampled record-breaking data. In this paper, we study the problems of constructing the nonparametric confidence intervals for quantiles and quantile intervals of the parent distribution based on record data. For a single record-breaking sample, the confidence coefficients of the confidence intervals for the pth quantile cannot exceed p and 1?p, on the basis of upper and lower records, respectively; hence, replication is required. So, we develop the procedure based on k independent record-breaking samples. Various cases have been studied and in each case, the optimal k and the exact nonparametric confidence intervals are obtained, and exact expressions for the confidence coefficients of these confidence intervals are derived. Finally, the results are illustrated by numerical computations.     

Maximum likelihood estimation in a Weibull regression model with type-1 censoring: a Monte Carlo study

Results of the Monte Carlo study of the performance of a maximum likelihood estimation in a Weibull parametric regression model with two explanatory variables are presented. One simulation run contained 1000 samples censored on the average by the amount of 0-30%. Each simulatedsample was generated in a form of two-factor two-level balanced experiment. The confidence intervals were computed using the large-sample normal approximation via the matrix of observed information. For small sample sizes the estimates of the scale parameter b of the loglifetime were significantly negatively biased, which resulted in a poor quality of confidence intervals for b and the low-level quantiles. All estimators improved their quality when the nominal value of b decreased. A moderate amount of censoring improved the quality of point and confidence estimation. The reparametrization b 7 produced rather accurate confidence intervals. Exact confidence intervals for b in case of non-censoring were obtained using the pivotal quantity b/b.     

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.     

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.     

Confidence interval estimation and hypothesis testing for a common correlation coefficient

Several procedures have been proposed for testing the hypothesis that all off-diagonal elements of the correlation matrix of a multivariate normal distribution are equal. If the hypothesis of equal correlation can be accepted, it is then of interest to estimate and perhaps test hypotheses for the common correlation. In this paper, two versions of five different test statistics are compared via simulation in terms of adequacy of the normal approximation, coverage probabilities of confidence intervals, control of Type I error, and power. The results indicate that two test statistics based on the average of the Fisher z-transforms of the sample correlations should be used in most cases. A statistic based on the sample eigenvalues also gives reasonable results for confidence intervals and lower-tailed tests.     

Comparison of Several Confidence Intervals for Median Residual Lifetime with Left-truncated and Right-censored Data

For left-truncated and right-censored data, the technique proposed by Brookmeyer and Crowley (1982) is extended to construct a point-wise confidence interval for median residual lifetime. This procedure is computationally simpler than the score type confidence interval in Jeong et al. (2008) and empirical likelihood ratio confidence interval in Zhou and Jeong (2011). Further, transformations of the estimator are applied to improve the approximation to the asymptotic distribution for small sample sizes. A simulation study is conducted to investigate the accuracy of these confidence intervals and the implementation of these confidence intervals to two real datasets is illustrated.