Moments of scale mixtures of skew-normal distributions and their quadratic forms

We obtain the first four moments of scale mixtures of skew-normal distributions allowing for scale parameters. The first two moments of their quadratic forms are obtained using those moments. Previous studies derived the moments, but all relevant results do not allow for scale parameters. In particular, it is shown that the mean squared error becomes an unbiased estimator of σ² with skewed and heavy-tailed errors. Two measures of multivariate skewness are calculated.     

A class of maximum-entropy multivariate distributions

This paper characterizes a class of multivariate distributions that includes the multinormal and is contained in the exponential family. The wide range of possible applications of these distributions is suggested by some of hte characteristics germane to them: First, they maximize Shannon's entropy among all distributions that have finite moments of given orders. As such, they constitute a class of distributions that includes the multinormal and some likely alternatives. Second, they can exhibit several modes, and, further-more, they do so with a relatively small number of parameters (compared to mixtures of multinormals). Third, they are the stationary distributions of certain diffusion processes. Fourth, they approximate, near the multinormal, the multivariate Pearson family. And fifth, the maximum likelihood estimators of their population moments are the sample moments. Two possible methods of estimating the distributions are studied in this paper: maximum likelihood estimation, and a fast procedure that can be used to find consistent estimators of the parameters via sample moments. A FORTTAN subroutine that implements the latter method is also provided.     

Pairwise likelihood estimation for multivariate mixed Poisson models generated by Gamma intensities

Estimating the parameters of multivariate mixed Poisson models is an important problem in image processing applications, especially for active imaging or astronomy. The classical maximum likelihood approach cannot be used for these models since the corresponding masses cannot be expressed in a simple closed form. This paper studies a maximum pairwise likelihood approach to estimate the parameters of multivariate mixed Poisson models when the mixing distribution is a multivariate Gamma distribution. The consistency and asymptotic normality of this estimator are derived. Simulations conducted on synthetic data illustrate these results and show that the proposed estimator outperforms classical estimators based on the method of moments. An application to change detection in low-flux images is also investigated.     

Asymptotic normality of estimators for parameters of a multivariate skew-normal distribution

In this paper, asymptotic normality is established for the parameters of the multivariate skew-normal distribution under two parametrizations. Also, an analytic expression and an asymptotic normal law are derived for the skewness vector of the skew-normal distribution. The estimates are derived using the method of moments. Convergence to the asymptotic distributions is examined both computationally and in a simulation experiment.     

A Method to Generate Multivariate Data with the Desired Moments

We show how it is possible to generate multivariate data which has moments arbitrary close to the desired ones. They are generated as linear combinations of variables with known theoretical moments. It is shown how to derive the weights of the linear combinations in both the univariate and the multivariate setting. The use in bootstrapping is discussed and the method is exemplified with a Monte Carlo simulation where the importance of the ability of generating data with control of higher moments is shown.     

A New Spatial Skew-Normal Random Field Model

Skewness is often present in a wide range of geostatistical problems, and modeling it in the spatial context remains a challenging problem. In this article, we propose and study a new class of spatial skew-normal random fields, defined in terms of the closed multivariate skew-normal distribution. Such fields can be written as the sum of two independent fields: one Gaussian and the other truncated Gaussian. We derive theoretical expressions for the first- and second-order moments, and use them within a method of moments based procedure to estimate the parameters of the model. Data simulated from the model are used to illustrate the methodology developed.     

Multivariate lagrange distributions and a subfamiliy

Second order moments about its means, i.e. the variances and covari-ances for multivariate Lagrange distributions are derived in a matrix form. A subfamily of multivariate Lagrange distributions which can be characterized as the distributions of customers served in a busy period in queues with some conditions are considered. Theorems about their probability functions, one of which is a multivariate generalization of a formula by Takà cs(1989). are given and the means and second order moments about its means are considered. As an example, a multivariate Borel-Tanner distribution is derived.     

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.     

Initializing EM using the properties of its trajectories in Gaussian mixtures

A strategy is proposed to initialize the EM algorithm in the multivariate Gaussian mixture context. It consists in randomly drawing, with a low computational cost in many situations, initial mixture parameters in an appropriate space including all possible EM trajectories. This space is simply defined by two relations between the two first empirical moments and the mixture parameters satisfied by any EM iteration. An experimental study on simulated and real data sets clearly shows that this strategy outperforms classical methods, since it has the nice property to widely explore local maxima of the likelihood function.     

Maximum Likelihood Estimation and Inference in Multivariate Conditionally Heteroscedastic Dynamic Regression Models With Student t Innovations

We provide numerically reliable analytical expressions for the score, Hessian, and information matrix of conditionally heteroscedastic dynamic regression models when the conditional distribution is multivariatet. We also derive one-sided and two-sided Lagrange multiplier tests for multivariate normality versus multivariate t based on the first two moments of the squared norm of the standardized innovations evaluated at the Gaussian pseudo-maximum likelihood estimators of the conditional mean and variance parameters. Finally, we illustrate our techniques through both Monte Carlo simulations and an empirical application to 26 U.K. sectorial stock returns that confirms that their conditional distribution has fat tails.     

Regularized parameter estimation of high dimensional t distribution

We propose penalized-likelihood methods for parameter estimation of high dimensional t distribution. First, we show that a general class of commonly used shrinkage covariance matrix estimators for multivariate normal can be obtained as penalized-likelihood estimator with a penalty that is proportional to the entropy loss between the estimate and an appropriately chosen shrinkage target. Motivated by this fact, we then consider applying this penalty to multivariate t distribution. The penalized estimate can be computed efficiently using EM algorithm for given tuning parameters. It can also be viewed as an empirical Bayes estimator. Taking advantage of its Bayesian interpretation, we propose a variant of the method of moments to effectively elicit the tuning parameters. Simulations and real data analysis demonstrate the competitive performance of the new methods.     

Higher order moments of the estimated tangency portfolio weights

In this paper, we consider the estimated weights of the tangency portfolio. We derive analytical expressions for the higher order non-central and central moments of these weights when the returns are assumed to be independently and multivariate normally distributed. Moreover, the expressions for mean, variance, skewness and kurtosis of the estimated weights are obtained in closed forms. Later, we complement our results with a simulation study where data from the multivariate normal and t-distributions are simulated, and the first four moments of estimated weights are computed by using the Monte Carlo experiment. It is noteworthy to mention that the distributional assumption of returns is found to be important, especially for the first two moments. Finally, through an empirical illustration utilizing returns of four financial indices listed in NASDAQ stock exchange, we observe the presence of time dynamics in higher moments.     

Inverse moments of discrete distributions

In this paper an expression for the inverse moment of order r is given for the truncated binomial and Poisson distributions. This enables one to obtain inverse moments in a finite series. Some applications and multivariate generalizations are also given. The method also enables one to obtain relations between inverse moments and factorial moments and distributions of sums of variables.     

A fractionally integrated Wishart stochastic volatility model

There has recently been growing interest in modeling and estimating alternative continuous time multivariate stochastic volatility models. We propose a continuous time fractionally integrated Wishart stochastic volatility (FIWSV) process, and derive the conditional Laplace transform of the FIWSV model in order to obtain a closed form expression of moments. A two-step procedure is used, namely estimating the parameter of fractional integration via the local Whittle estimator in the first step, and estimating the remaining parameters via the generalized method of moments in the second step. Monte Carlo results for the procedure show a reasonable performance in finite samples. The empirical results for the S&P 500 and FTSE 100 indexes show that the data favor the new FIWSV process rather than the one-factor and two-factor models of the Wishart autoregressive process for the covariance structure.     

Statistical Inference for the Multidimensional Mixed Rasch Model

Inference in generalized linear mixed models with multivariate random effects is often made cumbersome by the high-dimensional intractable integrals involved in the marginal likelihood. This article presents an inferential methodology based on the GEE approach. This method involves the approximations of the marginal likelihood and joint moments of the variables. It is also proposed an approximate Akaike and Bayesian information criterions based on the approximate marginal likelihood using the estimation of the parameters by the GEE approach. The different results are illustrated with a simulation study and with an analysis of real data from health-related quality of life.     

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.     

A new family of multivariate skew slash distribution

In this article, the new family of multivariate skew slash distribution is defined. According to the definition, a stochastic representation of the multivariate skew slash distribution is derived. The first four moments and measures of skewness and kurtosis of a random vector with the multivariate skew slash distribution are obtained. The distribution of quadratic forms for the multivariate skew slash distribution and the non central skew slash χ² distribution are studied. Maximum likelihood inference and real data illustration are discussed. In the end, the potential extension of multivariate skew slash distribution is discussed.     

Some results for multidimensional stationary independent increment processes

The characteristic function, cumulants and moments of vector-valued multidimensional processes, satisfying properties similar to stationary independent increments, are derived. By considering a set of additional postulates for such processes, it is shown that the marginal distribution of such processes is multivariate Poisson. Some of the results in this paper are extensions of the properties of the first two moments of a univariate one-dimensional process with stationary independent increments.     

Population moments and moments of the sample mean

This paper presents new formulae which simultaneously express and estimate moments of the sample mean and estimate population moments, from a simple random sample drawn without replcement from a finite population. By avoiding the generality of the multivariate case, these two problems are not only unified but are made significantly more tractable. Explicit solution are given up to eighth moments. Asymptotic results for infinite populations are also given.     

On the Moments of Multivariate Discrete Distributions Using Finite Difference Operators

This article discusses a general approach to finding the moments of two classes of multivariate discrete distributions, which include those widely used in applied and theoretical statistics. The two classes of multivariate discrete distributions are the multivariate generalized power series distributions (GPSD) and the unified multivariate hypergeometric (UMH) Distributions. The results of Link (1981) follow as special cases.