Bootstrap confidence regions for correspondence analysis

A bootstrap-based method for constructing confidence regions (CRs) around row or column points projected onto a pair of axes from the correspondence analysis (CA) of a two-way contingency table is presented. These regions deal with the specific question of the sampling variation of sample row and column profile points around population row and column profile points when both are projected onto the observed axes, rather than the decomposition of the χ²-test of independence or the general question of the stability of the observed CA display which has been considered in previous work. The method therefore constructs the regions in a different way to what has been proposed before. A simulation experiment shows that the method performs well in most of the situations in which it might be used, with a few exceptions being noted. An example illustrates that the method produces conclusions which are consistent with those from detailed parametric modelling.     

On confidence regions for the mean of a multivariate time series

We consider the problem of setting up a confidence region for the mean of amultivariate timeseries ont he basis of a part-realisation of that series.A procedure for setting up a confidence interval for the mean of a univariate time series Is implicitin Jones(1976).We present an analogous procedure for setting up a confidence region for the mean of a multivariatet ime series.This procedure is base donastatistic which is an analogue of Hotelling'sT'.Our results are applied to a comparison of climate means obtained from experiments with a General Circulation Model of the earth's atmosphere.     

Two‐Sample Bootstrap Hypothesis Tests for Three‐Dimensional Labelled Landmark Data

Abstract. We investigate resampling methodologies for testing the null hypothesis that two samples of labelled landmark data in three dimensions come from populations with a common mean reflection shape or mean reflection size‐and‐shape. The investigation includes comparisons between (i) two different test statistics that are functions of the projection onto tangent space of the data, namely the James statistic and an empirical likelihood statistic; (ii) bootstrap and permutation procedures; and (iii) three methods for resampling under the null hypothesis, namely translating in tangent space, resampling using weights determined by empirical likelihood and using a novel method to transform the original sample entirely within refection shape space. We present results of extensive numerical simulations, on which basis we recommend a bootstrap test procedure that we expect will work well in practise. We demonstrate the procedure using a data set of human faces, to test whether humans in different age groups have a common mean face shape.     

Bootstrap confidence intervals for the pareto index

In the present paper we develop second-order theory using the subsample bootstrap in the context of Pareto index estimation. We show that the bootstrap is not second-order accurate, in the sense that it fails to correct the first term describing departure from the limit distribution. Worse than this, even when the subsample size is chosen optimally, the error between the subsample bootstrap approximation and the true distribution is often an order of magnitude larger than that oi tue asymptotic approximation. To overcome this deficiency, we show that an extrapolation method, based quite literally on a mixture of asymptotic and subsample bootstrap methods, can lead to second-order correct confidence intervals for the Pareto index.     

Exact confidence intervals for the shape parameter of the gamma distribution

In this paper exact confidence intervals (CIs) for the shape parameter of the gamma distribution are constructed using the method of Bølviken and Skovlund [Confidence intervals from Monte Carlo tests. J Amer Statist Assoc. 1996;91:1071–1078]. The CIs which are based on the maximum likelihood estimator or the moment estimator are compared to bootstrap CIs via a simulation study.     

A new test for the mean vector in large dimension and small samples

In this article, we consider the problem of testing the mean vector in the multivariate normal distribution, where the dimension p is greater than the sample size N. We propose a new test T_Block and obtain its asymptotic distribution. We also compare the proposed test with other two tests. The simulation results suggest that the performance of the new test is comparable to the existing two tests, and under some circumstances it may have higher power. Therefore, the new statistic can be employed in practice as an alternative choice.     

Generalized simultaneous confidence regions for regression coefficients and their ratios

This study constructs a simultaneous confidence region for two combinations of coefficients of linear models and their ratios based on the concept of generalized pivotal quantities. Many biological studies, such as those on genetics, assessment of drug effectiveness, and health economics, are interested in a comparison of several dose groups with a placebo group and the group ratios. The Bonferroni correction and the plug-in method based on the multivariate-t distribution have been proposed for the simultaneous region estimation. However, the two methods are asymptotic procedures, and their performance in finite sample sizes has not been thoroughly investigated. Based on the concept of generalized pivotal quantity, we propose a Bonferroni correction procedure and a generalized variable (GV) procedure to construct the simultaneous confidence regions. To address a genetic concern of the dominance ratio, we conduct a simulation study to empirically investigate the probability coverage and expected length of the methods for various combinations of sample sizes and values of the dominance ratio. The simulation results demonstrate that the simultaneous confidence region based on the GV procedure provides sufficient coverage probability and reasonable expected length. Thus, it can be recommended in practice. Numerical examples using published data sets illustrate the proposed methods.     

Bootstrap confidence bands for the CDF using ranked-set sampling

In ranked-set sampling (RSS), a stratification by ranks is used to obtain a sample that tends to be more informative than a simple random sample of the same size. Previous work has shown that if the rankings are perfect, then one can use RSS to obtain Kolmogorov–Smirnov type confidence bands for the CDF that are narrower than those obtained under simple random sampling. Here we develop Kolmogorov–Smirnov type confidence bands that work well whether the rankings are perfect or not. These confidence bands are obtained by using a smoothed bootstrap procedure that takes advantage of special features of RSS. We show through a simulation study that the coverage probabilities are close to nominal even for samples with just two or three observations. A new algorithm allows us to avoid the bootstrap simulation step when sample sizes are relatively small.     

Bootstrap confidence intervals for the coefficient of quartile variation

ABSTRACT

In non-normal populations, it is more convenient to use the coefficient of quartile variation rather than the coefficient of variation. This study compares the percentile and t-bootstrap confidence intervals with Bonett's confidence interval for the quartile variation. We show that empirical coverage of the bootstrap confidence intervals is closer to the nominal coverage (0.95) for small sample sizes (n = 5, 6, 7, 8, 9, 10 and 15) for most distributions studied. Bootstrap confidence intervals also have smaller average width. Thus, we propose using bootstrap confidence intervals for the coefficient of quartile variation when the sample size is small.     

Robust statistics for testing mean vectors of multivariate distributions

We develop a ‘robust’ statistic T² _R, based on Tiku's (1967, 1980) MML (modified maximum likelihood) estimators of location and scale parameters, for testing an assumed meam vector of a symmetric multivariate distribution. We show that T² _R is one the whole considerably more powerful than the prominenet Hotelling T² statistics. We also develop a robust statistic T² _D for testing that two multivariate distributions (skew or symmetric) are identical; T² _D seems to be usually more powerful than nonparametric statistics. The only assumption we make is that the marginal distributions are of the type (1/σ_k)f((x-μ_k)/σ_k) and the means and variances of these marginal distributions exist.     

A robust procedure for testing the equality of mean vectors of two bivariate populations with unequal covariance matrices

A robust test is developed for testing equality of the mean vectors of two bivariate (multivariate) populations when the variance-covariance matrices are not necessarily equal. The test is an extension of the univariate robust test given by Tiku and Singh (1981).     

Comparisons of approximate confidence intervals for the smallest extreme value distribution simple linear regression model under time censoring

Exact confidence interval estimation for accelerated life regression models with censored smallest extreme value (or Weibull) data is often impractical. This paper evaluates the accuracy of approximate confidence intervals based on the asymptotic normality of the maximum likelihood estimator, the asymptotic X²distribution of the likelihood ratio statistic, mean and variance correction to the likelihood ratio statistic, and the so-called Bartlett correction to the likelihood ratio statistic. The Monte Carlo evaluations under various degrees of time censoring show that uncorrected likelihood ratio intervals are very accurate in situations with heavy censoring. The benefits of mean and variance correction to the likelihood ratio statistic are only realized with light or no censoring. Bartlett correction tends to result in conservative intervals. Intervals based on the asymptotic normality of maximum likelihood estimators are anticonservative and should be used with much caution.     

A simulation study on the confidence interval procedures of some mean cumulative function estimators

Recurrence data are collected to study the recurrent events in biological, physical, and other systems. Quantities of interest include the mean cumulative number of events and the mean cumulative cost of the events. The mean cumulative function (MCF) can be estimated using non-parametric (NP) methods or by fitting parametric models, and many procedures have been suggested to construct the confidence intervals (CIs) for the MCF. This paper summarizes the results of a large simulation study that was designed to compare five CI procedures for both NP and parametric estimation. When performing parametric estimation, we assume the power law non-homogeneous Poisson process (NHPP) model. Our results include the evaluation of these procedures when they are used for window-observation recurrence data where recurrence histories of some systems are available only in observation windows with gaps in between.     

Tests of fit for exponentiality based on a characterization via the mean residual life function

We study two new omnibus goodness of fit tests for exponentiality, each based on a characterization of the exponential distribution via the mean residual life function. The limiting null distributions of the tests statistics are the same as the limiting null distributions of the Kolmogorov-Smirnov and Cramér-von Mises statistics proposed when testing the simple hypothesis that the distribution of the sample variables is uniform on the interval [0, 1]. Work supported by the Deutsche Forschungsgemeinschaft     

Globally applicable control chart for online monitoring of stability of process mean

The shape features of run chart patterns of the most recent m observations arising from stable and unstable processes are different. Using this fact, a new monitoring statistic is defined whose value for given m depends on the pattern parameters but not on the process parameters. A control chart for this statistic for given m, therefore, will be globally applicable to normal processes. The simulation study reveals that the proposed statistic approximately follows normal distribution. The performances of the globally applicable control chart in terms of average run lengths (ARLs) are evaluated and compared with the X chart. Both in-control ARL and out-of-control ARLs with respect to different abnormal process conditions are found to be larger than the X chart. However, the proposed concept is promising because it can eliminate the burden of designing separate control charts for different quality characteristics or processes in a manufacturing set-up.     

Bootstrap approach to test the homogeneity of order restricted mean vectors when the covariance matrices are unknown

Testing homogeneity of multivariate normal mean vectors under an order restriction when the covariance matrices are unknown, arbitrary positive definite and unequal are considered. This problem of testing has been studied to some extent, for example, by Kulatunga and Sasabuchi (1984 Kulatunga, D. D. S., Sasabuchi, S. (1984). A test of homogeneity of mean vectors against multivariate isotonic alternatives. Mem Fac Sci, Kyushu Univ Ser A Mathemat 38:151–161. [Google Scholar]) when the covariance matrices are known and also Sasabuchi et al. (2003 Sasabuchi, S., Tanaka, K., Tsukamodo, T. (2003). Testing homogeneity of multivariate normal mean vectors under an order restriction when the covariance matrices are common but unknown. Annals of Statistics. 31(5):1517–1536.[Web of Science ®] , [Google Scholar]) and Sasabuchi (2007 Sasabuchi, S. (2007). More powerful tests for homogeneity of multivariate normal mean vectors under an order restriction. Sankhya 69(4):700–716. [Google Scholar]) when the covariance matrices are unknown but common. In this paper, a test statistic is proposed and because of the main advantage of the bootstrap test is that it avoids the derivation of the complex null distribution analytically, a bootstrap test statistic is derived and since the proposed test statistic is location invariance the bootstrap p-value defined logical and some steps are presented to estimate it. Our numerical studies via Monte Carlo simulation show that the proposed bootstrap test can correctly control the type I error rates. The power of the test for some of the p-dimensional normal distributions is computed by Monte Carlo simulation. Also, the null distribution of test statistic is estimated using kernel density. Finally, the bootstrap test is illustrated using a real data.     

One and two sample confidence intervals for estimating the mean of skewed populations: an empirical comparative study

In this paper we consider and propose some confidence intervals for estimating the mean or difference of means of skewed populations. We extend the median t interval to the two sample problem. Further, we suggest using the bootstrap to find the critical points for use in the calculation of median t intervals. A simulation study has been made to compare the performance of the intervals and a real life example has been considered to illustrate the application of the methods.     

Improved confidence intervals for the exponential mean via tail functions

Abstract

The method of tail functions is applied to confidence estimation of the exponential mean in the presence of prior information. It is shown how the “ordinary” confidence interval can be generalized using a class of tail functions and then engineered for optimality, in the sense of minimizing prior expected length over that class, whilst preserving frequentist coverage. It is also shown how to derive the globally optimal interval, and how to improve on this using tail functions when criteria other than length are taken into consideration. Probabilities of false coverage are reported for some of the intervals under study, and the theory is illustrated by application to confidence estimation of a reliability coefficient based on some survival data.