Geographically weighted bivariate zero inflated generalized Poisson regression model and its application.
Ontology highlight
ABSTRACT:
SUBMITTER: Purhadi
PROVIDER: S-EPMC8319482 | biostudies-literature |
REPOSITORIES: biostudies-literature
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: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:Geographically weighted regression (GWR) models handle geographical dependence through a spatially varying coefficient model and have been widely used in applied science, but its general Bayesian extension is unclear because it involves a weighted log-likelihood which does not imply a probability distribution on data. We present a Bayesian GWR model and show that its essence is dealing with partial misspecification of the model. Current modularized Bayesian inference models accommodate partial misspecification from a single component of the model. We extend these models to handle partial misspecification in more than one component of the model, as required for our Bayesian GWR model. Information from the various spatial locations is manipulated via a geographically weighted kernel and the optimal manipulation is chosen according to a Kullback–Leibler (KL) divergence. We justify the model via an information risk minimization approach and show the consistency of the proposed estimator in terms of a geographically weighted KL divergence.
| S-EPMC7614111 | biostudies-literature
Project description:There is heightened interest in using high-throughput sequencing technologies to quantify abundances of microbial taxa and linking the abundance to human diseases and traits. Proper modeling of multivariate taxon counts is essential to the power of detecting this association. Existing models are limited in handling excessive zero observations in taxon counts and in flexibly accommodating complex correlation structures and dispersion patterns among taxa. In this article, we develop a new probability distribution, zero-inflated generalized Dirichlet multinomial (ZIGDM), that overcomes these limitations in modeling multivariate taxon counts. Based on this distribution, we propose a ZIGDM regression model to link microbial abundances to covariates (e.g. disease status) and develop a fast expectation-maximization algorithm to efficiently estimate parameters in the model. The derived tests enable us to reveal rich patterns of variation in microbial compositions including differential mean and dispersion. The advantages of the proposed methods are demonstrated through simulation studies and an analysis of a gut microbiome dataset.
| S-EPMC7410344 | 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:Geographically Weighted Regression (GWR) is one of the local statistical models that can capture the effects of spatial heterogeneity. This model can be used for both univariate and multivariate responses. However, it should be noted that GWR models require the assumption of error normality. To overcome this problem, we propose a GWR model for generalized gamma distributed responses that can capture the phenomenon of some special continuous distributions. The proposed model is known as Geographically Weighted Multivariate Generalized Gamma Regression (GWMGGR). Parameter estimation is performed using the Maximum Likelihood Estimation (MLE) method optimized with the Bernt-Hall-Hall-Haussman (BHHH) algorithm. To determine the significance of the spatial heterogeneity effect, a hypothesis test was conducted using the Maximum Likelihood Ratio Test (MLRT) approach. We made a spatial cluster based on the estimated model parameters for each response using the k-means clustering method to interpret the obtained results. Some highlights of the proposed method are:•A new model for GWR with multivariate generalized gamma distributed responses to overcome the assumption of normally distributed errors.•Goodness of fit test to test the spatial effects in GWMGGR model.•Spatial clustering of districts/cities in Central Java based on three dimensions of educational indicators.
| S-EPMC11372798 | 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:Detecting the temporally and spatially varying correlations is important to understand the biological and disease systems. Here we proposed a geographically weighted temporally correlated logistic regression (GWTCLR) model to identify such dynamic correlation of predictors on binomial outcome data, by incorporating spatial and temporal information for joint inference. The local likelihood method is adopted to estimate the spatial relationship, while the smoothing method is employed to estimate the temporal variation. We present the construction and implementation of GWTCLR and the study of the asymptotic properties of the proposed estimator. Simulation studies were conducted to evaluate the robustness of the proposed model. GWTCLR was applied on real epidemiologic data to study the climatic determinants of human seasonal influenza epidemics. Our method obtained results largely consistent with previous studies but also revealed certain spatial and temporal varying patterns that were unobservable by previous models and methods.
| S-EPMC5780421 | biostudies-literature
Project description: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 covariates. The method, which is applicable to other missing data problems, is illustrated and compared with complete case analysis by using simulations and dental data on the caries preventive effects of a school-based fluoride mouthrinse program.
| S-EPMC6453121 | biostudies-literature
Project description:The zero-inflated negative binomial distribution has been widely used for count data analyses in various biomedical settings due to its capacity of modeling excess zeros and overdispersion. When there are correlated count variables, a bivariate model is essential for understanding their full distributional features. Examples include measuring correlation of two genes in sparse single-cell RNA sequencing data and modeling dental caries count indices on two different tooth surface types. For these purposes, we develop a richly parametrized bivariate zero-inflated negative binomial model that has a simple latent variable framework and eight free parameters with intuitive interpretations. In the scRNA-seq data example, the correlation is estimated after adjusting for the effects of dropout events represented by excess zeros. In the dental caries data, we analyze how the treatment with Xylitol lozenges affects the marginal mean and other patterns of response manifested in the two dental caries traits. An R package "bzinb" is available on Comprehensive R Archive Network.
| S-EPMC10500952 | biostudies-literature