Inference in the presence of ranking error in ranked set sampling

The author proposes inference techniques for ranked set sample data in the presence of judgment ranking errors. He bases his analysis on the models of Bohn & Wolfe (1994) and Frey (2007a, b), of which parameters are estimated by minimizing a distance measure. He then uses the fitted models to calibrate confidence intervals and tests. He shows the validity of his approach through simulation and illustrates its application through the construction of distribution‐free confidence intervals for the median area of apple tree leaves covered by a spray.     

On an asymptotic relative efficiency concept based on expected volumes of confidence regions

ABSTRACT

The paper deals with an asymptotic relative efficiency concept for confidence regions of multidimensional parameters that is based on the expected volumes of the confidence regions. Under standard conditions the asymptotic relative efficiencies of confidence regions are seen to be certain powers of the ratio of the limits of the expected volumes. These limits are explicitly derived for confidence regions associated with certain plugin estimators, likelihood ratio tests and Wald tests. Under regularity conditions, the asymptotic relative efficiency of each of these procedures with respect to each one of its competitors is equal to 1. The results are applied to multivariate normal distributions and multinomial distributions in a fairly general setting.     

Stein's two sample test of a general linear hypothesis and a noncentral f–type distribution

Stein's two–sample procedure for a general linear model is studied and derived in terms of matrices in which the error tems are distributed as multivatriate student t–error terms. Tests and confidence regions are constructed in a similar way to classical linear models which involves percentage points of student t and F distributions. The advantages of taking two samples are: the variance of the error terms is known, and the power of tests are size of confidence regions are controllable. A new distribution called noncentral F–type distribution different from the nencentral F is found when considerinf the power of the test of general linear hypothesis.     

A stein two-sample procedure for the general linear model with unequal error variances

The pronerties of the tests and confidence regions for the parameters in the classical general linear model depend upon the equality of the variances of the error terms. The level and power of tests and the confidence coefficients associated with confidence regions are vitiated when the assumption of equality is not true. Even when the error variances are equal the power of tests and the size of confidence regions depend upon the unknown common variance and hence are uncontrollable. This paper presents a two-stage procedure which yields tests and confidence regions which are completely independent of the variances of the errors and hence tests with controllable power and confidence regions of fixed controllable size are obtained.     

Hypothesis Testing and Confidence Regions for the Mean Sojourn Time of an M/M/1 Queueing System

This article discusses testing hypotheses and confidence regions with correct levels for the mean sojourn time of an M/M/1 queueing system. The uniformly most powerful unbiased tests for three usual hypothesis testing problems are obtained and the corresponding p values are provided. Based on the duality between hypothesis tests and confidence sets, the uniformly most accurate confidence bounds are derived. A confidence interval with correct level is proposed.     

Total least squares solution for compositional data using linear models

The restrictive properties of compositional data, that is multivariate data with positive parts that carry only relative information in their components, call for special care to be taken while performing standard statistical methods, for example, regression analysis. Among the special methods suitable for handling this problem is the total least squares procedure (TLS, orthogonal regression, regression with errors in variables, calibration problem), performed after an appropriate log-ratio transformation. The difficulty or even impossibility of deeper statistical analysis (confidence regions, hypotheses testing) using the standard TLS techniques can be overcome by calibration solution based on linear regression. This approach can be combined with standard statistical inference, for example, confidence and prediction regions and bounds, hypotheses testing, etc., suitable for interpretation of results. Here, we deal with the simplest TLS problem where we assume a linear relationship between two errorless measurements of the same object (substance, quantity). We propose an iterative algorithm for estimating the calibration line and also give confidence ellipses for the location of unknown errorless results of measurement. Moreover, illustrative examples from the fields of geology, geochemistry and medicine are included. It is shown that the iterative algorithm converges to the same values as those obtained using the standard TLS techniques. Fitted lines and confidence regions are presented for both original and transformed compositional data. The paper contains basic principles of linear models and addresses many related problems.     

Some properties of tests for parameters that can be arbitrarily close to being unidentified

Confidence intervals for parameters that can be arbitrarily close to being unidentified are unbounded with positive probability [e.g. Dufour, J.-M., 1997. Some impossibility theorems in econometrics with applications to instrumental variables and dynamic models. Econometrica 65, 1365–1388; Pfanzagl, J. 1998. The nonexistence of confidence sets for discontinuous functionals. Journal of Statistical Planning and Inference 75, 9–20], and the asymptotic risks of their estimators are unbounded [Pötscher, B.M., 2002. Lower risk bounds and properties of confidence sets for ill-posed estimation problems with applications to spectral density and persistence estimation, unit roots, and estimation of long memory parameters. Econometrica 70, 1035–1065]. We extend these “impossibility results” and show that all tests of size α concerning parameters that can be arbitrarily close to being unidentified have power that can be as small as α for any sample size even if the null and the alternative hypotheses are not adjacent. The results are proved for a very general framework that contains commonly used models.     

Linear bootstrap methods for vector autoregressive moving-average models

We provide the theoretical justification of bootstrapping stationary invertible echelon vector autoregressive moving-average (VARMA) models using linear methods. The asymptotic validity of the bootstrap is established with strong white noise under parametric and nonparametric assumptions. Our methods are practical and useful for building reliable simulation-based inference and forecasting without implementing nonlinear estimation techniques such as ML which is usually burdensome, time demanding or impractical, particularly in big or highly persistent systems. The relevance of our procedures is more pronounced in the context of dynamic simulation-based techniques such as maximized Monte Carlo (MMC) tests [see Dufour J-M. Monte Carlo tests with nuisance parameters: a general approach to finite-sample inference and nonstandard asymptotics in econometrics. J Econom. 2006;133(2):443–477 and Dufour J-M, Jouini T. Finite-sample simulation-based tests in VAR models with applications to Granger causality testing. J Econom. 2006;135(1–2):229–254 for the VAR case]. Simulation evidence shows that, compared with conventional asymptotics, our bootstrap methods have good finite-sample properties in approximating the actual distribution of the studentized echelon VARMA parameter estimates, and in providing echelon parameter confidence sets with satisfactory coverage.     

Empirical Likelihood for Single-Index Regression Models under Negatively Associated Errors

In this article, we use bockwise empirical likelihood technique to construct confidence regions for the parameter of the single-index models under negatively associated errors. It is shown that the blockwise empirical likelihood ratio statistic for the parameter of interest is asymptotically χ²-type distributed. The result can be used to obtain confidence regions for the parameter of interest.     

On confidence regions of semiparametric nonlinear reproductive dispersion models

This article investigates the confidence regions for semiparametric nonlinear reproductive dispersion models (SNRDMs), which is an extension of nonlinear regression models. Based on local linear estimate of nonparametric component and generalized profile likelihood estimate of parameter in SNRDMs, a modified geometric framework of Bates and Wattes is proposed. Within this geometric framework, we present three kinds of improved approximate confidence regions for the parameters and parameter subsets in terms of curvatures. The work extends the previous results of Hamilton et al. [in Accounting for intrinsic nonlinearity in nonlinear regression parameter inference regions, Ann. Statist. 10, pp. 386–393, 1982], Hamilton [in Confidence regions for parameter subset in nonlinear regression, Biometrika, 73, pp. 57–64, 1986], Wei [in On confidence regions of embedded models in regular parameter families (a geometric approch), Austral. J. Statist. 36, pp. 327–338, 1994], Tang et al. [in Confidence regions in quasi-likelihood nonlinear models: a geometric approach, J. Biomath. 15, pp. 55–64, 2000b] and Zhu et al. [in On confidence regions of semiparametric nonlinear regression models, Acta. Math. Scient. 20, pp. 68–75, 2000].     

Empirical Likelihood for Partial Linear Models under Negatively Associated Errors

In this paper, we use blockwise empirical likelihood (EL) technique to construct confidence regions for the parameter of the partial linear models under negatively associated errors. It is shown that the blockwise EL ratio statistic for the parameter of interest is asymptotically χ²-type distributed by employing the large-block and small-block arguments. The result can be used to obtain confidence regions for the parameter of interest.     

A comparison of multivariate signal discrimination techniques

This paper investigates several techniques to discriminate two multivariate stationary signals. The methods considered include Gaussian likelihood ratio tests for variance equality, a chi-squared time-domain test, and a spectral-based test. The latter two tests assess equality of the multivariate autocovariance function of the two signals over many different lags. The Gaussian likelihood ratio test is perhaps best viewed as principal component analyses (PCA) without dimension reduction aspects; it can be modified to consider covariance features other than variances via dimension augmentation tactics. A simulation study is constructed that shows how one can make inappropriate conclusions with PCA tests, even when dimension augmentation techniques are used to incorporate non-zero lag autocovariances into the analysis. The various discrimination methods are first discussed. A simulation study then illuminates the various properties of the methods. In this pursuit, calculations are needed to identify several multivariate time series models with specific autocovariance properties. To demonstrate the applicability of the methods, nine US and Canadian weather stations from three distinct regions are clustered. Here, the spectral clustering perfectly identified distinct regions, the chi-squared test performed marginally, and the PCA/likelihood ratio method did not perform well.     

Simultaneous inferences for variance ratios of some mixed linear models

Earlier investigations used a one-sided inequality to consltuct confidence regions for the variance ratios or balanced randoiu models. In this study, confidence regions are based on a two-sided generalisation of this inequality and the results are illustrated by estimating the parameters of some elementary random models.     

Empirical likelihood based confidence regions for first order parameters of heavy-tailed distributions

Let X₁,…,X_n be some i.i.d. observations from a heavy-tailed distribution F, i.e. the common distribution of the excesses over a high threshold u_n can be approximated by a generalized Pareto distribution G_γ,σn with γ>0. This paper deals with the problem of finding confidence regions for the couple (γ,σ_n): combining the empirical likelihood methodology with estimation equations (close but not identical to the likelihood equations) introduced by Zhang (2007), asymptotically valid confidence regions for (γ,σ_n) are obtained and proved to perform better than Wald-type confidence regions (especially those derived from the asymptotic normality of the maximum likelihood estimators). By profiling out the scale parameter, confidence intervals for the tail index are also derived.     

Exact tests and confidence regions in nonlinear regression

Nonlinear regression models with spherically symmetric error vectors and a single nonlinear parameter are considered. On the basis of a new geometric approach, exact one- and two-sided tests and confidence regions for the nonlinear parameter are derived in the cases of known and unknown error variances. A geometric measure representation formula is used to determine the power functions of the tests if the error variance is known and to derive different lower bounds for the power function of a one-sided test in the case of an unknown error variance. The latter can be done quite effectively by constructing and measuring several balls inside the critical region. A numerical study compares the results for different density generating functions of the error distribution.     

Conical tests: Observed levels of significance and confidence regions

Tests and confidence regions for a p-parameter nonnormal model can require integration in p dimensions, for example with a Bayesian or structural model. For small p, the computer integrations are manageable, see for example Fraser [3], Naylor and Smith [6], but for p beyond 5, 6, or 7, the integrations become unfeasible. This paper proposes conical tests of significance which involve manageable computer calculations. The conical tests also provide confidence regions, giving the confidence bound at a chosen level of significance in any direction from a central 0-confidence point.     

Robust confidence intervals for contrasts based upon a likelihood ratio test

This paper describes how to compute robust confidence intervals for differences of the effects using the likelihood ratio testF _M in the two-way analysis of variance. The probability for the α-error and the average length of the confidence intervals withF _m and the quadratic formQ _M are investigated and compared with the classical confidence intervals fort-distributed and lognormal errors. We also give a warning of building confidence intervals withF _M andQ _M in the presence of heterogeneous scale parameters, because these tests which do not regard heteroscedasticity are then much too liberal.     

Inference from PseudoLikelihoods with Plug‐In Estimates

Effective implementation of likelihood inference in models for high‐dimensional data often requires a simplified treatment of nuisance parameters, with these having to be replaced by handy estimates. In addition, the likelihood function may have been simplified by means of a partial specification of the model, as is the case when composite likelihood is used. In such circumstances tests and confidence regions for the parameter of interest may be constructed using Wald type and score type statistics, defined so as to account for nuisance parameter estimation or partial specification of the likelihood. In this paper a general analytical expression for the required asymptotic covariance matrices is derived, and suggestions for obtaining Monte Carlo approximations are presented. The same matrices are involved in a rescaling adjustment of the log likelihood ratio type statistic that we propose. This adjustment restores the usual chi‐squared asymptotic distribution, which is generally invalid after the simplifications considered. The practical implication is that, for a wide variety of likelihoods and nuisance parameter estimates, confidence regions for the parameters of interest are readily computable from the rescaled log likelihood ratio type statistic as well as from the Wald type and score type statistics. Two examples, a measurement error model with full likelihood and a spatial correlation model with pairwise likelihood, illustrate and compare the procedures. Wald type and score type statistics may give rise to confidence regions with unsatisfactory shape in small and moderate samples. In addition to having satisfactory shape, regions based on the rescaled log likelihood ratio type statistic show empirical coverage in reasonable agreement with nominal confidence levels.     

Second-order inference for generalized least squares

Confidence regions for generalized least squares are commonly derived from a measure of departure calculated on the tangent plane at the MLE or on the tangent plane at the true value; the first gives approximate confidence regions, the second exact. For surfaces with curvature, indeed with varying curvature, the exact regions typically are not likelihood regions and can include parameter values of highest and of lowest likelihood. This paper develops an alternative approach to deriving exact confidence regions and uses both surface curvature and distance from the surface as supporting ingredients. For this, conditionality is invoked in two ways beyond that supported by the usual conditionality principle. For the case of normal error the ordinary chi-squared departure is replaced by a Von Mises-type angular (or cosine) departure which is assessed using curvature properties in the data direction and radial distance of the data from the regression surface. For the usual linear model (constant curvature equal to zero) the method coincides with the ordinary tests and confidence regions; for the case of constant nonzero curvature, the method generalizes to spheres and sphere-cylinders the Fisher (Statistical Methods and Scientific Inference, 1956) analysis of a rotationally symmetric normal on ?² with mean constrained to a circle. The effects of conditioning are indicated by a computer plot for obtaining 95% confidence.     

Some results for type I censored sampling from geometric distributions

This paper gives the likelihood function for a Type I censored sample from the geometric distribution with parameter p, and the maximum likelihood estimator for p. Exact and asymptotic sampling distributions of joint sufficient statistics for p are derived. Such distributional results make it possible to develop tests or confidence intervals based on discrete censored data, which are not available now in the literature. Neyman-Pearson tests for p are developed. Examples are given to illustrate these results.