Treatment-seeking rates in malaria endemic countries
© Battle et al. 2016
Received: 29 August 2015
Accepted: 11 December 2015
Published: 11 January 2016
The proportion of individuals who seek treatment for fever is an important quantity in understanding access to and use of health systems, as well as for interpreting data on disease incidence from routine surveillance systems. For many malaria endemic countries (MECs), treatment-seeking information is available from national household surveys. The aim of this paper was to assemble sub-national estimates of treatment-seeking behaviours and to predict national treatment-seeking measures for all MECs lacking household survey data.
Data on treatment seeking for fever were obtained from Demographic and Health Surveys, Malaria Indicator Surveys and Multiple Cluster Indicator Surveys for every MEC and year that data were available. National-level social, economic and health-related variables were gathered from the World Bank as putative covariates of treatment-seeking rates. A generalized additive mixed model (GAMM) was used to estimate treatment-seeking behaviours for countries where survey data were unavailable. Two separate models were developed to predict the proportion of fever cases that would seek treatment at (1) a public health facility or (2) from any kind of treatment provider.
Treatment-seeking data were available for 74 MECs and modelled for the remaining 24. GAMMs found that the percentage of pregnant women receiving prenatal care, vaccination rates, education level, government health expenditure, and GDP growth were important predictors for both categories of treatment-seeking outcomes. Treatment-seeking rates, which varied both within and among regions, revealed that public facilities were not always the primary facility type used.
Estimates of treatment-seeking rates show how health services are utilized and help correct reported malaria case numbers to obtain more accurate measures of disease burden. The assembled and modelled data demonstrated that while treatment-seeking rates have overall increased over time, access remains low in some malaria endemic regions and utilization of government services is in some areas limited.
Although millions are affected by clinical malaria each year, the last 15 years have seen unprecedented gains from international efforts made to control this disease. The establishment of the Roll Back Malaria initiative and the Millennium Development Goals in the year 2000 were followed by a nearly 20-fold increase in international funding for malaria control . The scale-up of interventions that followed has resulted in a 40 % decline in Plasmodium falciparum clinical incidence in Africa while prevalence of the infection has nearly halved since the year 2000 . This and marked reductions in malaria-associated deaths  were largely attributed to increased coverage of the insecticide-treated bed nets (ITNs). Second to ITNs, which were the most widespread intervention, access to artemisinin-based combination therapy (ACTs) has been found to greatly impact the incidence of disease .
Prompt diagnosis and treatment of clinical malaria is the mainstay of all control or elimination programs [4–6]. The primary aim of treatment with ACT is to curtail clinical disease in patients, though access to effective treatment also impacts onward transmission to the wider community by reducing the infectious reservoir while at the same time containing the spread of drug resistance . Treatment coverage has been assessed across Africa , but understanding how access to treatment varies throughout the malaria endemic world is essential to evaluating the true public health impact of treatment coverage. This is all the more significant during the current transition from control of clinical disease towards regional focuses on malaria elimination and the post-2015 future is shaped as progress towards the Millennium Development Goals is evaluated .
In assessing care-seeking behaviours, it is also important to consider what proportion of care is sought at government-based facilities. Government facilities are more likely to comply with recommended diagnostic and treatment schedules [1, 9], and to have their routine records integrated into the national health data management system. Even in areas with strong surveillance systems, reports of passively detected malaria cases will capture only a certain fraction of all malaria cases [10–13] and so must be adjusted by a number of parameters before use as official burden estimates [1, 14]. These include (1) treatment-seeking behaviour (representing the proportion of cases not attending health facilities and thus being omitted from aggregated case reports [15–22] as well as the proportion seeking treatment outside the public health system); (2) malaria diagnoses made presumptively without parasitological confirmation (leading to reported case numbers including non-malaria illnesses [23, 24]); and (3) incomplete reporting (which leads to cases being lost from reported data) [1, 25]. Here we aim to improve the evidence-base of the first parameter to enable refined estimates of the true clinical burden of malaria disease.
Treatment-seeking rates vary widely between countries and greatly affect final burden estimates. Where available, these parameters are drawn from nationally representative, cross-sectional, household surveys such as Demographic and Health Surveys (DHS) , Malaria Indicator Surveys (MIS)  and UNICEF Multiple Indicator Cluster Surveys (MICS) . However, not all malaria-endemic countries (MECs) have such survey data available, resulting in these important parameters being either assumed or omitted and making comparisons of access to treatment across all malaria endemic regions and effective estimation of the burden of disease difficult.
The aim of this study was an exhaustive assembly of the available data on treatment-seeking for all MECs, sub-nationally where possible. For MECs lacking national survey data, predictive models were built to estimate use of public facility treatment as well as treatment of any kind. These estimates, including measurements of their uncertainty, will allow for improved understanding of how health services in endemic countries are routinely accessed and facilitate more accurate disease burden estimation.
Many DHS, MIS and MICS ask questions to determine the prevalence and treatment of fever in children less than 5 years of age. Here it was assumed that the treatment-seeking rates observed in children would be similar in older age groups, as supported by evidence from India, Indonesia and Ethiopia [16, 25]. For a handful of DHS surveys (n = 4) and earlier MICS rounds (MICS 3 and 4, n = 45), data on treatment-seeking for fever were not available, so data on treatment-seeking for acute respiratory infection (ARI) were used, which have previously been shown to correlate strongly with treatment-seeking for fever . For surveys that contained treatment-seeking data, the survey codes for questions regarding where treatment was sought were reviewed and categorized into public/government facilities (likely to have been captured by reporting systems) or ‘any’ medical treatment, which included private or NGO facilities, but excluded non-medical categories of care such as homeopathic doctors or ‘healers’. To generate a comprehensive database of treatment-seeking information, no time restriction was placed on the year of the survey. However, only MICS conducted in the third survey round or later (from 2005) were included due to inconsistency in interview questions on where treatment was sought within earlier surveys.
From each survey containing data on fever (or ARI) treatment-seeking, the total number of children reported to have fever was summarized by cluster and region. The numbers seeking treatment at either category of facility were also totalled to obtain the proportion of those who seek public/government facility-based treatment or any treatment. These data extractions from the DHS and MICS websites [26, 28] and summaries were all automated using Feature Manipulation Engine (FME) version 2015 by Safe Software . MIS data were extracted from the DHS platform and will be discussed as part of the DHS survey from here onward. The total number of children, fever cases and cases that sought treatment were then summarized nationally in order to generate predictions for those countries where survey data were unavailable. National survey data require the incorporation of sampling weights, which are provided in the survey results and help to adjust for differences in the probability of selection to produce a more accurate representation of population-based metrics. Sampling weights were extracted and applied at the individual level following DHS guidelines , such that the totals and proportions reported here were consistent with the same metrics (percentage of fever or illness) provided in the official survey reports.
Potential covariates were identified for inclusion in a predictive model using a review of the literature for ‘treatment-seeking’ and ‘care-seeking’ for both ‘malaria’ and ‘fever’ in PubMed on 18 July 2015 . The following have been reported as determinants of care-seeking rates: household wealth [15, 17, 19, 20, 22, 32–38], care-giver education [15, 19, 32, 34, 39] and household location (rural or urban) or access to health facilities [15, 20–22, 32, 33, 40–45]. The World Bank provides freely available national-level indicator data , and several indicators that were in keeping with the themes identified by the literature were downloaded: access to electricity (as a proxy for wealth and access to health facilities), gross domestic product (GDP; current US$), GDP per capita, GDP growth (annual %), gross national income (GNI) per capita (current US$), total health expenditure (% of GDP), public health expenditure (% of total), primary education completion rate (% of relevant age group), and rural population (% of total population). In addition to these metrics, the number of health workers was thought to influence treatment-seeking rates. However, there were not sufficient data on the number of community health workers per 1000 people available from the World Bank, so the number of nurses and midwives per 1000 population was used instead. Finally, explicit care-seeking variables were also included: the percentage of pregnant women receiving prenatal care and the percentage of children aged 12–23 months who were immunized against diphtheria, pertussis and tetanus (DPT).
Matching covariate data to national-level survey data was also achieved using FME. When available, covariate data were used from the same year as the survey, and otherwise were matched to the closest year. Covariate data were also matched to those countries without treatment data available. For these countries, covariate data was matched to 2013, which was the most recent year of indicator data available, or the closest year to that in order to generate the most up to date estimates of treatment-seeking rates.
Two generalized additive mixed models (GAMMs) were developed to predict the proportion of patients that sought treatment from (1) a facility covered by the government reporting system, and, (2) any type of medical care provider [47, 48]. Where multiple surveys were used from the same country, this was accounted for by adding a country-level random effects term. Both year and WHO region were included as fixed effect terms, along with the suite of country-level covariates. The list of the covariates described above was limited after testing for collinearity among all the potential variables. The variables included in GAMMs were: year of survey, WHO region, GDP growth, health expenditure, prenatal care rates, primary education, DPT immunization rates, nurses and midwives per population and proportion rural population. All statistical analyses were performed in the R statistical computing environment . Full details of the model development are provided in Additional file 1.
Existing treatment-seeking data were mapped at the regional level (boundaries shown in Fig. 1) for the 76 countries with available survey data. Although some of these data were available at a point level, only regional maps were generated for consistency with the other adjustment parameters. For the few countries with DHS cluster data and no regional data available, cluster level data were summarized and mapped to the first- (Colombia) or second-level (Zambia) administrative units . Predicted and observed values were then combined into a single, geographically complete, map of treatment-seeking for all MECs showing the most recent available data in each country.
Treatment-seeking data were collected from DHS (n = 195, 13 of which were MIS) and MICS (n = 59) for 76 countries. The number and year of surveys available from MECs are shown in Additional file 1: Figure S1. Data from DHS and MICS round 5 collected data on treatment-seeking for fever, whereas MICS rounds 3 and 4 assess care-seeking for cough (respiratory infection). Additional file 1: Figure S2 confirms that rates observed from the earlier MICS rounds are comparable to the treatment-seeking for subsequent surveys. Indicator data were available for all but nine surveys (there was no health expenditure data available for Zimbabwe and Somalia); these nine surveys could not therefore be included in the analysis. DHS surveys conducted prior to 1990 were excluded as they reported that no individuals sought treatment, and were therefore considered to be non-representative. Two MICS surveys, Yemen (2006) and Burkina Faso (2006), also reported zero treatment-seeking and were excluded. The GAMMs were thus fitted to 228 records from 72 countries. This left 22 MECs for which predictions of treatment-seeking outcomes were required. However, Brazil and Paraguay were also added to the prediction list since Brazil only had national level data available (DHS 1996) and Paraguay had large areas of missing data (DHS 1990). This increased the number of prediction countries to 24.
All observed and predicted values, as well as covariate data, are shown in Additional file 2. There were five countries for which insufficient covariate data precluded predictions. These countries were therefore assigned estimates from a similar and ideally neighbouring country in the same region (rounded to two decimal places to show that they are less precise). Myanmar (Burma) was matched to Timor-Leste, which had the most similar health indicator variables (DPT and pregnant women care). In the absence of indicator data or neighbouring countries in the same WHO region the Democratic People’s Republic of Korea (DPRK) was given the same estimate as Myanmar (Burma). Eritrea was given the same value of neighbouring Ethiopia, French Guiana was given Suriname’s values, and the Republic of Korea was given China’s estimate, which was closest to Korea for the indicators that were available.
Despite variability between geographic regions, treatment-seeking rates have generally increased over time. The non-linear relationship between treatment-seeking and time is shown in Additional file 1: Figure S3. The lowest rates of access to government treatment (<10 %) were observed prior to the year 2000, with the exception of findings from DHS surveys from Chad (2004), Bangladesh (2004 and 2007) and Pakistan (2007) and a MICS from Somalia in 2006 which were all also <10 %. The lowest reported rate of seeking any treatment was 24 % in Togo in 2006. All observations of treatment-seeking rates less than 40 % (from any health facility) were from the African continent (including Somalia, from EMRO). The highest measures of access to public facilities (>65 %) were observed in several regions: AFRO-C (Nigeria and Sierra Leone), AFRO-E (Zambia and Mozambique), PAHO (Costa Rica, Dominica Republic and Peru) and SEARO (Bangladesh (in 2011, but not in the earlier years referenced above), Bhutan, Indonesia and Timor-Leste). Bhutan had the highest reported measure of 74 % of individuals seeking care at government-based facilities (MICS 2010). High rates of access to any treatment type likewise spanned various regions, but values were far greater than public-facility data only. Regions with observations of access rates to any treatment greater than 80 % include: AFRO-C (Nigeria), AFRO-E (Tanzania and Uganda), EMRO (Pakistan), PAHO (Argentina, Belize, Costa Rica, Dominican Republic and Panama), SEARO (Indonesia (since 1991) and Thailand) and WPRO (Cambodia and Viet Nam).
Weighted means of treatment-seeking rates by WHO region
Information regarding the treatment-seeking behaviours of fever cases in MECs is essential to assessing the feasibility and success of malaria control and elimination programs. Malaria is a treatable disease and while effective therapies exist, population health-seeking behaviours may limit the extent to which they are utilized. Understanding care-seeking rates in MECs also helps to quantify the scale of the malaria burden. Routine surveillance data are used to measure the burden of clinical disease and only those cases who seek care at a government-based facility are likely to be included in regional estimates. Quantification of the proportion of cases that are severe enough to seek care (clinical), but are missed from passive surveillance, is therefore necessary to more accurately estimate ‘true’ case numbers.
Data obtained from national surveys (DHS and MICS) could be assembled for the majority of the MECs (78 %, n = 76). Treatment seeking rates specific to fevers were available from all DHS surveys and MICS 5, while only cough-based treatment rates were available from the earlier MICS surveys. Precedent exists for using care-seeking for respiratory infection as a proxy for fever treatment [14, 25], and this was evidenced by comparable treatment-seeking rates for cough and fever in the observed data gathered here (Additional file 1: Figure S2). Therefore, all available treatment-seeking data were used to inform the model. However, because fever is a primary symptom of malaria and because DHS data are georeferenced and could be mapped sub-nationally, DHS data were given precedence when generating final mapped outputs and summary estimates of the data gathered. Sub-national predictions and mapping of the treatment-seeking outcomes using cluster-level data to produce smooth surfaces like those produced by the Malaria Atlas Project for prevalence was also explored [2, 53, 54]. However, there were not sufficient covariate data available at smooth resolutions at the time of this analysis to support this. The quantity and quality of higher resolution sub-national covariate data that can be used in geostatistical analyses continue to improve with time and there may be greater potential for this type of analysis in the future .
Gathering treatment-seeking data for all MECs was the primary aim of this analysis, with the intent to also show sub-national rates where possible. Data that corresponded to the covariate variables identified in the literature review of key factors determining treatment seeking behaviour were readily available from the World Bank and produced models to show that government treatment-seeking and any treatment-seeking could be predicted at the national level from a limited set of covariate variables: year, WHO region, percent of pregnant women that receive prenatal care, immunization rates, primary education completion rate, GDP growth, and national health expenditure (public and total). The percentage of women receiving prenatal care was a strong indicator of fever treatment-seeking that likely drove the low model error values observed. The differing drivers of seeking care from either source were evidenced through the difference in the best models for each treatment-seeking outcome. Government-based treatment seeking was predicted primarily by other health-seeking indicators such as the childhood immunization (DPT) and prenatal care. Access to any treatment, on the other hand, was also influenced by country wealth and education. Educated individuals in more economically stable countries are therefore more likely to spend money on health care outside of the government system.
From the assembled observed data and modelled missing data, there emerged geographic patterns in both the outcomes and the certainty of the predictions. There were areas with low access or use of public treatment facilities in all regions of the malaria-endemic world. Figure 3 highlights areas such as Central Africa and Indian sub-continent. Accessing treatment of any kind was inherently higher in all countries because government facilities are included in that metric. However, the any treatment data revealed that treatment-seeking in some endemic areas, such as India, Pakistan and Afghanistan, was largely pursued outside the public sector. The CI ranges (Fig. 2) show that outcomes were well predicted in the Americas, but less so in Asia and least accurate in the Eastern Mediterranean countries. This implies that models and indicator variables were better suited to the Central and South American countries. Future predictions of this nature may be improved by including additional covariates or, following further research into treatment-seeking indicators, parameters that are tailored to each region.
The comparison of the two treatment-seeking outcomes and the measures of uncertainty provide valuable information as countries define control and elimination goals. The predicted measurements resulted in the greatest uncertainty, and signal the need for treatment-seeking to be formally assessed in those regions. Most notably, this data assembly and analysis reveals parts of the malaria endemic world where treatment is primarily sought outside of government programs. In light of concerns regarding the spread of antimalarial resistance [56, 57], it is essential for countries to ensure that cases are being diagnosed and treated properly using approved and legitimate drugs . If treatment is most commonly sought outside of government-based facilities, control programs must consider how best to monitor treatment safety and efficacy as well as numbers of cases presenting for treatment.
Information on treatment-seeking behaviours in malaria endemic countries can be readily assembled from national survey data. Where data on treatment-seeking behaviours were not available from national surveys, modelling techniques using freely available data were applied to fill data gaps. Both the results and methods presented have potential application beyond those described here and may inform the control and burden of other febrile diseases. However, in this context, data on treatment-seeking for fever are essential to understanding the efficacy with which malaria cases are treated and detected. Gathering and visualizing these data for all MECs, sub-nationally when possible, is of use to estimate the burden in areas of low endemicity where passive surveillance is the primary tool through which cases are monitored. These results will facilitate downstream efforts to produce a hybridized burden estimation approach that employs both surveillance-based and cartographic techniques in an effort to more accurately quantify the global burden of falciparum and vivax malarias and provide immediate feedback regarding parts of the malaria endemic world where treatment for malaria is not readily accessed or is more commonly sought beyond the government or control programme sectors.
KEB and PWG conceived the study and oversaw its implementation, with input from SIH and REH. KEB wrote the first draft of the manuscript and assembled treatment-seeking and indicator data with assistance from HG, DW, BM, and UD. DB assisted the design of the modelling framework with input from SB and EC. All authors read and approved the final manuscript.
KEB and REH are supported by a Senior Research Fellowship from the Wellcome Trust (095066) awarded to SIH. PWG is a Career Development Fellow (#K00669X) jointly funded by the UK Medical Research Council (MRC) and the UK Department for International Development (DFID) under the MRC/DFID Concordat agreement and receives support from the Bill & Melinda Gates Foundation (#OPP1068048, #OPP1106023), which also supports SB, EC, DW, HG, BM, and UD.
The authors declare that they have no competing interests.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 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.
- WHO. World malaria report 2014. Geneva: World Health Organization; 2014. p. 226.Google Scholar
- Bhatt S, Weiss DJ, Cameron E, Bisanzio D, Mappin B, Dalrymple U, et al. The impact of malaria control on Plasmodium falciparum in Africa, 2000–2015. Nature. 2015;526:207–11.View ArticlePubMedGoogle Scholar
- Eisele TP, Larsen DA, Walker N, Cibulskis RE, Yukich JO, Zikusooka CM. Estimates of child deaths prevented from malaria prevention scale-up in Africa 2001–2010. Malar J. 2012;11:93.View ArticlePubMedPubMed CentralGoogle Scholar
- Greenwood BM. Control to elimination: implications for malaria research. Trends Parasitol. 2008;24:449–54.View ArticlePubMedGoogle Scholar
- Tanner M, Greenwood B, Whitty CJ, Ansah EK, Price RN, Dondorp AM. Malaria eradication and elimination: views on how to translate a vision into reality. BMC Med. 2015;13:167.View ArticlePubMedPubMed CentralGoogle Scholar
- Tanner M, Vlassoff C. Treatment-seeking behaviour for malaria: a typology based on endemicity and gender. Soc Sci Med. 1998;46:523–32.View ArticlePubMedGoogle Scholar
- Galactionova K, Tediosi F, de Savigny D, Smith T, Tanner M. Effective coverage and systems effectiveness for malaria case management in sub-Saharan African countries. PLoS One. 2015;10:e0127818.View ArticlePubMedPubMed CentralGoogle Scholar
- Murray CJ, Ortblad KF, Guinovart C, Lim SS, Wolock TM, Roberts DA, et al. Global, regional, and national incidence and mortality for HIV, tuberculosis, and malaria during 1990–2013: a systematic analysis for the Global Burden of Disease Study 2013. Lancet. 2014;384:1005–70.View ArticlePubMedPubMed CentralGoogle Scholar
- WHO. Guidelines for the treatment of malaria. Third ed. Geneva: World Health Organization; 2015.Google Scholar
- Gething PW, Noor AM, Gikandi PW, Ogara EA, Hay SI, Nixon MS, et al. Improving imperfect data from health management information systems in Africa using space-time geostatistics. PLoS Med. 2006;3:e271.View ArticlePubMedPubMed CentralGoogle Scholar
- Hay SI, Okiro EA, Gething PW, Patil AP, Tatem AJ, Guerra CA, et al. Estimating the global clinical burden of Plasmodium falciparum malaria in 2007. PLoS Med. 2010;7:e1000290.View ArticlePubMedPubMed CentralGoogle Scholar
- Murray CJ, Lopez AD, Wibulpolprasert S. Monitoring global health: time for new solutions. BMJ. 2004;329:1096–100.View ArticlePubMedPubMed CentralGoogle Scholar
- Snow RW. Sixty years trying to define the malaria burden in Africa: have we made any progress? BMC Med. 2014;12:227.View ArticlePubMedPubMed CentralGoogle Scholar
- Cibulskis RE, Aregawi M, Williams R, Otten M, Dye C. Worldwide incidence of malaria in 2009: estimates, time trends, and a critique of methods. PLoS Med. 2011;8:e1001142.View ArticlePubMedPubMed CentralGoogle Scholar
- Chaturvedi HK, Mahanta J, Pandey A. Treatment-seeking for febrile illness in north-east India: an epidemiological study in the malaria endemic zone. Malar J. 2009;8:301.View ArticlePubMedPubMed CentralGoogle Scholar
- Deressa W. Treatment-seeking behaviour for febrile illness in an area of seasonal malaria transmission in rural Ethiopia. Malar J. 2007;6:49.View ArticlePubMedPubMed CentralGoogle Scholar
- Deressa W, Ali A, Berhane Y. Household and socioeconomic factors associated with childhood febrile illnesses and treatment seeking behaviour in an area of epidemic malaria in rural Ethiopia. Trans R Soc Trop Med Hyg. 2007;101:939–47.View ArticlePubMedGoogle Scholar
- Espino F, Manderson L. Treatment seeking for malaria in Morong, Bataan, the Philippines. Soc Sci Med. 2000;50:1309–16.View ArticlePubMedGoogle Scholar
- Fenny AP, Asante FA, Enemark U, Hansen KS. Malaria care seeking behavior of individuals in Ghana under the NHIS: are we back to the use of informal care? BMC Public Health. 2015;15:370.View ArticlePubMedPubMed CentralGoogle Scholar
- Ilunga-Ilunga F, Leveque A, Ngongo LO, Laokri S, Dramaix M. Treatment-seeking paths in the management of severe malaria in children under 15 years of age treated in reference hospitals of Kinshasa, Democratic Republic of Congo. Trop Med Health. 2015;43:11–9.View ArticlePubMedPubMed CentralGoogle Scholar
- Thandar MM, Kyaw MP, Jimba M, Yasuoka J. Caregivers’ treatment-seeking behaviour for children under age five in malaria-endemic areas of rural Myanmar: a cross-sectional study. Malar J. 2015;14:1.View ArticlePubMedPubMed CentralGoogle Scholar
- Xu JW, Xu QZ, Liu H, Zeng YR. Malaria treatment-seeking behaviour and related factors of Wa ethnic minority in Myanmar: a cross-sectional study. Malar J. 2012;11:417.View ArticlePubMedPubMed CentralGoogle Scholar
- Jensen TP, Bukirwa H, Njama-Meya D, Francis D, Kamya MR, Rosenthal PJ, et al. Use of the slide positivity rate to estimate changes in malaria incidence in a cohort of Ugandan children. Malar J. 2009;8:213.View ArticlePubMedPubMed CentralGoogle Scholar
- Petti CA, Polage CR, Quinn TC, Ronald AR, Sande MA. Laboratory medicine in Africa: a barrier to effective health care. Clin Infect Dis. 2006;42:377–82.View ArticlePubMedGoogle Scholar
- WHO. World malaria report 2008. Geneva: World Health Organization; 2008. p. 215.Google Scholar
- The DHS Program. http://www.dhsprogram.com/. Accessed 22 June 2015.
- Malaria Indicator Surveys. http://www.malariasurveys.org/. Accessed 22 June 2015.
- Multiple Indicator Cluster Surveys. http://mics.unicef.org/. Accessed 22 June 2015.
- Safe Software: Feature Manipulation Engine. Surrey, BC Canada; 2015.Google Scholar
- Rutstein SO, Rojas G. Guide to DHS statistics. DHS Toolkit of methodology for the MEASURE DHS Phase III project, implemented from 2008–2013. Calverton: USAID; 2006. p. 168.Google Scholar
- PubMed. http://0-www.ncbi.nlm.nih.gov.brum.beds.ac.uk/pubmed/. Accessed 18 July 2015.
- Hossain SM, Bhuiya A, Rasheed S. Correlates of perceived malarial episodes and treatment-seeking behavior in a malaria-endemic rural area in Bangladesh. Southeast Asian J Trop Med Public Health. 2001;32:707–19.PubMedGoogle Scholar
- Mustafa HS, Malik EM, Tuok HT, Mohamed AA, Julla AI, Bassili A. Malaria preventive measures, health care seeking behaviour and malaria burden in different epidemiological settings in Sudan. Trop Med Int Health. 2009;14:1488–95.View ArticlePubMedGoogle Scholar
- Nabyonga Orem J, Mugisha F, Okui AP, Musango L, Kirigia JM. Health care seeking patterns and determinants of out-of-pocket expenditure for malaria for the children under-five in Uganda. Malar J. 2013;12:175.Google Scholar
- Najnin N, Bennett CM, Luby SP. Inequalities in care-seeking for febrile illness of under-five children in urban Dhaka, Bangladesh. J Health Popul Nutr. 2011;29:523–31.View ArticlePubMedPubMed CentralGoogle Scholar
- Rutebemberwa E, Kallander K, Tomson G, Peterson S, Pariyo G. Determinants of delay in care-seeking for febrile children in eastern Uganda. Trop Med Int Health. 2009;14:472–9.View ArticlePubMedGoogle Scholar
- Smith LA, Bruce J, Gueye L, Helou A, Diallo R, Gueye B, et al. From fever to anti-malarial: the treatment-seeking process in rural Senegal. Malar J. 2010;9:333.View ArticlePubMedPubMed CentralGoogle Scholar
- Sumba PO, Wong SL, Kanzaria HK, Johnson KA, John CC. Malaria treatment-seeking behaviour and recovery from malaria in a highland area of Kenya. Malar J. 2008;7:245.View ArticlePubMedPubMed CentralGoogle Scholar
- Klein EY, Lewis IA, Jung C, Llinas M, Levin SA. Relationship between treatment-seeking behaviour and artemisinin drug quality in Ghana. Malar J. 2012;11:110.View ArticlePubMedPubMed CentralGoogle Scholar
- Ewing VL, Lalloo DG, Phiri KS, Roca-Feltrer A, Mangham LJ, SanJoaquin MA. Seasonal and geographic differences in treatment-seeking and household cost of febrile illness among children in Malawi. Malar J. 2011;10:32.View ArticlePubMedPubMed CentralGoogle Scholar
- Holtz TH, Kachur SP, Marum LH, Mkandala C, Chizani N, Roberts JM. Care seeking behaviour and treatment of febrile illness in children aged less than five years: a household survey in Blantyre District, Malawi. Trans R Soc Trop Med Hyg. 2003;97:491–7.View ArticlePubMedGoogle Scholar
- Kassile T, Lokina R, Mujinja P, Mmbando BP. Determinants of delay in care seeking among children under five with fever in Dodoma region, central Tanzania: a cross-sectional study. Malar J. 2014;13:348.View ArticlePubMedPubMed CentralGoogle Scholar
- Matovu F, Nanyiti A, Rutebemberwa E. Household health care-seeking costs: experiences from a randomized, controlled trial of community-based malaria and pneumonia treatment among under-fives in eastern Uganda. Malar J. 2014;13:222.View ArticlePubMedPubMed CentralGoogle Scholar
- Muller O, Traore C, Becher H, Kouyate B. Malaria morbidity, treatment-seeking behaviour, and mortality in a cohort of young children in rural Burkina Faso. Trop Med Int Health. 2003;8:290–6.View ArticlePubMedGoogle Scholar
- Tynan A, Atkinson JA, Toaliu H, Taleo G, Fitzgerald L, Whittaker M et al. Community participation for malaria elimination in Tafea Province, Vanuatu: part II. Social and cultural aspects of treatment-seeking behaviour. Malar J. 2011;10:204.Google Scholar
- The World Bank. Data: indicators. http://data.worldbank.org/indicator/all. Accessed 23 June 2015.
- Wood SN. Stable and efficient multiple smoothing parameter estimation for generalized additive models. J Am Stat Assoc. 2004;99:673–86.View ArticleGoogle Scholar
- Lin XH, Zhang DW. Inference in generalized additive mixed models by using smoothing splines. J Roy Stat Soc B. 1999;61:381–400.View ArticleGoogle Scholar
- R Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2013.Google Scholar
- FAO. The global administrative unit layers (GAUL): technical aspects. Rome: Food and Agriculture Organization of the United Nations, EC-FAO Food Security Programme (ESTG); 2013.Google Scholar
- NASA/CIESIN. Gridded Population of the World, Version 3 (GPWv3): Population Density Grid. Palisades, NY: NASA Socioeconomic Data and Applications Center (SEDAC); 2015.Google Scholar
- ESRI. ArcGIS Desktop 10.3. Redlands: Environmental Systems Resource Institute; 2014.Google Scholar
- Gething PW, Elyazar IR, Moyes CL, Smith DL, Battle KE, Guerra CA, et al. A long neglected world malaria map: Plasmodium vivax endemicity in 2010. PLoS Negl Trop Dis. 2012;6:e1814.View ArticlePubMedPubMed CentralGoogle Scholar
- Gething PW, Patil AP, Smith DL, Guerra CA, Elyazar IR, Johnston GL, et al. A new world malaria map: Plasmodium falciparum endemicity in 2010. Malar J. 2011;10:378.View ArticlePubMedPubMed CentralGoogle Scholar
- Weiss DJ, Mappin B, Dalrymple U, Bhatt S, Cameron E, Hay SI. Re-examining environmental correlates of Plasmodium falciparum malaria endemicity: a data-intensive variable selection approach. Malar J. 2015;14:68.View ArticlePubMedPubMed CentralGoogle Scholar
- Dondorp AM, Yeung S, White L, Nguon C, Day NP, Socheat D, et al. Artemisinin resistance: current status and scenarios for containment. Nat Rev Microbiol. 2010;8:272–80.View ArticlePubMedGoogle Scholar
- Price RN, von Seidlein L, Valecha N, Nosten F, Baird JK, White NJ. Global extent of chloroquine-resistant Plasmodium vivax: a systematic review and meta-analysis. Lancet Infect Dis. 2014;14:982–91.View ArticlePubMedPubMed CentralGoogle Scholar
- Tabernero P, Fernandez FM, Green M, Guerin PJ, Newton PN. Mind the gaps—the epidemiology of poor-quality anti-malarials in the malarious world–analysis of the WorldWide Antimalarial Resistance Network database. Malar J. 2014;13:139.View ArticlePubMedPubMed CentralGoogle Scholar