A Recursive Formulation of the Old Problem of Obtaining Moments from Cumulants and Vice Versa

It seems difficult to find a formula in the literature that relates moments to cumulants (and vice versa) and is useful in computational work rather than in an algebraic approach. Hence I present four very simple recursive formulas that translate moments to cumulants and vice versa in the univariate and multivariate situations.     

Moments for Some Kumaraswamy Generalized Distributions

Explicit expansions for the moments of some Kumaraswamy generalized (Kw-G) distributions (Cordeiro and de Castro, 2011 Cordeiro, G.M., de Castro, M. (2011). A new family of generalized distributions. J. Statist. Computat. Simul. 81:883–898.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) are derived using special functions. We explore the Kw-normal, Kw-gamma, Kw-beta, Kw-t, and Kw-F distributions. These expressions are given as infinite weighted linear combinations of well-known special functions for which numerical routines are readily available.     

A Method for Simulating Multivariate Non Normal Distributions with Specified Standardized Cumulants and Intraclass Correlation Coefficients

Intraclass correlation coefficients (ICCs) are commonly used indices in subject areas such as biometrics, longitudinal data analysis, measurement theory, quality control, and survey research. The properties of the ICCs most often used are derived under the assumption of normality. However, real-world data often violate the normality assumption. In view of this, a computationally efficient procedure is developed for simulating multivariate non normal continuous distributions with specified (a) standardized cumulants, (b) Pearson intercorrelations, and (c) ICCs. The linear model specified is a two-factor design with either fixed or random effects. A numerical example is worked and the results of a Monte Carlo simulation are provided to demonstrate and confirm the methodology.     

Discrete Gamma Distributions: Properties and Parameter Estimations

A two-parameter discrete gamma distribution is derived corresponding to the continuous two parameters gamma distribution using the general approach for discretization of continuous probability distributions. One parameter discrete gamma distribution is obtained as a particular case. A few important distributional and reliability properties of the proposed distribution are examined. Parameter estimation by different methods is discussed. Performance of different estimation methods are compared through simulation. Data fitting is carried out to investigate the suitability of the proposed distribution in modeling discrete failure time data and other count data.     

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.     

Characterizations of Distributions Based on Moments of Residual Life

Let T be a random variable having an absolutely continuous distribution function. It is known that linearity of E(T | T > t) can be used to characterize distributions such as exponential, power and Pareto distribution. In this work, we will extend the above results. More precisely, we characterize the distribution of T by using certain relationships of conditional moments of T. Our results can also be used to obtain new characterization of distributions based on adjacent order statistics or record values.     

Explicit formulas for the cumulants and the vector-valued odd moments of the multivariate linearly skewed elliptical distributions

In this letter explicit expressions are derived for the cumulants and the vector-valued odd moments of the multivariate linearly skewed elliptical family of distributions. The general calculations of such moments are described by multivariate integrals which complicate the calculations. We show how such multivariate computations can be projected into a univariate framework, which extremely simplifies the computations.     

On the Weibull-Gamma frailty model, its infinite moments, and its connection to generalized log-logistic, logistic, Cauchy, and extreme-value distributions

It is shown that the commonly used Weibull-Gamma frailty model has a finite number of finite moments only and that its marginal distribution generalizes the log-logistic distribution. In some cases there is not even a finite variance, and there are cases without a single finite moment. Upon transformation to the entire real line, generalized logistic and generalized Cauchy distributions are introduced and their connection with the previous ones established, as well as with the extreme-value distribution. Apart from intrinsic and classroom value, the family can be of use when formulating non-informative priors in Bayesian data analysis. Also, gauging the amount of finite moments is important when checking regularity conditions in the Weibull-Gamma model. Our findings are illustrated using data from survival in cancer patients.     

Small Sample Tests for Shape Parameters of Gamma Distributions

The introduction of shape parameters into statistical distributions provided flexible models that produced better fit to experimental data. The Weibull and gamma families are prime examples wherein shape parameters produce more reliable statistical models than standard exponential models in lifetime studies. In the presence of many independent gamma populations, one may test equality (or homogeneity) of shape parameters. In this article, we develop two tests for testing shape parameters of gamma distributions using chi-square distributions, stochastic majorization, and Schur convexity. The first one tests hypotheses on the shape parameter of a single gamma distribution. We numerically examine the performance of this test and find that it controls Type I error rate for small samples. To compare shape parameters of a set of independent gamma populations, we develop a test that is unbiased in the sense of Schur convexity. These tests are motivated by the need to have simple, easy to use tests and accurate procedures in case of small samples. We illustrate the new tests using three real datasets taken from engineering and environmental science. In addition, we investigate the Bayes’ factor in this context and conclude that for small samples, the frequentist approach performs better than the Bayesian approach.     

Cumulant Plots for Assessing the Gamma Distribution

This article introduces graphical procedures for assessing the fit of the gamma distribution. The procedures are based on a standardized version of the cumulant generating function. Plots with bands of 95% simultaneous confidence level are developed by utilizing asymptotic and finite-sample results. The plots have linear scales and do not rely on the use of tables or values of special functions. Further, it is found through simulation, that the goodness-of-fit test implied by these plots compares favorably with respect to power to other known tests for the gamma distribution in samples drawn from lognormal and inverse Gaussian distributions.     

Approximate Identifiability, Moments and Censored Data

It is well known that the joint distribution of a pair of random variables ( X,Y ) is not identifiable on the basis of the joint distribution of the function (min ( X,Y ), 1[ X < Y ]). This paper introduces the concept of approximate identifiability and studies its relevance to the function (min ( X,Y ), Y ). It shows that the distribution of ( X,Y ) is approximately identifiable on the basis of the distribution of (min ( X,Y ), Y ). The identification is explicitly executed by a method of moments. The method is applied to the analysis of censored distributions arising in the theory of clinical trials and is compared to the standard method of Kaplan and Meier.     

Gamma-Generalized Inverse Gaussian Class of Distributions with Applications

In this article, a new family of probability distributions with domain in ?⁺ is introduced. This class can be considered as a natural extension of the exponential-inverse Gaussian distribution in Bhattacharya and Kumar (1986 Bhattacharya , S. K. , Kumar , S. ( 1986 ). E-IG model in life testing . Calcutta Statist. Assoc. Bull. 35 : 85 – 90 . [Google Scholar]) and Frangos and Karlis (2004 Frangos , N. , Karlis , D. ( 2004 ). Modelling losses using an exponential-inverse Gaussian distribution . Insur. Math. Econo. 35 : 53 – 67 .[Crossref], [Web of Science ®] , [Google Scholar]). This new family is obtained through the mixture of gamma distribution with generalized inverse Gaussian distribution. We also show some important features such as expressions of probability density function, moments, etc. Special attention is paid to the mixture with the inverse Gaussian distribution, as a particular case of the generalized inverse Gaussian distribution. From the exponential-inverse Gaussian distribution two one-parameter family of distributions are obtained to derive risk measures and credibility expressions. The versatility of this family has been proven in numerical examples.     

Approximations of the Distributions of Test Statistics for Homogeneity of a Product Multinomial Model

Statistics R ^a based on power divergence can be used for testing the homogeneity of a product multinomial model. All R ^a have the same chi-square limiting distribution under the null hypothesis of homogeneity. R ⁰ is the log likelihood ratio statistic and R ¹ is Pearson's X ² statistic. In this article, we consider improvement of approximation of the distribution of R ^a under the homogeneity hypothesis. The expression of the asymptotic expansion of distribution of R ^a under the homogeneity hypothesis is investigated. The expression consists of continuous and discontinuous terms. Using the continuous term of the expression, a new approximation of the distribution of R ^a is proposed. A moment-corrected type of chi-square approximation is also derived. By numerical comparison, we show that both of the approximations perform much better than that of usual chi-square approximation for the statistics R ^a when a ≤ 0, which include the log likelihood ratio statistic.     

Some Relationships Between Skew-Normal Distributions and Order Statistics from Exchangeable Normal Random Vectors

The article explores the relationship between distributions of order statistics from random vectors with exchangeable normal distributions and several skewed generalizations of the normal distribution. In particular, we show that the order statistics of exchangeable normal observations have closed skew-normal distributions, and that the corresponding density function is log-concave when the order statistic is extreme. Special attention is given to the bivariate case, which is related to the univariate skew-normal distribution. The applications discussed focus on the lifetimes of coherent systems.     

On a Relation Between a Family of Distributions Attaining the Bhattacharyya Bound and That of Linear Combinations of the Distributions from an Exponential Family

Abstract

An unbiased estimation problem of a function g(θ) of a real parameter is considered. A relation between a family of distributions for which an unbiased estimator of a function g(θ) attains the general order Bhattacharyya lower bound and that of linear combinations of the distributions from an exponential family is discussed. An example on a family of distributions involving an exponential and a double exponential distributions with a scale parameter is given. An example on a normal distribution with a location parameter is also given.     

Goodness-of-Fit Tests for the Gamma Distribution Based on the Empirical Laplace Transform

We propose a class of goodness-of-fit tests for the gamma distribution that utilizes the empirical Laplace transform. The consistency of the tests as well as their asymptotic distribution under the null hypothesis are investigated. As the decay of the weight function tends to infinity, the test statistics approach limit values related to the first non zero component of Neyman's smooth test for the gamma law. The new tests are compared with other omnibus tests for the gamma distribution.     

Estimation of the Lognormal-Pareto Distribution Using Probability Weighted Moments and Maximum Likelihood

This article deals with the estimation of the lognormal-Pareto and the lognormal-generalized Pareto distributions, for which a general result concerning asymptotic optimality of maximum likelihood estimation cannot be proved. We develop a method based on probability weighted moments, showing that it can be applied straightforwardly to the first distribution only. In the lognormal-generalized Pareto case, we propose a mixed approach combining maximum likelihood and probability weighted moments. Extensive simulations analyze the relative efficiencies of the methods in various setups. Finally, the techniques are applied to two real datasets in the actuarial and operational risk management fields.     

Mixture Distributions Based Methods of Calibration for the Empirical Log-Likelihood Ratio

Empirical likelihood ratio confidence regions based on the chi-square calibration suffer from an undercoverage problem in that their actual coverage levels tend to be lower than the nominal levels. The finite sample distribution of the empirical log-likelihood ratio is recognized to have a mixture structure with a continuous component on [0, + ∞) and a point mass at + ∞. The undercoverage problem of the Chi-square calibration is partly due to its use of the continuous Chi-square distribution to approximate the mixture distribution of the empirical log-likelihood ratio. In this article, we propose two new methods of calibration which will take advantage of the mixture structure; we construct two new mixture distributions by using the F and chi-square distributions and use these to approximate the mixture distributions of the empirical log-likelihood ratio. The new methods of calibration are asymptotically equivalent to the chi-square calibration. But the new methods, in particular the F mixture based method, can be substantially more accurate than the chi-square calibration for small and moderately large sample sizes. The new methods are also as easy to use as the chi-square calibration.     

Testing the Fit of Gamma Distributions Using the Empirical Moment Generating Function

ABSTRACT

This article presents goodness-of-fit tests for two and three-parameter gamma distributions that are based on minimum quadratic forms of standardized logarithmic differences of values of the moment generating function and its empirical counterpart. The test statistics can be computed without reliance to special functions and have asymptotic chi-squared distributions. Monte Carlo simulations are used to compare the proposed test for the two-parameter gamma distribution with goodness-of-fit tests employing empirical distribution function or spacing statistics. Two data sets are used to illustrate the various tests.     

Moments of the trace of a noncentral wishart matrix

Several articles have appeared on the moments of the trace of a noncentral Wishart matrix. Partial results are given in these papers but the representations of these partial results are often too cumbersome for any practical use. In this paper general results are given in compact form which are readily usable in practical problems.