Marginalized zero-inflated Poisson models with missing covariates.
Ontology highlight
ABSTRACT:
SUBMITTER: Benecha HK
PROVIDER: S-EPMC6453121 | biostudies-literature | 2018 Jul
REPOSITORIES: biostudies-literature
Publications
Marginalized zero-inflated Poisson models with missing covariates.
Benecha Habtamu K HK Preisser John S JS Divaris Kimon K Herring Amy H AH Das Kalyan K
Biometrical journal. Biometrische Zeitschrift 20180511 4
Unlike zero-inflated Poisson regression, marginalized zero-inflated Poisson (MZIP) models for counts with excess zeros provide estimates with direct interpretations for the overall effects of covariates on the marginal mean. In the presence of missing covariates, MZIP and many other count data models are ordinarily fitted using complete case analysis methods due to lack of appropriate statistical methods and software. This article presents an estimation method for MZIP models with missing covari ...[more]
Similar Datasets
Project description:The zero-inflated Poisson (ZIP) regression model is often employed in public health research to examine the relationships between exposures of interest and a count outcome exhibiting many zeros, in excess of the amount expected under sampling from a Poisson distribution. The regression coefficients of the ZIP model have latent class interpretations, which correspond to a susceptible subpopulation at risk for the condition with counts generated from a Poisson distribution and a non-susceptible subpopulation that provides the extra or excess zeros. The ZIP model parameters, however, are not well suited for inference targeted at marginal means, specifically, in quantifying the effect of an explanatory variable in the overall mixture population. We develop a marginalized ZIP model approach for independent responses to model the population mean count directly, allowing straightforward inference for overall exposure effects and empirical robust variance estimation for overall log-incidence density ratios. Through simulation studies, the performance of maximum likelihood estimation of the marginalized ZIP model is assessed and compared with other methods of estimating overall exposure effects. The marginalized ZIP model is applied to a recent study of a motivational interviewing-based safer sex counseling intervention, designed to reduce unprotected sexual act counts.
| S-EPMC4227945 | biostudies-literature
Project description:BackgroundTsetse flies are the major vectors of human trypanosomiasis of the form Trypanosoma brucei rhodesiense and T.b.gambiense. They are widely spread across the sub-Saharan Africa and rendering a lot of challenges to both human and animal health. This stresses effective agricultural production and productivity in Africa. Delimiting the extent and magnitude of tsetse coverage has been a challenge over decades due to limited resources and unsatisfactory technology. In a bid to overcome these limitations, this study attempted to explore modelling skills that can be applied to spatially estimate tsetse abundance in the country using limited tsetse data and a set of remote-sensed environmental variables.MethodologyEntomological data for the period 2008-2018 as used in the model were obtained from various sources and systematically assembled using a structured protocol. Data harmonisation for the purposes of responsiveness and matching was carried out. The key tool for tsetse trapping was itemized as pyramidal trap in many instances and biconical trap in others. Based on the spatially explicit assembled data, we ran two regression models; standard Poisson and Zero-Inflated Poisson (ZIP), to explore the associations between tsetse abundance in Uganda and several environmental and climatic covariates. The covariate data were constituted largely by satellite sensor data in form of meteorological and vegetation surrogates in association with elevation and land cover data. We finally used the Zero-Inflated Poisson (ZIP) regression model to predict tsetse abundance due to its superiority over the standard Poisson after model fitting and testing using the Vuong Non-Nested statistic.ResultsA total of 1,187 tsetse sampling points were identified and considered as representative for the country. The model results indicated the significance and level of responsiveness of each covariate in influencing tsetse abundance across the study area. Woodland vegetation, elevation, temperature, rainfall, and dry season normalised difference vegetation index (NDVI) were important in determining tsetse abundance and spatial distribution at varied scales. The resultant prediction map shows scaled tsetse abundance with estimated fitted numbers ranging from 0 to 59 flies per trap per day (FTD). Tsetse abundance was found to be largest at low elevations, in areas of high vegetative activity, in game parks, forests and shrubs during the dry season. There was very limited responsiveness of selected predictors to tsetse abundance during the wet season, matching the known fact that tsetse disperse most significantly during wet season.ConclusionsA methodology was advanced to enable compilation of entomological data for 10 years, which supported the generation of tsetse abundance maps for Uganda through modelling. Our findings indicate the spatial distribution of the G. f. fuscipes as; low 0-5 FTD (48%), medium 5.1-35 FTD (18%) and high 35.1-60 FTD (34%) grounded on seasonality. This approach, amidst entomological data shortages due to limited resources and absence of expertise, can be adopted to enable mapping of the vector to provide better decision support towards designing and implementing targeted tsetse and tsetse-transmitted African trypanosomiasis control strategies.
| S-EPMC8648107 | biostudies-literature
Project description:Dimension reduction of high-dimensional microbiome data facilitates subsequent analysis such as regression and clustering. Most existing reduction methods cannot fully accommodate the special features of the data such as count-valued and excessive zero reads. We propose a zero-inflated Poisson factor analysis model in this paper. The model assumes that microbiome read counts follow zero-inflated Poisson distributions with library size as offset and Poisson rates negatively related to the inflated zero occurrences. The latent parameters of the model form a low-rank matrix consisting of interpretable loadings and low-dimensional scores that can be used for further analyses. We develop an efficient and robust expectation-maximization algorithm for parameter estimation. We demonstrate the efficacy of the proposed method using comprehensive simulation studies. The application to the Oral Infections, Glucose Intolerance, and Insulin Resistance Study provides valuable insights into the relation between subgingival microbiome and periodontal disease.
| S-EPMC7594114 | biostudies-literature
Project description:Community water fluoridation is an important public health measure to prevent dental caries, but it continues to be somewhat controversial. The Iowa Fluoride Study (IFS) is a longitudinal study on a cohort of Iowa children that began in 1991. The main purposes of this study (http://www.dentistry.uiowa.edu/preventive-fluoride-study) were to quantify fluoride exposures from both dietary and nondietary sources and to associate longitudinal fluoride exposures with dental fluorosis (spots on teeth) and dental caries (cavities). We analyze a subset of the IFS data by a marginal regression model with a zero-inflated version of the Conway-Maxwell-Poisson distribution for count data exhibiting excessive zeros and a wide range of dispersion patterns. In general, we introduce two estimation methods for fitting a ZICMP marginal regression model. Finite sample behaviors of the estimators and the resulting confidence intervals are studied using extensive simulation studies. We apply our methodologies to the dental caries data. Our novel modeling incorporating zero inflation, clustering, and overdispersion sheds some new light on the effect of community water fluoridation and other factors. We also include a second application of our methodology to a genomic (next-generation sequencing) dataset that exhibits underdispersion.
| S-EPMC4948193 | biostudies-literature
Project description:This study discusses the development of Zero Inflated Generalized Poisson Regression (ZIGPR) with two response variables, that is Bivariate ZIGPR (BZIGPR). The extension of the ZIGPR model by considering spatial factor called Geographically Weighted Zero Inflated Generalized Poisson Regression (GWBZIGPR). The GWBZIGPR produces a local parameter estimator for each location of observation. The parameter estimation using the Maximum Likelihood Estimation (MLE) method obtained an equation that did not closed-form so that the numerical iteration of Berndt Hall Hall Hausman (BHHH) is used. The data used in this study are the number of pregnant maternal mortality and postpartum maternal mortality data in 91 sub-districts in Pekalongan Residency, Central Java Province. The results showed that the Akaike Information Criterion Corrected (AICc) value in the GWBZIGPR model is smaller than BZIGPR, so it means that the GWBZIGPR is better than the BZIGPR for modeling the number of pregnant maternal mortality and postpartum maternal mortality in Pekalongan Residency. The results of this study will assist local governments in anticipating the causes of maternal mortality.
| S-EPMC8319482 | biostudies-literature
Project description:RNA sequencing (RNA-Seq) has been frequently used in genomic studies and has generated a vast amount of data. The RNA-Seq data are composed of two parts: (a) a sequence of nucleotides of the genome; and (b) a corresponding sequence of counts, standing for the number of short reads whose mapped positions start at each position of the genome. One common feature of these count data is that they are typically nonuniform; recent studies have revealed that the nonuniformity is partially owing to a systematic bias resulted from the sequencing preference. Existing works in the literature model the nonuniformity with a single component Poisson linear model that incorporates the effects of the sequencing preference. However, we observe consistently that the short reads mapped to a gene may have a mixture structure and can be zero-inflated. A single component model may not suffice to model the complexity of such data. In this paper, we propose a zero-inflated mixture Poisson linear model for the RNA-Seq count data and derive a fast expectation-maximisation-based algorithm for estimating the unknown parameters. Numerical studies are conducted to illustrate the effectiveness of our method.
| S-EPMC6763381 | biostudies-literature
Project description:Count data often exhibit more zeros than predicted by common count distributions like the Poisson or negative binomial. In recent years, there has been considerable interest in methods for analyzing zero-inflated count data in longitudinal or other correlated data settings. A common approach has been to extend zero-inflated Poisson models to include random effects that account for correlation among observations. However, these models have been shown to have a few drawbacks, including interpretability of regression coefficients and numerical instability of fitting algorithms even when the data arise from the assumed model. To address these issues, we propose a model that parameterizes the marginal associations between the count outcome and the covariates as easily interpretable log relative rates, while including random effects to account for correlation among observations. One of the main advantages of this marginal model is that it allows a basis upon which we can directly compare the performance of standard methods that ignore zero inflation with that of a method that explicitly takes zero inflation into account. We present simulations of these various model formulations in terms of bias and variance estimation. Finally, we apply the proposed approach to analyze toxicological data of the effect of emissions on cardiac arrhythmias.
| S-EPMC5111636 | biostudies-literature
Project description:Zero-inflated regression models have emerged as a popular tool within the parametric framework to characterize count data with excess zeros. Despite their increasing popularity, much of the literature on real applications of these models has centered around the latent class formulation where the mean response of the so-called at-risk or susceptible population and the susceptibility probability are both related to covariates. While this formulation in some instances provides an interesting representation of the data, it often fails to produce easily interpretable covariate effects on the overall mean response. In this article, we propose two approaches that circumvent this limitation. The first approach consists of estimating the effect of covariates on the overall mean from the assumed latent class models, while the second approach formulates a model that directly relates the overall mean to covariates. Our results are illustrated by extensive numerical simulations and an application to an oral health study on low income African-American children, where the overall mean model is used to evaluate the effect of sugar consumption on caries indices.
| S-EPMC4988952 | biostudies-literature
Project description:Recently developed accelerometer devices have been used in large epidemiological studies for continuous and objective monitoring of physical activities. Typically, physical movements are summarized as minutes in light, moderate, and vigorous physical activities in each wearing day. Because of preponderance of zeros, zero-inflated distributions have been used for modeling the daily moderate or higher levels of physical activity. Yet, these models do not fully account for variations in daily physical activity and cannot be extended to model weekly physical activity explicitly, while the weekly physical activity is considered as an indicator for a subject's average level of physical activity. To overcome these limitations, we propose to use a zero-inflated Poisson mixture distribution that can model daily and weekly physical activity in same family of mixture distributions. Under this method, the likelihood of an inactive day and the amount of exercise in an active day are simultaneously modeled by a joint random effects model to incorporate heterogeneity across participants. If needed, the method has the flexibility to include an additional random effect to address extra variations in daily physical activity. Maximum likelihood estimation can be obtained through Gaussian quadrature technique, which is implemented conveniently in an R package GLMMadaptive. Method performances are examined using simulation studies. The method is applied to data from the Hispanic Community Health Study/Study of Latinos to examine the relationship between physical activity and BMI groups and within a participant the difference in physical activity between weekends and weekdays.
| S-EPMC8521567 | biostudies-literature
Project description:Bulk gene expression experiments relied on aggregations of thousands of cells to measure the average expression in an organism. Advances in microfluidic and droplet sequencing now permit expression profiling in single cells. This study of cell-to-cell variation reveals that individual cells lack detectable expression of transcripts that appear abundant on a population level, giving rise to zero-inflated expression patterns. To infer gene co-regulatory networks from such data, we propose a multivariate Hurdle model. It is comprised of a mixture of singular Gaussian distributions. We employ neighborhood selection with the pseudo-likelihood and a group lasso penalty to select and fit undirected graphical models that capture conditional independences between genes. The proposed method is more sensitive than existing approaches in simulations, even under departures from our Hurdle model. The method is applied to data for T follicular helper cells, and a high-dimensional profile of mouse dendritic cells. It infers network structure not revealed by other methods; or in bulk data sets. An R implementation is available at https://github.com/amcdavid/HurdleNormal.
| S-EPMC6684253 | biostudies-literature