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 Birnbaum–Saunders distribution: properties and associated inference

The univariate fatigue life distribution proposed by Birnbaum and Saunders [A new family of life distributions. J Appl Probab. 1969;6:319–327] has been used quite effectively to model times to failure for materials subject to fatigue and for modelling lifetime data and reliability problems. In this article, we introduce a Birnbaum–Saunders (BS) distribution in the multivariate setting. The new multivariate model arises in the context of conditionally specified distributions. The proposed multivariate model is an absolutely continuous distribution whose marginals are univariate BS distributions. General properties of the multivariate BS distribution are derived and the estimation of the unknown parameters by maximum likelihood is discussed. Further, the Fisher's information matrix is determined. Applications to real data of the proposed multivariate distribution are provided for illustrative purposes.     

Copules archimédiennes et families de lois bidimensionnelles dont les marges sont données

Every bivariate distribution function with continuous marginals can be represented in terms of a unique copula, that is, in terms of a distribution function on the unit square with uniform marginals. This paper is concerned with a special class of copulas called Archimedean, which includes the uniform representation of many standard bivariate distributions. Conditions are given under which these copulas are stochastically ordered and pointwise limits of sequences of Archimedean copulas are examined. We also provide two new one-parameter families of bivariate distributions which include as limiting cases the Frechet bounds and the independence distribution.     

A family of bivariate distributions with non-elliptical contours

In this paper, a family of copulas with two parameters is proposed and its dependence analysis is performed. The corresponding family of bivariate distributions with specified marginals is constructed. For normal marginals, the new distributions are non-elliptical and can be applied in data analysis. They provide various alternative hypotheses for testing normality. Finally, an example is given.     

Multivariate tail conditional expectation for scale mixtures of skew-normal distribution

In practice, a financial or actuarial data set may be a skewed or heavy-tailed and this motivates us to study a class of distribution functions in risk management theory that provide more information about these characteristics resulting in a more accurate risk analysis. In this paper, we consider a multivariate tail conditional expectation (MTCE) for multivariate scale mixtures of skew-normal (SMSN) distributions. This class of distributions contains skewed distributions and some members of this class can be used to analyse heavy-tailed data sets. We also provide a closed form for TCE in a univariate skew-normal distribution framework. Numerical examples are also provided for illustration.     

Multivariate distributions with generalized inverse gaussian marginals,and associated poisson mixtures

Several types of multivariate extensions of the inverse Gaussian (IG) distribution and the reciprocal inverse Gaussian (RIG) distribution are proposed. Some of these types are obtained as random-additive-effect models by means of well-known convolution properties of the IG and RIG distributions, and they have one-dimensional IG or RIG marginals. They are used to define a flexible class of multivariate Poisson mixtures.     

On the matrix-variate generalized hyperbolic distribution and its Bayesian applications

In the first part of the paper, we introduce the matrix-variate generalized hyperbolic distribution by mixing the matrix normal distribution with the matrix generalized inverse Gaussian density. The p-dimensional generalized hyperbolic distribution of [Barndorff-Nielsen, O. (1978). Hyperbolic distributions and distributions on hyperbolae. Scand. J. Stat., 5, 151–157], the matrix-T distribution and many well-known distributions are shown to be special cases of the new distribution. Some properties of the distribution are also studied. The second part of the paper deals with the application of the distribution in the Bayesian analysis of the normal multivariate linear model.     

Bivariate Kumaraswamy distribution: properties and a new method to generate bivariate classes

In this paper, we introduce a bivariate Kumaraswamy (BVK) distribution whose marginals are Kumaraswamy distributions. The cumulative distribution function of this bivariate model has absolutely continuous and singular parts. Representations for the cumulative and density functions are presented and properties such as marginal and conditional distributions, product moments and conditional moments are obtained. We show that the BVK model can be obtained from the Marshall and Olkin survival copula and obtain a tail dependence measure. The estimation of the parameters by maximum likelihood is discussed and the Fisher information matrix is determined. We propose an EM algorithm to estimate the parameters. Some simulations are presented to verify the performance of the direct maximum-likelihood estimation and the proposed EM algorithm. We also present a method to generate bivariate distributions from our proposed BVK distribution. Furthermore, we introduce a BVK distribution which has only an absolutely continuous part and discuss some of its properties. Finally, a real data set is analysed for illustrative purposes.     

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.     

Complex elliptically symmetric distributions

In this paper, the authors introduce a class of distributions known as complex elliptically symmetric distributions. The complex multivariate normal and complex multivariate t distributions are members of this class. Various properties of the complex elliptically symmetric distributions are studied. Finally, the robustness of certain test procedures are discussed when the assumption of complex multivariate normality is violated but the underlying distribution still belongs to the class of elliptically symmetric distributions.     

A Bivariate Distribution Family with Specified Correlation Coefficient and Given Marginals

For given continuous distribution functions F(x) and G(y) and a Pearson correlation coefficient ρ, an algorithm is provided to construct a sequence of continuous bivariate distributions with marginals equal to F(x) and G(y) and the corresponding correlation coefficient converges to ρ. The algorithm can be easily implemented using S-Plus or R. Applications are given to generate bivariate random variables with marginals including Gamma, Beta, Weibull, and uniform distributions.     

Multivariate Poisson hidden Markov models with a case study of modelling seismicity

While the literature on multivariate models for continuous data flourishes, there is a lack of models for multivariate counts. We aim to contribute to this framework by extending the well known class of univariate hidden Markov models to the multidimensional case, by introducing multivariate Poisson hidden Markov models. Each state of the extended model is associated with a different multivariate discrete distribution. We consider different distributions with Poisson marginals, starting from the multivariate Poisson distribution and then extending to copula based distributions to allow flexible dependence structures. An EM type algorithm is developed for maximum likelihood estimation. A real data application is presented to illustrate the usefulness of the proposed models. In particular, we apply the models to the occurrence of strong earthquakes (surface wave magnitude ≥5), in three seismogenic subregions in the broad region of the North Aegean Sea for the time period from 1 January 1981 to 31 December 2008. Earthquakes occurring in one subregion may trigger events in adjacent ones and hence the observed time series of events are cross‐correlated. It is evident from the results that the three subregions interact with each other at times differing by up to a few months. This migration of seismic activity is captured by the model as a transition to a state of higher seismicity.     

Effect of frailty on marginal regression estimates in survival analysis

Unexplained heterogeneity in univariate survival data and association in multivariate survival can both be modelled by the inclusion of frailty effects. This paper investigates the consequences of ignoring frailty in analysis, fitting misspecified Cox proportional hazards models to the marginal distributions. Regression coefficients are biased towards 0 by an amount which depends in magnitude on the variability of the frailty terms and the form of frailty distribution. The bias is reduced when censoring is present. Fitted marginal survival curves can also differ substantially from the true marginals.     

Form-invariance under weighted sampling

In this paper, form-invariant weighted distributions are considered in an exponential family. The class of bivariate distribution with invariant property is identified under exponential weight function. The class includes some of the custom bivariate models. The form-invariant multivariate normal distributions are obtained under a quadratic weight function.     

Statistical applications of the multivariate skew normal distribution

Azzalini and Dalla Valle have recently discussed the multivariate skew normal distribution which extends the class of normal distributions by the addition of a shape parameter. The first part of the present paper examines further probabilistic properties of the distribution, with special emphasis on aspects of statistical relevance. Inferential and other statistical issues are discussed in the following part, with applications to some multivariate statistics problems, illustrated by numerical examples. Finally, a further extension is described which introduces a skewing factor of an elliptical density.     

Bivariate inverse Weibull distribution

ABSTRACT

Recently it is observed that the inverse Weibull (IW) distribution can be used quite effectively to analyse lifetime data in one dimension. The main aim of this paper is to define a bivariate inverse Weibull (BIW) distribution so that the marginals have IW distributions. It is observed that the joint probability density function and the joint cumulative distribution function can be expressed in compact forms. Several properties of this distribution such as marginals, conditional distributions and product moments have been discussed. We obtained the maximum likelihood estimates for the unknown parameters of this distribution and their approximate variance– covariance matrix. We perform some simulations to see the performances of the maximum likelihood estimators. One data set has been re-analysed and it is observed that the bivariate IW distribution provides a better fit than the bivariate exponential distribution.     

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.     

Weighted Marshall–Olkin bivariate exponential distribution

Recently, Gupta and Kundu [R.D. Gupta and D. Kundu, A new class of weighted exponential distributions, Statistics 43 (2009), pp. 621–634] have introduced a new class of weighted exponential (WE) distributions, and this can be used quite effectively to model lifetime data. In this paper, we introduce a new class of weighted Marshall–Olkin bivariate exponential distributions. This new singular distribution has univariate WE marginals. We study different properties of the proposed model. There are four parameters in this model and the maximum-likelihood estimators (MLEs) of the unknown parameters cannot be obtained in explicit forms. We need to solve a four-dimensional optimization problem to compute the MLEs. One data set has been analysed for illustrative purposes and finally we propose some generalization of the proposed model.     

Power-normal distribution

Recently, Gupta and Gupta [Analyzing skewed data by power-normal model, Test 17 (2008), pp. 197–210] proposed the power-normal distribution for which normal distribution is a special case. The power-normal distribution is a skewed distribution, whose support is the whole real line. Our main aim of this paper is to consider bivariate power-normal distribution, whose marginals are power-normal distributions. We obtain the proposed bivariate power-normal distribution from Clayton copula, and by making a suitable transformation in both the marginals. Lindley–Singpurwalla distribution also can be used to obtain the same distribution. Different properties of this new distribution have been investigated in detail. Two different estimators are proposed. One data analysis has been performed for illustrative purposes. Finally, we propose some generalizations to multivariate case also along the same line and discuss some of its properties.     

Algorithms

Abstract

The main reason for the limited use of multivariate discrete models is the difficulty in calculating the required probabilities. The task is usually undertaken via recursive relationships which become quite computationally demanding for high dimensions and large values. The present paper discusses efficient algorithms that make use of the recurrence relationships in a manner that reduces the computational effort and thus allow for easy and cheap calculation of the probabilities. The most common multivariate discrete distribution, the multivariate Poisson distribution is treated. Real data problems are provided to motivate the use of the proposed strategies. Extensions of our results are discussed. It is shown that probabilities, for a large family of multivariate distributions, can be computed efficiently via our algorithms.