A spatial beta-binomial model for clustered count data on dental caries.

The Poisson regression model is widely used for count data. This model assumes equidispersion. In practice, equidispersion is seldom reflected in data. However, in real-life data, the variance usually exceeds the mean. This situation is known as overdispersion. Negative binomial distribution and other Poisson mix models are often used to model overdispersion count data. Another extension of the negative binomial distribution in another model for count data is the univariate generalized Waring. In addition, the model developed by Famoye can be used in the analysis of count data. When the count data contains a large number of zeros, it is necessary to use zero-inflated models. In this study, different generalized regression models are emphasized for the analysis of excessive zeros count data. For this purpose, a real data set was analysed with the generalized Poisson model, generalized negative binomial model, generalized negative binomial Famoye, generalized Waring model, and the foregoing zero-inflated models. Log-likelihood, Akaike information criterion, Bayes information criterion, Vuong statistics were used for model comparisons.

BackgroundThe performance of the Beta Binomial (BB) model is compared with several existing models for mapping the EORTC QLQ-C30 (QLQ-C30) on to the EQ-5D-3L using data from lung cancer trials.MethodsData from 2 separate non small cell lung cancer clinical trials (TOPICAL and SOCCAR) are used to develop and validate the BB model. Comparisons with Linear, TOBIT, Quantile, Quadratic and CLAD models are carried out. The mean prediction error, R2, proportion predicted outside the valid range, clinical interpretation of coefficients, model fit and estimation of Quality Adjusted Life Years (QALY) are reported and compared. Monte-Carlo simulation is also used.ResultsThe Beta-Binomial regression model performed ‘best’ among all models. For TOPICAL and SOCCAR trials, respectively, residual mean square error (RMSE) was 0.09 and 0.11; R2 was 0.75 and 0.71; observed vs. predicted means were 0.612 vs. 0.608 and 0.750 vs. 0.749. Mean difference in QALY’s (observed vs. predicted) were 0.051 vs. 0.053 and 0.164 vs. 0.162 for TOPICAL and SOCCAR respectively. Models tested on independent data show simulated 95% confidence from the BB model containing the observed mean more often (77% and 59% for TOPICAL and SOCCAR respectively) compared to the other models. All algorithms over-predict at poorer health states but the BB model was relatively better, particularly for the SOCCAR data.ConclusionThe BB model may offer superior predictive properties amongst mapping algorithms considered and may be more useful when predicting EQ-5D-3L at poorer health states. We recommend the algorithm derived from the TOPICAL data due to better predictive properties and less uncertainty.Electronic supplementary materialThe online version of this article (doi:10.1186/s12955-014-0163-7) contains supplementary material, which is available to authorized users.

This paper discusses the formation of one temporal spatial model, the Generalized Space Time Autoregressive (GSTAR) model, if the data of model is count data. The development of the GSTAR model is an update or problem completion of the count data which tends to be stationary and non-normal / Gaussian data, because GSTAR model is assumed normal distributed and stationarity. GSTAR modelling for count data refers to the Time series model for count data, which are the Generalized Autoregressive Moving Average (GARMA) model and modeling Count Time Series which has the Generalized Linear Model (GLM) concept. The model formed is called the Poisson GSTAR model

Paradigms for substance abuse cue-reactivity research involve pharmacological or stressful stimulation designed to elicit stress and craving responses in cocaine-dependent subjects. It is unclear as to whether stress induced from participation in such studies increases drug-seeking behavior. We propose a 2-state Hidden Markov model to model the number of cocaine abuses per week before and after participation in a stress-and cue-reactivity study. The hypothesized latent state corresponds to 'high' or 'low' use. To account for a preponderance of zeros, we assume a zero-inflated Poisson model for the count data. Transition probabilities depend on the prior week's state, fixed demographic variables, and time-varying covariates. We adopt a Bayesian approach to model fitting, and use the conditional predictive ordinate statistic to demonstrate that the zero-inflated Poisson hidden Markov model outperforms other models for longitudinal count data.

Microorganisms play critical roles in human health and disease. They live in diverse communities in which they interact synergistically or antagonistically. Thus for estimating microbial associations with clinical covariates, such as treatment effects, joint (multivariate) statistical models are preferred. Multivariate models allow one to estimate and exploit complex interdependencies among multiple taxa, yielding more powerful tests of exposure or treatment effects than application of taxon-specific univariate analyses. Analysis of microbial count data also requires special attention because data commonly exhibit zero inflation, i.e., more zeros than expected from a standard count distribution. To meet these needs, we developed a Bayesian variable selection model for multivariate count data with excess zeros that incorporates information on the covariance structure of the outcomes (counts for multiple taxa), while estimating associations with the mean levels of these outcomes. Though there has been much work on zero-inflated models for longitudinal data, little attention has been given to high-dimensional multivariate zero-inflated data modeled via a general correlation structure. Through simulation, we compared performance of the proposed method to that of existing univariate approaches, for both the binary ("excess zero") and count parts of the model. When outcomes were correlated the proposed variable selection method maintained type I error while boosting the ability to identify true associations in the binary component of the model. For the count part of the model, in some scenarios the univariate method had higher power than the multivariate approach. This higher power was at a cost of a highly inflated false discovery rate not observed with the proposed multivariate method. We applied the approach to oral microbiome data from the Pediatric HIV/AIDS Cohort Oral Health Study and identified five (of 44) species associated with HIV infection.

The Poisson regression model is often used as a first model for count data with covariates. Since this model is a GLM with canonical link, regression parameters can be easily fitted using standard software. However the model requires equidispersion, which might not be valid for the data set under consideration. There have been many models proposed in the literature to allow for overdispersion. One such model is the negative binomial regression model. In addition, score tests have been commonly used to detect overdispersion in the data. However these tests do not allow to quantify the effects of overdispersion. In this paper we propose easily interpretable discrepancy measures which allow to quantify the overdispersion effects when comparing a negative binomial regression to Poisson regression. We propose asymptotic $\alpha$-level tests for testing the size of overdispersion effects in terms of the developed discrepancy measures. A graphical display of p-values curves can then be used to allow for an exact quantification of the overdispersion effects. This can lead to a validation of the Poisson regression or a discrimination of the Poisson regression with respect to the negative binomial regression. The proposed asymptotic tests are investigated in small samples using simulation and applied to two examples.

Clustering with variable selection is a challenging yet critical task for modern small-n-large-p data. Existing methods based on sparse Gaussian mixture models or sparse $K$-means provide solutions to continuous data. With the prevalence of RNA-seq technology and lack of count data modeling for clustering, the current practice is to normalize count expression data into continuous measures and apply existing models with a Gaussian assumption. In this article, we develop a negative binomial mixture model with lasso or fused lasso gene regularization to cluster samples (small $n$) with high-dimensional gene features (large $p$). A modified EM algorithm and Bayesian information criterion are used for inference and determining tuning parameters. The method is compared with existing methods using extensive simulations and two real transcriptomic applications in rat brain and breast cancer studies. The result shows the superior performance of the proposed count data model in clustering accuracy, feature selection, and biological interpretation in pathways.

Panel count data occur in many fields including clinical, demographical and industrial studies and an extensive literature has been established for their regression analysis. However, most of the existing methods apply only to the situations where both covariates and their effects are constant or one of them may be time-dependent. In the first part of this dissertation, we consider a situation where both covariates and their effects may be time-dependent and an estimating equation-based approach is developed for estimating those time-varying effects. In the method, B-spline functions are employed to approximate time-dependent coefficients and the asymptotic properties of the proposed estimators are established. To assess the performance of the proposed approach, an extensive simulation study is conducted and suggests that it works well in practical situations. An application to the China Health and Nutrition Survey (CHNS) study is provided. In practice, there could exist more than one type of event of interest, such as two types of tumor recurrence, leading to multivariate panel count data. The second part of this dissertation considers marginal mean model for multivariate panel count data with time-dependent coefficient and covariate effects, which has limited previous research. Based on the conditional estimating equation method developed for time-dependent covariates, we approximate the coefficients by B-splines, hence allow both coefficients and covariates to be time-dependent. Simulation studies show that the proposed estimation procedures work well for practical situations. The methodology is again applied to the China Health and Nutrition Survey (CHNS) study. When we consider time-varying covariates and coefficients effects, most of the previous study focused on the proportional mean model because the likelihood function under the rate model involves intractable integration. However, the rate model is more realistic and efficient. Hence, in the third part of this dissertation, we propose a semi-parametric MLE method under the rate model for panel count data with time-dependent covariates and time-varying effects. B-spline functions are employed again to approximate time-dependent coefficients and an efficient Expectation-Maximization-type algorithm is developed to overcome the computational difficulty. The resulting estimators are shown to be consistent and asymptotically efficient. Monte Carlo simulation studies demonstrate that the proposed method enjoys desirable finite-sample properties. An application to The Young Women's Project (YWP) is provided.

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.

This article applies growth curve models to longitudinal count data characterized by an excess of zero counts. We discuss a zero-inflated Poisson regression model for longitudinal data in which the impact of covariates on the initial counts and the rate of change in counts over time is the focus of inference. Basic growth curve models using a single outcome are described, as well as a model in which two linked outcomes constitute a dual-trajectory growth process. This model is applied to assess the impact of changes in financial stress on longitudinal change in smoking.

Piping Plovers (Charadrius melodus) undertake comparatively short migrations for a shorebird and were not previously thought to congregate in large numbers during migration. Superpopulation size (individuals occurring at the study site during the study period) and stopover duration were estimated for migratory Piping Plovers on South Point, Ocracoke Island, North Carolina, USA, from 3 July–28 October 2016 by integrating a Jolly-Seber model and a binomial model of resighting and count data. We estimated 569 Piping Plovers (95% CI: 502–651) used South Point during fall migration in 2016, approximately 10% of the global population and approximately 15% of the Atlantic Coast population. Stopover duration differed between Piping Plovers that bred on or near our study site (Southern USA, 46 days, SD = 1.7) compared to individuals that bred in the northern area of the breeding range (Atlantic Canada, 41 days, SD = 2.0; New England States, USA, 44 days, SD = 1.8) and the central area of the breeding range (New York and New Jersey, USA, 26 days, SD = 1.4). South Point may be unique in that no other areas are known to host similar numbers of Piping Plovers during fall migration.

ABSTRACTPanel count data frequently arise in various applications such as medical research, social sciences and so on. In this paper, a partially varying coefficient model of the panel count data with informative observation times is developed to accommodate the nonlinear interact effects between covariates. For statistical inference of the unknown parameters, quantile regression approaches are proposed, in which the baseline function and the varying coefficients are approximated by B-spline functions. Moreover, asymptotic properties for the estimators are established. Some numerical studies are performed to confirm and evaluate the finite-sample behaviours of the proposed approaches. Finally, the proposed model is applied to the bladder cancer tumour data as an application.

The identification of accident hot spots is a central task of road safety management. Bayesian count data models have emerged as the workhorse method for producing probabilistic rankings of hazardous sites in road networks. Typically, these methods assume simple linear link function specifications, which, however, limit the predictive power of a model. Furthermore, extensive specification searches are precluded by complex model structures arising from the need to account for unobserved heterogeneity and spatial correlations. Modern machine learning (ML) methods offer ways to automate the specification of the link function. However, these methods do not capture estimation uncertainty, and it is also difficult to incorporate spatial correlations. In light of these gaps in the literature, this paper proposes a new spatial negative binomial model which uses Bayesian additive regression trees to endogenously select the specification of the link function. Posterior inference in the proposed model is made feasible with the help of the Pólya-Gamma data augmentation technique. We test the performance of this new model on a crash count data set from a metropolitan highway network. The empirical results show that the proposed model performs at least as well as a baseline spatial count data model with random parameters in terms of goodness of fit and site ranking ability.

In medical and health studies, heterogeneities in clustered count data have been traditionally modeled by positive random effects in Poisson mixed models; however, excessive zeros often occur in clustered medical and health count data. In this paper, we consider a three-level random effects zero-inflated Poisson model for health-care utilization data where data are clustered by both subjects and families. To accommodate zero and positive components in the count response compatibly, we model the subject level random effects by a compound Poisson distribution. Our model displays a variance components decomposition which clearly reflects the hierarchical structure of clustered data. A quasi-likelihood approach has been developed in the estimation of our model. We illustrate the method with analysis of the health-care utilization data. The performance of our method is also evaluated through simulation studies.

Maximum likelihood estimation of an extended latent Markov model for clustered binary panel data