On a class of bivariate mixed Sarmanov distributions

Multivariate distributions are more and more used to model the dependence encountered in many fields. However, classical multivariate distributions can be restrictive by their nature, while Sarmanov's multivariate distribution, by joining different marginals in a flexible and tractable dependence structure, often provides a valuable alternative. In this paper, we introduce some bivariate mixed Sarmanov distributions with the purpose to extend the class of bivariate Sarmanov distributions and to obtain new dependency structures. Special attention is paid to the bivariate mixed Sarmanov distribution with Poisson marginals and, in particular, to the resulting bivariate Sarmanov distributions with negative binomial and with Poisson‐inverse Gaussian marginals; these particular types of mixed distributions have possible applications in, for example modelling bivariate count data. The extension to higher dimensions is also discussed. Moreover, concerning the dependency structure, we also present some correlation formulas.     

Multivariate Dispersion Models Generated From Gaussian Copula

In this paper a class of multivariate dispersion models generated from the multivariate Gaussian copula is presented. Being a multivariate extension of Jørgensen's (1987a) dispersion models, this class of multivariate models is parametrized by marginal position, dispersion and dependence parameters, producing a large variety of multivariate discrete and continuous models including the multivariate normal as a special case. Properties of the multivariate distributions are investigated, some of which are similar to those of the multivariate normal distribution, which makes these models potentially useful for the analysis of correlated non-normal data in a way analogous to that of multivariate normal data. As an example, we illustrate an application of the models to the regression analysis of longitudinal data, and establish an asymptotic relationship between the likelihood equation and the generalized estimating equation of Liang & Zeger (1986).     

On principal points for location mixtures of spherically symmetric distributions

In this paper, we investigate some properties of 2-principal points for location mixtures of spherically symmetric distributions with focus on a linear subspace in which a set of 2-principal points must lie. Our results can be viewed as an extension of those of Yamamoto and Shinozaki [2000. Two principal points for multivariate location mixtures of spherically symmetric distributions. J. Japan Statist. Soc. 30, 53–63], where a finite location mixture of spherically symmetric distributions is treated. As an extension of their paper, this paper defines a wider class of distributions, and derives a linear subspace in which a set of 2-principal points must exist. A theorem useful for comparing the mean squared distances is also established.     

A note on the Cauchy-type mixture distributions

An interesting class of continuous distributions, called Cauchy-type mixture, with potential applications in modelling erratic phenomena is introduced by Soltani and Tafakori [A class of continuous kernels and Cauchy type heavy tail distributions. Statist Probab Lett. 2013;83:1018–1027]. In this work, we provide more insights into the Cauchy-type mixture distributions, involving certain characterizations, connections with the generalized Linnik distributions and the class of discrete distributions induced by stable laws. We also prove that the Laplace transform of Cauchy-type mixture distributions when normalized by constant terms become as a density functions in terms of distributional conjugate property.     

A family of repeated measurements models

A general family of multivariate distributions for repeated measures can be obtained by applying the Laplace transform of a gamma distribution to the integrated intensity function of any continuous distribution on the positive real line. Both clustering and serial dependence can be handled. The response variable may be counts, durations between events, or any continuous positive-valued measurements.     

Efficient time/space algorithm to compute rectangular probabilities of multinomial, multivariate hypergeometric and multivariate Pólya distributions

The computation of rectangular probabilities of multivariate discrete integer distributions such as the multinomial, multivariate hypergeometric or multivariate Pólya distributions is of great interest both for statistical applications and for probabilistic modeling purpose. All these distributions are members of a broader family of multivariate discrete integer distributions for which computationaly efficient approximate methods have been proposed for the evaluation of such probabilities, but with no control over their accuracy. Recently, exact algorithms have been proposed for computing such probabilities, but they are either dedicated to a specific distribution or to very specific rectangular probabilities. We propose a new algorithm that allows to perform the computation of arbitrary rectangular probabilities in the most general case. Its accuracy matches or even outperforms the accuracy exact algorithms when the rounding errors are taken into account. In the worst case, its computational cost is the same as the most efficient exact method published so far, and is much lower in many situations of interest. It does not need any additional storage than the one for the parameters of the distribution, which allows to deal with large dimension/large counting parameter applications at no extra memory cost and with an acceptable computation time, which is a major difference with respect to the methods published so far.     

Some bivariate uniform distributions

Bivariate uniform distributions with dependent components are readily derived by distribution function transformations of the components of non-uniform dependent continuous bivariate random variables (X,Y). Contour plots of joint density functions show the various, and varying, forms of dependence which can arise from different distributional forms for (X,Y) and aids the choice of bivariate uniform distributions as empirical models.     

A new class of multivariate skew distributions with applications to bayesian regression models

Abstract: The authors develop a new class of distributions by introducing skewness in multivariate elliptically symmetric distributions. The class, which is obtained by using transformation and conditioning, contains many standard families including the multivariate skew‐normal and t distributions. The authors obtain analytical forms of the densities and study distributional properties. They give practical applications in Bayesian regression models and results on the existence of the posterior distributions and moments under improper priors for the regression coefficients. They illustrate their methods using practical examples.     

Parameterizations and Fitting of Bi-directed Graph Models to Categorical Data

Abstract.  We discuss two parameterizations of models for marginal independencies for discrete distributions which are representable by bi-directed graph models, under the global Markov property. Such models are useful data analytic tools especially if used in combination with other graphical models. The first parameterization, in the saturated case, is also known as thenation multivariate logistic transformation, the second is a variant that allows, in some (but not all) cases, variation-independent parameters. An algorithm for maximum likelihood fitting is proposed, based on an extension of the Aitchison and Silvey method.     

Duality of the Multivariate Distributions of Marshall–Olkin Type and Tail Dependence

A dual class of the multivariate distributions of Marshall–Olkin type is introduced, and their copulas are presented and utilized to derive explicit expressions of the distributional tail dependencies, which describe the amount of dependence in the upper-orthant tail or lower-orthant tail of a multivariate distribution and can be used in the study of dependence among extreme values. A sufficient condition under which tail dependencies of two such distributions can be compared are obtained. Some examples are also presented to illustrate our results.     

Empirical likelihood ratio test for or against a set of inequality constraints

We use the empirical likelihood ratio approach introduced by Owen (Biometrika 75 (1988), 237–249) to test for or against a set of inequality constraints when the parameters are defined by estimating functions. Our objective in this paper is to show that under fairly general conditions, the limiting distributions of the empirical likelihood ratio test statistics are of chi-bar square type (as in the parametric case) and give the expression of the weighting values. The results obtained here are similar to those in El Barmi and Dykstra (1995) where a full distributional model is assumed. This work presents also an extension of the results in Qin and Lawless (1995).     

Simulation of some multivariate distributions related to the Dirichlet distribution with application to Monte Carlo simulations

Simulation of multivariate distributions is important in many applications but remains computationally challenging in practice. In this article, we introduce three classes of multivariate distributions from which simulation can be conducted by means of their stochastic representations related to the Dirichlet random vector. More emphasis is made to simulation from the class of uniform distributions over a polyhedron, which is useful for solving some constrained optimization problems and ha`s many applications in sampling and Monte Carlo simulations. Numerical evidences show that, by utilizing state-of-the-art Dirichlet generation algorithms, the introduced methods become superior to other approaches in terms of computational efficiency.     

Positive dependence of random variables with a common marginal distribution

In this paper we review some notions of positive dependence of random variables with a common univariate marginal distribution and describe the related moment and probability inequalities. We first present a comparison between i.i.d. random variables and exchangeable random variables via an application of de Finetti's theorem, then describe some useful probability inequalities via partial orderings of the strength of their positive dependence. Finally, we state a result for random variables which are not necessarily exchangeable. Special applications to the multivariate normal distribution will be discussed, and the results involve only the correlation matrix of the distribution.     

A general algorithm fob,univariate and multivariate goodness of pit tests based on graphical representation

This paper presents a general algorithm tor assessing the distributional assumptions. Empirical distributions of the corresponding test statistics are obtained and examples are given to illustrate various applications of the proposed test. By using the squared radii and angles, it is shown that the problem of assessing multivariate normality can be reduced to that of testing for a univariate distribution. A limited comparison is made to investigate the power of the proposed test. This work was supported in part by the National Science Foundation under Grant NO.G88135. Support from the Computer Applications ami Software Engineering (CASE) Center of Syracuse University is also gratefully acknowledged     

Two multivariate generalized beta families

A multivariate generalized beta distribution is introduced that extends the univariate generalized beta distribution and includes many multivariate distributions, such as the multivariate beta of the first and second kind, the generalized gamma, and the Burr and Dirichlet distributions as special and limiting cases. These interrelationships can be illustrated using a distributional family tree. The corresponding marginal distributions are univariate generalized beta distributions and their special cases. Selected expressions for the moments are reported, and an application to the joint distribution of income and wealth is presented. A simple transformation of the multivariate generalized beta distribution leads to what will be referred to as a multivariate exponential generalized beta distribution, which includes a multivariate form of the logistics and Burr distributions as special cases.     

A Multivariate Extension of Inverse Gaussian Distribution Derived from Inverse Relationship

Abstract

We propose a new multivariate extension of the inverse Gaussian distribution derived from a certain multivariate inverse relationship. First we define a multivariate extension of the inverse relationship between two sets of multivariate distributions, then define a reduced inverse relationship between two multivariate distributions. We derive the multivariate continuous distribution that has the reduced multivariate inverse relationship with a multivariate normal distribution and call it a multivariate inverse Gaussian distribution. This distribution is also characterized as the distribution of the location of a multivariate Brownian motion at some stopping time. The marginal distribution in one direction is the inverse Gaussian distribution, and the conditional distribution in the space perpendicular to this direction is a multivariate normal distribution. Mean, variance, and higher order cumulants are derived from the multivariate inverse relationship with a multivariate normal distribution. Other properties such as reproductivity and infinite divisibility are also given.     

Simultaneous inference in structured additive conditional copula regression models: a unifying Bayesian approach

While most regression models focus on explaining distributional aspects of one single response variable alone, interest in modern statistical applications has recently shifted towards simultaneously studying multiple response variables as well as their dependence structure. A particularly useful tool for pursuing such an analysis are copula-based regression models since they enable the separation of the marginal response distributions and the dependence structure summarised in a specific copula model. However, so far copula-based regression models have mostly been relying on two-step approaches where the marginal distributions are determined first whereas the copula structure is studied in a second step after plugging in the estimated marginal distributions. Moreover, the parameters of the copula are mostly treated as a constant not related to covariates and most regression specifications for the marginals are restricted to purely linear predictors. We therefore propose simultaneous Bayesian inference for both the marginal distributions and the copula using computationally efficient Markov chain Monte Carlo simulation techniques. In addition, we replace the commonly used linear predictor by a generic structured additive predictor comprising for example nonlinear effects of continuous covariates, spatial effects or random effects and furthermore allow to make the copula parameters covariate-dependent. To facilitate Bayesian inference, we construct proposal densities for a Metropolis–Hastings algorithm relying on quadratic approximations to the full conditionals of regression coefficients avoiding manual tuning. The performance of the resulting Bayesian estimates is evaluated in simulations comparing our approach with penalised likelihood inference, studying the choice of a specific copula model based on the deviance information criterion, and comparing a simultaneous approach with a two-step procedure. Furthermore, the flexibility of Bayesian conditional copula regression models is illustrated in two applications on childhood undernutrition and macroecology.     

Empirical Study of Six Tests for Equality of Populations with Zero-Inflated Continuous Distributions

We evaluated the properties of six statistical methods for testing equality among populations with zero-inflated continuous distributions. These tests are based on likelihood ratio (LR), Wald, central limit theorem (CLT), modified CLT (MCLT), parametric jackknife (PJ), and nonparametric jackknife (NPJ) statistics. We investigated their statistical properties using simulated data from mixed distributions with an unknown portion of non zero observations that have an underlying gamma, exponential, or log-normal density function and the remaining portion that are excessive zeros. The 6 statistical tests are compared in terms of their empirical Type I errors and powers estimated through 10,000 repeated simulated samples for carefully selected configurations of parameters. The LR, Wald, and PJ tests are preferred tests since their empirical Type I errors were close to the preset nominal 0.05 level and each demonstrated good power for rejecting null hypotheses when the sample sizes are at least 125 in each group. The NPJ test had unacceptable empirical Type I errors because it rejected far too often while the CLT and MCLT tests had low testing powers in some cases. Therefore, these three tests are not recommended for general use but the LR, Wald, and PJ tests all performed well in large sample applications.     

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.