Bridging the Gap: A Generalized Stochastic Process for Count Data

The Bernoulli and Poisson processes are two popular discrete count processes; however, both rely on strict assumptions. We instead propose a generalized homogenous count process (which we name the Conway–Maxwell–Poisson or COM-Poisson process) that not only includes the Bernoulli and Poisson processes as special cases, but also serves as a flexible mechanism to describe count processes that approximate data with over- or under-dispersion. We introduce the process and an associated generalized waiting time distribution with several real-data applications to illustrate its flexibility for a variety of data structures. We consider model estimation under different scenarios of data availability, and assess performance through simulated and real datasets. This new generalized process will enable analysts to better model count processes where data dispersion exists in a more accommodating and flexible manner.     

A (Semi)Parametric Functional Coefficient Logarithmic Autoregressive Conditional Duration Model

In this article, we propose a class of logarithmic autoregressive conditional duration (ACD)-type models that accommodates overdispersion, intermittent dynamics, multiple regimes, and asymmetries in financial durations. In particular, our functional coefficient logarithmic autoregressive conditional duration (FC-LACD) model relies on a smooth-transition autoregressive specification. The motivation lies on the fact that the latter yields a universal approximation if one lets the number of regimes grows without bound. After establishing sufficient conditions for strict stationarity, we address model identifiability as well as the asymptotic properties of the quasi-maximum likelihood (QML) estimator for the FC-LACD model with a fixed number of regimes. In addition, we also discuss how to consistently estimate a semiparametric variant of the FC-LACD model that takes the number of regimes to infinity. An empirical illustration indicates that our functional coefficient model is flexible enough to model IBM price durations.     

The Steady-State Performance of Cumulative Count of a Conforming Control Chart Based On Runs Rules

Cumulative count of conforming control chart is usually used to monitor fraction nonconforming in high-yield processes. In this article, we propose m-of-m control chart based on cumulative count of conforming units for high-yield processes. The steady-state properties of the m-of-m control chart are investigated. We compare performance of the m-of-m control chart with control chart based on cumulative count of conforming units. We present Markov chain model of the m-of-m control chart to evaluate average run length, standard deviation of run length and quartiles.     

A study of economic design of control charts for cumulative count of conforming items

In this paper the economic design of Cumulative Count of Conforming (CCC) control charts to maintain the current control of fraction nonconforming of a process is studied. CCC chart is an attribute chart for monitoring high quality processes by plotting the cumulative count of conforming items between two nonconforming ones on a suitable chart. A process model is proposed to obtain an appropriate loss function. An alogorithm to search for the optimal setting of the sampling and control parameters is derived. Numerical illustrations of the method and some properties of the optimal economic design are provided.     

A new economic scheme for CCC charts with run rules based on average number of inspected items

New statistical techniques and procedures have been developed to control high-yield processes along with looking for process improvement opportunities and minimizing production cost. Cumulative count of conforming control chart is generally a technique for high-quality processes, when nonconforming items are rarely produced. The objective of this study is to design control chart based on cumulative count of conforming items and run rules that develops an economic model based on the average number of inspected items to design m-of-m CCC chart in order to facilitate minimum average cost per item produced. The optimal design parameters for different values of nonconforming fraction and different cost parameters in each scenario are determined. Finally, to analyze the behavior of optimal economic solutions, sensitivity analysis of the model parameters is performed.     

On Baseline Conditions for Zero-Inflated Longitudinal Count Data

We describe a mixed-effect hurdle model for zero-inflated longitudinal count data, where a baseline variable is included in the model specification. Association between the count data process and the endogenous baseline variable is modeled through a latent structure, assumed to be dependent across equations. We show how model parameters can be estimated in a finite mixture context, allowing for overdispersion, multivariate association and endogeneity of the baseline variable. The model behavior is investigated through a large-scale simulation experiment. An empirical example on health care utilization data is provided.     

On the assessment of various factors effecting the improvement in CD4 count of aids patients undergoing antiretroviral therapy using generalized Poisson regression

An important marker for identifying the progression of human immunodeficiency virus (HIV) infection in an individual is the CD4 cell count. Antiretroviral therapy (ART) is a treatment for HIV/AIDS (AIDS, acquired immune-deficiency syndrome) which prolongs and improves the lives of patients by improving the CD4 cell count and strengthen the immune system. This strengthening of the immune system in terms of CD4 count, not only depends on various biological factors, but also other behavioral factors. Previous studies have shown the effect of CD4 count on the mortality, but nobody has attempted to study the factors which are likely to influence the improvement in CD4 count of patients diagnosed of AIDS and undergoing ART. In this paper, we use Poisson regression model (GPR) for exploring the effect of various socio-demographic covariates such as age, gender, geographical location, and drug usage on the improvement in the CD4 count of AIDS patients. However, if the CD4 count data suffers from under or overdispersion, we use GPR model and compare it with negative binomial distribution. Finally, the model is applied for the analysis of data on patients undergoing the ART in the Ram Manohar Lohia Hospital, Delhi, India. The data exhibited overdispersion and hence, GPR model provided the best fit.     

Statistical estimation of time-trends in two-course crop rotations

An assessment of time-trends in yield parameters is essential to the utilization of data from long-term field trials for the comparison of different crop rotations and input regimes, and the identification of sustainable production systems. The barley-vetch rotation established at Breda in northern Syria has provided the basis for estimation of the time-trends in yield data from selected treatments in a two-course crop rotation trial. The model used for the estimation accounts for the effect of rainfall, a major determinant of each annual yield value, and the first-order autocorrelation structure in the errors arising from the same plot over time. An expression for the minimum number of cycles required to detect a significant time-trend has been obtained. Results from the barley-vetch rotation under two fertilizer regimes have been discussed.     

Comparisons of some bivariate regression models

The bivariate negative binomial regression (BNBR) and the bivariate Poisson log-normal regression (BPLR) models have been used to describe count data that are over-dispersed. In this paper, a new bivariate generalized Poisson regression (BGPR) model is defined. An advantage of the new regression model over the BNBR and BPLR models is that the BGPR can be used to model bivariate count data with either over-dispersion or under-dispersion. In this paper, we carry out a simulation study to compare the three regression models when the true data-generating process exhibits over-dispersion. In the simulation experiment, we observe that the bivariate generalized Poisson regression model performs better than the bivariate negative binomial regression model and the BPLR model.     

A joint model for hierarchical continuous and zero-inflated overdispersed count data

Many applications in public health, medical and biomedical or other studies demand modelling of two or more longitudinal outcomes jointly to get better insight into their joint evolution. In this regard, a joint model for a longitudinal continuous and a count sequence, the latter possibly overdispersed and zero-inflated (ZI), will be specified that assembles aspects coming from each one of them into one single model. Further, a subject-specific random effect is included to account for the correlation in the continuous outcome. For the count outcome, clustering and overdispersion are accommodated through two distinct sets of random effects in a generalized linear model as proposed by Molenberghs et al. [A family of generalized linear models for repeated measures with normal and conjugate random effects. Stat Sci. 2010;25:325–347]; one is normally distributed, the other conjugate to the outcome distribution. The association among the two sequences is captured by correlating the normal random effects describing the continuous and count outcome sequences, respectively. An excessive number of zero counts is often accounted for by using a so-called ZI or hurdle model. ZI models combine either a Poisson or negative-binomial model with an atom at zero as a mixture, while the hurdle model separately handles the zero observations and the positive counts. This paper proposes a general joint modelling framework in which all these features can appear together. We illustrate the proposed method with a case study and examine it further with simulations.     

Integer autoregressive models with structural breaks

Even though integer-valued time series are common in practice, the methods for their analysis have been developed only in recent past. Several models for stationary processes with discrete marginal distributions have been proposed in the literature. Such processes assume the parameters of the model to remain constant throughout the time period. However, this need not be true in practice. In this paper, we introduce non-stationary integer-valued autoregressive (INAR) models with structural breaks to model a situation, where the parameters of the INAR process do not remain constant over time. Such models are useful while modelling count data time series with structural breaks. The Bayesian and Markov Chain Monte Carlo (MCMC) procedures for the estimation of the parameters and break points of such models are discussed. We illustrate the model and estimation procedure with the help of a simulation study. The proposed model is applied to the two real biometrical data sets.     

On the Correlation Structure of Gaussian Copula Models for Geostatistical Count Data

We describe a class of random field models for geostatistical count data based on Gaussian copulas. Unlike hierarchical Poisson models often used to describe this type of data, Gaussian copula models allow a more direct modelling of the marginal distributions and association structure of the count data. We study in detail the correlation structure of these random fields when the family of marginal distributions is either negative binomial or zero‐inflated Poisson; these represent two types of overdispersion often encountered in geostatistical count data. We also contrast the correlation structure of one of these Gaussian copula models with that of a hierarchical Poisson model having the same family of marginal distributions, and show that the former is more flexible than the latter in terms of range of feasible correlation, sensitivity to the mean function and modelling of isotropy. An exploratory analysis of a dataset of Japanese beetle larvae counts illustrate some of the findings. All of these investigations show that Gaussian copula models are useful alternatives to hierarchical Poisson models, specially for geostatistical count data that display substantial correlation and small overdispersion.     

Integer‐valued Trawl Processes: A Class of Stationary Infinitely Divisible Processes

This paper introduces a new continuous‐time framework for modelling serially correlated count and integer‐valued data. The key component in our new model is the class of integer‐valued trawl processes, which are serially correlated, stationary, infinitely divisible processes. We analyse the probabilistic properties of such processes in detail and, in addition, study volatility modulation and multivariate extensions within the new modelling framework. Moreover, we describe how the parameters of a trawl process can be estimated and obtain promising estimation results in our simulation study. Finally, we apply our new modelling framework to high‐frequency financial data.     

Flexible Bivariate INAR(1) Processes Using Copulas

Multivariate count time series data occur in many different disciplines. The class of INteger-valued AutoRegressive (INAR) processes has the great advantage to consider explicitly both the discreteness and autocorrelation characterizing this type of data. Moreover, extensions of the simple INAR(1) model to the multi-dimensional space make it possible to model more than one series simultaneously. However, existing models do not offer great flexibility for dependence modelling, allowing only for positive correlation. In this work, we consider a bivariate INAR(1) (BINAR(1)) process where cross-correlation is introduced through the use of copulas for the specification of the joint distribution of the innovations. We mainly emphasize on the parametric case that arises under the assumption of Poisson marginals. Other marginal distributions are also considered. A short application on a bivariate financial count series illustrates the model.     

Estimation and inference for exponential smooth transition nonlinear volatility models

A family of threshold nonlinear generalised autoregressive conditionally heteroscedastic models is considered, that allows smooth transitions between regimes, capturing size asymmetry via an exponential smooth transition function. A Bayesian approach is taken and an efficient adaptive sampling scheme is employed for inference, including a novel extension to a recently proposed prior for the smoothing parameter that solves a likelihood identification problem. A simulation study illustrates that the sampling scheme performs well, with the chosen prior kept close to uninformative, while successfully ensuring identification of model parameters and accurate inference for the smoothing parameter. An empirical study confirms the potential suitability of the model, highlighting the presence of both mean and volatility (size) asymmetry; while the model is favoured over modern, popular model competitors, including those with sign asymmetry, via the deviance information criterion.     

Finite mixtures of multivariate Poisson distributions with application

In the present paper we examine finite mixtures of multivariate Poisson distributions as an alternative class of models for multivariate count data. The proposed models allow for both overdispersion in the marginal distributions and negative correlation, while they are computationally tractable using standard ideas from finite mixture modelling. An EM type algorithm for maximum likelihood (ML) estimation of the parameters is developed. The identifiability of this class of mixtures is proved. Properties of ML estimators are derived. A real data application concerning model based clustering for multivariate count data related to different types of crime is presented to illustrate the practical potential of the proposed class of models.     

The effect of parameter estimation on the performance of one-sided Shewhart control charts for zero-inflated processes

ABSTRACT

Zero-inflated probability models are used to model count data that have an excessive number of zeros. Shewhart-type control charts have been proposed for the monitoring of zero-inflated processes. Usually their performance is evaluated under the assumption of known process parameters. However, in practice, their values are rarely known and they have to be estimated from an in-control historical Phase I sample. In the present paper, we investigate the performance of Shewhart-type control charts for zero-inflated processes with estimated parameters and propose practical guidelines for the statistical design of the examined charts, when the size of the preliminary sample is predetermined.     

Multilevel multivariate modelling of childhood growth, numbers of growth measurements and adult characteristics

Summary.  A general latent normal model for multilevel data with mixtures of response types is extended in the case of ordered responses to deal with variates having a large number of categories and including count data. An example is analysed by using repeated measures data on child growth and adult measures of body mass index and glucose. Applications are described that are concerned with the flexible prediction of adult measurements from collections of growth measurements and for studying the relationship between the number of measurement occasions and growth trajectories.     

Parameter estimation for binomial AR(1) models with applications in finance and industry

Methods for analyzing and modeling count data time series are used in various fields of practice, and they are particularly relevant for applications in finance and economy. We consider the binomial autoregressive (AR(1)) model for count data processes with a first-order AR dependence structure and a binomial marginal distribution. We present four approaches for estimating its model parameters based on given time series data, and we derive expressions for the asymptotic distribution of these estimators. Then we investigate the finite-sample performance of the estimators and of the respective asymptotic approximations in a simulation study, including a discussion of the 2-block jackknife. We illustrate our methods and findings with a real-data example about transactions at the Korea stock market. We conclude with an application of our results for obtaining reliable estimates for process capability indices.     

Modeling spatial overdispersion via the generalized Waring point process

Modeling spatial overdispersion requires point process models with finite‐dimensional distributions that are overdisperse relative to the Poisson distribution. Fitting such models usually heavily relies on the properties of stationarity, ergodicity, and orderliness. In addition, although processes based on negative binomial finite‐dimensional distributions have been widely considered, they typically fail to simultaneously satisfy the three required properties for fitting. Indeed, it has been conjectured by Diggle and Milne that no negative binomial model can satisfy all three properties. In light of this, we change perspective and construct a new process based on a different overdisperse count model, namely, the generalized Waring (GW) distribution. While comparably tractable and flexible to negative binomial processes, the GW process is shown to possess all required properties and additionally span the negative binomial and Poisson processes as limiting cases. In this sense, the GW process provides an approximate resolution to the conundrum highlighted by Diggle and Milne.