On finite mixtures of geometric and negative binomial distributions

Finite mixtures of distributions have been getting increasing use in the applied literature. In the continuous case, linear combinations of exponentials and gammas have been shown to be well suited for modeling purposes. In the discrete case, the focus has primarily been on continuous mixing, usually of Poisson distributions and typically using gammas to describe the random parameter, But many of these applications are forced, especially when a continuous mixing distribution is used. Instead, it is often prefe-rable to try finite mixtures of geometries or negative binomials, since these are the fundamental building blocks of all discrete random variables. To date, a major stumbling block to their use has been the lack of easy routines for estimating the parameters of such models. This problem has now been alleviated by the adaptation to the discrete case of numerical procedures recently developed for exponential, Weibull, and gamma mixtures. The new methods have been applied to four previously studied data sets, and significant improvements reported in goodness-of-fit, with resultant implications for each affected study.     

Marginal-frechet-bounds for multidimensional distribution functions

In this paper for an n-dimensional distribution function F and each natural number m, m<n, n-dimensional distribution functions F^l(m) and F^y(m) are constructed such that the m-dimensional marginal distributions are the same as for F and the following inequality holds F^l(m) ≦ F ≦ F^u(m)     

The inverse power Lindley distribution

Several probability distributions have been proposed in the literature, especially with the aim of obtaining models that are more flexible relative to the behaviors of the density and hazard rate functions. Recently, two generalizations of the Lindley distribution were proposed in the literature: the power Lindley distribution and the inverse Lindley distribution. In this article, a distribution is obtained from these two generalizations and named as inverse power Lindley distribution. Some properties of this distribution and study of the behavior of maximum likelihood estimators are presented and discussed. It is also applied considering two real datasets and compared with the fits obtained for already-known distributions. When applied, the inverse power Lindley distribution was found to be a good alternative for modeling survival data.     

A note on minimum variance unbiased estimators of power series distributions

We give a simple theorem which easily enables us to get the minimum variance unbiased estimators of manv useful parametric functions of the parmecer in a left cruncated power series distribution. The theorem can be used in both cases:when the truncation is know and (ii) when truncation point is unknown.     

On the type-B integral equation and the distribution of wilks statist for testing independence of several groups of variables

In this paper, WILKS'type-B integral equation is solved in the general form of a series of beta functions and a series of weighted gamma functions as proposed by WALD and BROOKNER 1941. The coefficients in both representations can be obtained by explicit recurrence relartions, therefore the results solve many distributional problems and have the fewest computational difficulties of any representation that has surfaced to date. The radius of convergence of the second series representation is given, whereas the convergence property of the first series representation is given, whereas the convergence property of the first series representation was studied by WALD and Brookner. The exact null distributions of WILKS' statistic A for testing the independence of several groups of variables and of V = -log A are given. The coefficients in all the series representration can be computed recursilvely and hence can be obtained easily with the help of modern computatinal facilities     

Testing the Stochastic Structure of Production: A Flexible Moment-Based Approach

Conventional production function specifications are shown to impose restrictions on the probability distribution of output that cannot be tested with the conventional models. These restrictions have important implications for firm behavior under uncertainty. A flexible representation of a firm's stochastic technology is developed based on the moments of the probability distribution of output. These moments are a unique representation of the technology and are functions of inputs. Large-sample estimators are developed for a linear moment model that is sufficiently flexible to test the restrictions implied by conventional production function specifications. The flexible moment-based approach is applied to milk production data. The first three moments of output are statistically significant functions of inputs. The cross-moment restrictions implied by conventional models are rejected.     

Bivariate generalized gamma distributions of Kibble's type

In this paper, a new type of bivariate generalized gamma (BGG) distribution derived from the bivariate gamma distribution of Kibble [Two-variate gamma-type distribution. Sankh?a 1941;5:137–150] by means of a power transformation is presented. The explicit expressions of statistical properties of the BGG distribution are presented. The estimation of marginal and dependence parameters using the method of moments and the method of inference functions for margins are discussed, and their performance through a Monte Carlo simulation study is assessed. Finally, an example is given to illustrate the applicability of the distributions introduced here.     

Testing the difference between spectral densities of two independent periodically correlated (cyclostationary) time series models

In this work, the asymptotic distribution for the discrete Fourier transform of periodically correlated (PC) processes is applied to test the equality of two PC time series. Then the performance of the proposed method is investigated through the Monte Carlo simulations.     

Extreme statistics from the cosine distribution

The cosine distribution serves as a very good teaching example for which exact moment properties of extreme order statistics can be expressed in terms of elementary functions. This article presents the limiting joint distribution of the extreme order statistics for the cosine distribution. It is shown that burrows'(1986) result is a special case of the result presented in this article.     

Saddlepoint density and distribution functions for the ratio of two linear functions and the product of generalized gamma variates

The generalized gamma distribution is a flexible and attractive distribution because it incorporates several well-known distributions, i.e., gamma, Weibull, Rayleigh, and Maxwell. This article derives saddlepoint density and distribution functions for the ratio of two linear functions of generalized gamma variables and the product of n independent generalized gamma variables. Simulation studies are used to evaluate the accuracy of the saddlepoint approximations. The saddlepoint approximations are fast, easy, and very accurate.     

A simple formula based on quantiles for the moments of beta generalized distributions

In this article, we derive explicit expansions for the moments of beta generalized distributions from power series expansions for the quantile functions of the baseline distributions. We apply our formula to the beta normal, beta Student t, beta gamma and beta beta generalized distributions. We propose a simple way to express the quantile function of any beta generalized distribution as a power series expansion with known coefficients.     

A taxonomy of mixing and outcome distributions based on conjugacy and bridging

Abstract

The generalized linear mixed model (GLMM) is commonly used for the analysis of hierarchical non Gaussian data. It combines an exponential family model formulation with normally distributed random effects. A drawback is the difficulty of deriving convenient marginal mean functions with straightforward parametric interpretations. Several solutions have been proposed, including the marginalized multilevel model (directly formulating the marginal mean, together with a hierarchical association structure) and the bridging approach (choosing the random-effects distribution such that marginal and hierarchical mean functions share functional forms). Another approach, useful in both a Bayesian and a maximum-likelihood setting, is to choose a random-effects distribution that is conjugate to the outcome distribution. In this paper, we contrast the bridging and conjugate approaches. For binary outcomes, using characteristic functions and cumulant generating functions, it is shown that the bridge distribution is unique. Self-bridging is introduced as the situation in which the outcome and random-effects distributions are the same. It is shown that only the Gaussian and degenerate distributions have well-defined cumulant generating functions for which self-bridging holds.     

On Some Properties of Autopersistence Functions and Autopersistence Graphs

The classical autocorrelation function might not be very informative when measuring a dependence in binary time series. Recently, alternative tools, namely the autopersistence functions (APF) and their sample counterparts, the autpersistence graphs (APG), have been proposed for the analysis of dependent dichotomous variables. In this article, we summarize properties of the autopersistence functions for general binary series as well as for some important particular cases. We suggest a normalized version of APF which might be more convenient for a practical use. The asymptotic properties of autopersistence graphs are investigated. The consistency and asymptotic normality is discussed. The theoretical results are illustrated by a simulation study.     

Distribution of the product of a pair of independent two-sided power variates

ABSTRACT

The distributions of algebraic functions of random variables are important in theory of probability and statistics and other areas such as engineering, reliability, and actuarial applications, and many results based on various distributions are available in the literature. The two-sided power distribution is defined on a bounded range, and it is a generalization of the uniform, triangular, and power-function probability distributions. This paper gives the exact distribution of the product of two independent two-sided power-distributed random variables in a computable representation. The percentiles of the product are then computed, and a real data application is given.     

Bayes Hilbert Spaces

A Bayes linear space is a linear space of equivalence classes of proportional σ‐finite measures, including probability measures. Measures are identified with their density functions. Addition is given by Bayes' rule and substraction by Radon–Nikodym derivatives. The present contribution shows the subspace of square‐log‐integrable densities to be a Hilbert space, which can include probability and infinite measures, measures on the whole real line or discrete measures. It extends the ideas from the Hilbert space of densities on a finite support towards Hilbert spaces on general measure spaces. It is also a generalisation of the Euclidean structure of the simplex, the sample space of random compositions. In this framework, basic notions of mathematical statistics get a simple algebraic interpretation. A key tool is the centred‐log‐ratio transformation, a generalization of that used in compositional data analysis, which maps the Hilbert space of measures into a subspace of square‐integrable functions. As a consequence of this structure, distances between densities, orthonormal bases, and Fourier series representing measures become available. As an application, Fourier series of normal distributions and distances between them are derived, and an example related to grain size distributions is presented. The geometry of the sample space of random compositions, known as Aitchison geometry of the simplex, is obtained as a particular case of the Hilbert space when the measures have discrete and finite support.     

Evaluation of Tweedie exponential dispersion model densities by Fourier inversion

The Tweedie family of distributions is a family of exponential dispersion models with power variance functions V(μ)=μ ^p for . These distributions do not generally have density functions that can be written in closed form. However, they have simple moment generating functions, so the densities can be evaluated numerically by Fourier inversion of the characteristic functions. This paper develops numerical methods to make this inversion fast and accurate. Acceleration techniques are used to handle oscillating integrands. A range of analytic results are used to ensure convergent computations and to reduce the complexity of the parameter space. The Fourier inversion method is compared to a series evaluation method and the two methods are found to be complementary in that they perform well in different regions of the parameter space.     

A Note On Beta Inverse-Weibull Distribution

The inverse Weibull distribution is one of the widely applied distribution for problems in reliability theory. In this article, we introduce a generalization—referred to as the Beta Inverse-Weibull distribution—generated from the logit of a beta random variable. We provide a comprehensive treatment of the mathematical properties of the Beta Inverse-Weibull distribution. The shapes of the corresponding probability density function and the hazard rate function have been obtained and graphical illustrations have been given. The distribution is found to be unimodal. Results for the non central moments are obtained. The relationship between the parameters and the mean, variance, skewness, and kurtosis are provided. The method of maximum likelihood is proposed for estimating the model parameters. We hope that this generalization will attract wider applicability to the problems in reliability theory and mechanical engineering.     

EM-based maximum likelihood parameter estimation for multivariate generalized hyperbolic distributions with fixed λ

Generalized Hyperbolic distribution (Barndorff-Nielsen 1977) is a variance-mean mixture of a normal distribution with the Generalized Inverse Gaussian distribution. Recently subclasses of these distributions (e.g., the hyperbolic distribution and the Normal Inverse Gaussian distribution) have been applied to construct stochastic processes in turbulence and particularly in finance, where multidimensional problems are of special interest. Parameter estimation for these distributions based on an i.i.d. sample is a difficult task even for a specified one-dimensional subclass (subclass being uniquely defined by ) and relies on numerical methods. For the hyperbolic subclass ( = 1), computer program hyp (Blæsild and Sørensen 1992) estimates parameters via ML when the dimensionality is less than or equal to three. To the best of the author's knowledge, no successful attempts have been made to fit any given subclass when the dimensionality is greater than three. This article proposes a simple EM-based (Dempster, Laird and Rubin 1977) ML estimation procedure to estimate parameters of the distribution when the subclass is known regardless of the dimensionality. Our method relies on the ability to numerically evaluate modified Bessel functions of the third kind and their logarithms, which is made possible by currently available software. The method is applied to fit the five dimensional Normal Inverse Gaussian distribution to a series of returns on foreign exchange rates.     

A General Quantile Function Model for Economic and Financial Time Series

This article proposed a general quantile function model that covers both one- and multiple-dimensional models and that takes several existing models in the literature as its special cases. This article also developed a new uniform Bayesian framework for quantile function modelling and illustrated the developed approach through different quantile function models. Many distributions are defined explicitly only via their quanitle functions as the corresponding distribution or density functions do not have an explicit mathematical expression. Such distributions are rarely used in economic and financial modelling in practice. The developed methodology makes it more convenient to use these distributions in analyzing economic and financial data. Empirical applications to economic and financial time series and comparisons with other types of models and methods show that the developed method can be very useful in practice.     

Posterior ranges of functions of parameters under priors with specified quantiles

Suppose some quantiles of the prior distribution of a nonnegative parameter θ are specified. Instead of eliciting just one prior density function, consider the class Γ of all the density functions compatible with the quantile specification. Given a likelihood function, find the posterior upper and lower bounds for the expected value of any real-valued function h(θ), as the density varies in Γ. Such a scheme agrees with a robust Bayesian viewpoint. Under mild regularity conditions about h(θ) and the likelihood, a procedure for finding bounds is derived and applied to an example, after transforming the given functional optimisation problems into finite-dimensional ones.