Background Globally, at least 4.7 million people die from being overweight or obese. In Ethiopia, the level of overweight and obesity among women grew from 3% to 8%. However, as far as my literature searching, studies concerning the spatial variation of overweight/obesity and factors associated are not researched in Ethiopia using geospatial techniques. Therefore, this study aimed to explore the spatial variation of overweight/obesity and factor associated among reproductive age group women in Ethiopia using geospatial techniques. Mothed A total weighted sample of 10,928 reproductive age women were included in the study. ArcGIS version10.7 was used to explore the spatial variation of overweight/obesity. Bernoulli based model was used to analyze the purely spatial cluster detection of overweight/obesity through SaTScan version 9.6.1 software. Ordinary Least Square analysis and geographically weighted regression analysis was employed to assess the association between an outcome variable and explanatory variables by using ArcGIS 10.7 software. P value of less than 0.05 was used to declare statically significant. Result The spatial distribution of overweight/obesity in Ethiopia was clustered. Statistically, a significant-high hot spot overweight/obesity was identified at Addis Ababa, harrari, Dire Dawa. SaTScan identified 66 primary spatial clusters (RR = 4.17, P < 0.001) located at Addis Ababa, southeast amhara, central part of oromia region and northern part of SNNP region. In geographically weighted regression, rich wealth index, women’s age (35–39 and 40–44 years), watching TV, internet use and not working were statistically significant that affecting spatial variation of overweight/obesity. Conclusion In Ethiopia, overweight/obesity varies across the region. Statistically, significant-high hot spots of overweight/obesity were detected in Addis Ababa, Harari, Dire Dawa, some parts of Amhara and afar region, most of the Oromia and Somalia region, and the South Nation Nationality and People region of Ethiopia. Therefore, the ministry of health and the Ethiopian public health institute, try to initiate policies and practices that could include providing funding for physical education as well as recreational centers in communities most in need. In addition, public and private mass media create awareness of healthy lifestyles is promoted by health education regarding increased physical activity and reduced sedentary behavior through various media platforms.
In EDHS 2016, a community-based cross-sectional study was conducted by the Central Statistical Agency (CSA) from January 18 to June 27, 2016, in Ethiopia [15]. It has nine Regional states (Afar, Amhara, Benishangul-Gumuz, Gambella, Harari, Oromia, Somali, Southern Nations, Nationalities, and People’s Region (SNNP) and Tigray) and two city Administrative (Addis Ababa and Dire-Dawa) (Fig 1). The source of population was all reproductive-age women within five years before the survey in Ethiopia. A weighted total of 10,928 reproductive age group women who had a complete answer to all variables of interest were included from the total 15,683 women aged 15–49 years were interviewed. All women those who had underweight and women for whom there was missing information on height and/or weight and/or women those who had pregnant and observation from enumeration areas with zero coordinate were excluded. The data was acquired from the Demographic and Health Surveys (DHS) Program by requesting for this work and accessed www.dhsprogram.com website. In the 2016 EDHS, a standardized and validated questionnaire were adapted from the DHS Program’s standard questionnaires in a way to reveal the population and health issues relevant to Ethiopia. A two stage stratified sampling technique was employed to select representative samples for the country as whole. The regions of the country were stratified into urban and rural areas. Then, samples of enumeration areas (EAs) were selected in each stratum in two stages. In the first stage, 645 EAs were selected with probability proportional to the EA size. The EA size is the number of residential households in the EA as determined in the 2007 Ethiopian Population and Housing Census. In the second stage, a fixed number of 28 households per cluster were selected randomly from the household listing [15]. For the purpose of this study, the women’s data (IR) from the 2016 EDHS was utilized. The detailed sampling procedure is available in the Ethiopian Demographic and Health Survey reports from Measure DHS website (www.dhsprogram.com). The outcome variable of this study was Body Mass Index (BMI). BMI is a measure for quantifying nutritional status in adults which is calculated by using the formula weight in kg divided by height in meter square. This study utilized the WHO-specific cut-offs value to categorize BMI of the reproductive age group women. Normal weight: 18.5–24.9 kg/m2, Overweight: 25.0–29.99 kg/m2 and Obese: ≥30.0 kg/m2 [1]. For the purpose of this study, a BMI ≥25.0 kg/m2 was categorized as overweight/obesity and coded as 1, while a BMI 18.5–24.99 kg/m2 was categorized as normal and coded as 0. Socio-demographic factor marital status (labelled as single, married, divorced/separated and widowed), Women’s education status (labelled as no education, primary, secondary, and higher education, maternal occupation (labelled as working or not working), parity(labelled as 0 children, 1–3 children and 4 and above children), wealth index (labelled as poor, middle and rich) and maternal age. Behavioral factors such as alcohol use (labeled as yes and No), cigarette smoking (labeled as yes and No) and contraceptive use (labeled as yes and No). Exposure to mass media (whether the woman read any newspapers/magazines, listened to the radio, or watched television and internet use) were labeled as (less than once a week and at least once a week) and not at all). The data was clean by STATA version 14.1 software and Microsoft Excel, for data analysis we used Arc GIS 10.7 and SaTScan 9.6. Spatial autocorrelation (Global Moran’s I) statistic measure is used to assess whether overweight/obesity is dispersed, clustered, or randomly distributed in Ethiopia. Moran’s I value close to − 1 shows dispersed overweight/obesity, close to + 1 shows clustered, and if Moran’s I value zero shows randomly distributed and a statistically significant Moran’s I (p < 0.05) leads to rejection of the null hypothesis [16]. Hot Spot Analysis (the Getis-Ord Gi * statistic) of the z-scores and significant p-values tells the features with either hot spot or cold spot values for the clusters spatially [17]. For measuring spatial autocorrelation for a set of distances, a line graph of those distances and their related Z-scores was generated. The level of spatial clustering and statistical significance are represented by Z-scores. Peak Z -scores show distances where spatial processes promoting clustering are most pronounced. These peak distances are often proper values to use for tools with a Distance Band or Distance Radius parameter. This tool can help to select an appropriate distance threshold or radius for tools that have these parameters, such as hot spot analysis [18]. The spatial interpolation technique is used to predict overweight/obesity for unsampled areas based on sampled clusters [19]. We used deterministic and geostatistical interpolation methods. To compare the above interpolation method we employed geostatistical analysis based on result with lowest mean predicted error (MPE) and root mean square predicted error (RMSPE) was the best fitted interpolation technique for overweight/obesity. Those small values indicate that predicted values are close to the observed values and vice versa [20]. For this study ordinary kriging interpolation method was selected, since lowest mean predicted error (MPE) and root mean square predicted error (RMSPE). Bernoulli based model spatial Kuldorff’s Scan statistics was used to determine the geographical locations of statistically significant spatial window for overweight/obesity using SaTScan version 9.6.1 software [21]. The outcome variable has a Bernoulli distribution, so Bernoulli model was used by applying the Kuldorff’s method for purely spatial analysis. The scanning window that moves across the study area in which women give overweight/obesity was taken as case and those women who give normal body weight was taken as control to fit the Bernoulli model. The default maximum spatial cluster size of < 50% of the population was used as an upper limit and most likely clusters was identified by using p-values and likelihood ratio tests based on 999 Monte Carlo replications. To generate secondary clusters, we employed non-overlapping options by SaTScan version 9.6.1 and ArcGIS software version 10.7 was used to map the cluster and attribute of overweight produced by SaTScan™. Exploratory Regression was used to find a model that meets the OLS method’s assumptions, all while identifying models with a high Adjusted R2 value. Ordinary Least Square regression (OLS) model is a global model that estimates only one single coefficient per explanatory variable over the entire study region. We used to check assumptions of spatial regression using explanatory regression with the particular tests. The Jarque-Bera test was used to assess the normality assumption for residuals. As residuals are not spatially auto-correlated, the statistically significant Koenker (BP) statistic shows that the relationships modeled are not consistent (either due to non-stationarity or heteroscedasticity). Multicollinearity (Variance Inflation Factor) was used to check redundancy among predictor variables, coefficients have the expected sign and statically significant, and strong adjusted R2 values. A geographically weighted regression model is an extension of the OLS regression model. It gives local parameter estimates to reflect variations over space in the association between an outcome and predictor variables [22]. For geographically weighted regression analysis, the aggregated proportion of overweight/obesity among reproductive-age women and all the predictor variables are considered for each cluster. To determine the predictor variables for overweight/obesity among reproductive-age women, we used a geographically weighted regression model. The model structure of geographically weighted regression written as: Where Yi is the outcome variable, (ui, vi) represents the coordinates of the ith point in space, β0 is the intercept at the (ui, vi) coordinate, βk is the coefficient of the covariate X at the (ui, vi) coordinate, and i is the random error term. Geographical heterogeneity for each coefficient can be measured by comparing the AICc between the GWR model and the global OLS regression model. The corrected Akaike Information Criteria (AIC) and Adjusted R-squared was used for model comparison of OLS (global model) and GWR (local) model. A model with the lowest AICc value and a higher adjusted R-squared value was used to determine the best fit model for local parameter estimates. Ethical clearance was obtained from the ethical review board of the University Of Gondar Institute of Public Health, CMHS. The guidelines expressed in Ethiopia’s Declaration of Central Statistical Agency (DCSA) guided the countrywide survey. The survey was also approved by CSA’s Ethical Review Board (ERB), and everyone who decided to take part in it completed a consent form. Permission for data access was acquired from the measure demographic and health survey through an online request by written letter of objective and significance of the study from http://www.dhsprogram.com. Moreover, for Geographic information system coordinates, the coordinates are only for the enumeration area (EA) as a whole and the measured coordinates were randomly displaced within a large geographic area so that no particular enumeration areas can be identified.
N/A