Saddlepoint density and distribution functions for the ratio of two linear functions and the product of generalized gamma variates

The generalized gamma distribution is a flexible and attractive distribution because it incorporates several well-known distributions, i.e., gamma, Weibull, Rayleigh, and Maxwell. This article derives saddlepoint density and distribution functions for the ratio of two linear functions of generalized gamma variables and the product of n independent generalized gamma variables. Simulation studies are used to evaluate the accuracy of the saddlepoint approximations. The saddlepoint approximations are fast, easy, and very accurate.     

Moments of generalized quadratic forms and distributions of certain multivariate test statistics

Summary Moments and distributions of quadratic forms or quadratic expressions in normal variables are available in literature. Such quadratic expressions are shown to be equivalent to a linear function of independent central or noncentral chi-square variables. Some results on linear functions of generalized quadratic forms are also available in literature. Here we consider an arbitrary linear function of matrix-variate gamma variables. Moments of the determinant of such a linear function are evaluated when the matrix-variate gammas are independently distributed. By using these results, arbitrary non-null moments as well as the non-null distribution of the likelihood ratio criterion for testing the hypothesis of equality of covariance matrices in independent multivariate normal populations are derived. As a related result, the distribution of a linear function of independent matrix-variate gamma random variables, which includes linear functions of independent Wishart matrices, is also obtained. Some properties of generalized special functions of several matrix arguments are used in deriving these results.     

The effect of dependence on distribution of the functions of random variables

In this article, we study the effect of dependence on the distributional properties of functions of two random variables. Expressions for the cumulative distribution functions of the linear combinations, products, and ratios of two dependent random variables in terms of their associated copula are derived. We discuss the effect of dependence on quantities such as the variances of linear combinations of functions, the value-at-risk measure, and the stress–strength parameter. Several examples, a simulation study, and a real data analysis are provided to illustrate the result.     

On approximating the central and noncentral multivariate gamma distributions

In this paper a finite series approximation involving Laguerre polynomials is derived for central and noncentral multivariate gamma distributions. It is shown that if one approximates the density of any k nonnegative continuous random variables by a finite series of Laguerre polynomials up to the (n₁, …, n_k)th degree, then all the mixed moments up to the order (n₁, …, n_k) of the approximated distribution equal to the mixed moments up to the same order of the random variables. Some numerical results are given for the bivariate central and noncentral multivariate gamma distributions to indicate the usefulness of the approximations.     

Bayesian inference for the variance components in general mixed linear models

The general mixed linear model, containing both the fixed and random effects, is considered. Using gamma priors for the variance components, the conditional posterior distributions of the fixed effects and the variance components, conditional on the random effects, are obtained. Using the normal approximation for the multiple t distribution, approximations are obtained for the posterior distributions of the variance components in infinite series form. The same approximation Is used to obtain closed expressions for the moments of the variance components. An example is considered to illustrate the procedure and a numerical study examines the closeness of the approximations.     

The Use of Noncentral F Approximations for Calculation of Power and Sample Size

Approximations to the noncentral F distribution yield surprisingly accurate results for power and sample size problems arising from linear hypotheses about normal random variables. The approximations are easy to use with a desk (or hand-held) calculator that computes cumulative F probabilities. These approximations are particularly advantageous for testing the hypothesis that differences among the means are small against the alternative that the differences are large.     

LOG-LINEAR MODELS FOR MEAN AND DISPERSION IN MIXED POISSON REGRESSION MODELS

This paper is concerned with the analysis of repeated measures count data overdispersed relative to a Poisson distribution, with the overdispersion possibly heterogeneous. To accommodate the overdispersion, the Poisson random variable is compounded with a gamma random variable, and both the mean of the Poisson and the variance of the gamma are modelled using log linear models. Maximum likelihood estimates (MLE) are then obtained. The paper also gives extended quasi-likelihood estimates for a more general class of compounding distributions which are shown to be approximations to the MLEs obtained for the gamma case. The theory is illustrated by modelling the determination of asbestos fibre intensity on membrane filters mounted on microscope slides.     

Distribution of products of independently distributed pathway random variables

A lot of work has been done on products and ratios of random variables by Provost and his co-workers, see, for example, Provost [S.B. Provost, The exact distribution of the ratio of a linear combination of chi-square variables over the root of a product of chi-square variables, Canad. J. Statist. 14 (1986), pp. 61–67; S.B. Provost, The distribution function of a statistic for testing the equality of scale parameters in two gamma populations, Metrika 36 (1989), pp. 337–345]. Here, we extend this idea by introducing a pathway model. The exact density functions of the products of pathway random variables are obtained using the Mellin transform technique. Their computable series forms are derived. The particular cases of the derived results are shown to be associated with the thermonuclear functions and reaction rate probability integral in the theory of nuclear reaction rate, Krätzel integral in applied analyses and inverse Gaussian density in stochastic processes. Graphical representations of the density functions of the product of random variables for the different values of the pathway parameters are shown. The new probability model is fitted to revenue data.     

Testing the equality of correlation matrices

We present results that extend an existing test of equality of correlation matrices. A new test statistic is proposed and is shown to be asymptotically distributed as a linear combination of independent x ² random variables. This new formulation allows us to find the power of the existing test and our extensions by deriving the distribution under the alternative using a linear combination of independent non-central x ² random variables. We also investigate the null and the alternative distribution of two related statistics. The first one is a quadratic form in deviations from a control group with which the remaining k-1 groups are to be compared. The second test is designed for comparing adjacent groups. Several approximations for the null and the alternative distribution are considered and two illustrative examples are provided.     

Study of the quality of several asymptotic and near-exact approximations based on moments for the distribution of the Wilks Lambda statistic

In this paper a measure of proximity of distributions, when moments are known, is proposed. Based on cases where the exact distribution is known, evidence is given that the proposed measure is accurate to evaluate the proximity of quantiles (exact vs. approximated). The measure may be applied to compare asymptotic and near-exact approximations to distributions, in situations where although being known the exact moments, the exact distribution is not known or the expression for its probability density function is not known or too complicated to handle. In this paper the measure is applied to compare newly proposed asymptotic and near-exact approximations to the distribution of the Wilks Lambda statistic when both groups of variables have an odd number of variables. This measure is also applied to the study of several cases of telescopic near-exact approximations to the exact distribution of the Wilks Lambda statistic based on mixtures of generalized near-integer gamma distributions.     

Modified Normal-based Approximation to the Percentiles of Linear Combination of Independent Random Variables with Applications

A modified normal-based approximation for calculating the percentiles of a linear combination of independent random variables is proposed. This approximation is applicable in situations where expectations and percentiles of the individual random variables can be readily obtained. The merits of the approximation are evaluated for the chi-square and beta distributions using Monte Carlo simulation. An approximation to the percentiles of the ratio of two independent random variables is also given. Solutions based on the approximations are given for some classical problems such as interval estimation of the normal coefficient of variation, survival probability, the difference between or the ratio of two binomial proportions, and for some other problems. Furthermore, approximation to the percentiles of a doubly noncentral F distribution is also given. For all the problems considered, the approximation provides simple satisfactory solutions. Two examples are given to show applications of the approximation.     

Approximations to noncentral distributions

Two methods for approximating the distribution of a noncentral random variable by a central distribution in the same family are presented. The first consists of relating a stochastic expansion of a random variable to a corresponding asymptotic expansion for its distribution function. The second approximates the cumulant generating function and is used to provide central χ² and gamma approximations to the noncentral χ² and gamma distributions.     

Exact and Approximate Statistical Inference for Nonlinear Regression and the Estimating Equation Approach

The exact density distribution of the non‐linear least squares estimator in the one‐parameter regression model is derived in closed form and expressed through the cumulative distribution function of the standard normal variable. Several proposals to generalize this result are discussed. The exact density is extended to the estimating equation (EE) approach and the non‐linear regression with an arbitrary number of linear parameters and one intrinsically non‐linear parameter. For a very special non‐linear regression model, the derived density coincides with the distribution of the ratio of two normally distributed random variables previously obtained by Fieler almost a century ago, unlike other approximations previously suggested by other authors. Approximations to the density of the EE estimators are discussed in the multivariate case. Numerical complications associated with the non‐linear least squares are illustrated, such as non‐existence and/or multiple solutions, as major factors contributing to poor density approximation. The non‐linear Markov–Gauss theorem is formulated on the basis of the near exact EE density approximation.     

On the Exact and Near-Exact Distributions of the Product of Generalized Gamma Random Variables and the Generalized Variance

In this article, the authors first obtain the exact distribution of the logarithm of the product of independent generalized Gamma r.v.’s (random variables) in the form of a Generalized Integer Gamma distribution of infinite depth, where all the rate and shape parameters are well identified. Then, by a routine transformation, simple and manageable expressions for the exact distribution of the product of independent generalized Gamma r.v.’s are derived. The method used also enables us to obtain quite easily very accurate, manageable and simple near-exact distributions in the form of Generalized Near-Integer Gamma distributions. Numerical studies are carried out to assess the precision of different approximations to the exact distribution and they show the high accuracy of the approximations provided by the near-exact distributions. As particular cases of the exact distributions obtained we have the distribution of the product of independent Gamma, Weibull, Frechet, Maxwell-Boltzman, Half-Normal, Rayleigh, and Exponential distributions, as well as the exact distribution of the generalized variance, the exact distribution of discriminants or Vandermonde determinants and the exact distribution of any linear combination of generalized Gumbel distributions, as well as yet the distribution of the product of any power of the absolute value of independent Normal r.v.’s.     

The use of Fourier series in the evaluation of probability distribution functions

In statistics, Fourier series have been used extensively in such areas as time series and stochastic processes. These series; however, to a large degree have been neglected with regard to their use in statistical distribution theory. This omission appears quite striking when one considers that, after the elementary functions, the trigonometric functions are the most important functions in applied mathematics. In this paper a procedure is developed for utilizing Fourier series to represent distribution functions of finite range random variables as Fourier series with coefficients easily expressible (using Chebyshev polynomials) In terms of the moments of the distribution. This method allows the evaluation of probabilities for a wide class of distributions. It is applied to the     

Approximated distributions of the weighted sum of correlated chi-squared random variables

We study the asymptotic behavior of the weighted sum of correlated chi-squared random variables. Both chi-squared and normal distributions are proved to approximate the exact distribution. These two approximations are established by matching the first two cumulants. Simulation comparison is made to study the performance of two approximations numerically. We find that the chi-squared approximation performs better than the normal one in the study.     

Approximations to the chi-square distribution

The chi-square distribution arises frequently in applied statistics.Associated with the chi-square random variable with v degrees of freedom are two interdependent variables: the probability integral and the percentage point.Given one of these variables,the other can be obtained from chi-square tables for selected values.In order to overcome the inconvenience of statistical tables and interpolation,many approximations have been suggested.The computational difficulty and accuracy of various approximations is compared.     

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.     

On the Linear Combination of Laplace and Logistic Random Variables

The distribution of linear combinations of random variables arises explicitly in many areas of engineering. This has increased the need to have available the widest possible range of statistical results on linear combinations of random variables. In this note, the exact distribution of the linear combination α X+β Y is derived when X and Y are Laplace and logistic random variables distributed independently of each other. Extensive tabulations of the associated percentage points obtained by inverting the derived distribution are also given.     

Computation of bivariate gamma and inverted beta distribution functions

This paper considers the evaluation of the distribution functions of the bivariate gamma distribution of Wicksell and Kibble, and a bivariate inverted beta distribution. Simple expansions of the distribution functions in terms of marginal quantities and the negative binomial probabilities are derived.