 Research
 Open access
 Published:
Socioenvironmental predictors of diabetes incidence disparities in Tanzania mainland: a comparison of regression models for count data
BMC Medical Research Methodology volumeÂ 24, ArticleÂ number:Â 75 (2024)
Abstract
Background
Diabetes is one of the top four noncommunicable diseases that cause death and illness to many people around the world. This study aims to use an efficient count data model to estimate socioenvironmental factors associated with diabetes incidences in Tanzania mainland, addressing lack of evidence on the efficient count data model for estimating factors associated with disease incidences disparities.
Methods
This study analyzed diabetes counts in 184 Tanzania mainland councils collected in 2020. The study applied generalized Poisson, negative binomial, and Poisson count data models and evaluated their adequacy using information criteria and Pearson chisquare values.
Results
The data were overdispersed, as evidenced by the mean and variance values and the positively skewed histograms. The results revealed uneven distribution of diabetes incidence across geographical locations, with northern and urban councils having more cases. Factors like population, GDP, and hospital numbers were associated with diabetes counts. The GP model performed better than NB and Poisson models.
Conclusion
The occurrence of diabetes can be attributed to geographical locations. To address this public health issue, environmental interventions can be implemented. Additionally, the generalized Poisson model is an effective tool for analyzing health information system count data across different population subgroups.
Background
To date, noncommunicable, diseases including diabetes, are still a global health challenge affecting people of all ages; however, elderly people are at higher risk [1, 2]. In 2016, statistics showed that NonCommunicable Diseases (NCDs) were responsible for 80% of all deaths worldwide. The NCD death risk is notably higher in SubSaharan Africa, Central Asia, and Eastern Europe [3]. In Tanzania, as in SubSaharan Africa, there is evidence of a high prevalence of NCD cases, including diabetes [4, 5].
The emergence of NCDs in humans is influenced by a complex combination of various factors, which include environmental conditions, cultural beliefs, selfmanagement, sociodemographic factors, genetics, and biology [1]. These diseases are sometimes referred to as behavioural diseases because, apart from other factors, selfmanagement, which is linked by a personâ€™s behavioural practice in running his/her daily life can increase oneâ€™s likelihood of developing NCDs [1, 6]. Cultural norms and values can also influence human behaviour, resulting in regional and national variations in the prevalence of NCDs [7]. This study aims to investigate the impact of socioenvironmental factors, which are the councilâ€™s zone and residence, along with other factors on the total diabetes incidences in the council. The link between environment and human behaviours is well explained in some behavioural theories, including the reciprocal deterministic concept of social cognitive theory [8].
Currently, many scenarios in public health and official statistics include count data. Count data include specific disease cases reported in a particular geographical unit, the total number of fatalities occurring within a given timeframe, etc. The Poisson model is a wellknown method for modeling count data and has been applied in many situations [9,10,11,12,13,14]. However, it assumes that the subject occurs randomly and at a constant rate, resulting in equality mean and variance, which is often unrealistic in reallife situations. When data exhibit overdispersion, the negative binomial (NB) model is often used as an alternative to the Poisson model [11, 15, 16]. Occasionally, under dispersion also occurs among count data, especially for rare events. To tackle this issue, researchers have developed new models that can model count data that exhibits over, under, or equal dispersion. These models were obtained as a result of generalization or mixing with the Poisson model. Examples of these models are the Generalized Poisson (GP) [17], the Weighted Poisson, the ConwayMaxwellPoisson (CMP), the HyperPoisson (HP) [18], Extended Biparametric Waring (EBW) [19], and the Complex Triparametric Pearson (CTP) [20].
Many of the distributions mentioned above have complex functional forms, which can lead to significant computational challenges and make them difficult to use. For this reason, the GP model was selected for this study. This model has a welldefined functional form and allows easy parametric estimation [21,22,23,24]. The GP model is the best option to be used in health and behavioural studies, for many reasons, including the nonuniformity of the population being studied, where individuals tend to cluster or aggregate within a particular combination with similar characteristics; dependence among observations due to environmental factors, where there are high incidences of diabetes cases in the same geographic area due to similarity in sociocultural factors; which causes unequal dispersion which happens in the data [15, 21].
Numerous studies on NCDs, including diabetes, have been conducted [4, 5, 25,26,27,28,29]. However, none have utilized the GP model or quantified socioenvironmental factors (such as zone and council residence) in NCD occurrences in mainland Tanzania. This study aims to establish a model that can be adopted in modeling NCDs count incidences associated with socioenvironmental and other risk factors. Hence, it emphasizes environmentalbased approaches to eradicating NCDs in Tanzania, and the model can also be adopted in similar scenarios. Many research articles elaborate on the application of GP regression in modeling overdispersed data [24, 30,31,32,33,34,35]. However, the articles do not describe or quantify how overestimation of the standard error occurs when using standard Poisson in modeling overdispersed data as the current does.
Methods
Design and settings
This study utilized crosssectional reseach design. Secondary data collected by the District Health Information System (DHIS2) and the National Bureau of Statistics (NBS) in 2020 were used for analysis. The response variables represent the number of patients diagnosed with diabetes mellitus admitted to all health facilities within a council, except for regional referral and zonal hospitals. Information collected in 2020 and across 184 councils in Tanzania mainland is used in this study.
Models description
Generalized linear models (GLMs) extend linear models (LMs) when the response variable is not normally distributed, allowing for the representation of nonnormal response variables. In GLMs, the distribution of the response variable can be counted, categorical, discrete, ordinal, and many others as long as it belongs to the exponential family of distributions. This family has several wellknown distributions including the Poisson distribution and its generalization, the binomial distribution, Gamma distribution, and many others. GLMs can be described using the following equation:
And is mainly characterized by its three components: (1) A random component, which describes the outcome variable \({Y}_{i}\) of the \({i}^{th}\) observation by its probability density function. (2) A linear component \({X}_{i}^{T}\beta\), where \({X}_{i}^{T}\) is the vector of predictors and \(\beta\) is a column vector of model coefficients. (3) Differentiable link function \(g\left(\mu \right)\), which relates the mean of the response variable and the linear function of the predictor variables [9, 36].
The Poisson process is often used to explain the variations in count data compared to a predicted average [9]. However, this model has certain assumptions, including that the data must be equally distributed and that the mean must always equal the variance. Poisson regression is a wellknown model for modeling the means of \(n\)Â nonnegative count response variables \({y}_{1},{y}_{2},\dots ,{y}_{2}\). Let \({Y}_{i} i=\text{0,1},\dots .\), be the response variable which represents the number of diabetic patients admitted to specific council in 2020 and \({X}_{i}^{{\prime }}={X}_{1i},\dots ,{X}_{ki}\) represents a kdimension vector of linear predictors associated with the response variable \(Y\). A Poisson regression of the response variable given predictors is written as:
For Poisson distribution, we haveÂ \(E\left({Y}_{i}={y}_{i}/{X}_{i}={x}_{i}\right)={\mu }_{i}=Var\left({Y}_{i}={y}_{i}/{X}_{i}={x}_{i}\right)\)
The logarithm of the likelihood of the equation above can be written as:
By substituting \({\mu }_{i}={e}^{{{x}_{i}}^{{\prime }}\beta }\), we obtain the logarithm of the likelihood function in terms of \(\beta^{\prime }s\) which can be written as:
Maximum likelihood estimates of Î²â€™s can be obtained by differentiating the logarithm of the likelihood equation with respect to Î²â€™s and setting the results equal to zero.
Thus, the Poisson regression model of the mean parameter \({\mu }_{i}\) is written as
As data are collected from councils across various geographical locations, including areas with differing behavioural patterns, there is a high probability of unequal dispersion in the data. This suggests that the data may have under or overdispersion. If Poisson regression is used to model these data, it could lead to incorrect conclusions because the standard error may be overestimated [15, 37]. The negative binomial model, also known as the Poisson Gamma mixture, is considered a better alternative to Poisson regression when dealing with overdispersed count data. The modelâ€™s mean and variance have a quadratic relationship, resulting in its being named NB2 [11].
The NB model was formulated as an extension of the Poisson model by considering the idea that the modeled outcomes cannot happen at a constant rate, leading to heterogeneity in the outcomes. The extended model can be formulated as follows:
A negative binomial distribution is generated using a series of Bernoulli trials with a constant success probability p. Let Y be the number of attempts that failed before the \({k}^{th}\) success \((k>0)\), then, Y follows a negative Binomial distribution with probability mass function (pmf) written as follows:
The mean and variance of Y are \(\frac{pk}{\left(1p\right)}\) and \(\frac{pk}{{\left(1p\right)}^{2}}\) respectively. In the negative binomial regression model the interest in modeling the mean of the outcome variable Y with its realization \({y}_{1},{y}_{2},\dots ,{y}_{2}\), and \({X}_{i}^{{\prime }}={X}_{1i},\dots ,{X}_{ki}\) denotes the matrix of predictors. Parametrization of Eq.Â (2) above in terms of \(\mu\) and dispersion parameter \(\alpha\) yield NB regression model as described below:
Let \(p=\frac{\alpha }{\alpha +\mu }\), where \(\alpha =k\), furthermore it is known that \(\left(\begin{array}{c}\alpha +{y}_{i}1\\ {y}_{i}\end{array}\right)=\frac{\left(\alpha +{y}_{i}1\right)!}{{y}_{i}!\left(\alpha 1\right)!}=\frac{{\Gamma }\left(\alpha +{y}_{i}\right)}{{y}_{i}!{\Gamma }\left(\alpha \right)}\), then the pmf of Y in Eq.Â (2) can be written as:
where \({\Gamma }\) represents the gamma function and \(\alpha\) is a dispersion index that has been modified to take positive values only. The NB can also be obtained by using the Poisson mixture gamma formula. Then, \(Y\sim NB\left(\mu,\alpha\right)\) and the mean and variance of Y are \({\mu }_{i}\) and \({\mu }_{i}+\frac{{{\mu }_{i}}^{2}}{\alpha }\) respectively. When \(\alpha \to \infty\), the mean and variance of Y tend to be equal, which implies that the Poisson model is a special case of the negative binomial model [9, 11, 15, 36, 38].
The likelihood of Eq.Â (3) is proportional to:
It is known that:
It follows that:
and the loglikelihood is given by:
It is known that
Estimates of the regression coefficients \(\beta^{\prime }s\) and dispersion index \(\alpha\) are obtained by substituting into the above equation and differentiating it with respect to \(\beta^{\prime }s\) and \(\alpha\) and setting the result equal to zero.
Then, the negative binomial regression model can be written as:
NB model cannot be used to model equal and underdispersed data. The finding reveals that the NB model faces convergence issues if inappropriately used to model count data, which does not exhibit overdispersion [38].
Many articles use the latest count model generalizations; however, the GP model remains beneficial and userfriendly [39]. This model can model stochastic processes with count data that have equal, under, or overdispersion. Moreover, estimating the parameters of this model is simple compared to other generalized models. Due to the reasons mentioned above, this study employs the model introduced by Consul and Jain [17, 40]. Let \({Y}_{i}\) represent diabetes incidences for inpatient recorded in a certain council for 2020. Then, \({Y}_{i}\) represents the response variable having response values \(y_1,y_2,\dots\;\dots,y_n\) associated with several explanatory variables. Then, \({Y}_{i}\) follows a GP distribution, and its probability mass function can be written as:
The mean and variance of the GP distribution are \(\mu =\frac{\alpha }{\left(1\delta \right)}=\vartheta \alpha\) and \(Var\left(y\right)={\frac{\alpha }{{\left(1\delta \right)}^{3}}=\vartheta }^{2}\mu\) respectively where \(\delta\) is called the dispersion parameter expressed by the dispersion factor \(\vartheta =\frac{1}{\left(1\delta \right)}\).
If \(\delta =0\), GP distribution reduces to the standard Poisson distribution when \(\delta <0\), it represents underdispersion, and if \(\delta >0\), it represents overdispersion.
Suppose explanatory variables are represented by \(\left(K1\right)\) dimensional vector \({X}_{i}^{{\prime }}={X}_{1i},\dots ,{X}_{ki}\). The conditional distribution of \({Y}_{i}\) for a given value of \({x}_{i}\) follows a GP distribution with the mean value given by:
where \(f\left({x}_{i},\beta \right)>0\) represents a differentiable function, \({C}_{i}\) represents a measure function and \(\underline\beta\) is the Kdimensional vector of regression parameters.
From the mean of GPD, \(\mu =\frac{\alpha }{\left(1\delta \right)}\) and \(\vartheta =\frac{1}{\left(1\delta \right)}\) the dispersion factor, then the generalized Poisson regression can be deduced as
where, \(\mu =\mu \left(x\right)>0\), \(\vartheta \ge \text{max}({}^{1}\!\left/ \!{}_{2}\right.,1{}^{\mu }\!\left/ \!{}_{4}\right.)\), \(\sqrt{\frac{var(y/x)}{\mu (y/x)}}\) stands for the square root of the dispersion index, and m is the largest positive integer for which \(\mu +m\left(\vartheta 1\right)>0\) when \(\vartheta\) is nonnegative.
When \(\vartheta =1\), GP distribution is condensed to standard Poisson regression (proper in modeling equal dispersed data); when \(\vartheta >1\) GPR is appropriate in modeling overdispersed data, and when \(\vartheta <1\), GPR is used to fit underdispersed data [34, 40].
Similar to the standard Poisson regression model, GPR uses a log link to connect the mean of the response variable and explanatory variables, as shown below:
where, \(\mu =\mu \left(x\right)=\frac{\alpha }{\left(1\delta \right)}\) is the mean, \(\underline x_i^T\) represents the \((k1)\) dimensional vector of explanatory variables and \(\underline\beta\) is the kdimensional vector of regression parameters.
In this study, diabetes counts in the council in 2020 have been used as a response variable regressed to the following elaborated explanatory variables:
S/no  Variable name  Variable description 

1  \({X}_{1}\)  Percentage of people living with HIV (PPLHIV) in the year 2020 in a specific council 
2  \({X}_{2}\)  Councilâ€™s 2020 population projection 
3  \({X}_{3}\)  Councilâ€™s estimated gross domestic product (GDP) at market price 
4  \({X}_{4}\)  Councilâ€™s residence (Rural or urban) 
5  \({X}_{5}\)  Number of health facilities in the council 
6  \({X}_{6}\)  Percentage of males in total diabetes count 
7  \({X}_{7}\)  Zone (Lake, Southern, Southern Highlands, Central, Eastern, Northern, and Western zones) 
The model can be written as:
Then,
Estimation of model coefficients \(\underline\beta\) was performed through the maximum likelihood method. Additionally, the goodness of fit of the GP model over the NB and Poisson models is also evaluated using AIC, AICc, and BIC.
Results
The histogram in Fig.Â 1 describes the dispersion property of diabetes incidence across councils. The plot indicates a significant positive skew, with more small numbers, including zero, and few large counts, suggesting overdispersion among diabetic patients incidence between two age groups, namely, age 5 to 59 and 60 or older. This is common among disease incidence datasets since sometimes disease severity is triggered by behavioural patterns among subpopulations being sampled, which vary from one society to another, leading to unequal dispersion. Since the data used reveals unequal dispersion, the GP model may give a precise estimate with meaningful inference [19, 20, 24, 39, 40].
The beeswarm plots in Fig.Â 2 display diabetes records per geographical location. The plots indicate a concentration of diabetes in areas with similar traits. Categories are arranged in ascending order based on the number of diabetes cases reported. The councils in the northern zone have more diabetes cases than the other zones, while the councils in the southern zone have fewer counts than other zones. Furthermore, councils inside highcount zones record fewer zeros and low counts than councils within lowcount zones. Additionally, there is a significant difference in diabetes records between rural and urban councils, with rural councils record many zero and small incidents while urban councils record a substantially large number of diabetes cases.
TablesÂ 1 and 2 summarize the association between diabetes count categories and categorical predictors in Tanzania mainland for patients with 5â€“59 and 60 years and above age groups respectively. The dataset consists of diabetes records from 184 councils, with the minimum and maximum recorded numbers being 0 and 958, respectively, across two age groups. The chisquare test for dependence was used to measure the presence of a statistically significant association between diabetes count and two categorical predictors associated with environmental factors: council residence (rural or urban) and council zone (northern, eastern, lake, southern, southern highland, western, and central zones) among datasets from two distinct age groups. For both scenarios, the pvalue is less than 0.001, indicating the presence of an association among categories.
Moreover, TablesÂ 1 and 2 demonstrate larger counts of diabetes incidences recorded among councils located in urban areas than councils in rural areas in diabetes patients aged 5â€“59 years. The presence of large numbers of diabetes records among people aged 5â€“59 years indicates the high chance of premature mortality and morbidity due to diabetes contrary to sustainable development goal number 3.4. For 60 years and older age groups, 60.8% of councils located in rural areas recorded fewer than 50 diabetes patients whereas 8.7% of councils in urban areas recorded fewer than 50 patients. Moreover, in the southern zone, none of the councils recorded more than 200 diabetes cases among both age groups.
TableÂ 3 describes the log of expected diabetes counts as a function of selected predictor variables using the GP model (located at the top of the table), the negative binomial model (located in the middle of TableÂ 3), and the standard Poisson (at the bottom of TableÂ 3).
Based on the GP modelâ€™s results, the number of diabetes cases in each council is influenced by its population. Thus, more populated councils are anticipated to have more cases of diabetes than the less populated ones. The logs of expected diabetes count in a council would be expected to increase by \(0.2264\;(pvalue<0.0001)\) when the councilâ€™s population increases by one unit. The number of health facilities is significantly associated with the number of diabetes cases in the councils. This may be because the availability of health facilities accelerates disease tracking and recording. Increase in the expected number of health facilities in the councils leads to increase the log of expected diabetes counts by \(0.0132\;(pvalue<0.0001)\). Moreover, the results suggest positive association between percentage of peoples living with HIV and diabetes incidences in the council. Conversely, GDP per capita shows a significantly negative association with the log of expected diabetes counts in the councils. This implies that diabetes cases happen more in councils with less GDP per capita. On the other hand, there is no significant association between the percentage of male diabetes patients and the number of diabetes cases in the councils.
Predictors representing environmental factors are significantly associated with diabetes counts in the councils. It can be demonstrated that, when other model covariates are held constant, the difference in logs of diabetes counts is predicted to be \(1.172\;(pvalue< 2\text{e}16 )\) larger for councils located in urban areas than those in rural areas. Compared to the northern zone, councils located in the central, eastern, lake, southern, and southern highlands zones have a decreased log diabetes counts of \(0.8480\;(pvalue=7.72\text{e}05)\), \(0.6483\;(pvalue=0.00024\), \(0.7265\;(pvalue=1.17\text{e}06)\), \(0.8064\;(pvalue=3.51\text{e}05)\), and \(0.7467\;(pvalue=4.78\text{e}05)\) respectively. The reason is because councils in the northern zone contribute more diabetes cases than the other zones. Additionally, there is no significant difference in log diabetes cases in the western zone compared to councils in the northern zone.
There are slight differences between the estimates and standard errors (SEs) obtained by the GP and NB models, resulting in different inferences for the western zone category. The GP model shows that the categoryâ€™s contribution to the logs of diabetes did not differ from that in the northern zone. In contrast, the NB model shows a significantly decreased log diabetes count by \(0.6395\;(pvalue=0.00899)\), compared to the northern zone when other factors in the model are kept constant. Additionally, the results in TableÂ 3 indicate that SEs in the Poisson model were underestimated because the values were visually smaller than those obtained in the GP and NB models. This occurs because the Poisson model cannot handle the overdispersion present in the analyzed datasets. Underestimating SEs leads to incorrect inferences being drawn about some predictors and factors.
Although the GP model finds that one predictor variable (percentage of males hospitalized by diabetes) and one category (western zone) are not important, the NB model only finds the percentage of males hospitalized by diabetes to be insignificant. However, all predictors are deemed significant in the Poisson model. This shows how the GP model excels in controlling overdispersion and producing precise estimates compared to the NB and Poisson models.
Based on the results in TableÂ 3 from the GP regression models, we have provided prediction equations for the average diabetes count as follows:
The antilogarithm of the prediction equation above gives the expected number of diabetes cases as given below:
TableÂ 4 gives the goodness of fit results obtained using different information criteria. The results show that the GP model earns the smallest information criteria values, which means that it outperforms the NB and Poisson models in modeling the used data. Moreover, the results show a slight difference among values obtained by NB and GPD, which may indicate that these two models have slight differences when used to model overdispersed data. However, the major difference between them is that the GP model is appropriate for modeling equal, over, and underdispersed data, while the NB model is used for modeling overdispersed data.
In TableÂ 4, there are Pearson chisquare (Pearson\({\chi }^{2}\)) and Pearson\({\chi }^{2}/DF\) values for the GP, NB, and Poisson models. A value of Pearson\({\chi }^{2}/DF\) greater than one means there is overdispersion, and if it is exactly or close to one, it means overdispersion is well controlled. The GP model has a value closest to one compared to the other models, making it the best choice for modeling overdispersed diabetes count data.
Discussion
This paper suggests utilizing the GP model to model socioenvironmental and other risk factors associated with diabetes incidences in Tanzania mainland. The GP modelâ€™s performance was compared to that of NB and Poisson, as these three models are related. The NB model was obtained through a parametrization process called Poisson mixture gamma, which can model overdispersed data that the standard Poisson model cannot. Additionally, the model can be reduced to the Poisson model when the dispersion parameter tends to infinity. Similarly, the GP model was obtained as a limit of the NB model and can model over, under, and equally dispersed count data. Similar to the NB model, the GP model can also be reduced to the Poisson model when its dispersion parameter equals zero. These models belong to the GLM category and are widely used in analyzing the relationship between a response variable that follows exponential families of distributions and one or more predictor variables. Linear models are a specific type of GLM with an identity link function [9, 41, 42]. The link function transforms the response variable to conform to the linear model assumption, connecting the mean of the response variable to a linear combination of predictor variables.
This studyâ€™s findings reveal that the unequal dominance of diabetes cases is associated with the type of council residence. Both descriptive and inferential analyses show that urban areas have more diabetes cases than rural areas probably due to the lifestyles in the two areas. Urban areas showed a strong positive contribution to diabetes cases, supporting that environmental factors, including urbanization, are a significant risk factor for diabetes and other NCDs [5, 43]. The findings also show a significant difference in the predicted log of diabetes cases among various zones. This indicates heterogeneity of the burden across socioenvironmental attributes. The northern zone, the reference category, appears to have made a significant contribution, causing the projected log of diabetes counts in other zones to be adverse. The western zone was found to have a negligible association with the log of diabetes cases compared to the northern zone according to the GP model. This finding is related to those of Stanifer et al. [44], who observed that hypertension is environmentally clustered since people living together share socialcultural norms like eating habits, crops produced, and other behavioural patterns that affect NCDs.
The study also investigated the contribution of other factors in diabetes cases, and the findings revealed that an increased log of diabetes counts is also associated with the councilâ€™s population and the number of health facilities. On the contrary, GDP at market price is shown to be negatively associated with the log of diabetes counts. This indicates diabetes incidences are also more common in lowincome societies. Several researchers have observed a high NCD rate in low and middleincome countries (LMICs), which aligns with these findings [3]. On the other hand, the total number of patients who attended hospitals for HIV care is not associated with diabetes cases. This result differs from those obtained by Castilho et al. [45]. The percentage of male diabetes cases does not significantly relate to total diabetes cases in the councils. This factor is used to measure the contribution of sex to diabetes incidence, as other studies concluded that there is a higher prevalence of NCD cases among males than females in Africa [3]. Also, there is empirical evidence of a high economic burden among poor households in Tanzania caused by NCDs [46]. This study findings reveals dominance of diabetes incidences among councils with low GDP which may increase poverty contrary to Sustainable Development Goal 1.
The GP model performs better than both the NB and traditional Poisson regression models based on the loglikelihood value, AIC, BIC, AICc and Pearson\({\chi }^{2}\) values. The model achieves the lowest value among all information criteria, suggesting that GP is better at controlling overdispersion among diabetes counts than its competitors. To determine the dispersion value of the data, one can also divide the Pearson chisquare value by its degree of freedom. This value should be close to or equal to 1 for equally dispersed data in the Poisson model. In this study, the value of Pearson\({\chi }^{2}/DF\) for the Poisson model is far greater than 1, indicating the presence of overdispersion. The problem is well handled in the GP model.
Conclusion
Considering the variability of count data when conducting statistical modeling is crucial. Ignoring this factor can lead to false estimates of the standard error, affecting the test statistic and pvalue. It is crucial to examine the dispersion nature of the data to avoid incorrect inferences during statistical modeling of count data.
Hence, this study recommends the use of the GP model in modeling risk factors associated with disease count incidences, specifically in data collected among population subgroups with varying social and environmental characteristics. The model can accommodate count data collected in population subgroups with equal and unequal dispersion. The model is advantageous because it does not involve a difficult computation burden, it does not suffer from convergence issues and gives precise results compared to the most applied NB and Poisson models.
Limitations of the study
The data in DHIS2 are recorded for very broad age groups which hinders further comparison regarding disease incidences. Additionally, the system does not include important patient information, which also limits model variables.
Availability of data and materials
The dataset that has been analyzed cannot be accessed by the public due to regulations regarding the usage and publication of health data in Tanzania. This study has obtained permission to use DHIS2 data but is strictly prohibited from sharing it. However, interested parties can request access to the data by submitting a proper request to the Tanzania Ministry of Health.
Abbreviations
 AIC:

Akaike Information Criteria
 AICc:

Corrected Akaike Information Criteria
 BIC:

Bayesian Information Criteria
 DF:

Degrees of Freedom
 DHIS2:

District Health Information System 2
 GDP:

Gross Domestic Product
 GLM:

Generalized Linear Model
 GP:

Generalized Poisson
 HIV:

Human Immunodeficiency Virus
 NB:

Negative Binomial
 NBS:

National Bureau of Statistics
 NCD:

NonCommunicable Disease
 SE:

Standard Error
References
Budreviciute A, Damiati S, Sabir DK, Onder K, SchullerGoetzburg P, Plakys G, et al. Management and Prevention Strategies for noncommunicable diseases (NCDs) and their risk factors. Front Public Heal. 2020;8:574111.
Blangiardo M, Boulieri A, Diggle P, Piel FB, Shaddick G, Elliott P. Advances in spatiotemporal models for noncommunicable disease surveillance. Int J Epidemiol. 2020;49:I2637.
Bennett JE, Stevens GA, Mathers CD, Bonita R, Rehm J, Kruk ME, et al. NCD countdown 2030: worldwide trends in noncommunicable disease mortality and progress towards sustainable development goal target 3.4. Lancet. 2018;392(10152):1072â€“88.
Ploth DW, Mbwambo JK, Fonner VA, Horowitz B, Zager P, Schrader R, et al. Prevalence of CKD, Diabetes, and Hypertension in Rural Tanzania. Kidney Int Rep. 2018;3(4):905â€“15. https://doi.org/10.1016/j.ekir.2018.04.006.
Roman WP, Martin HD, Sauli E. Cardiovascular diseases in Tanzania: The burden of modifiable and intermediate risk factors.Â J Xiangya Med. 2019;4(6):1â€“14. https://doi.org/10.21037/jxym.2019.07.03.
Thandi P, Lungiswa T, David S, Whadiah P. Chronic noncommunicable diseases: Primary Health Careâ€¯: programme areas. South African Heal Rev. 2008;2008(1):73â€“87. https://doi.org/10.10520/EJC35511.
Wang T, Wang H, Zeng Y, Cai X, Xie L. Health beliefs associated with preventive behaviours against noncommunicable diseases. Patient Educ Couns. 2022;105(1):173â€“81. https://doi.org/10.1016/j.pec.2021.05.024.
Middleton L, Hall H, Raeside R. Applications and applicability of Social Cognitive Theory in information science research. J Librariansh Inf Sci. 2019;51(4):927â€“37.
Agresti A et al. (University of F. An Introduction to Categorical Data Analysis. 3rd Editio. Balding DJ, Cressie NAC, Fitzmaurice GM, Givens GH, Goldstein H, Molenberghs G, editors. 111 River Street, Hoboken; 2019. 390 p.
Mccullagh P. Generalized linear models I Peter McCullagh. 2007.
Yirga AA, Melesse SF, Mwambi HG, Ayele DG. Negative binomial mixed models for analyzing longitudinal CD4 count data. Sci Rep. 2020;10(1):1â€“15. https://doi.org/10.1038/s41598020738837.
Algamal ZY. Variable selection in count data regression model based on firefly algorithm. Stat Optim Inf Comput. 2019;7(2):520â€“9.
Abdulahad FN, Algamal ZY. Almost unbiased ridge estimator in the count data regression models. Electron J Appl Stat Anal. 2021;14(1):44â€“57.
Algamal ZY. Diagnostic in poisson regression models. Electron J Appl Stat Anal. 2012;5(2):178â€“86.
Stoklosa J, Blakey RV, Hui FKC. An overview of modern applications of negative binomial modelling in Ecology and Biodiversity. Diversity. 2022;14(5):1â€“25.
Feng C, Feng C, Li L, Sadeghpour A. A comparison of residual diagnosis tools for diagnosing regression models for count data. BMC Med Res Methodol. 2020;20(1):1â€“21.
Consul PC, Jain GC. A generalization of the poisson distribution. Technometrics. 1973;15(4):791â€“9.
Harris T, Yang Z, Hardin JW. Modeling underdispersed count data with generalized Poisson regression. Comput Stat Data Anal. 2012;12(4):736â€“47.
CuevaLÃ³pez V, OlmoJimÃ©nez MJ, RodrÃguezAvi J. An over and underdispersed biparametric extension of the waring distribution. Mathematics. 2021;9(2):1â€“15.
OlmoJimÃ©nez MJ, RodrÃguezAvi J, CuevaLÃ³pez V. A review of the CTP distribution: a comparison with other over and underdispersed count data models. J Stat Comput Simul. 2018;88(14):2684â€“706.
Wagh YS, Kamalja KK. Comparison of methods of estimation for parameters of generalized Poisson distribution through simulation study. Commun Stat Simul Comput. 2017;46(5):4098â€“112. https://doi.org/10.1080/03610918.2015.1105971.
Shoukri MM. Maximum likelihood estimation for the generalized poisson distribution. Commun Stat  Theory Methods. 1984;13(12):1533â€“47.
Consul PC, Famoye F. Maximum likelihood estimation for the generalized poisson distribution when sample mean is larger than sample variance. Commun Stat  Theory Methods. 1988;17(1):299â€“309.
AlEid M, Shoukri MM. Inference procedures on the generalized Poisson distribution from multiple samples: comparisons with Nonparametric Models for Analysis of Covariance (ANCOVA) of Count Data. Open J Stat. 2021;11(03):420â€“36.
Mayige M, Kagaruki G, Ramaiya K, Swai A. Non communicable diseases in Tanzania: a call for urgent action. Tanzan J Health Res. 2011;13(5 SUPPLISS):1â€“11.
Mpondo B, Mchonde G, Munisi D, Nyundo A, Fabian F. The rising burden of noncommunicable diseases in Tanzania: Countryâ€™s response and preparedness e. Tanzan J Heal Sci. 2018;2(1):11â€“4.
Kato I, Tumaini B, Pallangyo K. Prevalence of noncommunicable diseases among individuals with HIV infection by antiretroviral therapy status in Dar es Salaam, Tanzania. PLoS One. 2020;15:1â€“13. https://doi.org/10.1371/journal.pone.0235542.
Kavishe B, Biraro S, Baisley K, Vanobberghen F, Kapiga S, Munderi P, Smeeth L, et al. High prevalence of hypertension and of risk factors for noncommunicable diseases (NCDs): a population based crosssectional survey of NCDS and HIV infection in Northwestern Tanzania and Southern Uganda. BMC Med. 2015;13:126.
Galson SW, Staton CA, Karia F, Kilonzo K, Lunyera J, Patel UD, et al. Epidemiology of hypertension in Northern Tanzania: a communitybased mixedmethods study. BMJ Open. 2017;7(11):1â€“10.
Consul PC, Famoye F. Generalized poisson regression model. Commun Stat  Theory Methods. 1992;21(Issue 1):89â€“109.
Husain MM, Bagmar MSH. Modeling underdispersed Count Data using generalized Poisson Regression Approach. Glob J Quant Sci. 2015;2(4):22â€“9.
Mahmoud MM, Alderiny MM. On estimating parameters of censored generalized poisson regression model. Appl Math Sci. 2010;4(13â€“16):623â€“35.
Islam MM, Alam M, Tariquzaman M, Kabir MA, Pervin R, Begum M, et al. Predictors of the number of underfive malnourished children in Bangladesh: application of the generalized poisson regression model. BMC Public Health. 2013;13(1):11.
Cui Y, Kim DY, Zhu J. On the generalized poisson regression mixture model for mapping quantitative trait loci with count data. Genetics. 2006;174(4):2159â€“72.
DurmuÅŸ B, GÃ¼neri Ã–Ä°. An application of the generalized Poisson Model for Over Dispersion Data on the number of strikes between 1984 and 2017. Alphanumeric J. 2020;8(2):25060.
Dobson AJ, Barnett AG. An introduction to Generalizd Linear Models. Fourth Edi. London New York: Taylor & Francis; 2018.
Yang Z, Hardin JW, Addy CL, Vuong QH. Testing approaches for Overdispersion in Poisson Regression versus the generalized Poisson Model. Biom J. 2007;49:565â€“84.
Fernandez GA, Vatcheva KP. A comparison of statistical methods for modeling count data with an application to hospital length of stay. BMC Med Res Methodol. 2022;22(1):1â€“21. https://doi.org/10.1186/s12874022016858.
Consul PC, Jain GC. On some interesting properties of the generalized POISSON distribution. Biom Z. 1973;15(7):495â€“500.
Consul PC, Famoye F. Generalized poisson regression model. Commun Stat  Theory Methods. 1992;21(1):89â€“109.
McCullagh P, Nelder JA. Generalized Linear models. Second edi. Statistical models in S. London New York: Chapman and Hall; 1989.
Hosmer DW, Lemeshow S, Sturdivant RX. Applied Logistic Regression: Chap. 4. Wiley series in probability and statistics. 2013;95:89â€“152.
Angkurawaranon C. Urbanization and internal migration as risk factors for noncommunicable diseases in Thailand. 2015. Available from: https://researchonline.lshtm.ac.uk/id/eprint/2267958/.
Stanifer JW, Egger JR, Turner EL, Thielman N, Patel UD. Neighborhood clustering of noncommunicable diseases: Results from a communitybased study in Northern Tanzania. BMC Public Health. 2016;16(1):1â€“10. https://doi.org/10.1186/s1288901629125.
Castilho JL, Escuder MM, Veloso V, Gomes JO, Jayathilake K, Ribeiro S, et al. Trends and predictors of noncommunicable disease multimorbidity among adults living with HIV and receiving antiretroviral therapy in Brazil. J Int AIDS Soc. 2019;22(1):e25233.
Murphy A, Palafox B, Walli M, Jackson TP, Rangarajan S, Alhabib KF et al. The household economic burden of communicable diseases in 18 countries. 2020;5:1â€“13.
Acknowledgements
The authors would like to express their gratitude for the invaluable support provided by their host institutions, namely the University of Dodoma and the University of Dar es Salaam. They would also like to extend their appreciation to their colleagues for their insightful comments and engaging discussions. Lastly, they would like to give special recognition to those who facilitated the availability of data.
Funding
There was no financial support provided for this study by any funding entity.
Author information
Authors and Affiliations
Contributions
SHM developed the concept of the study and was responsible for designing it, collecting and analyzing the data, and writing the manuscript. MCM made important revisions to the manuscript and provided extra input. GSR arranged the manuscript and conducted thorough editing before its submission. All authors have read and approved the final manuscript.
Corresponding author
Ethics declarations
Ethics approval and consent to participate
Ethical approval to conduct this study and use of DHIS2 data was obtained from the Tanzania Ministry of Health. The study utilized diabetes incidence records among councils; no patient interviews or experiments on humans or human tissue samples were conducted. The secondary data were collected from the DHIS2 platform under the Tanzania Ministry of Health, no health facility or research institution was involved in the data collection. Tanzania Ministry of Health also waived the need for informed consent from individuals since the study utilized human data at the community level. All methods were conducted in accordance with applicable guidelines and regulations.
Consent for publication
Not applicable.
Competing interests
The authors declare no competing interests.
Additional information
Publisherâ€™s Note
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated in a credit line to the data.
About this article
Cite this article
Mbwambo, S.H., Mbago, M.C. & Rao, G.S. Socioenvironmental predictors of diabetes incidence disparities in Tanzania mainland: a comparison of regression models for count data. BMC Med Res Methodol 24, 75 (2024). https://doi.org/10.1186/s1287402402166w
Received:
Accepted:
Published:
DOI: https://doi.org/10.1186/s1287402402166w