Edgeworth expansion and the bootstrap for stratified sampling without replacement from a finite population

An Edgeworth expansion for a linear combination of stratum means in stratified sampling without replacement from a finite population is derived. The expansion is applied to a bootstrap proposed for this context to show that the bootstrap captures the second-order term of the expansion.     

On cornish-fisher expansions as a ratio of determinants without cumulants

A new method of approximating one quantile of a distribution function in terms of the corresponding quantile of another distribution function is introduced. The method utilizes the Cornish-Fisher expansion so as to eliminate the requirement for knowing the cumulants while at the same time retaining the desired simplicity as well as the property of not affecting the order of the error of the approximation.     

Confidence intervals for the ratio of two variances

In this paper we consider confidence intervals for the ratio of two population variances. We propose a confidence interval for the ratio of two variances based on the t-statistic by deriving its Edgeworth expansion and considering Hall's and Johnson's transformations. Then, we consider the coverage accuracy of suggested intervals and intervals based on the F-statistic for some distributions.     

Generalized inferences on the common mean vector of several multivariate normal populations

The hypothesis testing and confidence region are considered for the common mean vector of several multivariate normal populations when the covariance matrices are unknown and possibly unequal. A generalized confidence region is derived using the concepts of generalized method based on the generalized p$p$-value. The generalized confidence region is illustrated with two numerical examples. The merits of the proposed method are numerically compared with those of existing methods with respect to their expected area or expected d-dimensional volumes and coverage probabilities under different scenarios.     

The neyman accuracy and the wolfowitz accuracy of the stein type confidence interval for the disturbance variance

In this paper we consider the Neyman accuracy and the Wolfowitz accuracy of the Stein type improved confidence interval I?_S for the disturbance variance in a linear regression model. The Neyman accuracy is a measure related to the unbiasedness of a confidence interval, and the Wolfowitz accuracy is related to the closeness of the endpoints to the true parameter. We show that I?_S is not unbiased and give some numerical results for the Neyman accuracy. As for the Wolfowitz accuracy we derive the sufficient condition for I?_S to improve on the usual confidence interval under this criterion and show numerically that a large degree of improvement can be obtainted.     

A note on the mean residual life function of a coherent system with exchangeable or nonidentical components

This article investigates some properties of the mean residual life function of (n−k+1)-out-of-n systems, when the lifetimes of the system components are independent random variables but not necessarily identically distributed and when the joint distribution of the component lifetimes is exchangeable, extending the results of Asadi and Goliforushani (2008) [On the mean residual life function of coherent systems. IEEE Transactions on Reliability 57 (4) 574-580] for the case of independent and identically distributed components. The extension to a coherent system with exchangeable components is also given.     

Cochran's rule for simple random sampling

Cochran's rule for the minimum sample size to ensure adequate coverage of nominal 95% confidence intervals is derived by using the Edgeworth expansion for the distribution function of the standardized sample mean. The rule is extended for confidence intervals based on the Studentized sample mean. The performance of the rule and Edgeworth approximations for smaller sample sizes are examined by simulation.     

A bias-reduced estimator for the mean of a heavy-tailed distribution with an infinite second moment

We use bias-reduced estimators of high quantiles of heavy-tailed distributions, to introduce a new estimator for the mean in the case of infinite second moment. The asymptotic normality of the proposed estimator is established and checked in a simulation study, by four of the most popular goodness-of-fit tests. The accuracy of the resulting confidence intervals is evaluated as well. We also investigate the finite sample behavior and compare our estimator with some versions of Peng's estimator of the mean (namely those based on Hill, t-Hill and Huisman et al. extreme value index estimators). Moreover, we discuss the robustness of the tail index estimators used in this paper. Finally, our estimation procedure is applied to the well-known Danish fire insurance claims data set, to provide confidence bounds for the means of weekly and monthly maximum losses over a period of 10 years.     

A note on bias reduction of maximum likelihood estimates for the scalar skew t distribution

The skew t distribution is a flexible parametric family to fit data, because it includes parameters that let us regulate skewness and kurtosis. A problem with this distribution is that, for moderate sample sizes, the maximum likelihood estimator of the shape parameter is infinite with positive probability. In order to try to solve this problem, Sartori (2006) has proposed using a modified score function as an estimating equation for the shape parameter. In this note we prove that the resulting modified maximum likelihood estimator is always finite, considering the degrees of freedom as known and greater than or equal to 2.     

Hypotheses testing for a multidimensional parameter of inhomogeneous Poisson processes

We observe a realization of an inhomogeneous Poisson process whose intensity function depends on an unknown multidimensional parameter. We consider the asymptotic behaviour of the Rao score test for a simple null hypothesis against the multilateral alternative. By using the Edgeworth type expansion (under the null hypothesis) for a vector of stochastic integrals with respect to the Poisson process, we refine the (classic) threshold of the test (obtained by the central limit theorem), which improves the first type probability of error. The expansion allows us to describe the power of the test under the local alternative, i.e. a sequence of alternatives, which converge to the null hypothesis with a certain rate. The rates can be different for components of the parameter.     

The Cornish-Fisher expansion for a class of statistics in first order autoregression

ABSTRACT

This article is concerned with the derivation and study of the Cornish-Fisher expansion for a wide class of estimators of the parameter in the first order autoregressive process. Second and third order Cornish-Fisher approximations to the quantile of the distribution of the corresponding asymptotically normal standardized statistic are stated explicitly and their accuracy is examined, both theoretically and numerically, by comparing them with the exact value of the quantile obtained by Monte Carlo simulation.     

Robustness of a bayes estimator for the mean of a normal population with nonnormal prior

The posterior distribution, Bayes estimator and risk are studied for the mean of a normal population (known variance) with respect to a moderately noimormal prior distribution. Extensive numerical work has suggested limiting behaviour of the derived results with respect to the sample size and precision of the prior distribution.     

On the relative performance of the block bootstrap for dependent data

In the independent setting, both Efron's bootstrap and “empiricai Edgeworth expansion” (E.E-expansion) give second-order accurate approximations to distributions of standardized and studentized statistics in the smooth function model. As a result, Efron's bootstrap was often regarded as roughly equivalent to the one-term E.E-expansion. However, a more detailed analysis shows that Efron's bootstrap outperforms the E.E-expansion in terms of loss functions by Bhattacharya and Qumsiyeh (1989) and in terms of probabilities for large deviations by Hall (1990) and Jing et a1 (1994). in this paper, we shall study the performances of the block bootstrap and the E.E-expansion for the weakly dependent data. It turns out that similar properties hold:both perform equally well at the center of the distribution but the block bootstrap provides accurate approximations even in the tails of the distributions. The study is focued on the simple case of standardized and studentized sample mean, but the conclusions can be easily extended to the smooth function of multivariate means.     

The Bayesian and frequentist approaches to testing a one-sided hypothesis about a multivariate mean

This paper compares the Bayesian and frequentist approaches to testing a one-sided hypothesis about a multivariate mean. First, this paper proposes a simple way to assign a Bayesian posterior probability to one-sided hypotheses about a multivariate mean. The approach is to use (almost) the exact posterior probability under the assumption that the data has multivariate normal distribution, under either a conjugate prior in large samples or under a vague Jeffreys prior. This is also approximately the Bayesian posterior probability of the hypothesis based on a suitably flat Dirichlet process prior over an unknown distribution generating the data. Then, the Bayesian approach and a frequentist approach to testing the one-sided hypothesis are compared, with results that show a major difference between Bayesian reasoning and frequentist reasoning. The Bayesian posterior probability can be substantially smaller than the frequentist p-value. A class of example is given where the Bayesian posterior probability is basically 0, while the frequentist p-value is basically 1. The Bayesian posterior probability in these examples seems to be more reasonable. Other drawbacks of the frequentist p-value as a measure of whether the one-sided hypothesis is true are also discussed.     

Testing linear hypotheses of mean vectors for high-dimension data with unequal covariance matrices

We propose a new test procedure for testing linear hypothesis on the mean vectors of normal populations with unequal covariance matrices when the dimensionality, p exceeds the sample size N  , i.e. p/N→c<∞$p/N\to c<\infty$. Our procedure is based on the Dempster trace criterion and is shown to be consistent in high dimensions.     

Nonparametric estimator for mean residual life and vitality function

In this paper, we study the estimation of the vitality function(v(x)=E(X|X>x) and mean residual life function(e(x)=E(X-x|X>x) from a sample ofX using the empirical estimator and kernel estimator. Under suitable conditions of regularity, the asymptotic normality of the kernel estimator is obtained. Partially supported by Consejeria de Cultura y Ed. (C.A.R.M.), under Grant PIB 95/90.     

Randomized confidence intervals of a parameter for a family of discrete exponential type distributions

For a family of one-parameter discrete exponential type distributions, the higher order approximation of randomized confidence intervals derived from the optimum test is discussed. Indeed, it is shown that they can be asymptotically constructed by means of the Edgeworth expansion. The usefulness is seen from the numerical results in the case of Poisson and binomial distributions.     

Bayesian Analysis for a Redundant Repairable System with Imperfect Coverage

System characteristics of a redundant repairable system with two primary units and one standby are studied from a Bayesian viewpoint with different types of priors assumed for unknown parameters, in which the coverage factor is the same for an operating unit failure as that for a standby unit failure. Times to failure and times to repair of the operating and standby units are assumed to follow exponential distributions. When times to failure and times to repair with uncertain parameters, a Bayesian approach is adopted to evaluate system characteristics. Monte Carlo simulation is used to derive the posterior distribution for the mean time to system failure and the steady-state availability. Some numerical experiments are performed to illustrate the results derived in this paper.     

Bootstrap confidence regions computed from autoregressions of arbitrary order

Given a linear time series, e.g. an autoregression of infinite order, we may construct a finite order approximation and use that as the basis for confidence regions. The sieve or autoregressive bootstrap, as this method is often called, is generally seen as a competitor with the better-understood block bootstrap approach. However, in the present paper we argue that, for linear time series, the sieve bootstrap has significantly better performance than blocking methods and offers a wider range of opportunities. In particular, since it does not corrupt second-order properties then it may be used in a double-bootstrap form, with the second bootstrap application being employed to calibrate a basic percentile method confidence interval. This approach confers second-order accuracy without the need to estimate variance. That offers substantial benefits, since variances of statistics based on time series can be difficult to estimate reliably, and—partly because of the relatively small amount of information contained in a dependent process—are notorious for causing problems when used to Studentize. Other advantages of the sieve bootstrap include considerably greater robustness against variations in the choice of the tuning parameter, here equal to the autoregressive order, and the fact that, in contradistinction to the case of the block bootstrap, the percentile t version of the sieve bootstrap may be based on the 'raw' estimator of standard error. In the process of establishing these properties we show that the sieve bootstrap is second order correct.     

Regularized parameter estimation of high dimensional t distribution

We propose penalized-likelihood methods for parameter estimation of high dimensional t distribution. First, we show that a general class of commonly used shrinkage covariance matrix estimators for multivariate normal can be obtained as penalized-likelihood estimator with a penalty that is proportional to the entropy loss between the estimate and an appropriately chosen shrinkage target. Motivated by this fact, we then consider applying this penalty to multivariate t distribution. The penalized estimate can be computed efficiently using EM algorithm for given tuning parameters. It can also be viewed as an empirical Bayes estimator. Taking advantage of its Bayesian interpretation, we propose a variant of the method of moments to effectively elicit the tuning parameters. Simulations and real data analysis demonstrate the competitive performance of the new methods.