Testing hypotheses about the common mean of two normal populations

Two simple tests which allow for unequal sample sizes are considered for testing hypothesis for the common mean of two normal populations. The first test is an exact test of size a based on two available t-statistics based on single samples made exact through random allocation of α among the two available t-tests. The test statistic of the second test is a weighted average of two available t-statistics with random weights. It is shown that the first test is more efficient than the available two t-tests with respect to Bahadur asymptotic relative efficiency. It is also shown that the null distribution of the test statistic in the second test, which is similar to the one based on the normalized Graybill-Deal test statistic, converges to a standard normal distribution. Finally, we compare the small sample properties of these tests, those given in Zhou and Mat hew (1993), and some tests given in Cohen and Sackrowitz (1984) in a simulation study. In this study, we find that the second test performs better than the tests given in Zhou and Mathew (1993) and is comparable to the ones given in Cohen and Sackrowitz (1984) with respect to power..     

Testing the equality of correlation matrices

We present results that extend an existing test of equality of correlation matrices. A new test statistic is proposed and is shown to be asymptotically distributed as a linear combination of independent x ² random variables. This new formulation allows us to find the power of the existing test and our extensions by deriving the distribution under the alternative using a linear combination of independent non-central x ² random variables. We also investigate the null and the alternative distribution of two related statistics. The first one is a quadratic form in deviations from a control group with which the remaining k-1 groups are to be compared. The second test is designed for comparing adjacent groups. Several approximations for the null and the alternative distribution are considered and two illustrative examples are provided.     

Testing Hypotheses in the Functional Linear Model

The functional linear model with scalar response is a regression model where the predictor is a random function defined on some compact set of R and the response is scalar. The response is modelled as Y =Ψ( X )+ ɛ , where Ψ is some linear continuous operator defined on the space of square integrable functions and valued in R . The random input X is independent from the noise ɛ . In this paper, we are interested in testing the null hypothesis of no effect, that is, the nullity of Ψ restricted to the Hilbert space generated by the random variable X . We introduce two test statistics based on the norm of the empirical cross-covariance operator of ( X , Y ). The first test statistic relies on a χ ² approximation and we show the asymptotic normality of the second one under appropriate conditions on the covariance operator of X . The test procedures can be applied to check a given relationship between X and Y . The method is illustrated through a simulation study.     

LM Test of Neglected Correlated Random Effects and Its Application

This article aims at achieving two distinct goals. The first is to extend the existing LM test of overdispersion to the situation where the alternative hypothesis is characterized by the correlated random effects model. We obtain a result that the test against the random effects model has a certain max-min type optimality property. We will call such a test the LM test of overdispersion. The second goal of the article is to draw a connection between panel data analysis and the analysis of multiplicity of equilibrium in games. Because such multiplicity can be viewed as a particular form of neglected heterogeneity, we propose an intuitive specification test for a class of two-step game estimators.     

PROPERTIES OF THE MINIMUM POINT OF AN UNBALANCED TWO-SIDED RANDOM WALK

In this paper, properties of minimum point of a unbalanced two-sided random walk are investigated. Under the condition that the parameters at both sides tend to zero at the same order, probabilities that the minimum point is on which side, and the second order expansions for the first two moments of the minimum point are obtained. Applications of these results are very promising. First, they can be used to study the properties of the maximum likelihood estimator for the change point in the large sample case; second, they can be used to study inference problems after CUSUM test.     

A test for second order stochastic dominance

In this paper we propose a test for second order stochastic dominance (SSD), for the case where both distribution functions are unknown. This is a generalization of a test proposed by Deshpande and Singh (1985), who compare a new random prospect with a known distribution function. We then show that our test is based on comparing the mean minus one half of Gini's mean difference of the distributions, which is known to be a necessary condition for SSD, as developed in the economics literature (Yitzhaki, 1982).     

ON CHARACTERIZING SOME BIVARIATE DISCRETE DISTRIBUTIONS

The bivariate (correlated) Poisson, binomial, negative binomial and logarithmic distributions are characterized by the conditional distribution of one random vector on the other and a regression function of the second random vector on the first. Characterizations of the above distributions when their component random variables are independent are also included.     

Tests for noncorrelation of two multivariate ARMA time series

In many situations, we want to verify the existence of a relationship between multivariate time series. In this paper, we generalize the procedure developed by Haugh (1976) for univariate time series in order to test the hypothesis of noncorrelation between two multivariate stationary ARMA series. The test statistics are based on residual cross-correlation matrices. Under the null hypothesis of noncorrelation, we show that an arbitrary vector of residual cross-correlations asymptotically follows the same distribution as the corresponding vector of cross-correlations between the two innovation series. From this result, it follows that the test statistics considered are asymptotically distributed as chi-square random variables. Two test procedures are described. The first one is based on the residual cross-correlation matrix at a particular lag, whilst the second one is based on a portmanteau type statistic that generalizes Haugh's statistic. We also discuss how the procedures for testing noncorrelation can be adapted to determine the directions of causality in the sense of Granger (1969) between the two series. An advantage of the proposed procedures is that their application does not require the estimation of a global model for the two series. The finite-sample properties of the statistics introduced were studied by simulation under the null hypothesis. It led to modified statistics whose upper quantiles are much better approximated by those of the corresponding chi-square distribution. Finally, the procedures developed are applied to two different sets of economic data.     

Testing for departures from nominal dispersion in generalized nonlinear models with varying dispersion and/or additive random effects

This paper discusses the tests for departures from nominal dispersion in the framework of generalized nonlinear models with varying dispersion and/or additive random effects. We consider two classes of exponential family distributions. The first is discrete exponential family distributions, such as Poisson, binomial, and negative binomial distributions. The second is continuous exponential family distributions, such as normal, gamma, and inverse Gaussian distributions. Correspondingly, we develop a unifying approach and propose several tests for testing for departures from nominal dispersion in two classes of generalized nonlinear models. The score test statistics are constructed and expressed in simple, easy to use, matrix formulas, so that the tests can easily be implemented using existing statistical software. The properties of test statistics are investigated through Monte Carlo simulations.     

A statistical approach to the resolution of point sources

SUMMARY The Rayleigh criterion in optics states that two point sources of equal intensity are 'barely resolved' when the maximum of the diffraction pattern of one source overlaps the first minimum of the diffraction pattern of the second source. Although useful for rough comparisons of optical systems, such a criterion does not take into account the randomness in the detection process and does not tell whether sources can actually be distinguished. We present a statistical approach that addressed these issues. From quantum optics, the photon counts in the pixels are independent Poisson random variables with means that depend on the distance 2theta between the sources. Resolving the sources corresponds to testing H0: theta =0 vs Ha: theta >0, under conditions that make the information number zero at theta =0. We define resolution as the (asymptotic) power function of the likelihood ratio test rather than as a single number. The asymptotic distribution of the test statistic is derived under H0 and under contiguous alternatives. The results are illustrated by an application to a sky survey to detect binary stars using the Hubble space telescope.     

Data-driven stochastic inversion via functional quantization

In this paper, we propose a new methodology for solving stochastic inversion problems through computer experiments, the stochasticity being driven by a functional random variables. This study is motivated by an automotive application. In this context, the simulator code takes a double set of simulation inputs: deterministic control variables and functional uncertain variables. This framework is characterized by two features. The first one is the high computational cost of simulations. The second is that the probability distribution of the functional input is only known through a finite set of realizations. In our context, the inversion problem is formulated by considering the expectation over the functional random variable. We aim at solving this problem by evaluating the model on a design, whose adaptive construction combines the so-called stepwise uncertainty reduction methodology with a strategy for an efficient expectation estimation. Two greedy strategies are introduced to sequentially estimate the expectation over the functional uncertain variable by adaptively selecting curves from the initial set of realizations. Both of these strategies consider functional principal component analysis as a dimensionality reduction technique assuming that the realizations of the functional input are independent realizations of the same continuous stochastic process. The first strategy is based on a greedy approach for functional data-driven quantization, while the second one is linked to the notion of space-filling design. Functional PCA is used as an intermediate step. For each point of the design built in the reduced space, we select the corresponding curve from the sample of available curves, thus guaranteeing the robustness of the procedure to dimension reduction. The whole methodology is illustrated and calibrated on an analytical example. It is then applied on the automotive industrial test case where we aim at identifying the set of control parameters leading to meet the pollutant emission standards of a vehicle.     

Estimating reference ranges from stratified two-stage samples

We present methodology for estimating age-specific reference ranges by using data from two-stage samples. On the basis of the information obtained in the first stage, the initial sample is stratified and random subsamples are drawn from each stratum, where the selection probabilities in this second-stage sampling may be different across strata in the population. The variable for which the reference ranges are to be established is measured at the second phase. The approach involves maximum likelihood estimation of the parameters of the age-specific distributions and separate estimation of the population stratum probabilities. These are combined to yield estimates of the quantiles of interest. The issue of variance estimation for the estimated quantiles is also addressed. The methodology is applied to the estimation of reference ranges for a cognitive test score in a study of non-demented older Japanese-Americans.     

On Minimum Variance Unbiased Estimators

The problem of finding minimum variance unbiased estimators of various parameters for parametric distributions is an important one in statistics. This article gives analytical formulas for the minimum variance unbiased estimators of parametric functions, which are usually used in a classroom, for two types of densities. The first type is the one-parameter regular exponential family, and the second is a two-parameter family of a continuous random variable whose range depends on the unknown parameters.     

Limit theory of bivariate dual generalized order statistics with random index

The class of limit distribution functions of bivariate extreme, intermediate and central dual generalized order statistics from independent and identically distributed random variables with random sample size is fully characterized. Two cases are considered. The first case is when the random sample size is assumed to be independent of all basic random variables. The second case is when the interrelation of the random size and the basic random variables is not restricted.     

Using -splines to test the linearity of partially linear models

The nonparametric component in a partially linear model is estimated by a linear combination of fixed-knot cubic B-splines with a second-order difference penalty on the adjacent B-spline coefficients. The resulting penalized least-squares estimator is used to construct two Wald-type spline-based test statistics for the null hypothesis of the linearity of the nonparametric function. When the number of knots is fixed, the first test statistic asymptotically has the distribution of a linear combination of independent chi-squared random variables, each with one degree of freedom, under the null hypothesis. The smoothing parameter is determined by specifying a value for the asymptotically expected value of the test statistic under the null hypothesis. When the number of knots is fixed and under the null hypothesis, the second test statistic asymptotically has a chi-squared distribution with K=q+2 degrees of freedom, where q is the number of knots used for estimation. The power performances of the two proposed tests are investigated via simulation experiments, and the practicality of the proposed methodology is illustrated using a real-life data set.     

Analyzing multiple end points with a two-stage design in clinical trials

In many clinical trials, the assessment of the response to interventions can include a large variety of outcome variables which are generally correlated. The use of multiple significance tests is likely to increase the chance of detecting a difference in at least one of the outcomes between two treatments. Furthermore, univariate tests do not take into account the correlation structure. A new test is proposed that uses information from the interim analysis in a two-stage design to form the rejection region boundaries at the second stage. Initially, the test uses Hotelling’s T² at the end of the first stage allowing only, for early acceptance of the null hypothesis and an O’Brien ‘type’ procedure at the end of the second stage. This test allows one to ‘cheat’ and look at the data at the interim analysis to form rejection regions at the second stage, provided one uses the correct distribution of the final test statistic. This distribution is derived and the power of the new test is compared to the power of three common procedures for testing multiple outcomes: Bonferroni’s inequality, Hotelling’s T²and O’Brien’s test. O’Brien’s test has the best power to detect a difference when the outcomes are thought to be affected in exactly the same direction and the same magnitude or in exactly the same relative effects as those proposed prior to data collection. However, the statistic is not robust to deviations in the alternative parameters proposed a priori, especially for correlated outcomes. The proposed new statistic and the derivation of its distribution allows investigators to consider information from the first stage of a two-stage design and consequently base the final test on the direction observed at the first stage or modify the statistic if the direction differs significantly from what was expected a prior.     

Normal probability plots with confidence for the residuals in linear regression

Normal probability plots for a simple random sample and normal probability plots for residuals from linear regression are not treated differently in statistical text books. In the statistical literature, 1 ? α simultaneous probability intervals for augmenting a normal probability plot for a simple random sample are available. The first purpose of this article is to demonstrate that the tests associated with the 1 ? α simultaneous probability intervals for a simple random sample may have a size substantially different from α when applied to the residuals from linear regression. This leads to the second purpose of this article: construction of four normal probability plot-based tests for residuals, which have size α exactly. We then compare the powers of these four graphical tests and a non-graphical test for residuals in order to assess the power performances of the graphical tests and to identify the ones that have better power. Finally, an example is provided to illustrate the methods.     

Modified entropy estimators for testing normality

In this paper, we present two new estimators for the entropy of absolutely continuous random variables and consider some of their properties. Consistency of the first estimator is shown by Monte Carlo method, and the consistency of the second estimator is proved theoretically. Using these estimators, two new tests for normality are presented and their powers are compared with the other entropy-based tests. Simulation results show that the proposed estimators and test statistics perform very well. Finally, a real example is presented and analysed.     

Failure Time of Non Homogeneous Gamma Process

We consider the non homogeneous gamma process as a degradation model. The hitting time of a deterministic or random level is studied here. We provide its distribution (both cdf and pdf) explicitly in the first case and in the second case when the threshold is exponentially or gamma distributed. The general case for a random threshold can be approximated by considering mixtures of Erlang distributions. Aging properties are also discussed in this article.     

Power of analysis of covariance tests under conditional specification

The procedure-wise power functions of two strategies for balanced single-factor analysis of covariance in the presence of possibly unequal regression slopes are evaluated and illustrated. The strategies differ in the action to be taken following a re-^ jection by the preliminary test for equal slopes. The first strategy simply discards the covariate and respecifies the model as the one-way ANOVA model for testing factor effects. The second leaves the unequal slopes covariance model intact, but respecifies the factor effects hypothesis to address the factor level means adjusted to the sample average of the covariate. One additional strategy, that of testing factor effects only if the preliminary slopes test does not reject, is included for comparison purposes. Computation of the power functions requires extensive use of the results obtained in Hawkins and Han (1986) concerning the bivariate distributions of certain ratios of independent noncentral chi-square random variables.