Analysis of the human sex ratio by using overdispersion models

For study of the human sex ratio, one of the most important data sets was collected in Saxony in the 19th century by Geissler. The data contain the sizes of families, with the sex of all children, at the time of registration of the birth of a child. These data are reanalysed to determine how the probability for each sex changes with family size. Three models for overdispersion are fitted: the beta–binomial model of Skellam, the 'multiplicative' binomial model of Altham and the double-binomial model of Efron. For each distribution, both the probability and the dispersion parameters are allowed to vary simultaneously with family size according to two separate regression equations. A finite mixture model is also fitted. The models are fitted using non-linear Poisson regression. They are compared using direct likelihood methods based on the Akaike information criterion. The multiplicative and beta–binomial models provide similar fits, substantially better than that of the double-binomial model. All models show that both the probability that the child is a boy and the dispersion are greater in larger families. There is also some indication that a point probability mass is needed for families containing children uniquely of one sex.     

A new mixture model for capture heterogeneity

Summary.  We propose a mixture of binomial and beta–binomial distributions for estimating the size of closed populations. The new mixture model is applied to several real capture–recapture data sets and is shown to provide a convenient, objective framework for model selection. The new model is compared with three alternative models in a simulation study, and the results shed light on the general performance of models in this area. The new model provides a robust flexible analysis, which automatically deals with small capture probabilities.     

Discrete Linnik Weibull distribution

Abstract

We introduce a new family of distributions using truncated discrete Linnik distribution. This family is a rich family of distributions which includes many important families of distributions such as Marshall–Olkin family of distributions, family of distributions generated through truncated negative binomial distribution, family of distributions generated through truncated discrete Mittag–Leffler distribution etc. Some properties of the new family of distributions are derived. A particular case of the family, a five parameter generalization of Weibull distribution, namely discrete Linnik Weibull distribution is given special attention. This distribution is a generalization of many distributions, such as extended exponentiated Weibull, exponentiated Weibull, Weibull truncated negative binomial, generalized exponential truncated negative binomial, Marshall-Olkin extended Weibull, Marshall–Olkin generalized exponential, exponential truncated negative binomial, Marshall–Olkin exponential and generalized exponential. The shape properties, moments, median, distribution of order statistics, stochastic ordering and stress–strength properties of the new generalized Weibull distribution are derived. The unknown parameters of the distribution are estimated using maximum likelihood method. The discrete Linnik Weibull distribution is fitted to a survival time data set and it is shown that the distribution is more appropriate than other competitive models.     

A litter-based approach to risk assessment in developmental toxicity studies via a power family of completely monotone functions

Summary.  A new class of distributions for exchangeable binary data is proposed that originates from modelling the joint success probabilities of all orders by a power family of completely monotone functions. The distribution proposed allows flexible modelling of the dose–response relationship for both the marginal response probability and the pairwise odds ratio and is especially well suited for a litter-based approach to risk assessment. Specifically, the risk of at least one adverse response within a litter takes on a simple form under the distribution proposed and can be reduced further to a generalized linear model if a complementary log–log-link function is used. Existing distributions such as the beta–binomial or folded logistic functions have a tendency to assign too much probability to zero, leading to an underestimation of the risk that at least one foetus is affected and an overestimation of the safe dose. The distribution proposed does not suffer from this problem. With the aid of symbolic differentiation, the distribution proposed can be fitted easily and quickly via the method of scoring. The usefulness of the class of distributions proposed and its superiority over existing distributions are demonstrated in a series of examples involving developmental toxicology and teratology data.     

The two-sided power-binomial distribution for modeling binary count data

The two-sided power (TSP) distribution is a flexible two-parameter distribution having uniform, power function and triangular as sub-distributions, and it is a reasonable alternative to beta distribution in some cases. In this work, we introduce the TSP-binomial model which is defined as a mixture of binomial distributions, with the binomial parameter p having a TSP distribution. We study its distributional properties and demonstrate its use on some data. It is shown that the newly defined model is a useful candidate for overdispersed binomial data.     

A New Markov Binomial Distribution

In this article, we introduce a two-state homogeneous Markov chain and define a geometric distribution related to this Markov chain. We define also the negative binomial distribution similar to the classical case and call it NB related to interrupted Markov chain. The new binomial distribution is related to the interrupted Markov chain. Some characterization properties of the geometric distributions are given. Recursion formulas and probability mass functions for the NB distribution and the new binomial distribution are derived.     

Power binomial exponential distribution: Modeling,simulation and application

ABSTRACT

The binomial exponential 2 (BE2) distribution was proposed by Bakouch et al. as a distribution of a random sum of independent exponential random variables, when the sample size has a zero truncated binomial distribution. In this article, we introduce a generalization of BE2 distribution which offers a more flexible model for lifetime data than the BE2 distribution. The hazard rate function of the proposed distribution can be decreasing, increasing, decreasing–increasing–decreasing and unimodal, so it turns out to be quite flexible for analyzing non-negative real life data. Some statistical properties and parameters estimation of the distribution are investigated. Three different algorithms are proposed for generating random data from the new distribution. Two real data applications regarding the strength data and Proschan's air-conditioner data are used to show that the new distribution is better than the BE2 distribution and some other well-known distributions in modeling lifetime data.     

A family of bivariate binomial distributions generated by extreme bernoulli distributions

In this paper, bivariate binomial distributions generated by extreme bivariate Bernoulli distributions are obtained and studied. Representation of the bivariate binomial distribution generated by a convex combination of extreme bivariate Bernoulli distributions as a mixture of distributions in the class of bivariate binomial distribution generated by extreme bivariate Bernoulli distribution is obtained. A subfamily of bivariate binomial distributions exhibiting the property of positive and negative dependence is constructed. Some results on positive dependence notions as it relates to the bivariate binomial distribution generated by extreme bivariate Bernoulli distribution and a linear combination of such distributions are obtained.     

Bivariate Negative Binomial Generalized Linear Models for Environmental Count Data

We propose a new bivariate negative binomial model with constant correlation structure, which was derived from a contagious bivariate distribution of two independent Poisson mass functions, by mixing the proposed bivariate gamma type density with constantly correlated covariance structure (Iwasaki & Tsubaki, 2005), which satisfies the integrability condition of McCullagh & Nelder (1989, p. 334). The proposed bivariate gamma type density comes from a natural exponential family. Joe (1997) points out the necessity of a multivariate gamma distribution to derive a multivariate distribution with negative binomial margins, and the luck of a convenient form of multivariate gamma distribution to get a model with greater flexibility in a dependent structure with indices of dispersion. In this paper we first derive a new bivariate negative binomial distribution as well as the first two cumulants, and, secondly, formulate bivariate generalized linear models with a constantly correlated negative binomial covariance structure in addition to the moment estimator of the components of the matrix. We finally fit the bivariate negative binomial models to two correlated environmental data sets.     

A bivariate mixture of negative binomial distributions and its applications

Abstract

We construct a new bivariate mixture of negative binomial distributions which represents over-dispersed data more efficiently. This is an extension of a univariate mixture of beta and negative binomial distributions. Characteristics of this joint distribution are studied including conditional distributions. Some properties of the correlation coefficient are explored. We demonstrate the applicability of our proposed model by fitting to three real data sets with correlated count data. A comparison is made with some previously used models to show the effectiveness of the new model.     

A model with long-term survivors: negative binomial Birnbaum-Saunders

Abstract

We propose a cure rate survival model by assuming that the number of competing causes of the event of interest follows the negative binomial distribution and the time to the event of interest has the Birnbaum-Saunders distribution. Further, the new model includes as special cases some well-known cure rate models published recently. We consider a frequentist analysis for parameter estimation of the negative binomial Birnbaum-Saunders model with cure rate. Then, we derive the appropriate matrices for assessing local influence on the parameter estimates under different perturbation schemes. We illustrate the usefulness of the proposed model in the analysis of a real data set from the medical area.     

On some distributions arising in inverse cluster sampling

In this paper, distributions of items sampled inversely in clusters are derived. In particular, negative binomial type of distributions are obtained and their properties are studied. A logarithmic series type of distribution is also defined as a limiting form of the obtained generalized negative binomial distribution.     

Confidence interval estimation for negative binomial group distribution

Negative binomial group distribution was proposed in the literature which was motivated by inverse sampling when considering group inspection: products are inspected group by group, and the number of non-conforming items of a group is recorded only until the inspection of the whole group is finished. The non-conforming probability p of the population is thus the parameter of interest. In this paper, the confidence interval construction for this parameter is investigated. The common normal approximation and exact method are applied. To overcome the drawbacks of these commonly used methods, a composite method that is based on the confidence intervals of the negative binomial distribution is proposed, which benefits from the relationship between negative binomial distribution and negative binomial group distribution. Simulation studies are carried out to examine the performances of our methods. A real data example is also presented to illustrate the application of our method.     

An Extended Binomial Distribution with Applications

Based on Skellam (Poisson difference) distribution, an extended binomial distribution is introduced as a byproduct of extending Moran's characterization of Poisson distribution to the Skellam distribution. Basic properties of the distribution are investigated. Also, estimation of the distribution parameters is obtained. Applications with real data are also described.     

A note on an integer valued time series model with Poisson–negative binomial marginal distribution

Abstract

This paper proposes a new model for autoregressive time series of counts in terms of a convolution of Poisson and negative binomial random variables, known as Poisson–negative binomial (PNB) distribution. The corresponding first-order integer valued time series models are developed and their properties are discussed. The geometric PNB and the geometric semi PNB distributions are also introduced and studied.     

The inverse binomial distribution as a statistical model

A particular case of Jain and Consul's (1971) generalized neg-ative binomial distribution is studied. The name inverse binomial is suggested because of its close relation with the inverse Gaussian distribution. We develop statistical properties including conditional inference of a parameter. An application using real data is given.     

A Family Which Integrates the Generalized Charlier Series and Extended Noncentral Negative Binomial Distributions

The generalized Charlier series distribution includes the binomial distribution, and the noncentral negative binomial distribution extends the negative binomial distribution. The present article proposes a family of counting distributions, which contains both the generalized Charlier series and extended noncentral negative binomial distributions. Compound and mixture formulations of the proposed distribution are given. The probability mass function is expressible in terms of the confluent hypergeometric function as well as the Gauss hypergeometric function. Recursive formulae for probability mass function have been studied by Panjer, Sundt and Jewell, Schröter, Sundt, and Kitano et al. in the context of insurance risk. This article explores horizontal, vertical, triangular, and diagonal recursions. Recursive formulae as well as exact expressions for descending factorial moments are studied. The proposed distribution allows overdispersion or underdispersion relative to a Poisson distribution. An illustrative example of data fitting is given.     

Identification of Salmonella high risk pig-herds in Belgium by using semiparametric quantile regression

Summary.  Consumption of pork that is contaminated with Salmonella is an important source of human salmonellosis world wide. To control and prevent salmonellosis, Belgian pig-herds with high Salmonella infection burden are encouraged to take part in a control programme supporting the implementation of control measures. The Belgian government decided that only the 10% of pig-herds with the highest Salmonella infection burden (denoted high risk herds) can participate. To identify these herds, serological data reported as sample-to-positive ratios (SP-ratios) are collected. However, SP-ratios have an extremely skewed distribution and are heavily subject to confounding seasonal and animal age effects. Therefore, we propose to identify the 10% high risk herds by using semiparametric quantile regression with P -splines. In particular, quantile curves of animal SP-ratios are estimated as a function of sampling time and animal age. Then, pigs are classified into low and high risk animals with high risk animals having an SP-ratio that is larger than the corresponding estimated upper quantile. Finally, for each herd, the number of high risk animals is calculated as well as the beta–binomial p -value reflecting the hypothesis that the Salmonella infection burden is higher in that herd compared with the other herds. The 10% pig-herds with the lowest p -values are then identified as high risk herds. In addition, since high risk herds are supported to implement control measures, a risk factor analysis is conducted by using binomial generalized linear mixed models to investigate factors that are associated with decreased or increased Salmonella infection burden. Finally, since the choice of a specific upper quantile is to a certain extent arbitrary, a sensitivity analysis is conducted comparing different choices of upper quantiles.     

A generalization of the binomial distribution

This paper presents a new departure in the generalization of the binomial distribution by adopting the assumption that the underlying Bernoulli trials take on the values α or β where α < β, rather than the conventional values 0 or 1. The adoption of this more general assumption renders the binomial distribution a four-parameter distribution of the form B(n,p,α,β), and requires the generalization of Romanovsky's (1923) reduction formula for central moments. This paper assesses the usefulness of B(n,p,α,β), and its reduction formula, in the numerical analysis of two problems of interest to decision theorists.     

A Statistical Reanalysis of the Classical Rutherford's Experiment

Modified chi-squared and some newly developed tests for the Poisson, binomial, and an approximated Feller's distribution are discussed. A reanalysis of the classical Rutherford's experimental data on alpha decay is done. Previous analyses of the data were not correct from the point of view of the theory of statistical testing. Tests used show that the data contradict to both Poisson and binomial distribution and do not contradict to a precise “binomial” approximation of Feller's distribution that takes into account a counter's dead time. This gives a plausible statistically correct confirmation of the well-established exponential law of radioactive decay.