Noncentral matrix quadratic forms of the skew elliptical variables

In this paper, the noncentral matrix quadratic forms of the skew elliptical variables are studied. A family of the matrix variate noncentral generalized Dirichlet distributions is introduced as the extension of the noncentral Wishart distributions, the Dirichlet distributions and the noncentral generalized Dirichlet distributions. Main distributional properties are investigated. These include probability density and closure property under linear transformation and marginalization, the joint distribution of the sub-matrices of the matrix quadratic forms in the skew elliptical variables and the moment generating functions and Bartlett's decomposition of the matrix quadratic forms in the skew normal variables. Two versions of the noncentral Cochran's Theorem for the matrix variate skew normal distributions are obtained, providing sufficient and necessary conditions for the quadratic forms in the skew normal variables to have the matrix variate noncentral generalized Dirichlet distributions. Applications include the properties of the least squares estimation in multivariate linear model and the robustness property of the Wilk's likelihood ratio statistic in the family of the matrix variate skew elliptical distributions.     

Normal theory likelihood ratio statistic for mean and covariance structure analysis under alternative hypotheses

The normal distribution based likelihood ratio (LR) statistic is widely used in structural equation modeling. Under a sequence of local alternative hypotheses, this statistic has been shown to asymptotically follow a noncentral chi-square distribution. In practice, the population mean vector and covariance matrix as well as the model and sample size are always fixed. It is hard to justify the validity of the noncentral chi-square distribution for the resulting LR statistic even when data are normally distributed and sample size is large. By extending results in the literature, this paper develops normal distributions to describe the behavior of the LR statistic for mean and covariance structure analysis. A sequence of local alternative hypotheses is not necessary for the proposed distributions to be asymptotically valid. When the effect size is medium and above or when the model is not trivially misspecified, empirical results indicate that a refined normal distribution describes the behavior of the LR statistic better than the commonly used noncentral chi-square distribution, as measured by the Kolmogorov-Smirnov distance. Quantile-quantile plots are also provided to better understand the different distributions.     

Inferences in linear mixed models with skew-normal random effects

For the linear mixed model with skew-normal random effects, this paper gives the density function, moment generating function and independence conditions. The noncentral skew chi-square distribution is defined and its density function is shown. The necessary and sufficient conditions under which a quadratic form is distributed as noncentral skew chi-square distribution are obtained. Also, a version of Cochran's theorem is given, which modifies the result of Wang et al. (2009) and is used to set up exact tests for fixed effects and variance components of the proposed model. For illustration, our main results are applied to a real data problem.     

Noncentral quadratic forms of the skew elliptical variables

In this paper the quadratic forms in the skew elliptical variables are studied. A family of the noncentral generalized Dirichlet distributions is introduced and their distribution functions and probability density functions are obtained. The moment generating functions of the quadratic forms in the skew normal variables are obtained. Sufficient and necessary conditions for the quadratic forms in the skew normal variables to have the noncentral generalized Dirichlet distributions are obtained. This leads to the noncentral Cochran's Theorem for the skew normal distribution.     

A parametric bootstrap approach for one-way classification model with skew-normal random effects

In this paper, several properties of one-way classification model with skew-normal random effects are obtained, such as moment generating function, density function and noncentral skew chi-square distribution, etc. Based on the EM algorithm, we discuss the maximum likelihood(ML) estimation of unknown parameters. For testing problem of fixed effect, a parametric bootstrap(PB) approach is developed. Finally, some simulation results on the Type I error rates and powers of the PB approach are obtained, which show that the PB approach provides satisfactory performances on the Type I error rates and powers, even for small samples.For illustration, our main results are applied to a real data problem.     

The bivariate noncentral chi-square distribution - A compound distribution approach

This paper proposes the bivariate noncentral chi-square (BNC) distribution by compounding the Poisson probabilities with the bivariate central chi-square distribution. The probability density and cumulative distribution functions of the joint distribution of the two noncentral chi-square variables are derived for arbitrary values of the correlation coefficient, degrees of freedom(s), and noncentrality parameters. Computational procedures to calculate the upper tail probabilities as well as the percentile points for selected values of the parameters, for both equal and unequal degrees of freedom, are discussed. The graphical representation of the distribution for different values of the parameters are provided. Some applications of the distribution are outlined.     

Nonparametric tests for conditional independence in two-way contingency tables

Testing for the independence between two categorical variables R and S forming a contingency table is a well-known problem: the classical chi-square and likelihood ratio tests are used. Suppose now that for each individual a set of p characteristics is also observed. Those explanatory variables, likely to be associated with R and S, can play a major role in their possible association, and it can therefore be interesting to test the independence between R and S conditionally on them. In this paper, we propose two nonparametric tests which generalise the chi-square and the likelihood ratio ideas to this case. The procedure is based on a kernel estimator of the conditional probabilities. The asymptotic law of the proposed test statistics under the conditional independence hypothesis is derived; the finite sample behaviour of the procedure is analysed through some Monte Carlo experiments and the approach is illustrated with a real data example.     

Equivalent conditions for noncentral generalized Laplacianness and independence of matrix quadratic forms

Let Y be an n×p multivariate normal random matrix with general covariance Σ_Y and W be a symmetric matrix. In the present article, the property that a matrix quadratic form Y^′WY is distributed as a difference of two independent (noncentral) Wishart random matrices is called the (noncentral) generalized Laplacianness (GL). Then a set of algebraic results are obtained which will give the necessary and sufficient conditions for the (noncentral) GL of a matrix quadratic form. Further, two extensions of Cochran’s theorem concerning the (noncentral) GL and independence of a family of matrix quadratic forms are developed.     

Nonnull distribution of likelihood ratio criterion for reality of covariance matrix

In this paper the distribution of the likelihood ratio test for testing the reality of the covariance matrix of a complex multivariate normal distribution is investigated. Some simplifications in the noncentral distribution are made and the noncentral distribution is derived for the special case where the rank of the noncentrality matrix is two. In the null case exact expressions for the distribution are given up to p = 6, and percentage points are tabulated. These percentage points were compared with percentage points derived from an asymptotic expansion of the distribution, and the accuracy of the approximation was found to be sufficient for several practical situations.     

On some testing problems for sparse contingency tables

A class of test statistics are proposed for sparse tables with ordered categories. It is shown that for different testing situations these test statistics asymptotically more powerful tests than Pearson's chi-square.     

Exact probabilities of correct classifications for uncorrelated repeated measurements from elliptically contoured distributions

Euclidean distance-based classification rules are derived within a certain nonclassical linear model approach and applied to elliptically contoured samples having a density generating function g. Then a geometric measure theoretical method to evaluate exact probabilities of correct classification for multivariate uncorrelated feature vectors is developed. When doing this one has to measure suitably defined sets with certain standardized measures. The geometric key point is that the intersection percentage functions of the areas under investigation coincide with those of certain parabolic cylinder type sets. The intersection percentage functions of the latter sets can be described as threefold integrals. It turns out that these intersection percentage functions yield simultaneously geometric representation formulae for the doubly noncentral g-generalized F-distributions. Hence, we get beyond new formulae for evaluating probabilities of correct classification new geometric representation formulae for the doubly noncentral g-generalized F-distributions. A numerical study concerning several aspects of evaluating both probabilities of correct classification and values of the doubly noncentral g-generalized F-distributions demonstrates the advantageous computational properties of the present new approach. This impression will be supported by comparison with the literature.It is shown that probabilities of correct classification depend on the parameters of the underlying sample distribution through a certain well-defined set of secondary parameters. If the underlying parameters are unknown, we propose to estimate probabilities of correct classification.     

Analytic,Computational, and Approximate Forms for Ratios of Noncentral and Central Gaussian Quadratic Forms

Many useful statistics equal the ratio of a possibly noncentral chi-square to a quadratic form in Gaussian variables with all positive weights. Expressing the density and distribution function as positively weighted sums of corresponding F functions has many advantages. The mixture forms have analytic value when embedded within a more complex problem. The mixture forms also have computational value. The expansions work well with quadratic forms having few components and small degrees of freedom. A more general algorithm from earlier literature can take longer or fail to converge in the same setting. Many approximations have been suggested for the problem. A positively weighted noncentral quadratic form can always have two moments matched to a noncentral chi-square. For a single quadratic form, the noncentral form performs neither uniformly more or less accurately than older approximations. The approach also gives a noncentral F approximation for any ratio of a positively weighted noncentral form to a positively weighted central quadratic form. The method provides better accuracy for noncentral ratios than approximations based on a single chi-square. The accuracy suffices for many practical applications, such as power analysis, even with few degrees of freedom. Naturally the approximation proves much faster and simpler to compute than any exact method. Embedding the approximation in analytic expressions provides simple forms which correctly guarantee only positive values have nonzero probabilities, and also automatically reduce to partially or fully exact results when either quadratic form has only one term.     

Graph presentations for moments of noncentral Wishart distributions and their applications

We provide formulas for the moments of the real and complex noncentral Wishart distributions of general degrees. The obtained formulas for the real and complex cases are described in terms of the undirected and directed graphs, respectively. By considering degenerate cases, we give explicit formulas for the moments of bivariate chi-square distributions and 2 × 2 Wishart distributions by enumerating the graphs. Noting that the Laguerre polynomials can be considered to be moments of a noncentral chi-square distributions formally, we demonstrate a combinatorial interpretation of the coefficients of the Laguerre polynomials.     

Trimmed portmanteau test for linear processes with infinite variance

In this article, we consider a model check test for linear processes with infinite variance. As a test statistic, we employ the portmanteau test with trimmed residuals. It is shown that the limiting null distribution of the test is a chi-square distribution. Simulation results are provided for illustration.     

Asymptotic distributions for the elementary symmetric functions of two matrices under the assumption of linearity

The asymptotic distributions of the elementary symmetric functions (esf's) of the characteristic roots of a noncentral multivariate beta matrix and of the generalized correlation matrix (noncentral under the assumption of linearity) are derived.     

Exact nonnull distribution of Wilks’ statistic: The ratio and product of independent components

The study of the noncentral matrix variate beta type distributions has been sidelined because the final expressions for the densities depend on an integral that has not been resolved in an explicit way. We derive an exact expression for the nonnull distribution of Wilks’ statistic and precise expressions for the densities of the ratio and product of two independent components of matrix variates where one matrix variate has the noncentral matrix variate beta type I distribution and the other has the matrix variate beta type I distribution. We provide the expressions for the densities of the determinant of the ratio and the product of these two components. These distributions play a fundamental role in various areas of statistics, for example in the criteria proposed by Wilks.     

Asymptotic optimal inference for multivariate branching-Markov processes via martingale estimating functions and mixed normality

Multivariate tree-indexed Markov processes are discussed with applications. A Galton-Watson super-critical branching process is used to model the random tree-indexed process. Martingale estimating functions are used as a basic framework to discuss asymptotic properties and optimality of estimators and tests. The limit distributions of the estimators turn out to be mixtures of normals rather than normal. Also, the non-null limit distributions of standard test statistics such as Wald, Rao’s score, and likelihood ratio statistics are shown to have mixtures of non-central chi-square distributions. The models discussed in this paper belong to the local asymptotic mixed normal family. Consequently, non-standard limit results are obtained.     

Bayesian Inference in the Noncentral Student-t Model

This article takes up Bayesian inference in linear models with disturbances from a noncentral Student-t distribution. The distribution is useful when both long tails and asymmetry are features of the data. The distribution can be expressed as a location-scale mixture of normals with inverse weights distributed according to a chi-square distribution. The computations are performed using Gibbs sampling with data augmentation. An empirical application to Standard and Poor's stock returns indicates that posterior odds strongly favor a noncentral Student-t specification over its symmetric counterpart.     

Estimation of the eigenvalues of noncentrality parameter matrix in noncentral Wishart distribution

We consider the problem of estimating the eigenvalues of noncentrality parameter matrix in noncentral Wishart distribution when the scale parameter is known. A decision theoretic approach is taken with squared error as the loss function. We propose two new estimators and show their superior performance to an usual estimator theoretically and numerically.     

The exact distribution of indefinite quadratic forms in noncentral normal vectors

The exact density of the difference of two linear combinations of independent noncentral chi-square variables is obtained in terms of Whittaker's function and expressed in closed forms. Two distinct representations are required in order to cover all the possible cases. The corresponding expressions for the exact distribution function are also given.