Lower dimensional marginal density functions of absolutely continuous truncated multivariate distributions

The lower dimensional marginal density functions of a truncated multivariate density function is derived in general, and shown that it is a function of untruncated marginal density function, appropriately defined conditional distribution function and size of the multivariate truncation region. As a special case, lower dimensional marginal density function of a truncated multivariate normal distribution is given.     

Diagonalization matrix and its application in distribution theory

Some matrix representations of diverse diagonal arrays are studied in this work; the results allow new definitions of classes of elliptical distributions indexed by kernels mixing Hadamard and usual products. A number of applications are derived in the setting of prior densities from the Bayesian multivariate regression model and families of non-elliptical distributions, such as the matrix multivariate generalized Birnbaum–Saunders density. The philosophy of the research about matrix representations of quadratic and inverse quadratic forms can be extended as a methodology for exploring possible new applications in non-standard distributions, matrix transformations and inference.     

Parameter estimation for multivariate diffusion processes with the time inhomogeneously positive semidefinite diffusion matrix

Statistical inference for the diffusion coefficients of multivariate diffusion processes has been well established in recent years; however, it is not the case for the drift coefficients. Furthermore, most existing estimation methods for the drift coefficients are proposed under the assumption that the diffusion matrix is positive definite and time homogeneous. In this article, we put forward two estimation approaches for estimating the drift coefficients of the multivariate diffusion models with the time inhomogeneously positive semidefinite diffusion matrix. They are maximum likelihood estimation methods based on both the martingale representation theorem and conditional characteristic functions and the generalized method of moments based on conditional characteristic functions, respectively. Consistency and asymptotic normality of the generalized method of moments estimation are also proved in this article. Simulation results demonstrate that these methods work well.     

Testing for lack of fit in linear multiresponse models based on exact or near replicates

Three procedures for testing the adequacy of a proposed linear multiresponse regression model against unspecified general alternatives are considered. The model has an error structure with a matrix normal distribution which allows the vector of responses for a particular run to have an unknown covariance matrix while the responses for different runs are uncorrelated. Furthermore, each response variable may be modeled by a separate design matrix. Multivariate statistics corresponding to the classical univariate lack of fit and pure error sums of squares are defined and used to determine the multivariate lack of fit tests. A simulation study was performed to compare the power functions of the test procedures in the case of replication. Generalizations of the tests for the case in which there are no independent replicates on all responses are also presented.     

Inequalities for tail probabilities for the multivariate normal distribution

Inequalities for tail probabilities of the multivariate normal distribution are obtained, as a generalization of those given by Feller (1966). Upper and lower bounds are given in the equi-correlated case. For an arbitrary correlation matrix R, an upper bound is obtained, using a result of Slepian (1962) which asserts that certain multivariate normal probabilities are a non-decreasing function of correlations.     

Bayesian multivariate normal analysis with a wishart prior

This paper considers the Bayesian analysis of the multivariate normal distribution when its covariance matrix has a Wishart prior density under the assumption of a multivariate quadratic loss function. New flexible marginal posterior distributions of the mean μ and of the covariance matrix Σ are developed and univariate cases with graphical representations are given.     

The exact noncentral distribution of the generalized multiple correlation matrix

Given p×n X N(βY, ∑?I), β, ∑ unknown, the noncentral multivariate beta density of the matrix L = [(YY′)^-1/2Y X′ (XX′)^-1XY′ (YY′)^-1/2] is desired. Khatri (1964) finds this density when β is of rank unity. The present paper derives the noncentral density of L and the density of the roots matrix of L for full rank β. The dual case density of L is also obtained. The derivations are based on generalized Sverdrup's lemma, Kabe (1965), and the relationship between primal and dual density of L is explicitly established.     

Multivariate lagrange distributions and a subfamiliy

Second order moments about its means, i.e. the variances and covari-ances for multivariate Lagrange distributions are derived in a matrix form. A subfamily of multivariate Lagrange distributions which can be characterized as the distributions of customers served in a busy period in queues with some conditions are considered. Theorems about their probability functions, one of which is a multivariate generalization of a formula by Takà cs(1989). are given and the means and second order moments about its means are considered. As an example, a multivariate Borel-Tanner distribution is derived.     

Identifiability of Finite Mixtures of Elliptical Distributions

Abstract.  We present general results on the identifiability of finite mixtures of elliptical distributions under conditions on the characteristic generators or density generators. Examples include the multivariate t -distribution, symmetric stable laws, exponential power and Kotz distributions. In each case, the shape parameter is allowed to vary in the mixture, in addition to the location vector and the scatter matrix. Furthermore, we discuss the identifiability of finite mixtures of elliptical densities with generators that correspond to scale mixtures of normal distributions.     

An application of linear feature selection to estimation of proportions

Let X be a random n-vector whose density function is given by a mixture of known multivariate normal density functions where the corresponding mixture proportions (a priori probabilities) are unknown. We present a numerically tractable method for obtaining estimates of the mixture proportions based on the linear feature selection technique of Guseman, Peters and Walker (1975).     

Maximum-likelihood estimates and likelihood-ratio criteria for multivariate elliptically contoured distributions

For a class of multivariate elliptically contoured distributions the maximum-likelihood estimators of the mean vector and covariance matrix are found under certain conditions. Likelihood-ratio criteria are obtained for a class of null hypotheses. These have the same form as in the normal case.     

A note on the inequalities for tail probability of the multivariate normal distribution

In this paper inequalities given by Harkness Godambe (1976) for the rail probabilities of the multivariate normal distribution in the equicorrelated case are improved by using the properties of the characteristic roots of a matrix and of the convex function.     

A class of weighted multivariate distributions related to doubly truncated multivariate t-distribution

This article introduces a class of weighted multivariate t-distributions, which includes the multivariate generalized Student t and multivariate skew t as its special members. This class is defined as the marginal distribution of a doubly truncated multivariate generalized Student t-distribution and studied from several aspects such as weighting of probability density functions, inequality constrained multivariate Student t-distributions, scale mixtures of multivariate normal and probabilistic representations. The relationships among these aspects are given, and various properties of the class are also discussed. Necessary theories and two applications are provided.     

A continuous general multivariate distribution and its properties

Starting from two known continuous univariate distributions, a bivariate distribution is constructed depending on a parameter which measures the degree of stochastic dependence between the two random variables. From the foregoing construction we then pass to a multivariate-type distribution, constructed using only univariate distributions and an association matrix. Some properties of the multivariate and bivariate case are studied.     

Cross-validation Bandwidth Matrices for Multivariate Kernel Density Estimation

Abstract.  The performance of multivariate kernel density estimates depends crucially on the choice of bandwidth matrix, but progress towards developing good bandwidth matrix selectors has been relatively slow. In particular, previous studies of cross-validation (CV) methods have been restricted to biased and unbiased CV selection of diagonal bandwidth matrices. However, for certain types of target density the use of full (i.e. unconstrained) bandwidth matrices offers the potential for significantly improved density estimation. In this paper, we generalize earlier work from diagonal to full bandwidth matrices, and develop a smooth cross-validation (SCV) methodology for multivariate data. We consider optimization of the SCV technique with respect to a pilot bandwidth matrix. All the CV methods are studied using asymptotic analysis, simulation experiments and real data analysis. The results suggest that SCV for full bandwidth matrices is the most reliable of the CV methods. We also observe that experience from the univariate setting can sometimes be a misleading guide for understanding bandwidth selection in the multivariate case.     

A predictive approach to the detection of additional information in a multivariate regression model

A predictive approach for the detection of additional information in a multivariate linear regression model is considered for the case of known and unknown error covariance matrices. The predictive density of future Observations on the additional variables under the model that they carry no information has been compared with the predictive density under the model that they do carry information. The Kullback-Leibler measure of divergence is used as a measure of comparison between the models.     

The uncertainty of forecasting: models with structural change versus those without changing parameters

The Bayesian predictive density is found for future observations of the unknown dependent variables for a multivariate linear model with a single shift in the regression matrix.

A numerical example shows that it is dangerous to predict future observations with an unchanging parameter model when the appropriate model should include structural change.     

Nonparametric Bayesian inference for multivariate density functions using Feller priors

Multivariate density estimation plays an important role in investigating the mechanism of high-dimensional data. This article describes a nonparametric Bayesian approach to the estimation of multivariate densities. A general procedure is proposed for constructing Feller priors for multivariate densities and their theoretical properties as nonparametric priors are established. A blocked Gibbs sampling algorithm is devised to sample from the posterior of the multivariate density. A simulation study is conducted to evaluate the performance of the procedure.     

TESTING SPHERICITY UNDER A MIXTURE MODEL

The robustness of Mauchly's sphericity test criterion when sampling from a mixture of two multivariate normal distributions is studied. The distribution of the sphericity test criterion when the sample covariance matrix has a non-central Wishart density of rank one is derived in terms of Meijer's G-functions; its distribution under the mixture model is then deduced. The robustness is studied by computing actual significance levels of the test under the mixture model using the critical values under the usual normal model.     

Methods for estimating proportions of convex combinations of normals using linear feature selection

Let X be a random n-vector whose density function is given by a mixtur.e of two density functions, h₁ and h₂ with unknown mixture proportions, Y₁ and Y₂ We assume that each of h₁ and h₂ is a convex combination of known multivariate normal density functions whose corresponding mixture proportions are also unknown. We present three numerically tractable methods for estimating Y₁ and Y₂ related to the technique of Guseman and Walton (1977), and based on the linear feature selection technique of Guseman, Peters and Walker (1975).