Numerical computation of rectangular bivariate and trivariate normal and t probabilities

Algorithms for the computation of bivariate and trivariate normal and t probabilities for rectangles are reviewed. The algorithms use numerical integration to approximate transformed probability distribution integrals. A generalization of Plackett's formula is derived for bivariate and trivariate t probabilities. New methods are described for the numerical computation of bivariate and trivariate t probabilities. Test results are provided, along with recommendations for the most efficient algorithms for single and double precision computations.     

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.     

The computation of standard normal distribution integral in any required precision based on reliability method

Based on reliability theory, the value of the standard normal distribution integral can be obtained by calculating the probability of the failure domain of the linear performance function. After the sample space is divided into some sub-sample spaces, a number of sub-failure domains are obtained. In the paper, the methods of computing the probabilities of sub-failure domains are discussed. All the formulae and the steps of computing the standard normal distribution integral which meet any required precision are given in the paper. Examples show that it is easy for the method to compute the standard normal distribution integral.     

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.     

Regression via order statistics and their concomitants

We consider a five-dimensional normal distribution and derive the exact joint distribution one variable, linear combinations of order statistics from two other variables, and linear combinations of the corresponding concomitants of these order statistics. We show that this joint distribution is a mixture of trivariate unified skew-normal distributions. This mixture representation enables us to predict one variable based on linear combinations of order statistics from two other variables and linear combinations of the corresponding concomitants. We finally illustrate the usefulness of these results by using a real data.     

A note on the inequalities for tail probability of the multivariate normal distribution

In this paper inequalities given by Harkness Godambe (1976) for the rail probabilities of the multivariate normal distribution in the equicorrelated case are improved by using the properties of the characteristic roots of a matrix and of the convex function.     

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.     

An algorithm to compute the cdf of the product of two normal random variables

This paper provides a fortran algorithm that can be used to compute the cdf of the product of two normal distribution random variables. We also give references that provide mathematical properties, tables, and applications of this distribution     

The distribution of the ratio of independent central wishart determinants

Ratios of independent central Wishart determinants are useful statistics in multivariate analyses, particularly in the study of multivariate linear models. A method based on the inversion of characteristic functions is outlined for deriving new experessions for the probability distribution functions of the logarithms of these statistics. Accurate tables of the percentiles of these distributions have been obtained covering many bivariate and trivariate cases which have been computed by approximating these expression.     

A Note on the Distribution of Linear Functions of Two Ordered Correlated Normal Random Variables

We consider the problem of finding the distribution of linear functions of two ordered correlated normal random variables. We derive some distributional properties for these linear statistics and briefly discuss the use of them in location estimation. The connection of the subject with the skew normal distribution is also noted.     

SOME SIMPLE APPROXIMATIONS FOR THE DOUBLY NONCENTRAL z DISTRIBUTION

We give two simple approximations for evaluating the cumulative probabilities of the doubly noncentral z distribution. These can easily be used for evaluating the cumulative probabilities of the doubly noncentral F distribution as well. We compare our results with those obtained by Tiku (1965) using series expansion. An industrial situation where a quality characteristic of interest follows the doubly noncentral z distribution is also cited. However, in this case the exact probabilities could be calculated using results on the ratio of two normal variables.     

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     

Hermite polynomials and truncated trivariate normal distributions

The moments of a trivariate and in general of a multivariate normal distribution, which is truncated with respect to a single variable, are obtained by using properties of Hermite polynomials. An expression for the truncated correlation coefficient is derived in terms of the true population correlation coefficient and the truncation point. The values of this truncated correlation coefficient are tabulated for given values of the true correlation coefficient and a few selected values of the truncation point. A listing of the computer program for this purpose is also given.     

Approximations to the multivariate normal integral

We propose a simple and efficient way to approximate multivariate normal probabilities using univariate and bivariate probabilities. The approximation is computationally tested for the trivariate and quadrivariate normal probabilities. A few problems of higher dimensions were also tested.     

Some properties and a characterization of trivariate and multivariate binomial distributions

By considering a trivariate binomial distribution, the regression equations are obtained and a set of necessary and sufficient conditions are given for the regression to be linear. Under the limiting conditions, the expressions approach those of trivariate Poisson distribution derived by Mahamunulu (1967). A characterization of the trivariate binomial distribution is also established based on the distribution of the sum of two trivariate random vectors. The results are shown to extend to the multidimensional case in a natural way.     

Maximum likelihood estimation of polychoric correlations in r×s ×t contingency tables

This paper discusses the maximum likelihood estimation of the polychoric correlation coefficient based on observed frequencies of three polytomous ordinal variables. The underlying latent variables are assumed to have a standardized trivariate normal distribution. The thresholds and correlations are estimated simultaneously via the scoring algorithm. Some practical applications of the method are discussed. An example is reported to illustrate the theory and some technical details are presented in the Appendix.     

Empirical Likelihood for Censored Linear Regression

In this paper we investigate the empirical likelihood method in a linear regression model when the observations are subject to random censoring. An empirical likelihood ratio for the slope parameter vector is defined and it is shown that its limiting distribution is a weighted sum of independent chi-square distributions. This reduces to the empirical likelihood to the linear regression model first studied by Owen (1991) if there is no censoring present. Some simulation studies are presented to compare the empirical likelihood method with the normal approximation based method proposed in Lai et al. (1995). It was found that the empirical likelihood method performs much better than the normal approximation method.     

BAYESIAN PREDICTION FOR SPATIAL GENERALISED LINEAR MIXED MODELS WITH CLOSED SKEW NORMAL LATENT VARIABLES

Spatial generalised linear mixed models are used commonly for modelling non‐Gaussian discrete spatial responses. In these models, the spatial correlation structure of data is modelled by spatial latent variables. Most users are satisfied with using a normal distribution for these variables, but in many applications it is unclear whether or not the normal assumption holds. This assumption is relaxed in the present work, using a closed skew normal distribution for the spatial latent variables, which is more flexible and includes normal and skew normal distributions. The parameter estimates and spatial predictions are calculated using the Markov Chain Monte Carlo method. Finally, the performance of the proposed model is analysed via two simulation studies, followed by a case study in which practical aspects are dealt with. The proposed model appears to give a smaller cross‐validation mean square error of the spatial prediction than the normal prior in modelling the temperature data set.     

Bayesian estimation methods for contingency tables

Summary A method of inputting prior opinion in contingency tables is described. The method can be used to incorporate beliefs of independence or symmetry but extensions are straightforward. Logistic normal distributions that express such beliefs are used as priors of the cell probabilities and posterior estimates are derived. Empirical Bayes methods are also discussed and approximate posterior variances are provided. The methods are illustrated by a numerical example.