Recurrence relations of coefficients of the generalized hypergeometric function in multivariate analysis

James(1960) defined the zonal polynomials and used it to represent the joint distributions of latent roots of VVisfiart matrix. The zonal polviionnals played an important role to define the generalized hypergeometric function of symmetric matrix argument Since then, many density functions and moments based on Wishart matrix have been expressed in terms of the generalized hy¬pergeometric Function. The purpose of this paper is to get the recurrence relations for the coefficients of it. In Section 1 we derive a partial differen¬tial equations having the generalized hypergeometric function as the unique solution. Then we ubtain the recurrence relations until order 7 in Section 2.     

Multivariate Matsumoto–Yor property is rather restrictive

Matsumoto and Yor [2001. An analogue of Pitman's 2M-X$2M-X$ theorem for exponential Wiener functionals. Part II: the role of the GIG laws. Nagoya Math. J. 162, 65–86] discovered an interesting invariance property of a product of the generalized inverse Gaussian (GIG) and the gamma distributions. For univariate random variables or symmetric positive definite random matrices it is a characteristic property for this pair of distributions. It appears that for random vectors the Matsumoto–Yor property characterizes only very special families of multivariate GIG and gamma distributions: components of the respective random vectors are grouped into independent subvectors, each subvector having linearly dependent components. This complements the version of the multivariate Matsumoto–Yor property on trees and related characterization obtained in Massam and Weso?owski [2004. The Matsumoto–Yor property on trees. Bernoulli 10, 685–700].     

On some distributions arising from certain generalized sampling schemes

With the notion of success in a series of trials extended tD refer to a run of like outcomes, several new distributions are obtained as the result of sampling from an urn without replacement. or with additional replacements., In this context, the hy-pergeometric, negative hypergeometric, logarithmic series, generalized Waring, Polya and inverse Polya distributions are extended and their properties are studied     

Properties and Inference on the Skew-Curved-Symmetric Family of Distributions

In this article, we study some results related to a specific class of distributions, called skew-curved-symmetric family of distributions that depends on a parameter controlling the skewness and kurtosis at the same time. Special elements of this family which are studied include symmetric and well-known asymmetric distributions. General results are given for the score function and the observed information matrix. It is shown that the observed information matrix is always singular for some special cases. We illustrate the flexibility of this class of distributions with an application to a real dataset on characteristics of Australian athletes.     

Robust finite mixture modeling of multivariate unrestricted skew-normal generalized hyperbolic distributions

In this paper, we introduce an unrestricted skew-normal generalized hyperbolic (SUNGH) distribution for use in finite mixture modeling or clustering problems. The SUNGH is a broad class of flexible distributions that includes various other well-known asymmetric and symmetric families such as the scale mixtures of skew-normal, the skew-normal generalized hyperbolic and its corresponding symmetric versions. The class of distributions provides a much needed unified framework where the choice of the best fitting distribution can proceed quite naturally through either parameter estimation or by placing constraints on specific parameters and assessing through model choice criteria. The class has several desirable properties, including an analytically tractable density and ease of computation for simulation and estimation of parameters. We illustrate the flexibility of the proposed class of distributions in a mixture modeling context using a Bayesian framework and assess the performance using simulated and real data.

     

Mixture of the Riesz distribution with respect to the generalized multivariate gamma distribution

Wishart natural exponential families (NEFs) characterized by Letac (1989) are extended to the Riesz NEFs on symmetric matrices. These families are characterized by their variance functions defined in Hassairi and Lajmi (2001). This work uses a particular basis of these NEFs to describe the class of the generalized multivariate gamma distributions and then to study the statistical model obtained by the mixture of this distribution with the Riesz one on the space of symmetric matrices.     

A generalization for two-sided power distributions and adjusted method of moments

A new rich class of generalized two-sided power (TSP) distributions, where their density functions are expressed in terms of the Gauss hypergeometric functions, is introduced and studied. In this class, the symmetric distributions are supported by finite intervals and have normal shape densities. Our study on TSP distributions also leads us to a new class of discrete distributions on {0, 1, …, k}. In addition, a new numerical method for parameter estimation using moments is given.     

On quadratic prediction problems in finite populations

For a finite population and its linear model, Liu and Rong proposed a notion of optimal invariant quadratic unbiased prediction (OIQUP) and offered two methods for studying this notion, in which the first is incomplete. In this note, we mainly aim at fulfilling the first approach used by Liu and Rong by considering a transformed matrix equation set through permutation matrix techniques. Solvability of the matrix equation set, optimality of the resulting predictor, and equivalence of the representations of OIQUP, derived in this note and by Liu and Rong, are investigated in detail. In addition, an application to predicting population variance is conducted based on a simulated population.     

GOODNESS OF FIT STATISTICS BASED ON THE SPACINGS OF COMPLETE OR CENSORED SAMPLES

A goodness-of-fit statistic Z is defined in terms of the spacings generated by the order statistics of a complete or a censored sample from a distribution of the type (l/)f((x-μ)/), μ and unknown. The distribution of Z is studied, mostly through Monte Carlo methods. The power properties of Z for testing Exponential, Uniform, Normal, Gamma and Logistic distributions are discussed; Z is shown to be more powerful than the Smith & Bain (1976) correlation statistic, except for testing Uniform, Normal and Logistic (symmetric distributions) against symmetric alternatives. The statistic Z is generalized to test the goodness-of-fit from κ 2 independent complete or censored samples.     

Antipodally symmetric distributions for orientation statistics

The conventional antipodally symmetric Bingham matrix distribution on the Stiefel manifold is generalized. Large sample maximum likelihood estimation and uniformity tests are discussed, and a parametric model for axial orientations (X-shapes) is suggested. A generalization of the Khatri-Mardia matrix distribution is developed to provide a model suitable for hybrids (T-shapes).     

On the matrix-variate generalized hyperbolic distribution and its Bayesian applications

In the first part of the paper, we introduce the matrix-variate generalized hyperbolic distribution by mixing the matrix normal distribution with the matrix generalized inverse Gaussian density. The p-dimensional generalized hyperbolic distribution of [Barndorff-Nielsen, O. (1978). Hyperbolic distributions and distributions on hyperbolae. Scand. J. Stat., 5, 151–157], the matrix-T distribution and many well-known distributions are shown to be special cases of the new distribution. Some properties of the distribution are also studied. The second part of the paper deals with the application of the distribution in the Bayesian analysis of the normal multivariate linear model.     

Diagonalization matrix and its application in distribution theory

Some matrix representations of diverse diagonal arrays are studied in this work; the results allow new definitions of classes of elliptical distributions indexed by kernels mixing Hadamard and usual products. A number of applications are derived in the setting of prior densities from the Bayesian multivariate regression model and families of non-elliptical distributions, such as the matrix multivariate generalized Birnbaum–Saunders density. The philosophy of the research about matrix representations of quadratic and inverse quadratic forms can be extended as a methodology for exploring possible new applications in non-standard distributions, matrix transformations and inference.     

Modified method on the means for several log-normal distributions

Among statistical inferences, one of the main interests is drawing the inferences about the log-normal means since the log-normal distribution is a well-known candidate model for analyzing positive and right-skewed data. In the past, the researchers only focused on one or two log-normal populations or used the large sample theory or quadratic procedure to deal with several log-normal distributions. In this article, we focus on making inferences on several log-normal means based on the modification of the quadratic method, in which the researchers often used the vector of the generalized variables to deal with the means of the symmetric distributions. Simulation studies show that the quadratic method performs well only for symmetric distributions. However, the modified procedure fits both symmetric and skew distribution. The numerical results show that the proposed modified procedure can provide the confidence interval with coverage probabilities close to the nominal level and the hypothesis testing performed with satisfactory results.     

An extension of the ranked set sampling theory

This paper is concerned with ranked set sampling theory which is useful to estimate the population mean when the order of a sample of small size can be found without measurements or with rough methods. Consider n sets of elements each set having size m. All elements of each set are ranked but only one is selected and quantified. The average of the quantified elements is adopted as the estimator. In this paper we introduce the notion of selective probability which is a generalization of a notion from Yanagawa and Shirahata (1976). Uniformly optimal unbiased procedures are found for some (n,m). Furthermore, procedures which are unbiased for all distributions and are good for symmetric distributions are studied for (n,m) which do not allow uniformly optimal unbiased procedures.     

On Properties of Predictive Priors in Linear Models

Utilizing the notion of matching predictives as in Berger and Pericchi, we show that for the conjugate family of prior distributions in the normal linear model, the symmetric Kullback-Leibler divergence between two particular predictive densities is minimized when the prior hyperparameters are taken to be those corresponding to the predictive priors proposed in Ibrahim and Laud and Laud and Ibrahim. The main application for this result is for Bayesian variable selection.     

On some formulas for weighted sum of invariant polynomials of two matrix arguments

Some properties of the generalized binomial coefficients and new coefficients related to symmetric functions are discussed. Several formulas of the weighted sum of invariant polynomials of two matrix arguments which are useful in multivariate distribution theory are presented.     

Dependent models for observations which include angular ones

This paper discusses some stochastic models for dependence of observations which include angular ones. First, we provide a theorem which constructs four-dimensional distributions with specified bivariate marginals on certain manifolds such as two tori, cylinders or discs. Some properties of the submodel of the proposed models are investigated. The theorem is also applicable to the construction of a related Markov process, models for incomplete observations, and distributions with specified marginals on the disc. Second, two maximum entropy distributions on the cylinder are discussed. The circular marginal of each model is distributed as the generalized von Mises distribution which represents a symmetric or asymmetric, unimodal or bimodal shape. The proposed cylindrical model is applied to two data sets.     

SCALE MIXTURES DISTRIBUTIONS IN STATISTICAL MODELLING

This paper presents two types of symmetric scale mixture probability distributions which include the normal, Student t, Pearson Type VII, variance gamma, exponential power, uniform power and generalized t (GT) distributions. Expressing a symmetric distribution into a scale mixture form enables efficient Bayesian Markov chain Monte Carlo (MCMC) algorithms in the implementation of complicated statistical models. Moreover, the mixing parameters, a by-product of the scale mixture representation, can be used to identify possible outliers. This paper also proposes a uniform scale mixture representation for the GT density, and demonstrates how this density representation alleviates the computational burden of the Gibbs sampler.     

Testing for departures from nominal dispersion in generalized nonlinear models with varying dispersion and/or additive random effects

This paper discusses the tests for departures from nominal dispersion in the framework of generalized nonlinear models with varying dispersion and/or additive random effects. We consider two classes of exponential family distributions. The first is discrete exponential family distributions, such as Poisson, binomial, and negative binomial distributions. The second is continuous exponential family distributions, such as normal, gamma, and inverse Gaussian distributions. Correspondingly, we develop a unifying approach and propose several tests for testing for departures from nominal dispersion in two classes of generalized nonlinear models. The score test statistics are constructed and expressed in simple, easy to use, matrix formulas, so that the tests can easily be implemented using existing statistical software. The properties of test statistics are investigated through Monte Carlo simulations.     

IG-symmetry and R-symmetry: Interrelations and applications to the inverse Gaussian theory

The two parameter inverse Gaussian (IG) distribution is often more appropriate and convenient for modelling and analysis of nonnegative right skewed data than the better known and now ubiquitous Gaussian distribution. Its convenience stems from its analytic simplicity and the striking similarities of its methodologies with those employed with the normal theory models. These, known as the G–IG analogies, include the concepts and measures of IG-symmetry, IG-skewness and IG-kurtosis, the IG-analogues of the corresponding classical notions and measures. The new IG-associated entities, although well defined and mathematically transparent, are intuitively and conceptually opaque. In this paper, we first elaborate the importance of the IG distribution and of the G–IG analogies. Then we consider the IG-related root-reciprocal IG (RRIG) distribution and introduce a physically transparent, conceptually clear notion of reciprocal symmetry (R-symmetry) and use it to explain the IG-symmetry. We study the moments and mixture properties of the R-symmetric distributions and the relationship of R-symmetry with IG-symmetry and note that RRIG distribution provides a link, in addition to Tweedie's Laplace transform link, between the Gaussian and inverse Gaussian distributions. We also give a structural characterization of the unimodal R-symmetric distributions. This work further expands the long list of G–IG analogies. Several applications including product convolution, monotonicity of power functions, peakedness and monotone limit theorems of R-symmetry are outlined.