Expedient universal robust likelihood method for general dispersed count data

ABSTRACT

We introduce a universal robust likelihood approach for regression analysis of general count data. The robust likelihood function is able to accommodate a wide range of dispersion and is insensitive to model failures. We use simulations and real data analysis to demonstrate the merit of the robust procedure.     

Negative moments of a modified power series distribution and bias of the maximum likelihood estimator

The class of Modified Power Series distributions (MPSD) containing Lagrangian Poisson (LPD) (Consul and Jain, 1973) and Lagrangian binomial distributions (LBD) (Jain and Consul, 1971) was studied by Gupta (1974). We investigate the problem of finding the negative momentsE[X^-r ], of displaced and decapitated Modified Power Series Distributions. We derive the relationship between rand (r-1) negative moments. The negative moments of the decapitated and displaced LPD are obtained. These results are, then, used to find the exact amount of bias in the ML estimators of the parameters in the LPD and the LBD. We have also given the variances of the ML estimator and the minimum variance unbiased estimator of the parameter in the LPD.     

Randomization. pro and contra 1

We discuss a stochastic urn model in which there are two urns A and B. B is originally empty and A contains some fixed number of white and black balls. A player selects integers n>O and b>O. Balls are drawn with replacement in A and balls of the same color are put in B as long as the number of white balls in B exceeds (b-1) times the number of black balls in B. Under this condition, the player stops after drawing n+bx balls and is declared to be a winnear if urn B has x black balls. This number of black balls, x, is shown to have the 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.     

“Faulty inspeciton” distributions

The distribution of the number of items observed to be defective in samples from a finite population. When detection of defectiveness is not certain, is obtained. The distribution of waiting time all a specified number of defective items is observed is also considered.

     

Exact distribution of the convolution of negative binomial random variables

In this note, we derive the exact distribution of S by using the method of generating function and BELL polynomials, where S = X₁ + X₂ + ??? + X_n, and each X_i follows the negative binomial distribution with arbitrary parameters. As a particular case, we also obtain the exact distribution of the convolution of geometric random variables.     

On the differences of two generalized negative binomial variates

The generalized negative binomial (GNB) distribution was defined by Jain and Consul (SIAM J. Appl. Math., 21 (1971)) and was obtained as a particular family of Lagrangian distributions by Consul and Shenton (SIAM J. Appl. Math., 23 (1973)). Consul and Shenton also gave the probability generating function (p.g.f.) and proved many properties of the GNBD. Consul and Gupta (SIAM J. Appl. Math., 39 (1980)) proved that the parameter β must be either zero or 1≤ β ≤ θ^-1 for the GNBD to be a true probability distribution and proved some other properties. Numerous applications and properties of this model have been studied by various researchers. Considering two independent GNB variates X and Y, with parameters (m,β,θ) and (n,β,θ) respectively, the probability distribuition of D = Y-X and its p.g.f. and cumulant generating function have been obtained. A recurrence relation between the cumulants has been established and the first four cumulants, β₁ and β₂ have been derived. Also some moments of the absolute difference |Y-X| have been obtained.     

Exact Short Confidence Intervals from Discrete Data

Suppose that X is a discrete random variable whose possible values are {0, 1, 2,⋯} and whose probability mass function belongs to a family indexed by the scalar parameter θ . This paper presents a new algorithm for finding a 1 − α confidence interval for θ based on X which possesses the following three properties: (i) the infimum over θ of the coverage probability is 1 − α ; (ii) the confidence interval cannot be shortened without violating the coverage requirement; (iii) the lower and upper endpoints of the confidence intervals are increasing functions of the observed value x . This algorithm is applied to the particular case that X has a negative binomial distribution.     

负二项分布的两种近似分布及其比较

负二项分布是一个重要的离散型随机变量的分布,可以用泊松分布和正态分布作为其近似分布,通过对两种近似分布进行比较分析,结果表明：在参数q很小时,泊松近似的精度好于正态近似,而且在参数q很小时,即便r不是很大,用泊松分布也能获得负二项分布较好的近似;当参数q较大时,泊松近似效果不好,相比之下,正态近似的结果不错。     

Approximation methods for lognormal characteristic functions

ABSTRACT

The characteristic function of the lognormal distribution is of interest in a number of scientific fields yet an analytic solution remains elusive, making reliable and efficient approximations necessary. In this article, we build on the results of N. C. Beaulieu and A. Saberali in ‘New approximations to the lognormal characteristic function’, by introducing a Taylor- and Bessel function-based partial expansion of the integrand and a Chebyshev quadrature approach. Through computer simulations we show that the Taylor expansion remains accurate and efficient for all commonly computed values, and specify the range of values for which the other two approaches show a significantly stronger performance.     

On the Bayes estimators of the parameters of size-biased generalized power series distributions

ABSTRACT

In this paper, we derive the Bayes estimators of functions of parameters of the size-biased generalized power series distribution under squared error loss function and weighted square error loss function. The results of size-biased GPSD are then used to obtain particular cases of the size-biased negative binomial, size-biased logarithmic series, and size-biased Poisson distributions. These estimators are better than the classical minimum variance unbiased estimators in the sense that they increase the range of the estimation. Finally, an example is provided to illustrate the results and a goodness of fit test is done using the maximum likelihood and Bayes estimators.     

Marshall–Olkin gamma–Weibull distribution with applications

Abstract

A Marshall–Olkin variant of the Provost type gamma–Weibull probability distribution is being introduced in this paper. Some of its statistical functions and numerical characteristics among others characteristics function, moment generalizing function, central moments of real order are derived in the computational series expansion form and various illustrative special cases are discussed. This density function is utilized to model two real data sets. The new distribution provides a better fit than related distributions as measured by the Anderson–Darling and Cramér–von Mises statistics. The proposed distribution could find applications for instance in the physical and biological sciences, hydrology, medicine, meteorology, engineering, etc.     

Efficient regression modeling for correlated and overdispersed count data

Abstract

The objective of this paper is to propose an efficient estimation procedure in a marginal mean regression model for longitudinal count data and to develop a hypothesis test for detecting the presence of overdispersion. We extend the matrix expansion idea of quadratic inference functions to the negative binomial regression framework that entails accommodating both the within-subject correlation and overdispersion issue. Theoretical and numerical results show that the proposed procedure yields a more efficient estimator asymptotically than the one ignoring either the within-subject correlation or overdispersion. When the overdispersion is absent in data, the proposed method might hinder the estimation efficiency in practice, yet the Poisson regression based regression model is fitted to the data sufficiently well. Therefore, we construct the hypothesis test that recommends an appropriate model for the analysis of the correlated count data. Extensive simulation studies indicate that the proposed test can identify the effective model consistently. The proposed procedure is also applied to a transportation safety study and recommends the proposed negative binomial regression model.     

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.     

An approximation to the convolution of gamma distributions

In general, the exact distribution of a convolution of independent gamma random variables is quite complicated and does not admit a closed form. Of all the distributions proposed, the gamma-series representation of Moschopoulos (1985 Moschopoulos, P. G. (1985). The distribution of the sum of independent gamma random variables. Annals of the Institute of Statistical Mathematics 37Part A:541–544. [Google Scholar]) is relatively simple to implement but for particular combinations of scale and/or shape parameters the computation of the weights of the series can result in complications with too much time consuming to allow a large-scale application. Recently, a compact random parameter representation of the convolution has been proposed by Vellaisamy and Upadhye (2009 Vellaisamy, P., Upadhye, N. S. (2009). On the sums of compound negative binomial and gamma random variables. Journal of Applied Probability 46:272–283.[Crossref], [Web of Science ®] , [Google Scholar]) and it allows to give an exact interpretation to the weights of the series. They describe an infinite discrete probability distribution. This result suggested to approximate Moschopoulos’s expression looking for an approximating theoretical discrete distribution for the weights of the series. More precisely, we propose a general negative binomial distribution. The result is an “excellent” approximation, fast and simple to implement for any parameter combination.     

Functional equations in characterizations of discrete distributions ev rao-rubin condition and its variants

The solutions of Cauchy functional equation and its special cases, studied by various authors while characterizing discrete distributions by Rao-Rubin condition and its variants, axe discussed in this survey paper. Certain new results on the solutions of generalizations of these functional equations are also mentioned     

Finite mixtures of natural exponential families

Let μ be a positive measure concentrated on R⁺ generating a natural exponential family (NEF) F with quadratic variance function V_F(m), m being the mean parameter of F. It is shown that v(dx) = (γ+x)μ(γ ≥ 0) (γ ≥ 0) generates a NEF G whose variance function is of the form l(m)Δ+cΔ(m), where l(m) is an affine function of m, Δ(m) is a polynomial in m (the mean of G) of degree 2, and c is a constant. The family G turns out to be a finite mixture of F and its length-biased family. We also examine the cases when F has cubic variance function and show that for suitable choices of γ the family G has variance function of the form P(m) + Q(m)m where P, Q are polynomials in m of degree m2 while Δ is an affine function of m. Finally we extend the idea to two dimensions by considering a bivariate Poisson and bivariate gamma mixture distribution.     

An EM algorithm for estimating negative binomial parameters

An EM algorithm is proposed for computing estimates of parameters of the negative bi-nomial distribution; the algorithm does not involve further iterations in the M-step, in contrast with the one given in Schader & Schmid (1985). The approach can be applied to the corresponding problem in the logarithmic series distribution. The convergence of the proposed scheme is investigated by simulation, the observed Fisher information is derivedand numerical examples based on real data are presented.     

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.     

Bayesian estimation and case influence diagnostics for the zero-inflated negative binomial regression model

In recent years, there has been considerable interest in regression models based on zero-inflated distributions. These models are commonly encountered in many disciplines, such as medicine, public health, and environmental sciences, among others. The zero-inflated Poisson (ZIP) model has been typically considered for these types of problems. However, the ZIP model can fail if the non-zero counts are overdispersed in relation to the Poisson distribution, hence the zero-inflated negative binomial (ZINB) model may be more appropriate. In this paper, we present a Bayesian approach for fitting the ZINB regression model. This model considers that an observed zero may come from a point mass distribution at zero or from the negative binomial model. The likelihood function is utilized to compute not only some Bayesian model selection measures, but also to develop Bayesian case-deletion influence diagnostics based on q-divergence measures. The approach can be easily implemented using standard Bayesian software, such as WinBUGS. The performance of the proposed method is evaluated with a simulation study. Further, a real data set is analyzed, where we show that ZINB regression models seems to fit the data better than the Poisson counterpart.