A note on the orthant probability of the elliptically contoured distribution

When a scale matrix satisfies certain conditions, the orthant probability of the elliptically contoured distribution with the scale matrix is expressed as the same probability of the equicorrelated normal distribution.     

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.     

The evaluation of general non-centred orthant probabilities

Summary. The evaluation of the cumulative distribution function of a multivariate normal distribution is considered. The multivariate normal distribution can have any positive definite correlation matrix and any mean vector. The approach taken has two stages. In the first stage, it is shown how non-centred orthoscheme probabilities can be evaluated by using a recursive integration method. In the second stage, some ideas of Schläfli and Abrahamson are extended to show that any non-centred orthant probability can be expressed as differences between at most ( m −1)! non-centred orthoscheme probabilities. This approach allows an accurate evaluation of many multivariate normal probabilities which have important applications in statistical practice.     

Computation of the quadrivariate and pentavariate normal cumulative distribution functions

This article provides explicit integration rules for the quadrivariate and the pentavariate normal distribution. By analytically reducing the dimension of the problem and simplifying the functions to be integrated, these rules form the basis for a numerical evaluation scheme yielding an observed maximum error in the order of 10^{? 7} and a computational time of less than 10^{? 6} s. The implementation is very straightforward as it is based on a classical Gauss–Legendre quadrature. Order statistics are also dealt with.     

Computation of multivariate normal probabilities with polar coordinate systems

We consider the problem of evaluation of the probability that all elements of a multivariate normally distributed vector have non-negative coordinates; this probability is called the non-centred orthant probability. The necessity for the evaluation of this probability arises frequently in statistics. The probability is defined by the integral of the probability density function. However, direct numerical integration is not practical. In this article, a method is proposed for the computation of the probability. The method involves the evaluation of a measure on a unit sphere surface in p-dimensional space that satisfies conditions derived from a covariance matrix. The required computational time for the p-dimensional problem is proportional to p²·2^p?1, and it increases at a rate that is lower than that in the case of the existing method.     

A Generalization of the Childs-Moran Result for Orthoscheme Probabilities

van der Vaart (1953, 1955) introduced the orthoscheme probability R_n (c ₁,..., c_n−1 ), meaning the orthant probability of an n -dimensional normal random vector with zero mean and tridiagonal correlation matrix with elements c ₁,..., c_n−1 on the upper diagonal. Childs (1967) conjectured and Moran (1983) proved that the generating function of { R_n (½,...,½)} equals tan z + sin z . This paper derives the generating function of { R_n (τ,½,...,½)}.     

A short table for computing normal orthant probabilities of dimensions 4 and 5

A short of the quadravariate normal orthant probability P ₄ is presented which, except for near singular correlation matrices, is generally accurate to 7+ significant digits. Since the quinvariate orthant probability P ₅is expressed in terms of five copies of P ₄, such a table is also useful for computing P ₅     

A general reduction method for n–variate normal orthant probability

Let X₁,X₂,…,X_n be n normal variates with zero means, unit variances and correlation matrix {p_ij). The orthant probability is the probability that all of the X₁'s are simultaneously positive. This paper presents a general reduction method by extending the method of Childs (1967), and shows that the probability can be represented by a linear combination of some multivariate integrals of order([n/2]?1). As illustrations, we apply the proposed method to the quadrivariate and six–variate cases. Some numerical results are also given.     

Convergence for weighted sums of widely orthant dependent random variables

ABSTRACT

In this article, a complete convergence result and a complete moment convergence result are obtained for the weighted sums of widely orthant dependent random variables under mild conditions. As corollaries, the corresponding results are also obtained under the extended negatively orthant dependent setup. In particular, the complete convergence result generalizes and improves the related known works in the literature.     

Monotonicity of certain multivariate normal probabilities

Let (X, Y) have a (p+q)-dimensional normal distribution and let C, K be convex symmetric sets of dimensions p, q respectively. Under certain restrictions on the mean vector it is shown that P(X ε C, Y ε K) is a monotonically increasing function of the first canonical correlation coefficient between X and Y, provided the remaining coefficients are zero.     

Parameter estimation under orthant restrictions

The authors consider the estimation of linear functions of a multivariate parameter under orthant restrictions. These restrictions are considered both for location models and for the Poisson distribution. For these models, situations are characterized for which the restricted maximum likelihood estimator dominates the unrestricted one for the estimation of any linear function of the parameter. The results obtained point directly to the importance of the dimension of the parameter space, the central direction of the cone and its vertex in these cases. Special attention is given to examples, such as the one‐way analysis of variance, where the estimation of individual interesting linear functions of the parameter, as the coordinates and the differences between them, is also treated.     

New approximation of the normal distribution function

Two new and simple expressions, one for 0 ≤ x ≤ 2 and another for x > 2, for the normal distribution function, are developed which can be easily computed on desk calculators. They are also comparable in accuracy to the one developed by Patry and Keller (1964).     

A simple algorithm for calculating values for folded normal distribution

Folded normal distribution originates from the modulus of normal distribution. In the present article, we have formulated the cumulative distribution function (cdf) of a folded normal distribution in terms of standard normal cdf and the parameters of the mother normal distribution. Although cdf values of folded normal distribution were earlier tabulated in the literature, we have shown that those values are valid for very particular situations. We have also provided a simple approach to obtain values of the parameters of the mother normal distribution from those of the folded normal distribution. These results find ample application in practice, for example, in obtaining the so-called upper and lower α-points of folded normal distribution, which, in turn, is useful in testing of the hypothesis relating to folded normal distribution and in designing process capability control chart of some process capability indices. A thorough study has been made to compare the performance of the newly developed theory to the existing ones. Some simulated as well as real-life examples have been discussed to supplement the theory developed in this article. Codes (generated by R software) for the theory developed in this article are also presented for the ease of application.     

A generalized two-piece skew normal distribution and some of its properties

In this paper, we develop a generalized version of the two-piece skew normal distribution of Kim [On a class of two-piece skew-normal distributions, Statistics 39(6) (2005), pp. 537–553] and derive explicit expressions for its distribution function and characteristic function and discuss some of its important properties. Further estimation of the parameters of the generalized distribution is carried out.     

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     

A generalized normal distribution

Undoubtedly, the normal distribution is the most popular distribution in statistics. In this paper, we introduce a natural generalization of the normal distribution and provide a comprehensive treatment of its mathematical properties. We derive expressions for the nth moment, the nth central moment, variance, skewness, kurtosis, mean deviation about the mean, mean deviation about the median, Rényi entropy, Shannon entropy, and the asymptotic distribution of the extreme order statistics. We also discuss estimation by the methods of moments and maximum likelihood and provide an expression for the Fisher information matrix.     

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.