The automatic percentile method: Accurate confidence limits in parametric models

The problem of constructing confidence limits for a scalar parameter is considered. Under weak conditions, Efron's accelerated bias-corrected bootstrap confidence limits are correct to second order in parametric familles. In this article, a new method, called the automatic percentile method, for setting approximate confidence limits is proposed as an attempt to alleviate two problems inherent in Efron's method. The accelerated bias-corrected method is not fully automatic, since it requires the calculation of an analytical adjustment; furthermore, it is typically not exact, though for many situations, particularly scalar-parameter familles, exact answers are available. In broader generality, the proposed method is exact when exact answers exist, and it is second-order accurate otherwise. The automatic percentile method is automatic, and for scalar parameter models it can be iterated to achieve higher accuracy, with the number of computations being linear in the number of iterations. However, when nuisance parameters are present, only second-order accuracy seems obtainable.     

Regeneration-based bootstrap for Markov chains

We consider a bootstrap method for Markov chains where the original chain is broken into a (random) number of cycles based on an atom (regeneration point) and the bootstrap scheme resamples from these cycles. We investigate the asymptotic accuracy of this method for the case of a sum (or a sample mean) related to the Markov chain. Under some standard moment conditions, the method is shown to be at least as good as the normal approximation, and better (second-order accurate) in the case of nonlattice summands. We give three examples to illustrate the applicability of our results.     

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.     

Simultaneous inferences for ordered exponential location parameters under unbalanced data and heteroscedasticity of scale parameters

A test procedure for testing homogeneity of location parameters against simple ordered alternative is proposed for k(k ≥ 2) members of two parameter exponential distribution under unbalanced data and heteroscedasticity of the scale parameters. The relevant one-sided and two-sided simultaneous confidence intervals (SCIs) for all k(k ? 1)/2 ordered pairwise differences of location parameters are also proposed. Simulation-based study revealed that the proposed procedure is better than the recently proposed procedure in terms of power, coverage probability, and average volume of SCIs. The implementation of proposed procedure is demonstrated through real life data.     

Parametric bootstrap inferences for unbalanced panel data models

This article presents parametric bootstrap (PB) approaches for hypothesis testing and interval estimation for the regression coefficients of panel data regression models with incomplete panels. Some simulation results are presented to compare the performance of the PB approaches with the approximate inferences. Our studies show that the PB approaches perform satisfactorily for various sample sizes and parameter configurations, and the performance of PB approaches is mostly better than the approximate methods with respect to the coverage probabilities and the Type I error rates. The PB inferences have almost exact coverage probabilities and Type I error rates. Furthermore, the PB procedure can be simply carried out by a few simulation steps, and the derivation is easier to understand and to be extended to the multi-way error component regression models with unbalanced panels. Finally, the proposed approaches are illustrated by using a real data example.     

Estimating precision,repeatability, and reproducibility from Gaussian and non- Gaussian data: a mixed models approach

Quality control relies heavily on the use of formal assessment metrics. In this paper, for the context of veterinary epidemiology, we review the main proposals, precision, repeatability, reproducibility, and intermediate precision, in agreement with ISO (international Organization for Standardization) practice, generalize these by placing them within the linear mixed model framework, which we then extend to the generalized linear mixed model setting, so that both Gaussian as well as non-Gaussian data can be employed. Similarities and differences are discussed between the classical ANOVA (analysis of variance) approach and the proposed mixed model settings, on the one hand, and between the Gaussian and non-Gaussian cases, on the other hand. The new proposals are applied to five studies in three diseases: Aujeszky's disease, enzootic bovine leucosis (EBL) and bovine brucellosis. The mixed-models proposals are also discussed in the light of their computational requirements.     

Second-order properties of improved two-stage procedure for a multivariate normal distribution

We consider the classical two problems: (I) fixed-size confidence region estimation problem, and (II) bounded risk point estimation problem, for the mean vector μ in the N_p(μ, ∑) population when ∑ is unknown and positive definite. Healy's (1956) solution given by a two-stage procedure for the prob¬lem (I) is improved in terms of the sample size. The second-order properties of the solution to each problem are verified under the assumptions that the maximum latent root which is simple of ∑ is bounded below for (I) and tr(∑) is bounded below for (II), by respective known and positive number.     

Adjusted empirical likelihood inference for additive hazards regression

ABSTRACT

This article develops an adjusted empirical likelihood (EL) method for the additive hazards model. The adjusted EL ratio is shown to have a central chi-squared limiting distribution under the null hypothesis. We also evaluate its asymptotic distribution as a non central chi-squared distribution under the local alternatives of order n^{? 1/2}, deriving the expression for the asymptotic power function. Simulation studies and a real example are conducted to evaluate the finite sample performance of the proposed method. Compared with the normal approximation-based method, the proposed method tends to have more larger empirical power and smaller confidence regions with comparable coverage probabilities.     

The relationship between the improvement on the point estimation and the improvement on the interval estimation for the disturbance variance in a linear regression model

In this paper the relationship between the improvement on the point estimation and the improvement on the interval estimation for the disturbance variance in a linear regression model is discussed It is shown that substituting the Stein-type estimatoi for the usual estimatoi in the confidence interval leads to the improvement on the interval estimation The equal-tailed and the shoitest unbiased intervals are dealt with Some appealing relationship is also found in the unbiased case.     

Tests and confidence intervals for the shape parameter of a gamma distribution

Inference based on the Central Limit Theorem has only first order accuracy. We give tests and confidence intervals (CIs) of second orderaccuracy for the shape parameter ρ of a gamma distribution for both the unscaled and scaled cases.

Tests and CIs based on moment and cumulant estimates are considered as well as those based on the maximum likelihood estimate (MLE).

For the unscaled case the MLE is the moment estimate of order zero; the most efficient moment estimate of integral order is the sample mean, having asymptotic relative efficiency (ARE) .61 when ρ= 1.

For the scaled case the most efficient moment estimate is a functionof the mean and variance. Its ARE is .39 when ρ = 1.

Our motivation for constructing these tests of ρ = 1 and CIs forρ is to provide a simple and convenient method for testing whether a distribution is exponential in situations such as rainfall models where such an assumption is commonly made.