Conditional distributions of generalized order statistics and some characterizations

Generalized order statistics have been introduced in Kamps (1995a). They enable a unified approach to several models of ordered random variables, e.g. (ordinary) order statistics, record values, sequential order statistics, record values from non-identical distributions. The purpose of this paper is to develop conditional distributions of one generalized order statistic given another and to characterize the underlying continuous distribution by different conditional expectations. Well-known results for ordinary order statistics and record values are extended to generalized order statistics. Received: July 1997     

ON THE DISTRIBUTIONS OF ORDER STATISTICS FOR A RANDOM SAMPLE SIZE*

The probability distribution of the i –th and j–th order statistics and of the range R of a sample of size n, taken from a population with probability density function f (x) have been obtained when the sample size n is a random variable N and has: (i) a generalized Poisson distribution; and (ii) a generalized negative bonimial distribution. Specific results are then obtained; (a) when f (x) is uniform over (0,1); and (b) when f(x) is exponential. All the results for N, being a Poisson, binomial and negative binomial rv follow as special cases.     

Fisher information under Gaussian quadrature models

This paper develops formulae to compute the Fisher information matrix for the regression parameters of generalized linear models with Gaussian random effects. The Fisher information matrix relies on the estimation of the response variance under the model assumptions. We propose two approaches to estimate the response variance: the first is based on an analytic formula (or a Taylor expansion for cases where we cannot obtain the closed form), and the second is an empirical approximation using the model estimates via the expectation–maximization process. Further, simulations under several response distributions and a real data application involving a factorial experiment are presented and discussed. In terms of standard errors and coverage probabilities for model parameters, the proposed methods turn out to behave more reliably than does the ‘disparity rule’ or direct extraction of results from the generalized linear model fitted in the last expectation–maximization iteration.     

On some properties of the linear function Poisson binomial and negative binomial distributions

Summary A general model in fluctuations of sums of random variables leading, under certain assumptions, to each of the generalized and linear function Poisson, binomial and negative binomial distributions is presented. Moreover the generating functions and the factorial moments of the linear function Poisson, binomial and negative binomial distributions are obtained in close forms and certain distributional properties are discussed.     

Generalized P‐values and Bayesian evidence in the one‐sided testing problems under exponential distributions

The problem of comparing the frequentist evidence and the Bayesian evidence in the one‐sided testing problems has been widely treated and many researches revealed that these two methods can reach an agreement approximately. However, most of the previous work dealt mainly with situations without nuisance parameters. Since the presence of nuisance parameters is very common in practice, whether these two kinds of evidence still reach an agreement is a problem worthy of study. In this article, we establish in a systematic way under the exponential distributions the agreement of the Bayesian evidence and the generalized frequentist evidence (the generalized P‐value) for a variety of one‐sided testing problems where the nuisance parameters are involved.     

Generalized densities of order statistics

Let X ₁, ... , X _n be independent identically distributed random variables with distribution F . We derive expressions for generalized joint 'densities' of order statistics of X ₁, ... , X _n , for arbitrary distributions F , in terms of Radon–Nikodym derivatives with respect to product measures based on F . We then give formulae for conditional distributions of order statistics and use them to derive results concerning Markov properties of order statistics, formulae for distributions of trimmed sums, and other useful representations. Our approach leads to simple and natural expressions which appear not to have been given before.     

Sharp exponential and entropy bounds on expectations of generalized order statistics

Sharp lower and upper bounds on expected values of generalized order statistics are proven by the use of Moriguti's inequality combined with the Young inequality. The bounds are expressed in terms of exponential moments or entropy. They are attainable providing new characterizations of some nontrivial distributions. Received October 2001/Revised May 2002     

On decompositions of some random variables

The aim of this paper is to present several stochastic analogs of classical formulas for the gamma function. The obtained results provide representation of some random variables as finite or infinite products of independent random variables. Examples include generalized gamma, normal, beta and other distributions.     

Inference on Reliability in Two-parameter Exponential Stress–strength Model

The problem of hypothesis testing and interval estimation of the reliability parameter in a stress–strength model involving two-parameter exponential distributions is considered. Test and interval estimation procedures based on the generalized variable approach are given. Statistical properties of the generalized variable approach and an asymptotic method are evaluated by Monte Carlo simulation. Simulation studies show that the proposed generalized variable approach is satisfactory for practical applications while the asymptotic approach is not satisfactory even for large samples. The results are illustrated using simulated data.     

On a class of discrete distributions arising from the birth-death-with-immigration process

This paper considers a class of distributions which may be regarded as the convolution of a negative binomial and a stopped-sum generalized hypergeometric factorial-moment random variables. Some properties are derived and it is shown that this class of distributions is a subset of distributions for the birth-and-death process with immigration (also reversible counter system). Formulations by mixing, limiting distributions and maximum likelihood equations are also discussed.     

The structure of strategy-proof random social choice functions over product domains and lexicographically separable preferences

We characterize the class of dominant-strategy incentive-compatible (or strategy-proof) random social choice functions in the standard multi-dimensional voting model where voter preferences over the various dimensions (or components) are lexicographically separable. We show that these social choice functions (which we call generalized random dictatorships) are induced by probability distributions on voter sequences of length equal to the number of components. They induce a fixed probability distribution on the product set of voter peaks. The marginal probability distribution over every component is a random dictatorship. Our results generalize the classic random dictatorship result in Gibbard (1977) and the decomposability results for strategy-proof deterministic social choice functions for multi-dimensional models with separable preferences obtained in LeBreton and Sen (1999).     

An application of the generalized Poisson difference distribution to the Bayesian modelling of football scores

The analysis of sports data, in particular football match outcomes, has always produced an immense interest among the statisticians. In this paper, we adopt the generalized Poisson difference distribution (GPDD) to model the goal difference of football matches. We discuss the advantages of the proposed model over the Poisson difference (PD) model, which was also used for the same purpose. The GPDD model, like the PD model, is based on the goal difference in each game that allows us to account for the correlation without explicitly modelling it. The main advantage of the GPDD model is its flexibility in the tails by considering shorter as well as longer tails than the PD distribution. We carry out the analysis in a Bayesian framework in order to incorporate external information, such as historical knowledge or data, through the prior distributions. We model both the mean and the variance of the goal difference and show that such a model performs considerably better than a model with a fixed variance. Finally, the proposed model is fitted to the 2012–2013 Italian Serie A football data, and various model diagnostics are carried out to evaluate the performance of the model.     

Top-<Emphasis Type="Italic">k</Emphasis>-lists

Top-k-lists are introduced as sequences of k-dimensional random vectors with ordered components being k largest observations from a sequence of independent identically distributed random variables. Such lists changing in time are natural stochastic models of ranking tables which appear in many situations in real life, when one wants to keep a track of several best results in a given field. Here we study basic properties of top-k-lists as joint distributions, conditional structures, representations, driving examples of top-k-lists from exponential and uniform distributions, asymptotics and a relation to generalized order statistics.     

Some results on discrete a–monotonicity

Some characterizations and preservation properties of a–monotonicity are investigated. In particular it is shown that a mixture of binomial or Poisson distributions preserves the a–monotonicity of the mixing distribution. Also, a characterization of a class of infinitely divisible distributions, within the discrete α–unimodal class, is given. This characterization is the discrete version of that of ALAMATSAZ (1985) and KOTZ and STEUTEL (1988) in the continuous case.     

Reparameterized inverse gamma regression models with varying precision

In this article, we propose a mean linear regression model where the response variable is inverse gamma distributed using a new parameterization of this distribution that is indexed by mean and precision parameters. The main advantage of our new parametrization is the straightforward interpretation of the regression coefficients in terms of the expectation of the positive response variable, as usual in the context of generalized linear models. The variance function of the proposed model has a quadratic form. The inverse gamma distribution is a member of the exponential family of distributions and has some distributions commonly used for parametric models in survival analysis as special cases. We compare the proposed model to several alternatives and illustrate its advantages and usefulness. With a generalized linear model approach that takes advantage of exponential family properties, we discuss model estimation (by maximum likelihood), black further inferential quantities and diagnostic tools. A Monte Carlo experiment is conducted to evaluate the performances of these estimators in finite samples with a discussion of the obtained results. A real application using minerals data set collected by Department of Mines of the University of Atacama, Chile, is considered to demonstrate the practical potential of the proposed model.     

The analysis of proportions in agricultural experiments by a generalized linear mixed model

This paper is concerned with the statistical analysis of proportions involving extra-binomial variation. Extra-binomial variation is inherent to experimental situations where experimental units are subject to some source of variation, e.g. biological or environmental variation. A generalized linear model for proportions does not account for random variation between experimental units. In this paper an extended version of the generalized linear model is discussed with special reference to experiments in agricultural research. In this model it is assumed that both treatment effects and random contributions of plots are part of the linear predictor. The methods are applied to results from two agricultural experiments.     

Tales from tails: On the empirical distributions of forecasting errors and their implication to risk

When evaluating the performances of time series extrapolation methods, both researchers and practitioners typically focus on the average or median performance according to some specific error metric, such as the absolute error or the absolute percentage error. However, from a risk-assessment point of view, it is far more important to evaluate the distributions of such errors, and especially their tails. For instance, a lack of normality and symmetry in error distributions can have significant implications for decision making, such as in stock control. Moreover, frequently these distributions can only be constructed empirically, as they may be the result of a computationally-intensive non-parametric approach, such as an artificial neural network. This study proposes an approach for evaluating the empirical distributions of forecasting methods and uses it to assess eleven popular time series extrapolation approaches across two different datasets (M3 and ForeDeCk). The results highlight some very interesting tales from the tails.     

On bivariate generalized binomial and negative binomial distributions

The present paper deals with two types of generalized general binomial (binomial or negative binomial) distributions: (i) a univariate general binomial generalized by a bivariate distribution and (ii) a bivariate general binomial generalized by two independent univariate distributions. The probabilities, moments, conditional distributions and regression functions for these distributions are obtained in terms of bipartitional polynomials. Moreover recurrence relations for the probabilities and moments, independent of the bipartitional polynomials, are given. Finally these general results are applied to the (i) Binomial-Bivariate Poisson and (ii) Bivariate Binomial-Poissons distributions.     

Variance-Gamma and Normal-Inverse Gaussian models: Goodness-of-fit to Chinese high-frequency index returns

In this study Variance-Gamma (VG) and Normal-Inverse Gaussian (NIG) distributions are compared with the benchmark of generalized hyperbolic distribution in terms of their fit to the empirical distribution of high-frequency stock market index returns in China. First, we estimate the considered models in a Markov regime switching framework for the identification of different volatility regimes. Second, the goodness-of-fit results are compared at different time scales of log-returns. Third, the goodness-of-fit results are validated through bootstrapping experiments. Our results show that as the time scale of log-returns decrease NIG model outperforms the VG model consistently and the difference between the goodness-of-fit statistics increase. For high-frequency Chinese index returns, NIG model is more robust and provides a better fit to the empirical distributions of returns at different time scales.