- Open Access
Spatial models for the rational allocation of routinely distributed bed nets to public health facilities in Western Kenya
Malaria Journalvolume 16, Article number: 367 (2017)
In high to moderate malaria transmission areas of Kenya, long-lasting insecticidal nets (LLINs) are provided free of charge to pregnant women and infants during routine antenatal care (ANC) and immunization respectively. Quantities of LLINs distributed to clinics are quantified based on a combination of monthly consumption data and population size of target counties. However, this approach has been shown to lead to stock-outs in targeted clinics. In this study, a novel LLINs need quantification approach for clinics in the routine distribution system was developed. The estimated need was then compared to the actual allocation to identify potential areas of LLIN over- or under-allocation in the high malaria transmission areas of Western Kenya.
A geocoded database of public health facilities was developed and linked to monthly LLIN allocation. A network analysis approach was implemented using the location of all public clinics and topographic layers to model travel time. Estimated travel time, socio-economic and ANC attendance data were used to model clinic catchment areas and the probability of ANC service use within these catchments. These were used to define the number of catchment population who were likely to use these clinics for the year 2015 equivalent to LLIN need. Actual LLIN allocation was compared with the estimated need. Clinics were then classified based on whether allocation matched with the need, and if not, whether they were over or under-allocated.
888 (70%) public health facilities were allocated 591,880 LLINs in 2015. Approximately 682,377 (93%) pregnant women and infants were likely to have attended an LLIN clinic. 36% of the clinics had more LLIN than was needed (over-allocated) while 43% had received less (under-allocated). Increasing efficiency of allocation by diverting over supply of LLIN to clinics with less stock and fully covering 43 clinics that did not receive nets in 2015 would allow for complete matching of need with distribution.
The proposed spatial modelling framework presents a rationale for equitable allocation of routine LLINs and could be used for quantification of other maternal and child health commodities applicable in different settings. Western Kenya region received adequate LLINs for routine distribution in line with government of Kenya targets, however, the model shows important inefficiencies in the allocation of the LLINs at clinic level.
Long-lasting insecticidal nets (LLINs) are distributed as the main tool of vector control for malaria prevention in Africa, through mass campaigns and routine distribution systems . Mass campaigns, target population at community level in a single time-limited operation every 3 years. Routine systems deliver LLINs through public health sector all year round to sustain coverage in the interval between mass campaigns to pregnant women and infants who are the most vulnerable. The World Health Organization (WHO) recommends both channels for countries in order to achieve and maintain universal coverage [1, 2].
At the launch of Kenya’s national malaria strategy in 2001, the aim was to achieve 60% insecticide-treated nets (ITNs) coverage of at risk children and pregnant women by 2006 . Prior to this, access to ITNs was limited to the private sector, a few research projects and non-governmental organizations. Different channels of ITNs delivery were explored including commercial retail sector ITNs and subsidized ITNs through public health sector. These approaches did not succeed in reaching the rural poor and achieving maximum coverage of the population at risk . In 2004, free routine LLIN distribution in antenatal care clinics (ANC) began and in 2006 the first free mass-campaign was conducted. Subsequent free mass campaigns were conducted in 2011, 2012, 2014 and 2015 in malaria endemic areas of Kenya .
Routine distribution of LLINs in Kenya are targeted to high and moderate malaria transmission areas, where pregnant women are provided with a free LLIN during their first ANC visit so that the mother and the unborn baby are protected at the earliest chance. Infants are provided with a free LLIN during expanded programme on immunization (EPI) [6, 7]. There is no specific WHO recommendation as to when LLINs should be distributed during EPI. However, a range of time-points is used in practice, from birth (tuberculosis vaccine) to 9 months (measles vaccine). Diphtheria-tetanus-pertussis-1 vaccination at 6 weeks is the most common distribution point [7, 8]. ANC clinics that distribute LLINs have increased rapidly from 1000 nationally in 2004 to over 4000 in 2015, with a total 23.3 million LLINs distributed .
Pregnant women receive a package of interventions during ANC visits which play an important role in ensuring a healthy mother and baby during pregnancy and after delivery [6, 9]. Healthcare utilization for ANC services is associated with intrapersonal, institutional, health systems and social demographic factors [10, 11]. Similarly, the ownership and use of LLINs during pregnancy is associated with household wealth, distance to health facility, residence, parity, marital status, and education among others . It is key, therefore, to include these spatial variables to account for the uneven healthcare utilization rates within facility catchment areas within a spatial modelling framework [13, 14].
The quantity of LLINs distributed to clinics is quantified at the central level using average monthly LLINs consumption data and the overall size of the population in targeted counties. However, this approach has resulted in frequent stock-outs [7, 15] with subsequent low population coverage of LLIN among pregnant women and infants [8, 15, 16]. In this study, a spatial modelling framework was proposed and used to quantify LLINs need by ANC clinic to allow more efficient targeting. A combination of modelled spatial access to clinics, household survey data on utilization and high resolution population grids were used to quantify the population of pregnant women and infants in 2015 in need of LLINs. The estimated need was compared to the actual allocation in 2015 and clinics with fewer (under-allocation), more (over-allocation) and enough (matched-allocation) LLIN relative to the estimated need identified. A fourth category is those facilities with zero distribution (non-allocation). Approaches to minimize targeting inefficiencies were then proposed.
Kenya is stratified into five malaria zones to address varied risks of malaria (Fig. 1a). They include Lake endemic (high, stable perennial transmission areas), Coast endemic (moderate, seasonal transmission areas), seasonal low transmission (acutely seasonal, low transmission areas), Highland epidemic (unstable, variable transmission areas) and Low risk (malaria free or extremely low transmission) [17, 18]. The analysis was restricted to the Lake endemic zone (Fig. 1b), because of its high, stable transmission throughout the year and presence of routine LLIN distribution to pregnant women and infants since 2004. In this zone, Plasmodium falciparum prevalence has historically been greater than 20%. The zone comprises of eight counties; Busia, Siaya, Bungoma, Homa Bay, Kakamega, Kisumu, Migori, and Vihiga (Fig. 1b) .
Previously mapped health facilities  were updated using the Ministry of Health  and the District Health Information System version 2 (DHIS 2) lists . Lists were compared to eliminate duplicates. Health facilities managed by the government, faith-based and non-governmental organizations were selected as public health facilities. Health facilities that did not have coordinates were geo-located using Google Earth . Monthly LLIN distribution data by clinic were obtained from Population Services Kenya and linked to the geocoded database.
Livebirths and pregnancies distribution maps for 2015 at 100 × 100 m spatial resolution were downloaded from WorldPop data portal . These maps were constructed using dasymetric spatial modelling techniques that redistributed the Kenya’s national census (2009) population counts within 6603 sub-locations [24, 25]. Population was shifted away from areas unlikely to be inhabited using weights assigned to each land cover with national estimates projected to match 2015 UN national estimates . Population surfaces were grouped by sex and 5 year-age groups and adjusted using age-specific fertility rates for females of childbearing age . The resultant live births surface was adjusted to match the Guttmacher Institute estimates  and converted to pregnancies by adjusting national totals to match national estimates . The livebirths and pregnancy surfaces was aggregated to 300 m spatial resolution while all the other surfaces were resampled to the same resolution in ArcGIS version 10.1 (ESRI Inc., Redlands, CA, USA) for ease of computation and compatibility with the travel time surface which was modelled at 300 m spatial resolution.
Socio-economic, demographic and ANC attendance data for pregnant women were obtained from the 2014 Kenya Demographic and Health Survey (KDHS). This survey employed a two-stage sampling design on a national sampling frame of 5360 clusters. In the first stage, 1612 clusters (995 in urban areas and 617 in rural areas) were selected with equal probability while in the second stage, 40,300 households were selected from the household listing of the 1612 clusters. More than 19% (2776) of all women in the reproductive age (15–49 years) who had a live birth 5 years preceding the survey in Kenya, were residents of 277 clusters (181 in rural and 96 in the urban areas) in the 8 study counties, all geo-located .
In order to protect and maintain the confidentiality of the respondents in the KDHS, Global Positioning System coordinates (latitude and longitude) were randomly displaced up to 2 km for urban clusters and 5 km for rural clusters with 1% of the rural clusters displaced up to 10 km . In this study, the scrambling effect was minimized by drawing 2 km (urban) and 5 km (rural) Euclidean buffers around the cluster points, which were then shifted to their most probable location within the buffer [32,33,34]. This was based on a point’s elevation (not scrambled) provided along with the coordinates and a digital elevation model (DEM). To supplement the elevation information, Google Earth  and population distribution surfaces  were used to relocate the point in populated or urban location within the buffer.
Road network data from OpenStreetMaps (OSM)  and Google Map Maker (GMM)  was updated via Google Earth . Duplicates and short sections of roads disconnected from the main network were removed in the resultant road network. Roads were then classified based on the inherent classification on OSM, GMM, Google Earth, the Kenya Roads Act  and Bill . The Bill and the Act outlines the criteria for classifying public roads in Kenya including primary, secondary, county and rural roads (Table 1). The road network was converted into a raster surface. Land cover at 30 m spatial resolution was downloaded from China’s global mapping project  generated from Landsat and Environmental Disaster Alleviation Satellites imagery. DEM at 30 m spatial resolution by Shuttle Radar Topographic Mission was obtained from Earth Explorer portal .
Estimating LLINs need per clinic
A flowchart showing the analytical process is presented in Fig. 2. Travel time and catchment areas for clinics were modeled. Within each catchment, the total population and proportion of pregnant women and infants in 2015 were estimated. Probability surface of ANC utilization was then modeled and used to adjust catchment population per ANC clinic to define those who are most likely to use it. The catchment population of pregnant women and infants for 2015 was then extracted per clinic and these were used as equivalent to the estimated LLIN need per clinic. This was then compared to the 2015 actual LLIN allocation.
To implement the travel time analysis, different speeds were assigned on the road network depending on the likelihood of using motorized, bicycling and walking travel modes (Table 1). Speeds assigned to road types and land cover classes were adopted from comparable studies [13, 41,42,43]. Slope derived from DEM was used to adjust walking  and bicycling speeds [45, 46]. The roads and land use grids were then combined into a single travel time raster grid by adding up the time needed to cross contiguous cells to the nearest facility via the shortest path in ArcView (version 3.2) using AccessMod (version 3.0) at 300 m spatial resolution for computation ease [13, 47].
For the analysis of the probability of utilization of ANC clinics, travel times were extracted for each KDHS 2014 cluster location. A distance decay curve was fitted to predict the probability of attendance with travel time  using a three parameter logistic regression . A cut-off time within which most pregnant women accessed a clinic was determined through the curve inflection point. The cut-off time and the travel time grid were used to define facility catchment areas  using Path Distance tool in ArcGIS.
To select factors associated with ANC attendance among those extracted from KDHS 2014 (wealth quintiles, maternal education, parity, maternal age, residence, marital status) and the modelled time, a hierarchical mixed effects logistic regression model (Eq. 1) was implemented using “lme4” package [32, 33, 50]. The parsimonious model was defined by the model with the lowest Akaike Information Criterion value .
where P- is the probability of ANC attendance, α-the intercept; β- the vector of unknown regression coefficients for the fixed effects; X- a matrix of known covariates (wealth, education, parity, age, residence, marital status and travel time. Z- a matrix of random effect (clusters nested within counties); γ- a vector of variance. A binary outcome was defined; 1 (1+ ANC visits) and 0 (0 ANC visit).
Spatial structure of the selected variables was assessed using variogram while interpolation was implemented using ordinary kriging in ‘geoR’ package, with a 10% hold out for validation . The travel time and the interpolated surfaces were used to define an ANC attendance probability surface (Eq. 2) in R software.
(β)- are the coefficients of the matrix of known covariates (X) and the intercept (α) from Eq. 1. N is a limiting factor on y axis. It defines the probability of attendance at zero minutes and in presence of ideal circumstances for other influencing variables (e.g. higher wealth quantile and higher maternal education). However, other variables including but not limited to the unmeasurable and unmeasured variables limit the probability to less than 1. All the other symbols have similar meaning as those in Eq. 1.
The gridded maps of pregnant women and live births in 2015 were summed together and the population per clinic catchment extracted. In addition, the summed population was multiplied with the ANC probability surface to get those likely to have used a clinic. Those likely to have used a clinic were extracted per each clinic catchment using Zonal statistics tool in ArcGIS.
LLIN allocation efficiency at each clinic
The difference between total population and those likely to have used a clinic per catchment area on average was a 100 people (both pregnant women and infants). Consequently, over-allocation was defined as >100 LLINs allocation higher than the need, under-allocation as >100 LLINs allocation lower than the need and matched-allocation as ±100 LLINs relative to the need. Non-allocated areas were defined as those with clinics where no LLINs were distributed in 2015. Those likely to have attended a clinic were compared with 2015 actual LLINs allocation in the corresponding clinic to define, under-, over-, and matched-allocation. LLINs in over-allocated clinics were rationalized to the under-allocated clinics. Non-covered areas with a population of ≤100 were assigned to the nearest clinic within the existing distribution chain. Those with population >100, were assigned to a public clinic within the pre-determined cut-off time of the facilities not allocated LLINs in 2015.
The final list contained 1271 public health facilities out of which 888 had distributed LLINs in 2015. The analysis of the KDHS 2014 data showed that all 277 cluster points were within approximately 1-h of the nearest public clinic distributing LLINs (Fig. 3). ANC utilization decreased with increasing travel time, with attendance declining rapidly for distances greater than 40 min away from the nearest clinic on the decay curve. Catchment areas defined to include all those locations within 40 min of a health facility for the 888 public clinics distributing LLINs (exclusive of 13 clinics in islands of Lake Victoria) are shown in Fig. 4 and contained over 97% of pregnant women and infants in 2015 in the study area.
In the study counties, among women who had at least one live birth in the 5 years preceding the survey, over 97% had at least one ANC visit. As expected and exhaustively documented in literature within the Kenyan context [53,54,55], the final model of variables associated with ANC utilization showed that women of higher education and wealth quintile, those who were married, or had low parity or lived closer to a clinic were associated with higher odds of attending ANC (Table 2). All these variables varied spatially and this structure was used in interpolating them across the study area with 81–97% of pregnant women attending ANC clinics.
In 2015, 730,947 pregnant women and infants resided in the eight study counties with 97% (706,450) residing within modelled clinics catchment areas. The estimated LLINs need per clinic for 2015 was between 15 and 16,900, with 48% having a need of less than 500 while 23% with a need of more than 1000 LLINs annually. Among pregnant women and infants in catchment areas, 97% (682,377) were likely to have attended a clinic distributing LLINs representing 93% of the total population of pregnant women and infants in the study area.
In 888 public health clinics, 591,880 LLINs were allocated in 2015. The number of LLINs allocated per clinic was between 10 and 6720. 53% of the clinics had less than 500 LLINs allocation per clinic while 19% had allocations greater than 1000 LLINs. Relative to an estimated overall need for 682,377 LLINs, the un-met need for 2015 was 90,497. The model suggested that 380 (43%) clinics were under-allocated by 255,628 LLINs (Fig. 4). The under-allocation per clinic was between 100 and >2000 LLINs (Fig. 5b). In 321 (36%) clinics (Fig. 4), 164,241 LLINs were over-allocated by up-to over 2000 LLINs (Fig. 5a) while 21% (187) of the clinics had a matched-allocation accounting for 86,990 LLINs (Fig. 4).
Using the spatial modelling framework to simulate an efficient LLIN allocation system, the 164, 241 over-allocated LLINs were redistributed reducing the shortfall from 255,628 to 91,387. Of the 95 areas that were outside 40 min and thus not covered by the LLIN distribution in 2015 (Fig. 4), 34 had a population of ≤100 and could be feasibly absorbed within the existing neighboring clinics resulting in increasing allocation of about 1606 LLINs to these clinics. 43 clinics not allocated LLINs in 2015, were within the determined 40 min cut-off time and would cater for 43 of remaining polygons equivalent to 17,703 extra LLINs. The remaining un-covered areas (18 polygons) all located between 40 and 60 min of the closest clinic distributing LLINs would require 5187 LLINs and accounts for less than 1% of those likely to attend a clinic.
A spatial modelling framework that can be adapted for the rational allocation of routinely distributed LLINs in Kenya was developed. LLIN need at each clinic was quantified by using a combination of modelled accessibility, population and ANC utilization and compared this to the actual allocations in 2015. Overall, the unmet need of LLINs for all pregnant women and infants and those likely to attend public clinics was 139,067 and 90,497 respectively. The current quantification approach of LLINs need covers approximately 88% of those likely to attend ANC clinics and about 81% of all pregnant women and infants in the study area. This shows that the routine distribution system in Western Kenya has a wide reach and is consistent with the Government of Kenya target of at least 80% of at risk population using appropriate malaria preventive intervention . However, the facility level analysis also shows that there are important inequities in the distribution of LLINs to pregnant women and infants. The model provides a means for reducing such inequities and can be upgraded as consumption data improves.
Currently, the Ministry of Health relies on consumption data reported from health facilities and adjusted for estimated census projections of pregnant women and infants to quantify routine LLIN distribution. At delivery points, mechanisms for rational distribution have been rolled out including verification that a client qualifies for an LLIN, stamping the recipients’ booklet to avoid repeat collection from the same or another clinic and improved overall training of health workers on the management of commodities, regular auditing at the clinics and monthly reports. Recently, with the introduction of DHIS 2 [21, 57] in Kenya, LLINs consumption data are now reported through this system. This could potentially improve rational allocation, especially if predefined catchment areas and population maps such as those developed in this analysis are used as additional information to check whether distributions match the need for LLINs.
The definition of catchment areas could either be theoretical and based on simple proximity of the population to a health facility [58, 59] or could be improved with the available empirical data on utilization of services derived from household survey data rates [13, 41]. Previous studies have used only travel time and its relationship with reported utilization to define the probability of service use [13, 32]. However, in this study, other determinants of utilization such as wealth and maternal education, were used to capture both the spatial and social, demographic and economic factors that influence the use of ANC clinics. Arguably, this approach offers a more complete framework for quantifying the catchment population (likely ANC clients) with a broader catchment area .
Some important limitations are worth noting in the modelling framework used. First, the use of the nearest clinic when modelling travel times was assumed. This is because household survey data does not allow for the identification of the actual facility used by women and infants. It is likely that a proportion of these populations often use facilities that are not the closest for a variety of reasons, including the quality of services on offer . In addition, all pregnant women were assumed to have an equal probability of getting an LLIN irrespective of the number of ANC visits (1+ ANC), however, pregnant women with 2 or more ANC visits may have an increased chance of receiving an LLIN in subsequent visits. Second, a threshold of 40 min was imposed in the definition of facility catchment areas but it is likely that this threshold is variable by facility and future models should account for them. Finally, there are uncertainties in the pregnancies and births density maps  used into compute the LLIN need. Live births were assumed to be are equivalent of infants, however, a neonatal mortality of 22 per 1000 live births in 2015 in Kenya  equivalent to 6429 deaths (2.2%) of all live births and 1.1% of all LLINs distributed in the study counties may empirically bias the estimation of need. However, the national policy planning makes allocation of live saving services on the basis that all avoidable deaths are averted and therefore services, such as LLINs, will not and should not be discounted on the basis on anticipated neonatal deaths. To ensure that all uncertainties are propagated within the modelling approach, Bayesian inference methods may provide a useful alternative, but these advantages must be balanced with the simplicity required in spatial modelling techniques in routine programme planning.
The spatial modelling approach presented here offers a rational basis for spatial allocation of routine LLINs that is applicable to malaria endemic settings and, potentially, other maternal and child related commodities. The analysis shows that with a few minor adjustments in distribution and better rationalization, equity in LLIN distribution through the routine system could be achieved in Western Kenya where the threat of malaria is the greatest.
long-lasting insecticidal nets
- DHIS 2:
District Health Information System (version 2)
- KDHS 2014:
Kenya Demographic and Health Survey for the year 2014
digital elevation model
World Health Organization (WHO). Methods for maintaining coverage with long-lasting insecticidal nets (LLIN): Vector Control Technical Expert Group Report to MPAC. 2013. http://www.who.int/malaria/mpac/mpac_sp13_vcteg_universal_llin_coverage_report.pdf. Accessed 2016 Feb 27.
World Health Organization (WHO). Achieving universal coverage with long-lasting insecticidal nets in malaria control. 2014. http://www.who.int/malaria/publications/atoz/who_recommendation_coverage_llin/en/. Accessed 2016 Jul 12.
Government of Kenya, Ministry of Health (GoK, MoH). National Malaria Strategy: 2001–2010. Nairobi, Kenya; 2001. http://www.eac.int/health/index.php?option=com_docman&task=doc_download&gid=40&Itemid=169. Accessed 2016 May 22.
Noor AM, Amin AA, Akhwale WS, Snow RW. Increasing coverage and decreasing inequity in insecticide-treated bed net use among rural Kenyan children. PLoS Med. 2007;4:1341–8.
Government of Kenya, Ministry of Health (GoK, MoH). The epidemiology and control profile of malaria in Kenya: reviewing the evidence to guide the future vector control. National Malaria Control Programme, Ministry of Health. Technical support provided by the LINK Project (London School of Hygiene and Tropical Medicine and the Information for Malaria (INFORM) Project, KEMRI-Wellcome Trust Research Programme). Nairobi, Kenya; 2016. http://virtual.lshtm.ac.uk/wpcontent/uploads/2016/11/Kenya-Epidemiological-Profile.pdf. Accessed 2016 Dec 13.
Government of Kenya, Ministry of Health (GoK, MoH). National guidelines for the diagnosis, treatment and prevention of malaria in Kenya; Nairobi, Kenya; 2010. http://www.thehealthcompass.org/sites/default/files/project_examples/Kenya_Malaria_Tx_Guideline_2010.pdf. Accessed 2016 Sep 8.
Roll Back Malaria (RBM). Country-to-country guide for implementers of LLIN keep up: A guide for continuous delivery of LLINs via ANC, EPI, and other routine Health Services. Geneva, Switzerland; 2011. http://www.rbm.who.int/mechanisms/vcwgWorkStream3.htmlAccessed 2016 Feb 27.
Theiss-Nyland K, Lynch M, Lines J. Assessing the availability of LLINs for continuous distribution through routine antenatal care and the Expanded Programme on Immunizations in sub-Saharan Africa. Malar J. 2016;15:255.
World Health Organization (WHO). Integrated Management of Pregnancy and Childbirth. WHO Recommended Interventions for Improving Maternal and Newborn Health. Geneva, Switzerland. 2009. http://apps.who.int/iris/bitstream/10665/69509/1/WHO_MPS_07.05_eng.pdf. 10. Accessed 8 Sep 2016.
Simkhada B, Van Teijlingen ER, Porter M, Simkhada P. Factors affecting the utilization of antenatal care in developing countries: systematic review of the literature. J Adv Nurs. 2008;61:244–60.
Magadi MA, Madise NJ, Rodrigues RN. Variations in antenatal care between women of different communities in Kenya. Soc Sci Med. 2000;51:551–6.
Singh M, Brown G, Rogerson SJ. Ownership and use of insecticide-treated nets during pregnancy in sub-Saharan Africa: a review. Malar J. 2013;12:268.
Alegana VA, Wright JA, Petrina U, Noor AM, Snow RW, Atkinson PM. Spatial modelling of healthcare utilisation for treatment of fever in Namibia. Int J Health Geogr. 2012;11:6.
Ouma PO, Agutu NO, Snow RW, Noor AM. Univariate and multivariate spatial models of health facility utilisation for childhood fevers in an area on the coast of Kenya (in press).
Theiss-Nyland K, Ejersa W, Karema C, Koné D, Koenker H, Cyaka Y, et al. Operational challenges to continuous LLIN distribution: a qualitative rapid assessment in four countries. Malar J. 2016;15:131.
O’Meara W, Smith N, Ekal E, Cole D, Ndege S. Spatial distribution of bednet coverage under routine distribution through the public health sector in a rural district in Kenya. PLoS ONE. 2011;6:e25949.
Government of Kenya, Ministry of Health (GoK, MoH). Kenya Malaria Monitoring and Evaluation Plan 2009–2017. July 2009. Nairobi, Kenya. 2009. https://www.medbox.org/kenya-malaria-monitoring-and-evaluation-plan-2009-2018/download.pdf. Accessed 5 May 2016.
Noor AM, Gething PW, Alegana VA, Patil AP, Hay SI, Muchiri E, et al. The risks of malaria infection in Kenya in 2009. BMC Infect Dis. 2009;9:180.
Noor AM, Alegana VA, Gething PW, Snow RW. A spatial national health facility database for public health sector planning in Kenya in 2008. Int J Health Geogr. 2009;8:13.
Government of Kenya, Ministry of Health (GoK, MoH). Kenya Master Health Facility Listing. 2016. http://kmhfl.health.go.ke/#/home. Accessed 2016 June 20.
Government of Kenya, Ministry of Health (GoK, MoH). DHIS 2: Kenya Health Information System. 2016. https://hiskenya.org/. Accessed 2016 Feb 20.
Google Earth. https://www.google.com/earth/. Accessed 2016 Feb 23.
Worldpop. http://www.worldpop.org.uk/data/. Accessed 2016 May 03.
Linard C, Gilbert M, Snow RW, Noor AM, Tatem AJ. Population distribution, settlement patterns and accessibility across Africa in 2010. PLoS ONE. 2012;7:e31743.
Mennis J. Dasymetric mapping for estimating population in small areas. Geogr Compass. 2009;3:727–45.
United Nations Population Division (UNPD). World Population Prospects: The 2015 Revision, Data Booklet. ST/ESA/SER.A/377. New York, USA: United Nations; 2015. https://esa.un.org/unpd/wpp/Publications/Files/WPP2015_DataBooklet.pdf.
Tatem AJ, Campbell J, Guerra-Arias M, de Bernis L, Moran A, Matthews Z. Mapping for maternal and newborn health: the distributions of women of childbearing age, pregnancies and births. Int J Health Geogr. 2014;13:2.
Guttmacher Institute. https://www.guttmacher.org/. Accessed 2016 May 13.
Singh S, Darroch JE, Ashford LS. Adding It Up. The need for and cost of maternal and newborn care -estimates for 2012. New York, USA; 2013. http://www.guttmacher.org/pubs/AIU-MNH-2012-estimates.pdf.
Kenya National Bureau of Statistics (KNBS), ICF International. Kenya Demographic and Health Survey 2014. Nairobi, Kenya; 2015. https://dhsprogram.com/pubs/pdf/FR308/FR308.pdf. Accessed 2016 Jan 20.
Burgert CR, Colston J, Roy T, Zachary B. Geographic displacement procedure and georeferenced data release policy for the Demographic and Health Surveys. DHS Spatial Analysis. Reports No. 7. 2013. http://dhsprogram.com/pubs/pdf/SAR7/SAR7.pdf. Accessed 2016 Jan 20.
Ruktanonchai CW, Ruktanonchai NW, Nove A, Lopes S, Pezzulo C, Bosco C, et al. Equality in maternal and newborn health: modelling geographic disparities in utilisation of care in five East African Countries. PLoS ONE. 2016;11:e0162006.
Masters SH, Burstein R, Amofah G, Abaogye P, Kumar S, Hanlon M. Travel time to maternity care and its effect on utilization in rural Ghana: a multilevel analysis. Soc Sci Med. 2013;93:147–54.
Warren JL, Perez-Heydrich C, Burgert CR, Emch ME. Influence of demographic and health survey point displacements on raster-based analyses. Spat Demogr. 2016;4:135–53.
OpenStreetMap. http://extract.bbbike.org/. Accessed 2015 May 9.
Google Map Maker Data Download. https://services.google.com/fb/forms/mapmakerdatadownload/. Accessed 2016 Apr 20.
Government of Kenya (GoK). The Kenya Roads Act,2007. Kenya; 2007. http://www.kenyalaw.org/lex/rest//db/kenyalex/Kenya/Legislation/English/AmendmentActs/No.2of2007.pdf. Accessed 2016 Feb 22.
Government of Kenya. Kenya Roads Bill 2015. 2015 Nairobi, Kenya. http://kenyalaw.org/kl/fileadmin/pdfdownloads/bills/2015/KenyaRoadsBill2015.pdf. Accessed 2016 Feb 22.
Globeland30. http://www.globallandcover.com/GLC30Download/index.aspx. Accessed 2016 Apr 20.
USGS Earth Explorer. https://earthexplorer.usgs.gov/. Accessed 2016 Apr 20.
Noor AM, Amin AA, Gething PW, Atkinson PM, Hay SI, Snow RW. Modelling distances travelled to government health services in Kenya. Trop Med Int Health. 2006;11:188–96.
Dixit A, Lee M-C, Goettsch B, Afrane Y, Githeko AK, Yan G. Discovering the cost of care: consumer, provider, and retailer surveys shed light on the determinants of malaria health-seeking behaviours. Malar J. 2016;15:179.
Macharia Peter, Ouma Paul, Gogo Ezekiel, Snow Robert, Noor Abdisalan. Spatial accessibility to basic public health services in South Sudan. Geospat Health. 2017;12:510.
Tobler W. Three presentations on geographical analysis and modeling: non-isotropic geographic modeling; speculations on the geometry of geography; global spatial analysis. Santa Barbara; 1993. http://www.geodyssey.com/papers/tobler93.html.
Zorn A. Bicycle Velocity and Power Calculator. 2008. http://www.kreuzotter.de/english/espeed.htm. Accessed 2016 Jul 4.
Austin C. Bike speed calculator. 2012. http://bikecalculator.com/. Accessed 2016 Jul 4.
Ray N, Ebener S. AccessMod 3.0: computing geographic coverage and accessibility to health care services using anisotropic movement of patients. Int J Health Geogr. 2008;7:63.
Stock R. Distance and the utilization of health facilities in rural Nigeria. Soc Sci Med. 1983;17:563–70.
Jones S, Wardlaw J, Crouch S, Carolan M. Modelling catchment areas for secondary care providers: a case study. Health Care Manag Sci. 2011;14:253–61.
Bates D, Machler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48.
Babyak MA. What you see may not be what you get: a brief, nontechnical introduction to overfitting in regression-type models. Psychosom Med. 2004;66:411–21.
Ribeiro PJ, Diggle PJ. geoR: A package for geostatistical data analysis using the R software. 2016. http://www.leg.ufpr.br/geoR/. Accessed 2016 Jul 20.
Van Eijk AM, Bles HM, Odhiambo F, Ayisi JG, Blokland IE, Rosen DH, et al. Use of antenatal services and delivery care among women in rural western Kenya: a community based survey. Reprod Health. 2006;3:2.
Brown CA, Sohani SB, Khan K, Lilford R, Mukhwana W. Antenatal care and perinatal outcomes in Kwale district, Kenya. BMC Pregnancy Childbirth. 2008;8:2.
Ochako R, Fotso JC, Ikamari L, Khasakhala A. Utilization of maternal health services among young women in Kenya: insights from the Kenya Demographic and Health Survey, 2003. Pregnancy Childbirth. 2011;11:4–9.
Government of Kenya, Ministry of Health (GoK MoH). Kenya Malaria Strategy 2009–2018 (Revised 2014). Nairobi, Kenya; 2015. http://www.nmcp.or.ke/index.php/resourcecentre/download-centre/category/4-programme-management?download=47:kenya-malaria-m-and-e-plan-2009-2018-revised-2014. Accessed 2016 Sep 15.
Githinji S, Oyando R, Malinga J, Ejersa W, Soti D, Rono J, et al. Completeness of malaria indicator data reporting via the District Health Information Software 2 in Kenya, 2011–2015. Malar J. 2017;16:344.
Schuurman N, Fiedler RS, Grzybowski SCW, Grund D. Defining rational hospital catchments for non-urban areas based on travel-time. Int J Health Geogr. 2006;5:43.
Noor AM, Zurovac D, Hay SI, Ochola Sa, Snow RW. Defining equity in physical access to clinical services using geographical information systems as part of malaria planning and monitoring in Kenya. Trop Med Int Heath. 2003;8:917–26.
Alford-Teaster J, Lange JM, Hubbard RA, Lee CI, Haas JS, Shi X, et al. Is the closest facility the one actually used? An assessment of travel time estimation based on mammography facilities. Int J Health Geogr. 2016;15:8.
Tatem AJ, Garcia AJ, Snow RW, Noor AM, Gaughan AE, Gilbert M, et al. Millennium development health metrics: where do Africa’s children and women of childbearing age live? Popul Health Metr. 2013;11:11.
PMM and AMN developed the study concept and design. PMM was responsible for data assembly and data analysis and wrote the first draft of the manuscript. AMN, RWS and PAO guided study design, model development and data interpretation and helped in drafting the manuscript. All authors reviewed the manuscript and contributed to the final submission. All authors read and approved the final manuscript.
The Population Services Kenya for providing LLIN distribution list.
The authors declare that they have no competing interests.
Availability of data and materials
The household survey datasets used in this study are publicly available to registered users on the Demographic and Health Surveys website. The LLIN database can be accessed from Population Services Kenya while all the other datasets are publically available through various online portals. The web links to the portals are provided in the manuscript.
Consent for publication
Not applicable. The manuscript does not contain any individual person’s data.
Ethics approval and consent to participate
This study used secondary data only. The procedures and questionnaires for DHS surveys have been reviewed and approved by the ICF International Institutional Review Board (IRB). The ICF International IRB ensures that the survey complies with the U.S. Department of Health and Human Services regulations for the protection of human subjects (45 CFR 46).
Funding for this project was provided through UK’s Department for International Development Support (# 203155) to a project on Strengthening the Use of Data for Malaria Decision Making in Africa first funded and piloted in 2013 and supported PMM. AMN is supported by the Wellcome Trust, UK as an Intermediate Fellow (#095127). RWS is supported by the Wellcome Trust, UK as a Principal Research Fellow (#103602). PMM, RWS and AMN acknowledge the support of the Wellcome Trust to the Kenya Major Overseas Programme (# 203077).
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.