An order-statistics-based method for constructing multivariate distributions with fixed marginals

A new system of multivariate distributions with fixed marginal distributions is introduced via the consideration of random variates that are randomly chosen pairs of order statistics of the marginal distributions. The distributions allow arbitrary positive or negative Pearson correlations between pairs of random variates and generalise the Farlie–Gumbel–Morgenstern distribution. It is shown that the copulas of these distributions are special cases of the Bernstein copula. Generation of random numbers from the distributions is described, and formulas for the Kendall and grade (Spearman) correlations are given. Procedures for data fitting are described and illustrated with examples.     

The matrix-t distribution and its applications in predictive inference

The predictive distributions of the future responses and regression matrix under the multivariate elliptically contoured distributions are derived using structural approach. The predictive distributions are obtained as matrix-t which are identical to those obtained under matrix normal and matrix-t distributions. This gives inference robustness with respect to departures from the reference case of independent sampling from the matrix normal or dependent but uncorrelated sampling from matrix-t distributions. Some successful applications of matrix-t distribution in the field of spatial prediction have been addressed.     

Generation of Over-Dispersed and Under-Dispersed Binomial Variates

Abstract

This article proposes an algorithm for generating over-dispersed and under-dispersed binomial variates with specified mean and variance. The over-dispersed/under-dispersed distributions are derived from correlated binary variables with an underlying continuous multivariate distribution. Different multivariate distributions or different correlation matrices result in different over-dispersed (or under-dispersed) distributions. The over-dispersed binomial distributions that are generated from three different correlation matrices of a multivariate normal are compared with the beta-binomial distribution for various mean and over-dispersion parameters by quantile-quantile (Q-Q) plots. The two distributions appear to be similar. The under-dispersed binomial distribution is simulated to model an example data set that exhibits under-dispersed binomial variation.     

Semi-self-decomposable distributions on Z<Subscript>+</Subscript>

We present a notion of semi-self-decomposability for distributions with support in Z ₊. We show that discrete semi-self-decomposable distributions are infinitely divisible and are characterized by the absolute monotonicity of a specific function. The class of discrete semi-self-decomposable distributions is shown to contain the discrete semistable distributions and the discrete geometric semistable distributions. We identify a proper subclass of semi-self-decomposable distributions that arise as weak limits of subsequences of binomially thinned sums of independent Z ₊-valued random variables. Multiple semi-self-decomposability on Z ₊ is also discussed.     

Computation of the quasi-stationary distributions inM(n)/GI/1/K andGI/M(n)/1/K queues

In this paper, we provide numerical means to compute the quasi-stationary (QS) distributions inM/GI/1/K queues with state-dependent arrivals andGI/M/1/K queues with state-dependent services. These queues are described as finite quasi-birth-death processes by approximating the general distributions in terms of phase-type distributions. Then, we reduce the problem of obtaining the QS distribution to determining the Perron-Frobenius eigenvalue of some Hessenberg matrix. Based on these arguments, we develop a numerical algorithm to compute the QS distributions. The doubly-limiting conditional distribution is also obtained by following this approach. Since the results obtained are free of phase-type representations, they are applicable for general distributions. Finally, numerical examples are given to demonstrate the power of our method.     

The nonexistence of regularization operators

We show that there are no continuous regularization procedures for the extension of distributions. We also show that there are no continuous projection operators from the spaces of distributions onto subspaces of distributions with support on a given closed set.     

Optimal Scaling of Random Walk Metropolis Algorithms with Non-Gaussian Proposals

The asymptotic optimal scaling of random walk Metropolis (RWM) algorithms with Gaussian proposal distributions is well understood for certain specific classes of target distributions. These asymptotic results easily extend to any light tailed proposal distribution with finite fourth moment. However, heavy tailed proposal distributions such as the Cauchy distribution are known to have a number of desirable properties, and in many situations are preferable to light tailed proposal distributions. Therefore we consider the question of scaling for Cauchy distributed proposals for a wide range of independent and identically distributed (iid) product densities. The results are somewhat surprising as to when and when not Cauchy distributed proposals are preferable to Gaussian proposal distributions. This provides motivation for finding proposal distributions which improve on both Gaussian and Cauchy proposals, both for finite dimensional target distributions and asymptotically as the dimension of the target density, d → ∞. Throughout we seek the scaling of the proposal distribution which maximizes the expected squared jumping distance (ESJD).     

Generalized hypergeometric,digamma and trigamma distributions

Summary In this paper we introduce and study new probability distributions named “digamma” and “trigamma” defined on the set of all positive integers. They are obtained as limits of the zero-truncated Type B3 generalized hypergeometric distributions (inverse Pólya-Eggenberger or negative binomial beta distributions), and also by compounding the logarithmic series distributions. The family of digamma distributions has the logarithmic series as a limit and the trigamma as another limit. The trigamma distributions are very close to the zeta (Zipf) distributions. Thus, our new distributions are useful as substitutes of the logarithmic series when the observed frequency data have such a long tail that cannot be fitted by the latter distributions. In the beginning sections we summarize properties of the Type B3 generalized hypergeometric distributions. It is emphasized that the distributions are obtained by compounding a Poisson distribution by “gamma product-ratio” distributions.     

Strictly operator-stable distributions

Strictly operator-stable distributions are defined and discussed. Characterization of strictly stable distributions with exponent 1 is generalized to strictly (α, Q)-stable distributions with α being an eigenvalue of Q.     

Discrete semi-stable distributions

The purpose of this paper is to introduce and study the concepts of discrete semi-stability and geometric semi-stability for distributions with support inZ ₊. We offer several properties, including characterizations, of discrete semi-stable distributions. We establish that these distributions posses the property of infinite divisibility and that their probability generating functions admit canonical representations that are analogous to those of their continuous counterparts. Properties of discrete geometric semi-stable distributions are deduced from the results obtained for discrete semi-stability. Several limit theorems are established and some examples are constructed.