Order statistics from trivariate normal and -distributions in terms of generalized skew-normal and skew- distributions

We consider here a generalization of the skew-normal distribution, GSN(λ₁,λ₂,ρ), defined through a standard bivariate normal distribution with correlation ρ, which is a special case of the unified multivariate skew-normal distribution studied recently by Arellano-Valle and Azzalini [2006. On the unification of families of skew-normal distributions. Scand. J. Statist. 33, 561–574]. We then present some simple and useful properties of this distribution and also derive its moment generating function in an explicit form. Next, we show that distributions of order statistics from the trivariate normal distribution are mixtures of these generalized skew-normal distributions; thence, using the established properties of the generalized skew-normal distribution, we derive the moment generating functions of order statistics, and also present expressions for means and variances of these order statistics.Next, we introduce a generalized skew-t_ν distribution, which is a special case of the unified multivariate skew-elliptical distribution presented by Arellano-Valle and Azzalini [2006. On the unification of families of skew-normal distributions. Scand. J. Statist. 33, 561–574] and is in fact a three-parameter generalization of Azzalini and Capitanio's [2003. Distributions generated by perturbation of symmetry with emphasis on a multivariate skew t distribution. J. Roy. Statist. Soc. Ser. B 65, 367–389] univariate skew-t_ν form. We then use the relationship between the generalized skew-normal and skew-t_ν distributions to discuss some properties of generalized skew-t_ν as well as distributions of order statistics from bivariate and trivariate t_ν distributions. We show that these distributions of order statistics are indeed mixtures of generalized skew-t_ν distributions, and then use this property to derive explicit expressions for means and variances of these order statistics.     

Nonlinear regression in a trivariate elliptical distribution via a linear combination of order statistics

In this paper, by assuming that (X, Y ₁, Y ₂)^T has a trivariate elliptical distribution, we derive the exact joint distribution of X and a linear combination of order statistics from (Y ₁, Y ₂)^T and show that it is a mixture of unified bivariate skew-elliptical distributions. We then derive the corresponding marginal and conditional distributions for the special case of t kernel. We also present these results for an exchangeable case with t kernel and illustrate the established results with an air-pollution data.     

An integral representation technique for calculating general multivariate probabilities with an application to multivariate x2

Recently in Dutt (1973, (1975), intgral representations over (0,A) were obtained for upper and lover multivariate normal and the probilities. It was pointed out that these integral representaitons when evaluated by Gauss-Hermite uadrature yield rapid and accurate numerical results.

Here integral representaitons, based on an integral formula due to Gurland (1948), are indicated for arbitrary multivariate probabilities. Application of this general representaion for computing multivariate x2 probabilities is discussed and numerical results using Gaussian quadrature are given for the bivariate and equicorre lated trivariate cases. Applications to the multivariate densities studied by Miller (1965) are also included     

The evaluation of trivariate normal probabilities defined by linear inequalities

This paper considers the evaluation of probabilities which are defined by a set of linear inequalities of a trivariate normal distribution. It is shown that these probabilities can be evaluated by a one-dimensional numerical integration. The trivariate normal distribution can have any covariance matrix and any mean vector, and the probability can be defined by any number of one-sided and two-sided linear inequalities. This affords a practical and efficient method for the calculation of these probabilities which is superior to basic simulation methods. An application of this method to the analysis of pairwise comparisons of four treatment effects is discussed.     

Robust mixture modeling using the skew <Emphasis Type="Italic">t</Emphasis> distribution

A finite mixture model using the Student's t distribution has been recognized as a robust extension of normal mixtures. Recently, a mixture of skew normal distributions has been found to be effective in the treatment of heterogeneous data involving asymmetric behaviors across subclasses. In this article, we propose a robust mixture framework based on the skew t distribution to efficiently deal with heavy-tailedness, extra skewness and multimodality in a wide range of settings. Statistical mixture modeling based on normal, Student's t and skew normal distributions can be viewed as special cases of the skew t mixture model. We present analytically simple EM-type algorithms for iteratively computing maximum likelihood estimates. The proposed methodology is illustrated by analyzing a real data example.     

Some useful integrals and their applications in correlation analysis

The distribution of the product moment correlation coefficient based on the bivariate normal distribution is well known. Recently in many business and economic data, fat tailed distributions especially some elliptical distributions have been considered as parent populations. The normal and t-distributions are well known special cases of elliptical distribution. In this paper we derive some theorems involving double integrals and apply them to derive the probability distribution of the correlation coefficient for some elliptical populations. The general nature of the theorems indicates their potential use in probability distribution theory.     

A self-validating numerical method for computation of central and non-central F probabilities and percentiles

A self-validating numerical method based on interval analysis for the computation of central and non-central F probabilities and percentiles is reported. The major advantage of this approach is that there are guaranteed error bounds associated with the computed values (or intervals), i.e. the computed values satisfy the user-specified accuracy requirements. The methodology reported in this paper can be adapted to approximate the probabilities and percentiles for other commonly used distribution functions.     

Bayesian analysis of mixture modelling using the multivariate t distribution

A finite mixture model using the multivariate t distribution has been shown as a robust extension of normal mixtures. In this paper, we present a Bayesian approach for inference about parameters of t-mixture models. The specifications of prior distributions are weakly informative to avoid causing nonintegrable posterior distributions. We present two efficient EM-type algorithms for computing the joint posterior mode with the observed data and an incomplete future vector as the sample. Markov chain Monte Carlo sampling schemes are also developed to obtain the target posterior distribution of parameters. The advantages of Bayesian approach over the maximum likelihood method are demonstrated via a set of real data.     

Approximating probability integrals of multivariate normal using association models

The bivariate plane is symmetrically partitioned into fine rectangular regions, and a symmetric uniform association model is used to represent the resulting discretized bivariate normal probabilities. A new algorithm is developed by utilizing a quadrature and the above association model to approximate the diagonal probabilities. The off-diagonal probabilities are then approximated using the model. This method is an alternative to Wang's (1987) approach, computationally advantageous and relatively easy to extend to higher dimensions. Bivariate and trivariate normal probabilities approximated by our method are observed to agree very closely with the corresponding known results.     

The product <Emphasis Type="Italic">t</Emphasis> density distribution arising from the product of two Student’s <Emphasis Type="Italic">t</Emphasis> PDFs

The Student’s t distribution has become increasingly prominent and is considered as a competitor to the normal distribution. Motivated by real examples in Physics, decision sciences and Bayesian statistics, a new t distribution is introduced by taking the product of two Student’s t pdfs. Various structural properties of this distribution are derived, including its cdf, moments, mean deviation about the mean, mean deviation about the median, entropy, asymptotic distribution of the extreme order statistics, maximum likelihood estimates and the Fisher information matrix. Finally, an application to a Bayesian testing problem is illustrated.     

Regression mean residual life of a parallel system of dependent components

In this paper, we consider a system consisting of two dependent components and we are interested in the average remaining life of the component that fails last when (i) the first failure occurs at time t and (ii) the first failure occurs after time t. For both the cases, expressions are derived in the case of general bivariate normal distribution and a class of bivariate exponential distribution including bivariate exponential distribution of Arnold and Strauss, absolutely continuous bivariate exponential distribution of Block and Basu, bivariate exponential distribution of Raftery, Freund's bivariate exponential distribution and Gumbel's bivariate exponential distribution.     

Robust mixture modelling using the t distribution

Normal mixture models are being increasingly used to model the distributions of a wide variety of random phenomena and to cluster sets of continuous multivariate data. However, for a set of data containing a group or groups of observations with longer than normal tails or atypical observations, the use of normal components may unduly affect the fit of the mixture model. In this paper, we consider a more robust approach by modelling the data by a mixture of t distributions. The use of the ECM algorithm to fit this t mixture model is described and examples of its use are given in the context of clustering multivariate data in the presence of atypical observations in the form of background noise.     

On some study of skew-t distributions

Abstract

In this note, through ratio of independent random variables, new families of univariate and bivariate skew-t distributions are introduced. Probability density function for each skew-t distribution will be given. We also derive explicit forms of moments of the univariate skew-t distribution and recurrence relations for its cumulative distribution function. Finally we illustrate the flexibility of this class of distributions with applications to a simulated data and the volcanos heights data.     

Self-validated Computations for the Probabilities of the Central Bivariate Chi-square Distribution and a Bivariate F Distribution. This work partially supported by National Science Foundation grant DMS-9500831

Self-validated computations using interval analysis produce results with a guaranteed error bound. This article presents methods for self-validated computation of probabilities and percentile points of the bivariate chi-square distribution and a bivariate F distribution. For the computation of critical points (c ₁,c ₂) in the equation P(Y ₁ @ c ₁, Y ₂ ≤ c ₂) = 1 ? α, the case c ₁ = c ₂ is considered. A combination of interval secant and bisection algorithms is developed for finding enclosures of the percentile points of the distribution. Results are compared to previously published tables.     

On the distribution of quotient of random variables conditioned to the positive quadrant

Abstract

Motivated by Caginalp and Caginalp [Physica A—Statistical Mechanics and Its Applications, 499, 2018, 457–471], we derive the exact distribution of X/Y conditioned on X?>?0, Y?>?0 for more than ten classes of distributions, including the bivariate t, bivariate Cauchy, bivariate Lomax, Arnold and Strauss’ bivariate exponential, Balakrishna and Shiji’s bivariate exponential, Mohsin et al.’s bivariate exponential, Morgenstern type bivariate exponential, bivariate gamma exponential and bivariate alpha skew normal distributions. The results can be useful in finance and other areas.     

A system for computing joint probabilities from jensen's bivariate f distribution

A system of subroutines is presented for efficient computation of joint probabilities from Jensen's bivariate F distribution. Any valid set of parameters is permitted, whereas previous work was limited to the special case of equal numerator degrees of freedom and equal canonical correlations in the underlying multinormal distribution. The use of joint Probabilities from Jensen's bivariate F distribution is demonstrated via an application to two-way ANOVA without interaction.     

Tail dependence of skew t-copulas

We examine tail behavior of skew t-copula in the bivariate case. The tail dependence coefficient is calculated for different skewing parameter values and compared with the corresponding coefficient for the t-copula. It is shown that depending on skewing parameter values, the tail dependence coefficient can differ considerably from the tail dependence of the t-copula. The speed of convergence of the estimator of tail dependence coefficient to its theoretical value is examined in a simulation experiment. Method of moments and maximum likelihood method are compared by simulation either. In the considered cases, maximum likelihood method converged faster to the theoretical value.     

On a Bivariate Pólya-Aeppli Distribution

In this article, we use the bivariate Poisson distribution obtained by the trivariate reduction method and compound it with a geometric distribution to derive a bivariate Pólya-Aeppli distribution. We then discuss a number of properties of this distribution including the probability generating function, correlation structure, probability mass function, recursive relations, and conditional distributions. The generating function of the tail probabilities is also obtained. Moment estimation of the parameters is then discussed and illustrated with a numerical example.