Skip to main content

Determinants of malaria among under-five children in Ethiopia: Bayesian multilevel analysis



In Ethiopia, malaria is one of the public health problems, and it is still among the ten top leading causes of morbidity and mortality among under-five children.

However, the studies conducted in the country have been inconclusive and inconsistent. Thus, this study aimed to assess factors associated with malaria among under-five children in Ethiopia.


We retrieved secondary data from the malaria indicator survey data collected from September 30 to December 10, 2015, in Ethiopia. A total of 8301 under-five-year-old children who had microscopy test results were included in the study. Bayesian multilevel logistic regression models were fitted and Markov chain Monte Carlo simulation was used to estimate the model parameters using Gibbs sampling. Adjusted Odd Ratio with 95% credible interval in the multivariable model was used to select variables that have a significant association with malaria.


In this study, sleeping under the insecticide-treated bed nets during bed time (ITN) [AOR 0.58,95% CI, 0.31–0.97)], having 2 and more ITN for the household [AOR 0.43, (95% CI, 0.17–0.88)], have radio [AOR 0.41, (95% CI, 0.19–0.78)], have television [AOR 0.19, (95% CI, 0.01–0.89)] and altitude [AOR 0.05, (95% CI, 0.01–0.13)] were the predictors of malaria among under-five children.


The study revealed that sleeping under ITN, having two and more ITN for the household, altitude, availability of radio, and television were the predictors of malaria among under-five children in Ethiopia. Thus, the government should strengthen the availability and utilization of ITN to halt under-five mortality due to malaria.

Peer Review reports


Malaria is one of the most public health problems worldwide, in which 300 to 500 million cases and more than one million deaths were reported in 2018. Of these deaths, 90% were reported from Sub-Saharan African countries [1, 2]. This is equivalent to one child in sub-Saharan Africa dying of malaria every 2 min, and it is one of the fourth leading cause of death of children under the age of 5 years in developing countries in general [1, 3].

In Ethiopia, 68% of the areas are endemic for malaria, and 60% of the country population are prone for infection of malaria [4]. Transmission of malaria mainly depends on the temperature, humidity, and rainfall [5]. In addition, its transmission is also determined by socioeconomic conditions, demographic factors, community-related factors, knowledge, and access to malaria prevention modalities [6].

After the increased utilization of indoor residual spray (IRS), Artemisinin-based combination therapy (ACT), and insecticide-treated bed nets (ITNs), deaths and admissions due to malaria among under-five children decreased by 81 and 73% respectively [7].

Nevertheless, still malaria is one of the ten top leading causes of morbidity and mortality among under-five children in Ethiopia [8]. Accordingly, morbidity and mortality rates of the disease blow up during epidemics time [9].

In the last decades, numerous strategies, policies, and malaria elimination programs had been implemented at the global and national levels [10,11,12]. As a result, over 6.2 million malaria deaths were prevented between 2000 and 2015 in sub-Saharan African countries [11]. Thus, despite the significant decline in the burden of malaria, the disease is still one of the major public health concerns in Ethiopia [13]. Moreover, the magnitude of malaria infection among under-five children is too high which ranges from 16 to 54% in Ethiopia [14,15,16,17,18].

Former studies conducted across the world have identified determinants of malaria among under-five children. Some of the determinants include; altitude of residence site, number of bed nets in the household, child sex, child age, living near to dam, housing conditions, presence of forest cover, and family size [15,16,17, 19,20,21]. Besides, there is still a need to investigate and identify the effect of these factors among under-five children in our country, Ethiopia, to allow an effective preparation of a national malaria prevention, control strategy, and intervention guidelines.

In Ethiopia, although some studies have been conducted to determine factors of malaria among under-five children, those previous studies [16, 17, 22,23,24,25,26,27] were limited (geographic area) and inconclusive (relatively small sample size) to show the determinants of malaria at the national level. Also, those previous studies did not include some community level factors like region, enumeration area, and altitude at national level. In this study, we have included those independent community-level variables, with a relatively large sample size that could help to make inference at a national level, besides we have also used a new statistical Bayesian multilevel analysis approach to assess the determinants of malaria infection. This Bayesian multilevel analysis approach is a powerful statistical method in medical and public health researches [28] that is done based on the data at hand and prior information that has been existed before.

Methods and materials

Study area

Ethiopia is an East African country with an estimated population of more than 100 million which makes it the second-most populous country in Africa [29]. Administratively, Ethiopia is divided into regions and regions are divided into zones, and zones, into administrative units called district (wereda). Each district is further subdivided into kebele that is the lowest administrative unit. Kebeles are also further subdivided into enumeration areas (EAs). Ethiopia, located within 3.30°–15°N, 33°–48°E in the northeastern part of Africa, and it has a total area of 1.1 million square kilometers. It’s topographic features ranges from mountains as high as Ras Dashen 4550 m (m) above sea level (ASL)—to 110 m below sea level in the Afar Depression [4].

Around two-thirds of the country’s territory is favorable for malaria transmission, with malaria primarily associated with altitude and rainfall. Approximately 60% of Ethiopia’s population lives in a malarious area. The highest period of malaria incidence occurs from September to December and from March to May in most parts of the country. The proportion of the population consisting of children under-5 years of age and pregnant women was estimated to be 14.6% and 3.3%, respectively [4].

Data source, study design, sampling procedure, and sample size

This study used cross-sectional survey data from a secondary source extracted from the Ethiopian Malaria Indicator Survey (EMIS), 2015. The EMIS 2015 was the third survey conducted in Ethiopia, a nationwide sample of 13,875 households from 555 EAs was selected.

The EMIS 2015 was used a two-stage cluster sampling methodology. 555EAs were selected in the first stage. Then a complete mapping and listing of all households in the selected EAs were conducted and 25 households were randomly selected for a total of 13, 875 households. Also, the survey involved testing for anemia and malaria among under-five children in all selected households [4].

Since the microscopic examination is the gold standard for the diagnosis of malaria, for this study, children were considered as malaria positive or negative based on theresult of this test only. In this study, we included all available relevant data for children under-5 years of age from the EMIS 2015. The sample size for this study was those all under-five children who were tested for malaria microscopy test. Thus, the number of children whose data were used in this study was 8301.

Study variables

Dependent variable

Malaria microscopy test result among under-five children (Yes/No). The dependent variable was the malaria microscopy test result which was dichotomized into Yes if the test is positive for Plasmodium falciparum, Plasmodium vivax, and mixed infection (both falciparum and vivax) and No if the test negative for all species.

Independent variable

The determinants of malaria among under-five children were grouped into individual-level and community-level determinants.

Individual-related predictors include; the age of a child, sex of a child, household insecticide-treated net (ITN) ownership, the household status of indoor residual spraying (IRS), utilization of nets, number of nets in the household, availability of electricity, housing conditions of the household (floor materials, roof materials, and wall materials), water sources for drinking, time to get water, availability of television, availability of radio, toilet facilities, community-related predictors include; enumeration area, altitude, and region.

Data management and statistical analysis

Sample allocation in the EMIS to different regions as well as urban and rural areas was not proportional. Thus, sample weights to the data were applied to estimate proportions and frequencies to adjust disproportionate sampling and non-response. Since the normality assumption was violated for continuous variables age and altitude, their median with interquartile range was reported in the descriptive analysis. Natural log transformation was applied for these variables before inclusion in the regression model to cure this problem.

The descriptive analysis was performed using both STATA (version 15) and R (version 3.5.2) statistical software and the inferential statistics were done by bayesian statistical software win BUGS (version1.4.3).

Bayesian multilevel logistic regression model

In the usual classical statistics, the analysis of multilevel logistic regression model is based on estimating parameters through Maximum Likelihood Estimation (MLE) and given the asymptotic properties [30]. However, the Bayesian approach has an advantage over the classical approach in the estimation of the model parameters, which is conducted based on their posterior distribution [31].

The EMIS data has hierarchical nature and clustering effect expected in this hierarchical data nature. Therefore, to account for this clustering effect and to get unbiased parameter estimates bayesian multilevel logistic regression analysis was applied to identify determinants of malaria among under-five children in Ethiopia.

In this study, the basic data structure of the two-level logistic regression is a collection of J groups (enumeration areas) and within-group j (j = 1,2,…, J), a random sample nj of level-one units (individual children). The outcome variable is denoted by;

Yij = 1 if ith children are in the jth enumeration area is positive for microscopy test.

0 if the ith children are in the jth enumeration area is negative for microscopy test.

With probabilities, Pij = 1/Xij, Uij) which is the probability of being positive for an ith child (i = 1,2, …nj) from the jth enumeration area. 1-Pij is the probability of being negative for the ith child (i = 1,2, …nj) from the jth enumeration area. Therefore, the model is

$$ Logit\ {\left({Y}_{ij}\right)}_{-}\left( Xij\beta \right)+{U}_{0j\kern1.25em } where\ {U}_{0j}\sim N\left(0,{\sigma^2}_u\right) $$

Xij is the observed value of the predictor variable for a child i in an enumeration area j and U0jis a random effect.

The convergence of the algorithm

In this study, Markov Chain Monte Carlo (MCMC) algorithm was carried out using the Bayesian statistical software Win BUGS version 1.4.3 [32]. The deviance information criterion (DIC) [33] was used to select the best fitted model. The empirical results from a given MCMC analysis are not deemed reliable until the chain has reached its stationary distribution. The term convergence of an MCMC algorithm denotes whether the algorithm has reached its equilibrium distribution or not. If the algorithm has reached its equilibrium distributions, then the generated sample derives from the true target distribution. Therefore, assessing the convergence of the algorithm is essential for producing results from the posterior distribution of interest.

To assess the convergence algorithm in our study, we used autocorrelation, time series plots, Gelman-Rubin statistic, and density plot. All the plots showed that the algorithm has reached its equilibrium (target) distribution for all parameters.

Summary statistics were carried out from the posterior distribution to describe the covariates and adjusted odds ratio (AOR) with corresponding 95% credible interval in the multilevel multivariable logistic regression model was used to select predictors of malaria among under-five children.


Descriptive statistics of the respondents

From the overall respondents, approximately 51.19% of the total study participants were males.

The median age of children was 32 (interquartile range: 16–48 months). The median altitude of residence site was 1703 (interquartile range: 1259–1938 m) in the sampled households. A total of 3362(40.50%) of households had at least one insecticide-treated net and 4215(50.78%) of the children slept under insecticide-treated net the night before the survey took place. Less than half of 1866(22.48%) the interior walls of houses were sprayed. Furthermore, the majority of 6551(78.92) of respondents have radio in their household and only 568(6.85%) of the respondents have a television in their household (Table 1).

Table 1 Weighted Socioeconomic and Demographic Characteristics Associated with Malaria among under-five children in Ethiopia, 2015

Malaria positivity in the study population

Among those who had a positive microscopy test result, the dominant plasmodium species were Plasmodium falciparum which accounts for 93.75% and followed by Plasmodium vivax which accounts 5.36%. In terms of regional variations, the highest prevalence of malaria was found among children from the Gambela region 5(7.14%) while the Benishangul gumuz region was the second-highest 41(3.04%) (Fig. 1).

Fig. 1

Prevalence of malaria under the age of 5 years old within each sampled region of Ethiopia, 2015

Result of empty (null) classical multilevel logistic regression model

In this null model, the random factor variance was 3.12 (95% CI: 2.18, 4.42). This variance value is greater than 0 and showed that there are EA differences in malaria status among under-five children in Ethiopia. The variance at the individual-level in logistic distributions equal to ∏2/3 (that is,3.29) [34,35,36].

Thus, the intra-enumeration area correlation coefficient (ICC) = 3.12/3.12 + 3.29 = 0.49, which showed that 49% of the overall variability in malaria status is through enumeration areas difference, with the rest 51% is due to the individual differences. Therefore, both the random factor variance and the ICC value suggested using multilevel logistic regression model to account for the EA differences in malaria status.

Result of empty (null) Bayesian multilevel logistic regression model

The Bayesian null model showed that precision of the random factor was 0.19 with 95% credible interval of (0.03–0.58), showing EA differences in malaria status in the country. Since the variance estimate which is reciprocal of the precision (1/0.19 = 5.26) is greater than zero it indicates that there are EA differences in malaria status among under-five children in Ethiopia.

Since the unobserved heterogeneity have a logistic distribution with a variance at the individual-level equivalent to ∏2/3 (that is, 3.29) [34,35,36]. So, the ICC = 5.26/5.26 + 3.29 = 0.62, which implied that 62% of the total variability in malaria status among under-five children is because of differences across enumeration areas and 38% of the variability is accounted by individual differences (Table 2). Both the random factor variance and the ICC value proposed to apply bayesian multilevel logistic regression model for additional analysis to handle the heterogeneity between EAs.

Table 2 Results of Null Bayesian Multilevel Logistic Regression Model on factors associated with malaria under the age of 5 years old, Ethiopia, 2015

When we compare the above two null models, even though, both of them suggested using multilevel analysis for additional analysis, the estimation of random effect shows a wide difference between the estimation of the classical analysis approach and the bayesian analysis approach. Based on classical estimation the random variance was 3.12 with an ICC of 0.49 whereas in Bayesian estimation it was 5.26 with an ICC of 0.62. Thus, the Bayesian multilevel random intercept model is more appropriate than the classical multilevel logistic regression in explaining the variation of malaria across enumeration areas in Ethiopia. Therefore, in this study multilevel analysis using the bayesian estimation technique was considered as a proper approach to estimate model parameters.

Bayesian multilevel logistic regression analysis

The gibbs sampler procedure was applied in three different chains with 123,790 total iterations. After10000 burn-in terms discarded a total of 22,758 samples generated from the full posterior distribution. Non-informative normal prior distribution with mean = 0 and precision = 0.001 for the fixed effect and gamma distribution with scale = 0.1, shape = 0.1 for the variance of random effect was used. The convergence produced from markov chains has been confirmed by convergence assessment plots, before taking any inference from the posterior distribution.

Assessment plots of model convergence

History plot

In this study, the three independently generated chains demonstrated good “chain mixture”, an indication of convergence (Fig. 2).

Fig. 2

History plot for some of the significant parameters in the model

Kernel density plot

The simulated samples from the posterior distribution for each regression coefficient are flat, the uni-modal shape of posterior marginal distribution indicating that simulated parameter value indicates convergence to the target distribution (Fig. 3).

Fig. 3

Density plot for some of the significant parameters in the model


The plots show that the three independent chains were mixed or overlapped to each other, which indicates the presence of convergence (Fig. 4).

Fig. 4

Autocorrelation plot for some of the significant parameters in the model

Gelman-Rubin statistic

On the plots, the green line represents the between variability, the blue line represents the within variance and the red line represents the ratio. Hence, in this plot, the red line seems exactly on 1 which provided the evidence for convergence of parameters (Fig. 5).

Fig. 5

Gelman-Rubin Statistics plot for some of the significant parameters in the model

Model comparison and selection

Models DIC values
Null model (Model with no covariates) 802.53
A model with only individual-related variables 703.84
A model with community-related variables 770.89
Full model (model with individual and community-related variables) 654.02

Among the four models fitted, the full model has the smallest (654.02) DIC value. Therefore, the full model is most likely to fit the data.

Bayesian multivariable multilevel logistic regression model

A Bayesian multivariable multilevel binary logistic regression model was fitted to identify the effect of explanatory variables. The parameters were estimated from the posterior distribution (Table 3).

Table 3 Posterior parameter estimates, SD, MC error, and odds ratios with their 95% credible intervals of the final model for factors associated with malaria under the age of 5 years old, Ethiopia, 2015

In the final model: sleeping under the insecticide-treated bed nets during bed time (ITN), number of ITN available in the household, altitude, availability of radio and television were statistically associated with malaria among under-five children at 95% credible interval.

Holding other covariates constant, the odds of having malaria for under-five-year-old children decreases by 95% with an increase in altitude (AOR 0.05, 95% CI 0.01–0.13).

After adjusted other covariates, children living in households with mosquito bed nets greater than two and more in their household were decreased the odds of having malaria by 57% than those in households without mosquito bed nets (AOR = 0.43, CI 0.17–0.88) and those children who were using mosquito bed nets were found to be decreased the risk of having malaria by 42% compared to those who were not using mosquito bed nets (AOR = 0.58, CI 0.31–0.97).

Keeping other covariates constant, children in households with television were decreased the risk of having malaria by 81% than those in households without television (AOR 0.19, 95% CI 0.01–0.89) and children in households with radio were decreased the odds of having malaria by 59% compared to those in households without radio (AOR 0.41, 95% CI 0.19–0.81).


This study aimed to investigate the relationship between the malaria status of children under the age of 5 years and socio-economic, demographic, and environmental factors based on the EMIS 2015 data using a bayesian analysis approach.

According to this study, the number of nets in the household, sleep under the net, altitude, availability of television and radio were statistically significant at 95% Bayesian credible interval.

Based on this study, among children residing in households with greater than two and more mosquito bed nets, the odds of having malaria was decreased by 57% compared to those in households without nets. This supports studies previously conducted in north-west Ethiopia, Nigeria, and north-east Tanzania [37,38,39].

According to the findings of this study, the use of a mosquito net has also been found as a protective factor against malaria cases among under-five children in Ethiopia. It revealed that those who were slept under mosquito bed net had a 42% decreased risk of malaria positivity as compared to those who did not use a bed net. This results supports also the study conducted by Hadya Zone, Southern Ethiopia which stated that those who were not using bed net were 4.67 times more likely to be infected [40] and another study conducted in Kenya also showed that malaria prevalence was observed to be more than two times higher among households that did not use mosquito nets compared to net user households [6]. This finding also consistent with studies done in different areas [21,22,23, 37, 38, 41]. This might be the fact that insecticides on the net repel mosquitoes, as a result, it decreases the number of mosquitos that enters into the house and attempts to feed on people inside. Besides this, ITN could also be a protective barrier and kill mosquitoes as well [42].

This study found that the availability of mass media played an important role in the prevention of malaria. Compared to children in households without television and radio, those in households with television and radio had 81 and 59% decreased risk of malaria positivity respectively. This finding is in agreement with the reports in Kenya, children of households with no access to television more likely to have malaria than their [6].This supports the positive findings of the influence of mass media for eliminating malaria in African settings [43]. This might be due to those who have watched television or listened to radio programs that might have information related to malaria prevention methods [6].

The association between malaria and altitude showed that the odds of malaria for children decreased with an increase in altitude. This is in agreement with the study conducted in southern Ruanda and Ghana [44, 45]. The finding also supports other studies conducted in Kenya, Uganda, and Northeastern Tanzania [6, 21, 37]. This might be due to the low altitude area is appropriate for mosquito breeding due to its flat topography that allows the collection of water following the rainy season [26].

Due to the fact this analysis was on a secondary data source, there is the possibility that some factors associated with malaria infection were not considered in this analysis because of dependence on variables collected only during the EMIS 2015 survey. Another limitation is that the survey used a cross-sectional design to collect data as such no causal inferences can be made between malaria infection and its determinants. Despite these limitations, the study used survey data collected from a nationally representative sample and has identified some of the predictors of malaria infection as such, the results can be generalized to determinants of malaria among under-five children in Ethiopia that can be used for planning future interventions for malaria prevention. Additional systematic review and meta-analysis study is recommended to have a pooled estimation of malaria prevalence and its determinants at the national level.


The study revealed that sleeping under ITN, availability of ITNs for the household, altitude, availability of radio, and television for the household were the predictors of malaria among under-five children in Ethiopia. Thus, the government should strengthen the availability and utilization of ITN to halt under-five mortality due to malaria. Moreover, the government should also emphasis for the community living in low altitude areas and those residing households with less availability of television and radio.

Availability of data and materials

The datasets used and/or analyzed during the current study are available from the corresponding author on reasonable request.



Addis Continental Institute of Public Health


Artemisinin-based Combination Therapy


Adjusted Odds Ratio


Above Sea Level


Credible Interval


Central Statistical Agency


Deviance Information Criteria


Enumeration Area


Ethiopian Demographic and Health Survey


Ethiopia National Malaria Indicator Survey


Ethiopian Public Health Institute


Federal Ministry of Health


Growth Domestic Product


Global Positioning System


Intra-cluster Correlation Coefficient


Indoor Residual Spraying


Insecticide Treated Net


Long-Lasting Insecticide-Treated Net


Markov Chain Monte Carlo


Medical Research Council (Cambridge)


National Strategic Plan


President’s Malaria Initiative


Rapid Diagnostic Test


Roll Back Malaria


Standard Deviation


Southern Nations, Nationalities, and Peoples’ Region


United Kingdom


World Health Organization


  1. 1.

    World Health Organization. World malaria report 2014. Geneva: WHO; 2014. Fecha de consulta. 2018; 23:247.

    Google Scholar 

  2. 2.

    Phillips RS. Current status of malaria and potential for control. Clin Microbiol Rev. 2001;14(1):208–26.

    CAS  PubMed  PubMed Central  Google Scholar 

  3. 3.

    World Health Organization. World malaria report 2011. Geneva: World Health Organization; 2011. Available at: Accessed 20 Mar 2020.

  4. 4.

    Ethiopian Public Health Institute. Ethiopia national malaria indicator survey 2015. Addis Ababa: Ethiopian Public Health Institute; 2016. Accessed 12 May 2020.

  5. 5.

    Abiodun GJ, Maharaj R, Witbooi P, et al. Modelling the influence of temperature and rainfall on the population dynamics of anopheles arabiensis. Malar J. 2016;15:364.

    Article  PubMed  PubMed Central  Google Scholar 

  6. 6.

    Sultana M, Sheikh N, Mahumud RA, Jahir T, Islam Z. Prevalence and associated determinants of malaria parasites among Kenyan children; 2017. p. 1–9.

    Google Scholar 

  7. 7.

    Aregawi M, Lynch M, Bekele W, Kebede H, Jima D, Taffese HS, et al. Time series analysis of trends in malaria cases and deaths at hospitals and the effect of antimalarial interventions, 2001–2011, Ethiopia. PLoS One. 2014;9:11.

    Google Scholar 

  8. 8.

    Deribew A, Tessema GA, Deribe K, et al. Trends, causes, and risk factors of mortality among children under 5 in Ethiopia, 1990–2013: findings from the global burden of disease study 2013. Popul Health Metrics. 2016;14:42.

    Article  Google Scholar 

  9. 9.

    Ministry of Health. National malaria guidelines. 3rd ed. Addis Ababa: Federal Ministry of Health of Ethiopia; 2012. Accessed 20 May 2020.

  10. 10.

    FMOH. Health sector transformation plan (HSTP 2016–2020). Addis Ababa: FMOH; 2015. p. 184.

    Google Scholar 

  11. 11.

    Nations U. The millennium development goals report; 2015. p. 75.

    Google Scholar 

  12. 12.

    United Nation (UN). Transforming our world: the 2030 agenda for sustainable development. Sustainabledevelopmentunorg; 2015.

    Google Scholar 

  13. 13.

    Jima D, Wondabeku M, Alemu A, Teferra A, Awel N, Deressa W, et al. Analysis of malaria surveillance data in Ethiopia: what can be learned from the integrated disease surveillance and response system? Malar J. 2012;11:330.

    PubMed  PubMed Central  Google Scholar 

  14. 14.

    Ashton RA, Kefyalew T, Tesfaye G, Pullan RL, Yadeta D, Reithinger R, et al. School-based surveys of malaria in Oromia regional state, Ethiopia: a rapid survey method for malaria in low transmission settings. Malar J. 2011;10:25.

    PubMed  PubMed Central  Google Scholar 

  15. 15.

    Ayalew S. The prevalence of malaria and the associated risk factors in Jiga area, Northwest Ethiopia. Addis Ababa; 2014.

    Google Scholar 

  16. 16.

    Debo GW, Kassa DH. Prevalence of malaria and associated factors in Benna Tsemay district of pastoralist community, southern Ethiopia. Trop Dis Travel Med Vaccines. 2016;2:16.

    PubMed  PubMed Central  Google Scholar 

  17. 17.

    Haji Y, Fogarty AW, Deressa W. Prevalence and associated factors of malaria among febrile children in Ethiopia: a cross-sectional health facility-based study. Acta Trop. 2016;155:63–70.

    PubMed  Google Scholar 

  18. 18.

    Molla E, Ayele B. Prevalence of malaria and associated factors in Dilla town and the surrounding rural areas, Gedeo zone, southern Ethiopia. J Bacteriol Parasitol. 2015;6:5.

    Google Scholar 

  19. 19.

    Tadesse F, Fogarty AW, Deressa W. Prevalence and associated risk factors of malaria among adults in east Shewa zone of Oromia regional state, Ethiopia: a cross-sectional study. BMC Public Health. 2018;18(1):25.

    Google Scholar 

  20. 20.

    Tassew A, Hopkins R, Deressa W. Factors influencing the ownership and utilization of long-lasting insecticidal nets for malaria prevention in Ethiopia. Malar J. 2017;16(1):262.

    PubMed  PubMed Central  Google Scholar 

  21. 21.

    Roberts D, Matthews G. Risk factors of malaria in children under the age of five years old in Uganda. Malar J. 2016;15(1):246.

    PubMed  PubMed Central  Google Scholar 

  22. 22.

    Alemu A, Abebe G, Tsegaye W, Golassa L. Climatic variables and malaria transmission dynamics in Jimma town, south West Ethiopia. Parasit Vectors. 2011;4(1):30.

    PubMed  PubMed Central  Google Scholar 

  23. 23.

    Ayele DG, Zewotir TT, Mwambi HG. Prevalence and risk factors of malaria in Ethiopia. Malar J. 2012;11(1):195.

    PubMed  PubMed Central  Google Scholar 

  24. 24.

    Tefera G. Prevalence of malaria and associated factors among patients attending at health center, southern Ethiopia. Immunol Infect Dis. 2014;2(3):25–9.

    Google Scholar 

  25. 25.

    Woday A, Mohammed A, Gebre A, Urmale K. Prevalence and associated factors of malaria among febrile children in Afar region, Ethiopia: a health facility based study. Ethiop J Health Sci. 2019;29:5.

    Google Scholar 

  26. 26.

    Woyessa A, Deressa W, Ali A, Lindtjørn B. Malaria risk factors in Butajira area, south-Central Ethiopia: a multilevel analysis. Malar J. 2013;12(1):273.

    PubMed  PubMed Central  Google Scholar 

  27. 27.

    Deribew A, Alemseged F, Tessema F, Sena L, Birhanu Z, Zeynudin A, et al. Malaria and under-nutrition: a community based study among under-five children at risk of malaria, south-West Ethiopia. PLoS One. 2010 May 21;5(5):e10775.

    PubMed  PubMed Central  Google Scholar 

  28. 28.

    Contreras LF, Brown ET, Ruest M. Bayesian data analysis to quantify the uncertainty of intact rock strength. J Rock Mech Geotech Eng. 2018;10(1):11–31.

    Google Scholar 

  29. 29.

    Central Statistical Agency (CSA) [Ethiopia] and ICF. Ethiopia demographic and health survey 2016. Addis Ababa and Rockville: CSA and ICF; 2016.

    Google Scholar 

  30. 30.

    Li M, Meeker WQ. Application of Bayesian methods in reliability data analyses. J Qual Technol. 2014;46(1):1–23.

    Google Scholar 

  31. 31.

    Lindley DV. The Bayesian Approach to Statistics. In: de Olivera TJ, editor. Some Recent Advances in Statistics. London: Academic Press; 1982. p. 65–87.

  32. 32.

    University of Cambridge, M.B.U., Cambridge. Bayesian statistical software win BUGS version 1.4.3; 2019.

    Google Scholar 

  33. 33.

    Spiegelhalter DJ, Best NG, Carlin BP, Van Der Linde A. Bayesian measures of model complexity and fit. J Roy Stat Soc. 2002;64(4):583–639.

    Google Scholar 

  34. 34.

    Goldstein H, Browne W, Rasbash J. Partitioning variation in multilevel models. Underst Stat. 2002;1(4):223–31.

    Google Scholar 

  35. 35.

    Rasbash J, Steele F, Browne W, Prosser B. Logistic models for binary and binomial responses. A user’s guide to MLwiN. Version, vol. 20; 2003.

    Google Scholar 

  36. 36.

    Snijders TAB, Bosker RJ. Multilevel analysis: an introduction to basic and advanced multilevel modeling. 1st ed. Thousand Oaks: Sage; 1999.

    Google Scholar 

  37. 37.

    Winskill P, Rowland M, Mtove G, Malima RC, Kirby MJ. Malaria risk factors in north-East Tanzania. Malar J. 2011;10(1):98.

    PubMed  PubMed Central  Google Scholar 

  38. 38.

    Yusuf OB, et al. Poverty and fever vulnerability in Nigeria: a multilevel analysis. Malar J. 2010;9(1):235.

    PubMed  PubMed Central  Google Scholar 

  39. 39.

    Aschale Y, Mengist A, Bitew A, Kassie B, Talie A. Prevalence of malaria and associated risk factors among asymptomatic migrant laborers in west Armachiho District, Northwest Ethiopia. Res Rep Trop Med. 2018;9:95.

    PubMed  PubMed Central  Google Scholar 

  40. 40.

    Delil RK, Dileba TK, Habtu YA, Gone TF, Leta TJ. Magnitude of malaria and factors among febrile cases in low transmission areas of Hadiya zone, Ethiopia: a facility based cross sectional study. PLoS One. 2016;11:5.

    Google Scholar 

  41. 41.

    Ayele DG, Zewotir T, Mwambi H. The risk factor indicators of malaria in Ethiopia. Int J Med Med Sci. 2013;5:335–47.

    Google Scholar 

  42. 42.

    Aderaw Z, Gedefaw M. Knowledge, attitude, and practice of the community towards malaria prevention and control options in anti-malaria association intervention zones of Amahara national regional state, Ethiopia. Global J Med Res. 2013;1:1.

    Google Scholar 

  43. 43.

    Yaya S, Uthman OA, Amouzou A, Bishwajit G. Mass media exposure and its impact on malaria prevention behaviour among adult women in sub-Saharan Africa: results from malaria indicator surveys. Global Health Res Policy. 2018;3(1):20.

    Google Scholar 

  44. 44.

    Nyirakanani C, Chibvongodze R, Habtu M, Masika M, Mukoko D, Njunwa KJ. Prevalence and risk factors of asymptomatic malaria among underfive children in Huye District, southern Rwanda. Tanzania J Health Res. 2018;20:1.

    Google Scholar 

  45. 45.

    Kweku M, Takramah W, Kudzo Axame W. Prevalence and risk factors of malaria among children under-five years in high and low altitude rural communities in the Hohoe municipality of Ghana. J Clin Immunol Res. 2017;1(1):1–8.

    Google Scholar 

Download references


Many thanks to the EPHI (Ethiopian Public Health Institute) for providing and granting permission for the use of the data in this study. A special thanks to the University of Gondar, Institute of Public Health for providing us Ethical clearance to conduct this study.


No specific funding source to conduct this study.

Author information




SB conceived the study, undertook statistical analysis and drafted the manuscript. KA, GA, DA, AW, AA and DD initiated the study and made major contributions to the study design and statistical analysis. All authors read and approved the final manuscript.

Corresponding author

Correspondence to Setognal Birara Aychiluhm.

Ethics declarations

Ethics approval and consent to participate

Ethical clearance was obtained from the Institutional Review Committee of Institute of Public Health, College of Medicine and Health Sciences, University of Gondar. Permission for data access was obtained from Ethiopian Public Health Institute (EPHI).

Consent for publication

Not Applicable.

Competing interests

The authors declare that they have 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 The Creative Commons Public Domain Dedication waiver ( applies to the data made available in this article, unless otherwise stated in a credit line to the data.

Reprints and Permissions

About this article

Verify currency and authenticity via CrossMark

Cite this article

Aychiluhm, S.B., Gelaye, K.A., Angaw, D.A. et al. Determinants of malaria among under-five children in Ethiopia: Bayesian multilevel analysis. BMC Public Health 20, 1468 (2020).

Download citation


  • Malaria microscopy test
  • Bayesian multilevel logistic regression
  • Ethiopia