On multivariate weighted distributions

The concept of weighted distributions is well-known in the literature concerning observational studies and surveys in research related to forestry, ecology, bio-medicine and many other areas (cf. Rao (1965)). This paper extends the idea of weighted distributions to multivariate case. A few multivariate orderings have been defined and some partial ordering results are presented. Some results regarding multivariate positive and negative dependence are also discussed. Multivariate weighted distributions - joint, marginal and conditional have been defined and some important results concerning them are presented along with an illustration. The Multivariate Poisson Negative Hypergeometric Distribution has been derived.     

Concomitants of multivariate order statistics from multivariate elliptical distributions

ABSTRACT

In this article, we consider a (k + 1)n-dimensional elliptically contoured random vector (X^T₁, X₂^T, …, X^T_k, Z^T)^T = (X₁₁, …, X_1n, …, X_k1, …, X_kn, Z₁, …, Z_n)^T and derive the distribution of concomitant of multivariate order statistics arising from X₁, X₂, …, X_k. Specially, we derive a mixture representation for concomitant of bivariate order statistics. The joint distribution of the concomitant of bivariate order statistics is also obtained. Finally, the usefulness of our result is illustrated by a real-life data.     

A class of maximum-entropy multivariate distributions

This paper characterizes a class of multivariate distributions that includes the multinormal and is contained in the exponential family. The wide range of possible applications of these distributions is suggested by some of hte characteristics germane to them: First, they maximize Shannon's entropy among all distributions that have finite moments of given orders. As such, they constitute a class of distributions that includes the multinormal and some likely alternatives. Second, they can exhibit several modes, and, further-more, they do so with a relatively small number of parameters (compared to mixtures of multinormals). Third, they are the stationary distributions of certain diffusion processes. Fourth, they approximate, near the multinormal, the multivariate Pearson family. And fifth, the maximum likelihood estimators of their population moments are the sample moments. Two possible methods of estimating the distributions are studied in this paper: maximum likelihood estimation, and a fast procedure that can be used to find consistent estimators of the parameters via sample moments. A FORTTAN subroutine that implements the latter method is also provided.     

Testing for symmetry in multivariate distributions

Using the empirical characteristic function, a Cramér–von Mises test for reflected symmetry about an unspecified point is derived for multivariate distributions. The test statistic is based on an empirical process for which the weak convergence is established. The null properties of the test are studied as well as its power and local power. Estimators for the unknown symmetric point are previously proposed. Their consistency and asymptotical normality are proved by studying the weak convergence of some multidimensional empirical process. A simulation experiment shows that the estimators of the symmetric point are good, and that the test performs well on the examples tested. The new test is compared to the one derived in [N. Henze, B. Klar, S.G. Meintanis, Invariant tests for symmetry about an unspecified point based on empirical characteristic function, J. Multivariate. Anal. 87 (2003) 275–297].     

Matrix variance inequalities for multivariate distributions

We derive generalizations of the Hoeffding identity for multivariate random vectors and study some measures of dependence in the multidimensional case. In addition, we derive bounds on the covariance matrix for some multivariate distributions in the sense of Loewner ordering.     

A strategy for constructing multivariate distributions

Given a random vector (X₁,…, X_n) for which the univariate and bivariate marginal distributions belong to some specified families of distributions, we present a procedure for constructing families of multivariate distributions with the specified univariate and bivariate margins. Some general properties of the resulting families of multivariate distributions are reviewed. This procedure is illustrated by generalizing the bivariate Plackett (1965) and Clayton (1978) distributions to three dimensions. In addition to providing rich families of models for data analysis, this method of construction provides a convenient way of simulating observations from multivariate distributions with specific types of univariate and bivariate marginal distributions. A general algorithm for simulating random observations from these families of multivariate distributions is presented     

Characterizing parameters of multivariate elliptical distributions*

This paper defines new parameters characterizing multivariate elliptical distributions. Mardia's coefficient of multivariate kurtosis is shown to be essentially one of these parameters. A simple relation is established between centered multivariate product moments and second moments of the variables. The general results are verified on the contaminated normal distribution as an example.     

On the distributions of multivariate sample skewness

In this paper, we consider the multivariate normality test based on measure of multivariate sample skewness defined by Srivastava (1984). Srivastava derived asymptotic expectation up to the order N⁻¹ for the multivariate sample skewness and approximate χ²${\chi }^2$ test statistic, where N   is sample size. Under normality, we derive another expectation and variance for Srivastava's multivariate sample skewness in order to obtain a better test statistic. From this result, improved approximate χ²${\chi }^2$ test statistic using the multivariate sample skewness is also given for assessing multivariate normality. Finally, the numerical result by Monte Carlo simulation is shown in order to evaluate accuracy of the obtained expectation, variance and improved approximate χ²${\chi }^2$ test statistic. Furthermore, upper and lower percentiles of χ²${\chi }^2$ test statistic derived in this paper are compared with those of χ²${\chi }^2$ test statistic derived by Mardia (1974) which is used multivariate sample skewness defined by Mardia (1970).     

Approximating discrete multivariate distributions prom known moments

To approximate the joint distribution of the two-colony stepping-stone model, a finite mixture approach is proposed for constructing discrete multi-variate distributions. This approach generahzes the classic method of linear combinations of independent variables. The stepping-stone model is approximated through matching known moments. Numerical examples from entomology are given. Comparisons are made with the work by Wehrly et al (1993).     

An empirical goodness-of-fit test for multivariate distributions

An empirical test is presented as a tool for assessing whether a specified multivariate probability model is suitable to describe the underlying distribution of a set of observations. This test is based on the premise that, given any probability distribution, the Mahalanobis distances corresponding to data generated from that distribution will likewise follow a distinct distribution that can be estimated well by means of a large sample. We demonstrate the effectiveness of the test for detecting departures from several multivariate distributions. We then apply the test to a real multivariate data set to confirm that it is consistent with a multivariate beta model.     

Regression for doubly inflated multivariate Poisson distributions

Dependent multivariate count data occur in several research studies. These data can be modelled by a multivariate Poisson or Negative binomial distribution constructed using copulas. However, when some of the counts are inflated, that is, the number of observations in some cells are much larger than other cells, then the copula-based multivariate Poisson (or Negative binomial) distribution may not fit well and it is not an appropriate statistical model for the data. There is a need to modify or adjust the multivariate distribution to account for the inflated frequencies. In this article, we consider the situation where the frequencies of two cells are higher compared to the other cells and develop a doubly inflated multivariate Poisson distribution function using multivariate Gaussian copula. We also discuss procedures for regression on covariates for the doubly inflated multivariate count data. For illustrating the proposed methodologies, we present real data containing bivariate count observations with inflations in two cells. Several models and linear predictors with log link functions are considered, and we discuss maximum likelihood estimation to estimate unknown parameters of the models.     

Generation of multivariate distributions by vertical density representation

Troutt (1991,1993) proposed the idea of the vertical density representation (VDR) based on Box-Millar method. Kotz, Fang and Liang (1997) provided a systematic study on the multivariate vertical density representation (MVDR). Suppose that we want to generate a random vector X[d]Rⁿthat has a density function ?(x). The key point of using the MVDR is to generate the uniform distribution on [D]_?(v) = {x :?(x) = v} for any v > 0 which is the surface in RⁿIn this paper we use the conditional distribution method to generate the uniform distribution on a domain or on some surface and based on it we proposed an alternative version of the MVDR(type 2 MVDR), by which one can transfer the problem of generating a random vector X with given density f to one of generating (X, X_n+i) that follows the uniform distribution on a region in Rⁿ⁺¹defined by ?. Several examples indicate that the proposed method is quite practical.     

Parametric estimation under a class of multivariate distributions

In this paper the parameters of some members of a class of multivariate distributions, which was constructed by AL-Hussaini and Ateya (2003), are estimated by using the maximum likelihood and Bayes methods.     

Conditional specifications of multivariate pareto and student distributions

A random vector has a multivariate Pareto distribution if one of its univariate conditional distribution is Pareto and some of its marginals are identically distributed.A general method developed in the course of the proof of this result is applied also to characterize the multivariate Student (Cauchy) measure by one univariate Student conditional distribution.     

The new unified representation of multivariate skewed distributions

There are so many proposals in construction skewed distributions, and it is worth finding an overall class which covers all of these proposals. We introduce a new unified representation of multivariate skewed distributions. We will show that this new unified multivariate form of skewed distributions includes all of the continuous multivariate skewed distributions in the literature. This new unified representation is based on the multivariate probability integral transformation and can be decomposed into one factor that is original multivariate symmetric probability density function (pdf) f on ? ^k and skewed factor defined by a pdf p on [0, 1] ^k . This decomposition leads us to prove some useful properties of this new unified form. Stochastic representations and basic properties of this new form are also investigated in this article. Our work is motivated by considering the different skewing mechanisms which lead to different skewed distributions and show that all of these common-used distributions can be viewed as a new unified form.     

Measures of centrality for multivariate and directional distributions

Measures of centrality that generalize the univariate median are studied and applied to multivariate and directional distributions. A standard example is developed for general multivariate settings, and the uniqueness of the median proved for distributions satisfying certain regularity conditions. In the presence of weaker regularity, this median is shown to be of codimension 2. Conditions are also provided for these measures of centrality to be equivariant under transformations on the sample space. The equivariance of the usual univariate median under monotone transformations is seen as a special case.     

Finite mixtures of multivariate Poisson distributions with application

In the present paper we examine finite mixtures of multivariate Poisson distributions as an alternative class of models for multivariate count data. The proposed models allow for both overdispersion in the marginal distributions and negative correlation, while they are computationally tractable using standard ideas from finite mixture modelling. An EM type algorithm for maximum likelihood (ML) estimation of the parameters is developed. The identifiability of this class of mixtures is proved. Properties of ML estimators are derived. A real data application concerning model based clustering for multivariate count data related to different types of crime is presented to illustrate the practical potential of the proposed class of models.     

Robust statistics for testing mean vectors of multivariate distributions

We develop a ‘robust’ statistic T² _R, based on Tiku's (1967, 1980) MML (modified maximum likelihood) estimators of location and scale parameters, for testing an assumed meam vector of a symmetric multivariate distribution. We show that T² _R is one the whole considerably more powerful than the prominenet Hotelling T² statistics. We also develop a robust statistic T² _D for testing that two multivariate distributions (skew or symmetric) are identical; T² _D seems to be usually more powerful than nonparametric statistics. The only assumption we make is that the marginal distributions are of the type (1/σ_k)f((x-μ_k)/σ_k) and the means and variances of these marginal distributions exist.     

Some multivariate singular unified skew-normal distributions and their application

Abstract

We introduce here the truncated version of the unified skew-normal (SUN) distributions. By considering a special truncations for both univariate and multivariate cases, we derive the joint distribution of consecutive order statistics X_{(r, ..., r + k)} = (X_(r), ..., X_{(r + K)})^T from an exchangeable n-dimensional normal random vector X. Further we show that the conditional distributions of X_{(r + j, ..., r + k)} given X_{(r, ..., r + j ? 1)}, X_{(r, ..., r + k)} given (X_(r) > t)?and X_{(r, ..., r + k)} given (X_{(r + k)} < t) are special types of singular SUN distributions. We use these results to determine some measures in the reliability theory such as the mean past life (MPL) function and mean residual life (MRL) function.