A unified approach to bivariate discrete distributions

Here we introduce a bivariate generalized hypergeometric factorial moment distribution (BGHFMD) through its probability generating function (p.g.f.) whose marginal distributions are the generalized hypergeometric factorial moment distributions introduced by Kemp and Kemp (Bull Int Stat Inst 43:336–338,1969). Well-known bivariate versions of distributions such as binomial, negative binomial and Poisson are special cases of this distribution. A genesis of the distribution and explicit closed form expressions for the probability mass function of the BGHFMD, its factorial moments and the p.g.f.’s of its conditional distributions are derived here. Certain recurrence relations for probabilities, moments and factorial moments of the bivariate distribution are also established.     

Bivariate non-central negative binomial distribution: Another generalisation

Summary Another bivariate generalisation (Type V) of the non-central negative binomial distribution is considered. This generalisation is constructed (i) as a latent structure model; (ii) as an extension of an accident proneness model investigated by Edwards/Gurland (1961); and (iii) as a reversible stochastic counter model. The third construction gives, as a result, an apparently new formulation of the Edwards/Gurland model. The probabilities, moments, recurrence formulas and some properties are given. An application to the data used by Holgate (1966) is considered.     

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.     

The bivariate non-central negative binomial distributions

Summary Four bivariate generalisations (Type I–IV) of the non-central negative binomial distribution (Ong/Lee) are considered. The Type I generalisation is constructed using the latent structure model scheme (Goodman) while the Type II generalisation arises from a variation of this scheme. The Type III generalisation is formed by using the method of random elements in common (Mardia). The Type IV is an extension of the Type I generalisation. Properties of these bivariate distributions including joint central and factorial moments are discussed; several recurrence formulae of the probabilities are given. An application to the childhood accident data of Mellinger et al. is considered with the precision of the Type I maximum likelihood estimates computed.     

A new bivariate generalized Poisson distribution

In this paper, a new bivariate generalized Poisson distribution (GPD) that allows any type of correlation is defined and studied. The marginal distributions of the bivariate model are the univariate GPDs. The parameters of the bivariate distribution are estimated by using the moment and maximum likelihood methods. Some test statistics are discussed and one numerical data set is used to illustrate the applications of the bivariate model.     

Highest posterior mass prediction intervals for binomial and poisson distributions

The problems of constructing prediction intervals (PIs) for the binomial and Poisson distributions are considered. New highest posterior mass (HPM) PIs based on fiducial approach are proposed. Other fiducial PIs, an exact PI and approximate PIs are reviewed and compared with the HPM-PIs. Exact coverage studies and expected widths of prediction intervals show that the new prediction intervals are less conservative than other fiducial PIs and comparable with the approximate one based on the joint sampling approach for the binomial case. For the Poisson case, the HPM-PIs are better than the other PIs in terms of coverage probabilities and precision. The methods are illustrated using some practical examples.     

Information loss for 2 × 2 tables with missing cell counts: binomial case

We formulate likelihood-based ecological inference for 2 × 2 tables with missing cell counts as an incomplete data problem and study Fisher information loss by comparing estimation from complete and incomplete data. In so doing, we consider maximum-likelihood (ML) estimators of probabilities governed by two independent binomial distributions and obtain simplified expressions for their covariance. These expressions reflect well the additional uncertainty arising from the unobserved data compared to complete data tables. We also discuss an approximation to the expected conditional variance of the unobserved entries and ML parameter bias correction. An empirical example is used to demonstrate the results.     

A sequential procedure to test three values of a binomial parameter

Sequential tests to decide among three binomial probabilities are needed in many situations, such as acceptance sampling used to determine the proportion of defective items and presence and absence sampling to decide whether pest species are causing economic damage to a crop such as corn. Approximate error probabilities associated with Armitage's (1950, JRSS B) method of simultaneously conducting three sequential probability ratio tests (SPRTs) are derived for the binomial distribution. These approximations provide a basis for adjusting the error rates used to establish the individual SPRTs so that the desired overall error rates are attained. Monte Carlo simulation is used to evaluate the revised procedure. Received: September 1998     

Estimating a generalized correlation coefficient for a generalized bivariate probit model

In this paper we consider semiparametric estimation of a generalized correlation coefficient in a generalized bivariate probit model. The generalized correlation coefficient provides a simple summary statistic measuring the relationship between the two binary decision processes in a general framework. Our semiparametric estimation procedure consists of two steps, combining semiparametric estimators for univariate binary choice models with the method of maximum likelihood for the bivariate probit model with nonparametrically generated regressors. The estimator is shown to be consistent and asymptotically normal. The estimator performs well in our simulation study.     

Some notes on the non-central negative binomial distribution

Summary In this paper we shall record some facts and further examine certain properties of the non-central negative binomial (NNB) distribution (Laguerre series distribution of Gurland, Chen and Hernandez 1983). We consider, among others, a stochastic formulation (birth-and-death process), the series expansion of the probability distribution and the corresponding series expansion of a generalized exponential distribution (Ong/Lee 1986), and the connection of the NNB distribution with the non-central beta, gamma and non-central gamma distributions. A four-parameter version of the NNB distribution is also presented.     

Maximum likelihood estimation of a binomial proportion using one‐sample misclassified binary data

In this article, we construct two likelihood‐based confidence intervals (CIs) for a binomial proportion parameter using a double‐sampling scheme with misclassified binary data. We utilize an easy‐to‐implement closed‐form algorithm to obtain maximum likelihood estimators of the model parameters by maximizing the full‐likelihood function. The two CIs are a naïve Wald interval and a modified Wald interval. Using simulations, we assess and compare the coverage probabilities and average widths of our two CIs. Finally, we conclude that the modified Wald interval, unlike the naïve Wald interval, produces close‐to‐nominal CIs under various simulations and, thus, is preferred in practice. Utilizing the expressions derived, we also illustrate our two CIs for a binomial proportion parameter using real‐data example.     

Model comparison of coordinate-free multivariate skewed distributions with an application to stochastic frontiers

We consider classes of multivariate distributions which can model skewness and are closed under orthogonal transformations. We review two classes of such distributions proposed in the literature and focus our attention on a particular, yet quite flexible, subclass of one of these classes. Members of this subclass are defined by affine transformations of univariate (skewed) distributions that ensure the existence of a set of coordinate axes along which there is independence and the marginals are known analytically. The choice of an appropriate m-dimensional skewed distribution is then restricted to the simpler problem of choosing m univariate skewed distributions. We introduce a Bayesian model comparison setup for selection of these univariate skewed distributions. The analysis does not rely on the existence of moments (allowing for any tail behaviour) and uses equivalent priors on the common characteristics of the different models. Finally, we apply this framework to multi-output stochastic frontiers using data from Dutch dairy farms.     

Efficient parameter estimation for independent and INAR(1) negative binomial samples

We consider moment based estimation methods for estimating parameters of the negative binomial distribution that are almost as efficient as maximum likelihood estimation and far superior to the celebrated zero term method and the standard method of moments estimator. Maximum likelihood estimators are difficult to compute for dependent samples such as samples generated from the negative binomial first-order autoregressive integer-valued processes. The power method of estimation is suggested as an alternative to maximum likelihood estimation for such samples and a comparison is made of the asymptotic normalized variance between the power method, method of moments and zero term method estimators.     

A negative binomial thinning-based bivariate INAR(1) process

This article considers a bivariate INAR(1) process based on an extension of the negative binomial thinning operator by prespecifying the distribution of the innovations. The dependence is introduced through the innovation components. The existence, uniqueness, strict stationarity, ergodicity, and some probabilistic properties of the process are derived. The estimation methods of conditional least squares and conditional maximum likelihood are considered. Some numerical results of the estimates are presented by simulation study. An application to crime data set is provided.     

Hodges–Lehmann optimality for testing moment conditions

This paper studies the Hodges and Lehmann (1956) optimality of tests in a general setup. The tests are compared by the exponential rates of growth to one of the power functions evaluated at a fixed alternative while keeping the asymptotic sizes bounded by some constant. We present two sets of sufficient conditions for a test to be Hodges–Lehmann optimal. These new conditions extend the scope of the Hodges–Lehmann optimality analysis to setups that cannot be covered by other conditions in the literature. The general result is illustrated by our applications of interest: testing for moment conditions and overidentifying restrictions. In particular, we show that (i) the empirical likelihood test does not necessarily satisfy existing conditions for optimality but does satisfy our new conditions; and (ii) the generalized method of moments (GMM) test and the generalized empirical likelihood (GEL) tests are Hodges–Lehmann optimal under mild primitive conditions. These results support the belief that the Hodges–Lehmann optimality is a weak asymptotic requirement.     

Sources of volatility persistence: A case study of the U.K. pound/U.S. dollar exchange rate returns

This paper utilizes a new approach to examine the inherent nonlinear dynamics of the exchange rate returns volatility. Specifically, we utilize a regime switching threshold (i) generalized autoregressive conditional heteroskedasticity (RS-TGARCH) and (ii) a fractional generalized autoregressive conditional heteroskedasticity (RS-TFIGARCH) model. The RS-TGARCH model is found to be adequate in analyzing the first two moments of the U.K. pound/U.S. dollar monthly exchange rate returns series. The RS-TFIGARCH is found to be adequate for the daily returns series. The volatility persistence and leverage effects associated with exchange rate returns series are jointly tested by means of a Wald Chi-square test.     

Characterizations and Modelling of Multivariate Lack of Memory Property

In this article we establish characterizations of multivariate lack of memory property in terms of the hazard gradient (whenever exists), the survival function and the cumulative hazard function. Based on one of these characterizations we establish a method of generating bivariate lifetime distributions possessing bivariate lack of memory property (BLMP) with specified marginals. It is observed that the marginal distributions have to satisfy certain conditions to be stated. The method generates absolutely continuous bivariate distributions as well as those containing a singular component. Bivariate exponential distributions due to Proschan and Sullo (Reliability and biometry, pp 423–440, 1974), Freund (in J Am Stat Assoc 56:971–977, 1961), Block and Basu (J Am Stat Assoc 89:1091–1097, 1974) and Marshall and Olkin (J Am Math Assoc 62:30–44, 1967) are generated as particular cases among others using the proposed method. Some other distributions generated using the method may be of practical importance. Shock models leading to bivariate distributions possessing BLMP are given. Some closure properties of a class of univariate failure rate functions that can generate distributions possessing BLMP and of the class of bivariate survival functions having BLMP are studied.     

Bounds for median and 50 percetage point of binomial and negative binomial distribution

Inequalities for comparing the binomial distribution function/operating characteristic function with value 0.5 have been established long ago. It is a trivial task to infer from these inequalities bounds for the 50 percentage point of the binomial operating characteristic function. However, it is also possible to use these inequalities to establish good bounds for the median of the binomial distribution function. By a relation of binomial distribution function/operating characteristic function and negative binomial distribution function/operating characteristic function we can get bounds for the respective parameters of negative binomial distribution, too.     

Sampling studies on a test against trend in binomial data

Summary A series of 375 “Monte-Carlo” realizations were performed on a number of groupings of 100 samples each fromt distinct binomial populations (t=3, 4, 5, 7, 9) in order to obtain estimates of the power function of theV andX ² index of dispersion tests in detection of a “trend” in the binomial probabilities.     

A characterization of the binomial and negative binomial distribution by regression of products

Summary Products of independent random variables with different shifts are considered. The binomial and negative binomial populations are characterized by proportionality of regressions of these products with respect to residuals.