On the distribution of quadratic forms in normal random variables

This paper provides necessary and sufficient conditions for a quadratic form in singular normal random variables to be distributed as a given linear combination of independent noncentral chi-square variables. Using this result, an extension of Cochran's theorem to quadratic forms of noncentral chi-square variables is derived.     

On the distribution of quadratic forms and applications in the analysis of repeated measurements

In this paper a unified approach is given to the distribution of scalar quadratic forms for dependent variables. Necessary and sufficient conditions are found for the sums of squares of the various hierarchical layers in ANOVA to be distributed as multiples of chi-square variables. Results concerning the usual univariate F-tests in ANOVA of repeated measurements are derived as a special case.     

Distribution of multivariate quadratic forms under certain covariance structures

Necessary and sufficient conditions are given for the covariance structure of all the observations in a multivariate factorial experiment under which certain multivariate quadratic forms are independent and distributed as a constant times a Wishart. It is also shown that exact multivariate test statistics can be formed for certain covariance structures of the observations when the assumption of equal covariance matrices for each normal population is relaxed. A characterization is given for the dependency structure between random vectors in which the sample mean and sample covariance matrix have certain properties.     

Minimum variance unbiased invariant estimation of variance components under normality

In a variance components model for normally distributed data, for a specified vector of linear combinations of the variance components, necessary and sufficient conditions are given under which the vector has a uniformly minimum variance unbiased translation-invariant estimator. The competing class of estimators is not restricted to those that are quadratic. For classification models, the conditions are translated into easy-to-check partial balance requirements on the incidence array.     

Moments of generalized quadratic forms and distributions of certain multivariate test statistics

Summary Moments and distributions of quadratic forms or quadratic expressions in normal variables are available in literature. Such quadratic expressions are shown to be equivalent to a linear function of independent central or noncentral chi-square variables. Some results on linear functions of generalized quadratic forms are also available in literature. Here we consider an arbitrary linear function of matrix-variate gamma variables. Moments of the determinant of such a linear function are evaluated when the matrix-variate gammas are independently distributed. By using these results, arbitrary non-null moments as well as the non-null distribution of the likelihood ratio criterion for testing the hypothesis of equality of covariance matrices in independent multivariate normal populations are derived. As a related result, the distribution of a linear function of independent matrix-variate gamma random variables, which includes linear functions of independent Wishart matrices, is also obtained. Some properties of generalized special functions of several matrix arguments are used in deriving these results.     

Convergence of the moment generating function of standardized quantiles

A sequence of independent, identically distributed random variables is considered. Given a simple local condition on the distribution of these random variables, we give necessary and sufficient conditions on the tails of the distribution for the moment generating function of a standardized quantile of the first n observations to converge to the moment generating function of an appropriate normal distribution as n →infinity;. This result is actually a special case of a more general result which can also be used to show convergence in distribution and convergence of moments of standardized quantiles.     

The n.s. conditions for the ration of two quadratic forms to have an f-distribution and its applications

Suppose that ξ and η be two random vectors and that (ξ^τ, η^τ have an elliptically contoured distribution or a multivariate normal distribution. In this article, we obtain some necessary and sufficient (N.S.) conditions such that the ratio of two quadratic forms, say ξ^τ Aξ and η^τ Bη(for some symmetric nonnegative matrices A and B), has an F-distribution. As applications, we extend the classical F-test to some dependent two group samples. Two cases are considered: elliptically contoured and normal distributions.     

A reformulation of the criteria for the independence of quadratic functions in normal variables

The Craig-Sakamoto theorem establishes a sufficient and necessary condition for the independence of two quadratic forms in normal variates, fascinating many statisticians and mathematicians, who continuously seek for simple and better proofs of the theorem and its extensions. In this article, we present a simple proof of a unified theorem on the independence of linear and quadratic functions in general normal variates.     

On the Admissibility of Linear Estimators in a Multivariate Normal Distribution Under LINEX Loss Function

Consider an estimation problem of a linear combination of population means in a multivariate normal distribution under LINEX loss function. Necessary and sufficient conditions for linear estimators to be admissible are given. Further, it is shown that the result is an extension of the quadratic loss case as well as the univariate normal case.     

A comparison of mixed and minimax estimators of linear models

In this paper the stochastic properties of two estimators of linear models, mixed and minimax, based on different types of prior information, are compared using quadratic risk as the criterion for superiority. A necessary and sufficient condition for the minimax estimator to be superior to the comparable mixed estimator is derived as well as a simpler necessary but not sufficient condition.     

A Matrix Variate Closed Skew-Normal Distribution with Applications to Stochastic Frontier Analysis

In this article, we introduce the matrix extension of the closed skew-normal distribution and give two constructions for it: a marginal one and another based on hidden truncation. Important basic properties of the distribution are presented such as its closure under linear transformation and moment generating function. We also give distributional results for quadratic forms involving random matrices distributed according to two particular cases of it. Using an additive construction, we derive a submodel which can be employed to describe the compound error structure of a very general multivariate stochastic frontier model. Finally, we consider the skew-elliptical extension of the proposed distribution.     

An Improved Result Relating Quadratic Forms and Chi-Square Distributions

Results from the theory of linear models establish a particular idempotency condition as being necessary and sufficient for a quadratic form in a nonsingular normal vector to follow a chi-square distribution. We give a theorem that is somewhat stronger than the standard ones, and provide a proof that is more accessible than those usually given, in that it uses only linear algebra and calculus.     

On approximating the distribution of indefinite quadratic forms

This paper provides a simple methodology for approximating the distribution of indefinite quadratic forms in normal random variables. It is shown that the density function of a positive definite quadratic form can be approximated in terms of the product of a gamma density function and a polynomial. An extension which makes use of a generalized gamma density function is also considered. Such representations are based on the moments of a quadratic form, which can be determined from its cumulants by means of a recursive formula. After expressing an indefinite quadratic form as the difference of two positive definite quadratic forms, one can obtain an approximation to its density function by means of the transformation of variable technique. An explicit representation of the resulting density approximant is given in terms of a degenerate hypergeometric function. An easily implementable algorithm is provided. The proposed approximants produce very accurate percentiles over the entire range of the distribution. Several numerical examples illustrate the results. In particular, the methodology is applied to the Durbin–Watson statistic which is expressible as the ratio of two quadratic forms in normal random variables. Quadratic forms being ubiquitous in statistics, the approximating technique introduced herewith has numerous potential applications. Some relevant computational considerations are also discussed.     

A simple condition for the multivariate CLT and the attraction to the Gaussian of Lévy processes at long and short times

We show that a necessary and sufficient condition for the sum of iid random vectors to converge (under appropriate centering and scaling) to a multivariate Gaussian distribution is that the truncated second moment matrix is slowly varying at infinity. This is more natural than the standard conditions, and allows for the possibility that the limiting Gaussian distribution is degenerate (so long as it is not concentrated at a point). We also give necessary and sufficient conditions for a d-dimensional Lévy process to converge (under appropriate centering and scaling) to a multivariate Gaussian distribution as time approaches zero or infinity.     

Admissibilities of matrix linear estimators multivariate linear models

This article respectively provides sufficient conditions and necessary conditions of matrix linear estimators of an estimable parameter matrix linear function in multivariate linear models with and without the assumption that the underlying distribution is a normal one with completely unknown covariance matrix. In the latter model, a necessary and sufficient condition is given for matrix linear estimators to be admissible in the space of all matrix linear estimators under each of three different kinds of quadratic matrix loss functions, respectively. In the former model, a sufficient condition is first provided for matrix linear estimators to be admissible in the space of all matrix estimators having finite risks under each of the same loss functions, respectively. Furthermore in the former model, one of these sufficient conditions, correspondingly under one of the loss functions, is also proved to be necessary, if additional conditions are assumed.     

The Whittle Estimator for Strongly Dependent Stationary Gaussian Fields

In this article we generalize results on the asymptotic behaviour of the Whittle estimator for certain stationary Gaussian long range dependent fields. These results have been established in the one-dimensional case under very general conditions. They require controlling the estimation bias and also giving convergence theorems for certain quadratic forms of the observations. In the multidimensional setting, our main interest will be controlling the bias. This can be done for d ≤ 3 using taper functions, and, depending on the shape of the singularity, also introducing certain regularizing functions. In this last case, however, the estimator will no longer be efficient. We also present certain partial results concerning the convergence to a limiting Gaussian distribution of the associated quadratic forms.     

Best unbiased estimation of fixed effects in mixed ANOVA models

In a linear model with an arbitrary variance–covariance matrix, Zyskind (Ann. Math. Statist. 38 (1967) 1092) provided necessary and sufficient conditions for when a given linear function of the fixed-effect parameters has a best linear unbiased estimator (BLUE). If these conditions hold uniformly for all possible variance–covariance parameters (i.e., there is a UBLUE) and if the data are assumed to be normally distributed, these conditions are also necessary and sufficient for the parametric function to have a uniformly minimum variance unbiased estimator (UMVUE). For mixed-effects ANOVA models, we show how these conditions can be translated in terms of the incidence array, which facilitates verification of the UBLUE and UMVUE properties and facilitates construction of designs having such properties.     

Characterization of discrete populations through conditional expectations of order statistics

In this paper we give some properties of the expected values of any order statistic when one of its adjacent order statistics is known (order mean function) from a sequence of sizen of independent and identically distributed random variables with discrete distribution. Furthermore, we obtain the explicit expressions of the distribution from these order mean functions, and finally, we show the necessary and sufficient conditions for any real function to be an order mean function. We also add some examples of characterization of discrete distributions from the order mean functions. Partially supported by Consejería de Cultura y Educación (C.A.R.M.), under Grant PIB 95/90.     

Characterization of weak convergence for smoothed empirical and quantile processes under ϕ-mixing

Let X_n, n ⩾ 1 be a sequence of ϕ-mixing random variables having a smooth common distribution function F. The smoothed empirical distribution function is obtained by integrating a kernel type density estimator. In this paper we provide necessary and sufficient conditions for the central limit theorem to hold for smoothed empirical distribution functions and smoothed sample quantiles. Also, necessary and sufficient conditions are given for weak convergence of the smoothed empirical process and the smoothed uniform quantile process.     

Estimation of the mean of a spherically symmetric distribution with constraints on the norm

The authors consider the problem of estimating, under quadratic loss, the mean of a spherically symmetric distribution when its norm is supposed to be known and when a residual vector is available. They give a necessary and sufficient condition for the optimal James‐Stein estimator to dominate the usual estimator. Various examples are given that are not necessarily variance mixtures of normal distributions. Consideration is also given to an alternative class of robust James‐Stein type estimators that take into account the residual vector. A more general domination condition is given for this class.