Asymptotic Expansions for i.i.d. Sums Via Lower-order Convolutions

In this article, we introduce new asymptotic expansions for probability functions of sums of independent and identically distributed random variables. Results are obtained by efficiently employing information provided by lower-order convolutions. In comparison with Edgeworth-type theorems, advantages include improved asymptotic results in the case of symmetric random variables and ease of computation of main error terms and asymptotic crossing points. The first-order estimate can perform quite well against the corresponding renormalized saddlepoint approximation and, pointwise, requires evaluation of only a single convolution integral. While the new expansions are fairly straightforward, the implications are fortuitous and may spur further related work.     

Distributions with Generalized Skewed Conditionals and Mixtures of Such Distributions

Mixtures of skewed distributions (univariate and bivariate) provide flexible models. An alternative modeling approach involves distributions with skewed conditional distributions and mixtures of such distributions. We consider the interrelationships between such models. Examples are provided to show that several skewed distributions already considered in the literature can be viewed as having been constructed via a combination of mixing and skewing.     

A modification of the Fisher-Cornish approximation for the student t percentiles

An asymptotic expansion of the Student t distribution is derived by expanding the standardized Student t distribution in terms of the normal distribution. This expansion is inverted to obtain corresponding asymptotic expansions for the Student t percentiles as functions of the standard normal percentiles0 Using the first two, three or four terms of these expansions, we get approximations of the Student t percentiles which are generally more accurate than the approximations given by Fisher and Cornish(1960) and Koehler (1983).An approximation of the distribution function obtained from this expansion is compared with the approximations discussed by Ling (1978) andfound to be more accurate for moderate degrees of freedom.     

A matrix variance inequality

In the course of solving a variational problem Chernoff (Ann. Probab. 9 (1981) 533) obtained what appears to be a specialized inequality for a variance, namely, that for a standard normal variable X, Var[g(X)]E[g^′(X)]². However, both the simplicity and usefulness of the inequality has generated a plethora of extensions, as well as alternative proofs. All previous papers have focused on a single function. We provide here an inequality for the covariance matrix of k functions, which leads to a matrix inequality in the sense of Loewner.     

Orthogonal polynomials generated by random vectors

Every random q-vector with finite moments generates a set of orthonormal polynomials. These are generated from the basis functions xn = xn11…xnqq using Gram–Schmidt orthogonalization. One can cycle through these basis functions using any number of ways. Here, we give results using minimum cycling. The polynomials look simpler when centered about the mean of X, and still simpler form when X is symmetric about zero. This leads to an extension of the multivariate Hermite polynomial for a general random vector symmetric about zero. As an example, the results are applied to the multivariate normal distribution.     

Series evaluation of Tweedie exponential dispersion model densities

Exponential dispersion models, which are linear exponential families with a dispersion parameter, are the prototype response distributions for generalized linear models. The Tweedie family comprises those exponential dispersion models with power mean-variance relationships. The normal, Poisson, gamma and inverse Gaussian distributions belong to theTweedie family. Apart from these special cases, Tweedie distributions do not have density functions which can be written in closed form. Instead, the densities can be represented as infinite summations derived from series expansions. This article describes how the series expansions can be summed in an numerically efficient fashion. The usefulness of the approach is demonstrated, but full machine accuracy is shown not to be obtainable using the series expansion method for all parameter values. Derivatives of the density with respect to the dispersion parameter are also derived to facilitate maximum likelihood estimation. The methods are demonstrated on two data examples and compared with with Box-Cox transformations and extended quasi-likelihoood.     

Statistical applications of the multivariate skew normal distribution

Azzalini and Dalla Valle have recently discussed the multivariate skew normal distribution which extends the class of normal distributions by the addition of a shape parameter. The first part of the present paper examines further probabilistic properties of the distribution, with special emphasis on aspects of statistical relevance. Inferential and other statistical issues are discussed in the following part, with applications to some multivariate statistics problems, illustrated by numerical examples. Finally, a further extension is described which introduces a skewing factor of an elliptical density.     

The Beta Product Distribution with Complex Parameters

The family consisting of the distributions of products of two independent beta variables is extended to include cases where some of the parameters are not positive but negative or complex. This “beta product” distribution is expressible as a Meijer G function. An example (from risk theory) where such a distribution arises is given: an infinite sum of products of independent random variables is shown to have a distribution that is the product convolution of a complex-parameter beta product and an independent exponential. The distribution of the infinite sum is a new explicit solution of the stochastic equation X = (in law) B(X + C). Characterizations of some G distributions are also proved.     

An approximation to the convolution of gamma distributions

In general, the exact distribution of a convolution of independent gamma random variables is quite complicated and does not admit a closed form. Of all the distributions proposed, the gamma-series representation of Moschopoulos (1985 Moschopoulos, P. G. (1985). The distribution of the sum of independent gamma random variables. Annals of the Institute of Statistical Mathematics 37Part A:541–544. [Google Scholar]) is relatively simple to implement but for particular combinations of scale and/or shape parameters the computation of the weights of the series can result in complications with too much time consuming to allow a large-scale application. Recently, a compact random parameter representation of the convolution has been proposed by Vellaisamy and Upadhye (2009 Vellaisamy, P., Upadhye, N. S. (2009). On the sums of compound negative binomial and gamma random variables. Journal of Applied Probability 46:272–283.[Crossref], [Web of Science ®] , [Google Scholar]) and it allows to give an exact interpretation to the weights of the series. They describe an infinite discrete probability distribution. This result suggested to approximate Moschopoulos’s expression looking for an approximating theoretical discrete distribution for the weights of the series. More precisely, we propose a general negative binomial distribution. The result is an “excellent” approximation, fast and simple to implement for any parameter combination.     

The Mean-shift Outlier Model under Skew Normal Distribution

Asymmetric models have been extensively studied in recent years, in situations where the normality assumption is not satisfied due to lack of symmetry of the data. Techniques for assessing the quality of fit and diagnostic analysis are important for model validation. This paper presents a study of the mean-shift method for detecting outliers in asymmetric normal regression models. Analytical solutions for the estimators of the parameters are obtained using the algorithm. Simulation studies and application to real data are presented, showing the efficiency of the method in detecting outliers.     

A Single Form for t-Distributions and Symmetric Beta Distributions

A single parametric form is given for the symmetric distributions in the Pearson system with finite variance. In effect, these are Student's t-distributions with ν > 2 and all centered symmetric beta distributions. A different parametrization allows the inclusion of the t-distributions with ν ≤2 at the expense of symmetric beta distributions with a low shape parameter.     

An approximation of the inverse moments of the positive hypergeometric distribution

The negative moments of the positive hyper geometric distribution are often approximated by the inverse of the positive moments of this distribution. In this paper, a suitable approximation to the positive hypergeometric distribution is used to obtain the negative moments.     

Tailoring the Gaussian Law for Excess Kurtosis and Skewness by Hermite Polynomials

This article aims at reshaping the normal law to account for tail-thickness and asymmetry, of which there is plenty of evidence in financial data. The inspiration to address the issue was provided by the orthogonality of Hermite polynomials with the Gaussian density as a weight function, with the Gram–Charlier expansion as background. A solution is then devised accordingly, by embodying skewness and excess-kurtosis in a normal kernel, via third- and forth-degree polynomial tune-up. Features of the densities so obtained are established in the main theorem of this article. In addition, a glance is cast at the issue of embodying between-squares correlation, and a solution is outlined.     

The flexible skew Laplace distribution

ABSTRACT

Skew-symmetric distributions have been discussed by several research-ers. In this article we construct a skew-symmetric Laplace distribution, which is the generalization of distribution given by Ali et al. (2009 Ali, M., Pal, M., Woo, J. (2009). Skewed reflected distributions generated by the Laplace kernel. Aust. J. Statist. 38:45–58. [Google Scholar]) and Nekoukhou and Alamatsaz (2012 Nekoukhou, V., Alamatsaz, M.H. (2012). A family of skew-symmetric-Laplace distributions. Statist. Papers. 53(3):685–696.[Crossref], [Web of Science ®] , [Google Scholar]). This new distribution contains more parameters, and this induces flexibility properties, such as unimodality or bimodality. We study on some properties of this distribution. In the last section we also provide an application with a real data. Concerning example has recently been discussed by Nekoukhou et al. (2013 Nekoukhou, V., Alamatsaz, M.H., Aghajani, A.H. (2013). A flexible skew-generalized normal distribution. Commun. Statist. Theory Methods. 42(13):2324–2334.[Taylor & Francis Online], [Web of Science ®] , [Google Scholar]) to apply to their model. We compare the behavior of our distribution to their distribution on this example.     

A spatial Bayesian semiparametric mixture model for positive definite matrices with applications in diffusion tensor imaging

Studies on diffusion tensor imaging (DTI) quantify the diffusion of water molecules in a brain voxel using an estimated 3 × 3 symmetric positive definite (p.d.) diffusion tensor matrix. Due to the challenges associated with modelling matrix‐variate responses, the voxel‐level DTI data are usually summarized by univariate quantities, such as fractional anisotropy. This approach leads to evident loss of information. Furthermore, DTI analyses often ignore the spatial association among neighbouring voxels, leading to imprecise estimates. Although the spatial modelling literature is rich, modelling spatially dependent p.d. matrices is challenging. To mitigate these issues, we propose a matrix‐variate Bayesian semiparametric mixture model, where the p.d. matrices are distributed as a mixture of inverse Wishart distributions, with the spatial dependence captured by a Markov model for the mixture component labels. Related Bayesian computing is facilitated by conjugacy results and use of the double Metropolis–Hastings algorithm. Our simulation study shows that the proposed method is more powerful than competing non‐spatial methods. We also apply our method to investigate the effect of cocaine use on brain microstructure. By extending spatial statistics to matrix‐variate data, we contribute to providing a novel and computationally tractable inferential tool for DTI analysis.     

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.     

Slepian's inequality via the central limit theorem

We define a partial ordering of distributions which is preserved under convolutions and scale transformations. Some properties of this partial ordering are developed and then used to give a new argument for Slepian's inequality (1962).     

A Polynomial Logistic Distribution and its Applications in Finance

The importance of Logistic distribution has been widely recognized in many applied areas such as, demography, population studies, finance, agriculture, etc. Since its introduction as a model, much attention has been paid to the study of several generalizations of it, which would offer additional flexibility when data fitting is chased. In the present paper we introduce and develop a natural generalization of the Logistic distribution by considering a probability model whose logit cumulative distribution function transformation is of polynomial type. The performance of the model's fitting to financial data, using different parameter estimation methods, is also investigated.     

The Beta Exponential-Geometric Distribution

A new four-parameter distribution with decreasing, increasing, and upside-down bathtub failure rate called the beta exponential-geometric distribution is proposed. The new distribution, generated from the logit of a beta random variable, extends the exponential-geometric distribution of Adamidis and Loukas (1998 Adamidis , K. , Loukas , S. ( 1998 ). A lifetime distribution with decreasing failure rate . Statistics and Probability Letters 39 : 35 – 42 .[Crossref], [Web of Science ®] , [Google Scholar]) and some other distributions. A comprehensive mathematical treatment of this distribution is provided. Some expressions for the moment generating function, moments, order statistics, and Rényi entropy of the new distribution are derived. Estimation of the stress-strength parameter is also obtained. The model parameters are estimated by the maximum likelihood method and Fisher information matrix is discussed. Finally, an application to a real data set is illustrated.