New bootstrap confidence intervals for means of positively skewed distributions

ABSTRACT

In this paper, we consider the problem of constructing non parametric confidence intervals for the mean of a positively skewed distribution. We suggest calibrated, smoothed bootstrap upper and lower percentile confidence intervals. For the theoretical properties, we show that the proposed one-sided confidence intervals have coverage probability α + O(n^{? 3/2}). This is an improvement upon the traditional bootstrap confidence intervals in terms of coverage probability. A version smoothed approach is also considered for constructing a two-sided confidence interval and its theoretical properties are also studied. A simulation study is performed to illustrate the performance of our confidence interval methods. We then apply the methods to a real data set.     

Confidence intervals for the mean of a normal distribution with restricted parameter space

For a normal distribution with known variance, the standard confidence interval of the location parameter is derived from the classical Neyman procedure. When the parameter space is known to be restricted, the standard confidence interval is arguably unsatisfactory. Recent articles have addressed this problem and proposed confidence intervals for the mean of a normal distribution where the parameter space is not less than zero. In this article, we propose a new confidence interval, rp interval, and derive the Bayesian credible interval and likelihood ratio interval for general restricted parameter space. We compare these intervals with the standard interval and the minimax interval. Simulation studies are undertaken to assess the performances of these confidence intervals.     

New Confidence Intervals for the Proportion of Interest in One-Sample Correlated Binary Data

Correlated binary data is obtained in many fields of biomedical research. When constructing a confidence interval for the proportion of interest, asymptotic confidence intervals have already been developed. However, such asymptotic confidence intervals are unreliable in small samples. To improve the performance of asymptotic confidence intervals in small samples, we obtain the Edgeworth expansion of the distribution of the studentized mean of beta-binomial random variables. Then, we propose new asymptotic confidence intervals by correcting the skewness in the Edgeworth expansion in one direct and two indirect ways. New confidence intervals are compared with the existing confidence intervals in simulation studies.     

Confidence Intervals for the Current Status Model

We discuss a new way of constructing pointwise confidence intervals for the distribution function in the current status model. The confidence intervals are based on the smoothed maximum likelihood estimator, using local smooth functional theory and normal limit distributions. Bootstrap methods for constructing these intervals are considered. Other methods to construct confidence intervals, using the non‐standard limit distribution of the (restricted) maximum likelihood estimator, are compared with our approach via simulations and real data applications.     

Confidence Intervals for Difference Between Two Poisson Rates

In this article, we develop four explicit asymptotic two-sided confidence intervals for the difference between two Poisson rates via a hybrid method. The basic idea of the proposed method is to estimate or recover the variances of the two Poisson rate estimates, which are required for constructing the confidence interval for the rate difference, from the confidence limits for the two individual Poisson rates. The basic building blocks of the approach are reliable confidence limits for the two individual Poisson rates. Four confidence interval estimators that have explicit solutions and good coverage levels are employed: the first normal with continuity correction, Rao score, Freeman and Tukey, and Jeffreys confidence intervals. Using simulation studies, we examine the performance of the four hybrid confidence intervals and compare them with three existing confidence intervals: the non-informative prior Bayes confidence interval, the t confidence interval based on Satterthwait's degrees of freedom, and the Bayes confidence interval based on Student's t confidence coefficient. Simulation results show that the proposed hybrid Freeman and Tukey, and the hybrid Jeffreys confidence intervals can be highly recommended because they outperform the others in terms of coverage probabilities and widths. The other methods tend to be too conservative and produce wider confidence intervals. The application of these confidence intervals are illustrated with three real data sets.     

Confidence Intervals for the Weighted Sum of Two Independent Binomial Proportions

Confidence intervals for the difference of two binomial proportions are well known, however, confidence intervals for the weighted sum of two binomial proportions are less studied. We develop and compare seven methods for constructing confidence intervals for the weighted sum of two independent binomial proportions. The interval estimates are constructed by inverting the Wald test, the score test and the Likelihood ratio test. The weights can be negative, so our results generalize those for the difference between two independent proportions. We provide a numerical study that shows that these confidence intervals based on large‐sample approximations perform very well, even when a relatively small amount of data is available. The intervals based on the inversion of the score test showed the best performance. Finally, we show that as for the difference of two binomial proportions, adding four pseudo‐outcomes to the Wald interval for the weighted sum of two binomial proportions improves its coverage significantly, and we provide a justification for this correction.     

Bootstrapping a Universal Pivot When Nuisance Parameters are Estimated

In complete samples from a continuous cumulative distribution with unknown parameters, it is known that various pivotal functions can be constructed by appealing to the probability integral transform. A pivotal function (or simply pivot) is a function of the data and parameters that has the property that its distribution is free of any unknown parameters. Pivotal functions play a key role in constructing confidence intervals and hypothesis tests. If there are nuisance parameters in addition to a parameter of interest, and consistent estimators of the nuisance parameters are available, then substituting them into the pivot can preserve the pivot property while altering the pivot distribution, or may instead create a function that is approximately a pivot in the sense that its asymptotic distribution is free of unknown parameters. In this latter case, bootstrapping has been shown to be an effective way of estimating its distribution accurately and constructing confidence intervals that have more accurate coverage probability in finite samples than those based on the asymptotic pivot distribution. In this article, one particular pivotal function based on the probability integral transform is considered when nuisance parameters are estimated, and the estimation of its distribution using parametric bootstrapping is examined. Applications to finding confidence intervals are emphasized. This material should be of interest to instructors of upper division and beginning graduate courses in mathematical statistics who wish to integrate bootstrapping into their lessons on interval estimation and the use of pivotal functions.

[Received November 2014. Revised August 2015.]     

A new confidence interval in errors-in-variables model with known error variance

This paper considers constructing a new confidence interval for the slope parameter in the structural errors-in-variables model with known error variance associated with the regressors. Existing confidence intervals are so severely affected by Gleser–Hwang effect that they are subject to have poor empirical coverage probabilities and unsatisfactory lengths. Moreover, these problems get worse with decreasing reliability ratio which also result in more frequent absence of some existing intervals. To ease these issues, this paper presents a fiducial generalized confidence interval which maintains the correct asymptotic coverage. Simulation results show that this fiducial interval is slightly conservative while often having average length comparable or shorter than the other methods. Finally, we illustrate these confidence intervals with two real data examples, and in the second example some existing intervals do not exist.     

Interval Estimation for the Correlation Coefficient

ABSTRACT

The correlation coefficient (CC) is a standard measure of a possible linear association between two continuous random variables. The CC plays a significant role in many scientific disciplines. For a bivariate normal distribution, there are many types of confidence intervals for the CC, such as z-transformation and maximum likelihood-based intervals. However, when the underlying bivariate distribution is unknown, the construction of confidence intervals for the CC is not well-developed. In this paper, we discuss various interval estimation methods for the CC. We propose a generalized confidence interval for the CC when the underlying bivariate distribution is a normal distribution, and two empirical likelihood-based intervals for the CC when the underlying bivariate distribution is unknown. We also conduct extensive simulation studies to compare the new intervals with existing intervals in terms of coverage probability and interval length. Finally, two real examples are used to demonstrate the application of the proposed methods.     

Confidence intervals for the difference between two median survival times for clustered survival data

Clustered survival data arise often in clinical trial design, where the correlated subunits from the same cluster are randomized to different treatment groups. Under such design, we consider the problem of constructing confidence interval for the difference of two median survival time given the covariates. We use Cox gamma frailty model to account for the within-cluster correlation. Based on the conditional confidence intervals, we can identify the possible range of covariates over which the two groups would provide different median survival times. The associated coverage probability and the expected length of the proposed interval are investigated via a simulation study. The implementation of the confidence intervals is illustrated using a real data set.     

Conservative confidence intervals for a single parameter

The method of constructing confidence intervals from hypothesis tests is studied in the case in which there is a single unknown parameter and is proved to provide confidence intervals with coverage probability that is at least the nominal level. The confidence intervals obtained by the method in several different contexts are seen to compare favorably with confidence intervals obtained by traditional methods. The traditional intervals are seen to have coverage probability less than the nominal level in several instances, This method can be applied to all confidence interval problems and reduces to the traditional method when an exact pivotal statistic is known.     

Interval estimation via tail functions

The authors describe a new method for constructing confidence intervals. Their idea consists in specifying the cutoff points in terms of a function of the target parameter rather than as constants. When it is suitably chosen, this so‐called tail function yields shorter confidence intervals in the presence of prior information. It can also be used to improve the coverage properties of approximate confidence intervals. The authors illustrate their technique by application to interval estimation of the mean of Bernoulli and normal populations. They further suggest guidelines for choosing the optimal tail function and discuss the relationship with Bayesian inference.     

Confidence intervals for endpoints of a c.d.f. via bootstrap

We propose methods of constructing confidence intervals for endpoints of a distribution. Under a mild condition on the tail of distribution, asymptotically correct confidence intervals are derived by bootstrapping Weissman's (Comm. Statist. Theory Method A 10 (1981) 549–557) statistics. It is also shown that a modification of this method works for type II censored data.     

Confidence Intervals for Quantiles of a Two-parameter Exponential Distribution under Progressive Type-II Censoring

Confidence intervals for the pth-quantile Q of a two-parameter exponential distribution provide useful information on the plausible range of Q, and only inefficient equal-tail confidence intervals have been discussed in the statistical literature so far. In this article, the construction of the shortest possible confidence interval within a family of two-sided confidence intervals is addressed. This shortest confidence interval is always shorter, and can be substantially shorter, than the corresponding equal-tail confidence interval. Furthermore, the computational intensity of both methodologies is similar, and therefore it is advantageous to use the shortest confidence interval. It is shown how the results provided in this paper can apply to data obtained from progressive Type II censoring, with standard Type II censoring as a special case. The applications of more complex confidence interval constructions through acceptance set inversions that can employ prior information are also discussed.     

Statistical inference for multiple change-point models

In this article, we propose a new technique for constructing confidence intervals for the mean of a noisy sequence with multiple change-points. We use the weighted bootstrap to generalize the bootstrap aggregating or bagging estimator. A standard deviation formula for the bagging estimator is introduced, based on which smoothed confidence intervals are constructed. To further improve the performance of the smoothed interval for weak signals, we suggest a strategy of adaptively choosing between the percentile intervals and the smoothed intervals. A new intensity plot is proposed to visualize the pattern of the change-points. We also propose a new change-point estimator based on the intensity plot, which has superior performance in comparison with the state-of-the-art segmentation methods. The finite sample performance of the confidence intervals and the change-point estimator are evaluated through Monte Carlo studies and illustrated with a real data example.     

Shortest confidence and prediction intervals for the log-normal

We provide the shortest prediction interval for X, and the shortest confidence interval for the median of X, when X has the log-normal distribution for both the case σ², the variance of log X, known and unknown. Tables are given to assist the practitioner in constructing these intervals. A real-life example is provided to illustrate the results.     

Interval estimation for a binomial proportion: a bootstrap approach

This paper discusses the classic but still current problem of interval estimation of a binomial proportion. Bootstrap methods are presented for constructing such confidence intervals in a routine, automatic way. Three confidence intervals for a binomial proportion are compared and studied by means of a simulation study, namely: the Wald confidence interval, the Agresti–Coull interval and the bootstrap-t interval. A new confidence interval, the Agresti–Coull interval with bootstrap critical values, is also introduced and its good behaviour related to the average coverage probability is established by means of simulations.     

Assessing and comparing the accuracy of various bootstrap methods

We evaluate the performance of various bootstrap methods for constructing confidence intervals for mean and median of several common distributions. Using Monte Carlo simulation, we assessed performance by looking at coverage percentages and average confidence interval lengths. Poor performance is characterized by coverage deviating from 0.95 and large confidence interval lengths. Undercoverage is of greater concern than overcoverage. We also assess the performance of bootstrap methods in estimating the parameters of the Cox Proportional Hazard model and Accelerated Failure Time model.     

Estimating kurtosis and confidence intervals for the variance under nonnormality

Exact confidence intervals for variances rely on normal distribution assumptions. Alternatively, large-sample confidence intervals for the variance can be attained if one estimates the kurtosis of the underlying distribution. The method used to estimate the kurtosis has a direct impact on the performance of the interval and thus the quality of statistical inferences. In this paper the author considers a number of kurtosis estimators combined with large-sample theory to construct approximate confidence intervals for the variance. In addition, a nonparametric bootstrap resampling procedure is used to build bootstrap confidence intervals for the variance. Simulated coverage probabilities using different confidence interval methods are computed for a variety of sample sizes and distributions. A modification to a conventional estimator of the kurtosis, in conjunction with adjustments to the mean and variance of the asymptotic distribution of a function of the sample variance, improves the resulting coverage values for leptokurtically distributed populations.     

Coverage accuracy for a mean without third moment

For constructing a confidence interval for the mean of a random variable with a known variance, one may prefer the sample mean standardized by the true standard deviation to the Student's t-statistic since the information of knowing the variance is used in the former way. In this paper, by comparing the leading error term in the expansion of the coverage probability, we show that the above statement is not true when the third moment is infinite. Our theory prefers the Student's t-statistic either when one-sided confidence intervals are considered for a heavier tail distribution or when two-sided confidence intervals are considered. Unlike other existing expansions for the Student's t-statistic, the derived explicit expansion for the case of infinite third moment can be used to estimate the coverage error so that bias correction becomes possible.