Bimodal skew-symmetric normal distribution

ABSTRACT

We introduce a new parsimonious bimodal distribution, referred to as the bimodal skew-symmetric Normal (BSSN) distribution, which is potentially effective in capturing bimodality, excess kurtosis, and skewness. Explicit expressions for the moment-generating function, mean, variance, skewness, and excess kurtosis were derived. The shape properties of the proposed distribution were investigated in regard to skewness, kurtosis, and bimodality. Maximum likelihood estimation was considered and an expression for the observed information matrix was provided. Illustrative examples using medical and financial data as well as simulated data from a mixture of normal distributions were worked.     

Gammaization and wishartness of dependent quadratic forms

In this paper two equivalent sets of necessary and sufficient conditions are derived for dependent quadratic forms to be distributed as multivariate gamma distribution. The procedure also gives a set of necessary and sufficient conditions for principal minors of generalized quadratic forms to be jointly distributed as the joint distribution of principal minors of a Kishart matrix.     

On the choice of co-ordinates in simultaneous estimation of normal means under misspecification of normal priors

The problem of choice of coordinates in Stein-type estimators,when simultaneously estimating normal means, is considered. The question of deciding whether to use all coordinates in one combined shrinkage estimators or to separate into groups and use separate shrinkage estimators on each group is considered in the situation in which part of the prior information may be " misspecified". It is observed that the amount of misspecification determines whether to use the combined shrinkage estimator the separate shrinkage estimator.     

Optimal minimax squared error risk estimation of the mean of a multivariate normal distribution

Minimax squared error risk estimators of the mean of a multivariate normal distribution are characterized which have smallest Bayes risk with respect to a spherically symmetric prior distribution for (i) squared error loss, and (ii) zero-one loss depending on whether or not estimates are consistent with the hypothesis that the mean is null. In (i), the optimal estimators are the usual Bayes estimators for prior distributions with special structure. In (ii), preliminary test estimators are optimal. The results are obtained by applying the theory of minimax-Bayes-compromise decision problems.     

A test for variance-covarianch parameters in normal linear models

This paper derives a test statistic for the variance-covariance parameters which is a quadratic function of their MINQUE (Minimum Norm Quadratic Unbiased Estimation) estimates. The test is a Wald-type test, and its development closely parallels the theory used to derive a similar test for the coefficients in linear models. In fact, the derivation proceeds by first setting up the estimation problem in a derived linear model in which the dispersion parameters are the coefficients. The test statistic is shown to be the sum of the squares of independent standardized x² variables.     

Some further aspects of the löwner-ordering antitonicity of the moore-penrose inverse

A survey is given of known proofs of the antitonicity of the inverse matrix function for positive definite matrices w.r.t. the Lowner partial ordering, and of the corresponding result for the Moore-Penrose inverse of nonnegative definite matrices [the theorem of Milliken and Akdeniz (1977)]. A short new proof of the latter result is obtained by employing an extremal representation of a nonnegative definite quadratic form. Another proof of this result involving Schur complements is also given, and is seen to be extendable to the case of symmetric (not necessarily nonnegative definite) matrices. A geometrical interpretation of Milliken and Akdeniz's theorem is presented. As an application, the relationship between the concepts of greater (maximum) concentration and smaller (minimum) dispersion is considered for a pair (class) of vector-valued statistics with possibly degenerate distributions.     

Quadratic perturbation expansions of certain functions of eigenvalues and eigenvectors and their application to sensitivity analysis in multivariate methods

Tanaka (1988) lias derived the influence functions, which are equivalent to the perturbation expansions up to linear terms, of two functions of eigenvalues and eigenvectors of a real symmetric matrix, and applied them to principal component analysis. The present paper deals with the perturbation expansions up to quadratic terms of the same functions and discusses their application to sensitivity analysis in multivariate methods, in particular, principal component analysis and principal factor analysis. Numerical examples are given to show how the approximation improves with the quadratic terms.     

Discrimination of AR,MA and ARMA time series models

The problem of discrimination between two stationary ARMA time series models is considered, and in particular AR(p), MA(p), ARMA(1,1) models. The discriminant based on the likelihood ration leads to a quadratic form that is generally too complicated to evaluated explicitly. The discriminant can be expressed approximately as a linear combination of independent chi–squared random varianles each with one degree of freedom, the coefficients, of which are eigenvalues of cumbersome matrices. An analytical solution which gives the coefficients approximately is suggested.