The international limits and population at risk of Plasmodium vivax transmission in 2009

BACKGROUND

A research priority for Plasmodium vivax malaria is to improve our understanding of the spatial distribution of risk and its relationship with the burden of P. vivax disease in human populations. The aim of the research outlined in this article is to provide a contemporary evidence-based map of the global spatial extent of P. vivax malaria, together with estimates of the human population at risk (PAR) of any level of transmission in 2009.

METHODOLOGY

The most recent P. vivax case-reporting data that could be obtained for all malaria endemic countries were used to classify risk into three classes: malaria free, unstable (<0.1 case per 1,000 people per annum (p.a.)) and stable (> or =0.1 case per 1,000 p.a.) P. vivax malaria transmission. Risk areas were further constrained using temperature and aridity data based upon their relationship with parasite and vector bionomics. Medical intelligence was used to refine the spatial extent of risk in specific areas where transmission was reported to be absent (e.g., large urban areas and malaria-free islands). The PAR under each level of transmission was then derived by combining the categorical risk map with a high resolution population surface adjusted to 2009. The exclusion of large Duffy negative populations in Africa from the PAR totals was achieved using independent modelling of the gene frequency of this genetic trait. It was estimated that 2.85 billion people were exposed to some risk of P. vivax transmission in 2009, with 57.1% of them living in areas of unstable transmission. The vast majority (2.59 billion, 91.0%) were located in Central and South East (CSE) Asia, whilst the remainder were located in America (0.16 billion, 5.5%) and in the Africa+ region (0.10 billion, 3.5%). Despite evidence of ubiquitous risk of P. vivax infection in Africa, the very high prevalence of Duffy negativity throughout Central and West Africa reduced the PAR estimates substantially.

CONCLUSIONS

After more than a century of development and control, P. vivax remains more widely distributed than P. falciparum and is a potential cause of morbidity and mortality amongst the 2.85 billion people living at risk of infection, the majority of whom are in the tropical belt of CSE Asia. The probability of infection is reduced massively across Africa by the frequency of the Duffy negative trait, but transmission does occur on the continent and is a concern for Duffy positive locals and travellers. The final map provides the spatial limits on which the endemicity of P. vivax transmission can be mapped to support future cartographic-based burden estimations.

Estimating the number of paediatric fevers associated with malaria infection presenting to Africa's public health sector in 2007

BACKGROUND

As international efforts to increase the coverage of artemisinin-based combination therapy in public health sectors gather pace, concerns have been raised regarding their continued indiscriminate presumptive use for treating all childhood fevers. The availability of rapid-diagnostic tests to support practical and reliable parasitological diagnosis provides an opportunity to improve the rational treatment of febrile children across Africa. However, the cost effectiveness of diagnosis-based treatment polices will depend on the presumed numbers of fevers harbouring infection. Here we compute the number of fevers likely to present to public health facilities in Africa and the estimated number of these fevers likely to be infected with Plasmodium falciparum malaria parasites.

METHODS AND FINDINGS

We assembled first administrative-unit level data on paediatric fever prevalence, treatment-seeking rates, and child populations. These data were combined in a geographical information system model that also incorporated an adjustment procedure for urban versus rural areas to produce spatially distributed estimates of fever burden amongst African children and the subset likely to present to public sector clinics. A second data assembly was used to estimate plausible ranges for the proportion of paediatric fevers seen at clinics positive for P. falciparum in different endemicity settings. We estimated that, of the 656 million fevers in African 0-4 y olds in 2007, 182 million (28%) were likely to have sought treatment in a public sector clinic of which 78 million (43%) were likely to have been infected with P. falciparum (range 60-103 million).

CONCLUSIONS

Spatial estimates of childhood fevers and care-seeking rates can be combined with a relational risk model of infection prevalence in the community to estimate the degree of parasitemia in those fevers reaching public health facilities. This quantification provides an important baseline comparison of malarial and nonmalarial fevers in different endemicity settings that can contribute to ongoing scientific and policy debates about optimum clinical and financial strategies for the introduction of new diagnostics. These models are made publicly available with the publication of this paper.

Climate change and the global malaria recession

The current and potential future impact of climate change on malaria is of major public health interest1,2. The proposed effects of rising global temperatures on the future spread and intensification of the disease3-5, and on existing malaria morbidity and mortality rates3, substantively influence global health policy6,7. The contemporary spatial limits of Plasmodium falciparum malaria and its endemicity within this range8, when compared with comparable historical maps, offer unique insights into the changing global epidemiology of malaria over the last century. It has long been known that the range of malaria has contracted through a century of economic development and disease control9. Here, for the first time, we quantify this contraction and the global decreases in malaria endemicity since c. 1900. We compare the magnitude of these changes to the size of effects on malaria endemicity hypothesised under future climate scenarios and associated with widely used public health interventions. Our findings have two key and often ignored implications with respect to climate change and malaria. First, widespread claims that rising mean temperatures have already led to increases in worldwide malaria morbidity and mortality are largely at odds with observed decreasing global trends in both its endemicity and geographic extent. Second, the proposed future effects of rising temperatures on endemicity are at least one order of magnitude smaller than changes observed since c. 1900 and up to two orders of magnitude smaller than those that can be achieved by the effective scale-up of key control measures. Predictions of an intensification of malaria in a warmer world, based on extrapolated empirical relationships or biological mechanisms, must be set against a context of a century of warming that has seen dramatic global declines in the disease and a substantial weakening of the global correlation between malaria endemicity and climate.

Optimal survey designs for targeting chemotherapy against soil-transmitted helminths: effect of spatial heterogeneity and cost-efficiency of sampling

Implementation of helminth control programs requires information on the distribution and prevalence of infection to target mass treatment to areas of greatest need. In the absence of data, the question of how many schools/communities should be surveyed depends on the spatial heterogeneity of infection and the cost efficiency of surveys. We used geostatistical techniques to quantify the spatial heterogeneity of soil-transmitted helminths in multiple settings in eastern Africa, and using the example of Kenya, conducted conditional simulation to explore the implications of alternative sampling strategies in identifying districts requiring mass treatment. Cost analysis is included in the simulations using data from actual field surveys and control programs. The analysis suggests that sampling four or five schools in each district provides a cost-efficient strategy in identifying districts requiring mass treatment, and that efficiency of sampling was relatively insensitive to the number of children sampled per school.

Estimating the global clinical burden of Plasmodium falciparum malaria in 2007

BACKGROUND

The epidemiology of malaria makes surveillance-based methods of estimating its disease burden problematic. Cartographic approaches have provided alternative malaria burden estimates, but there remains widespread misunderstanding about their derivation and fidelity. The aims of this study are to present a new cartographic technique and its application for deriving global clinical burden estimates of Plasmodium falciparum malaria for 2007, and to compare these estimates and their likely precision with those derived under existing surveillance-based approaches.

METHODS AND FINDINGS

In seven of the 87 countries endemic for P. falciparum malaria, the health reporting infrastructure was deemed sufficiently rigorous for case reports to be used verbatim. In the remaining countries, the mapped extent of unstable and stable P. falciparum malaria transmission was first determined. Estimates of the plausible incidence range of clinical cases were then calculated within the spatial limits of unstable transmission. A modelled relationship between clinical incidence and prevalence was used, together with new maps of P. falciparum malaria endemicity, to estimate incidence in areas of stable transmission, and geostatistical joint simulation was used to quantify uncertainty in these estimates at national, regional, and global scales. Combining these estimates for all areas of transmission risk resulted in 451 million (95% credible interval 349-552 million) clinical cases of P. falciparum malaria in 2007. Almost all of this burden of morbidity occurred in areas of stable transmission. More than half of all estimated P. falciparum clinical cases and associated uncertainty occurred in India, Nigeria, the Democratic Republic of the Congo (DRC), and Myanmar (Burma), where 1.405 billion people are at risk. Recent surveillance-based methods of burden estimation were then reviewed and discrepancies in national estimates explored. When these cartographically derived national estimates were ranked according to their relative uncertainty and replaced by surveillance-based estimates in the least certain half, 98% of the global clinical burden continued to be estimated by cartographic techniques.

CONCLUSIONS AND SIGNIFICANCE

Cartographic approaches to burden estimation provide a globally consistent measure of malaria morbidity of known fidelity, and they represent the only plausible method in those malaria-endemic countries with nonfunctional national surveillance. Unacceptable uncertainty in the clinical burden of malaria in only four countries confounds our ability to evaluate needs and monitor progress toward international targets for malaria control at the global scale. National prevalence surveys in each nation would reduce this uncertainty profoundly. Opportunities for further reducing uncertainty in clinical burden estimates by hybridizing alternative burden estimation procedures are also evaluated.

Quantifying aggregated uncertainty in Plasmodium falciparum malaria prevalence and populations at risk via efficient space-time geostatistical joint simulation

Risk maps estimating the spatial distribution of infectious diseases are required to guide public health policy from local to global scales. The advent of model-based geostatistics (MBG) has allowed these maps to be generated in a formal statistical framework, providing robust metrics of map uncertainty that enhances their utility for decision-makers. In many settings, decision-makers require spatially aggregated measures over large regions such as the mean prevalence within a country or administrative region, or national populations living under different levels of risk. Existing MBG mapping approaches provide suitable metrics of local uncertainty—the fidelity of predictions at each mapped pixel—but have not been adapted for measuring uncertainty over large areas, due largely to a series of fundamental computational constraints. Here the authors present a new efficient approximating algorithm that can generate for the first time the necessary joint simulation of prevalence values across the very large prediction spaces needed for global scale mapping. This new approach is implemented in conjunction with an established model for P. falciparum allowing robust estimates of mean prevalence at any specified level of spatial aggregation. The model is used to provide estimates of national populations at risk under three policy-relevant prevalence thresholds, along with accompanying model-based measures of uncertainty. By overcoming previously unchallenged computational barriers, this study illustrates how MBG approaches, already at the forefront of infectious disease mapping, can be extended to provide large-scale aggregate measures appropriate for decision-makers.

Reliable disease maps can support rational decision making. These maps are often made by interpolation: taking disease data from field studies and predicting values for the gaps between the data to make a complete map. Such maps always contain uncertainty, however, and measuring this uncertainty is vital so that the reliability of risk maps can be determined. A modern approach called model-based geostatistics (MBG) has led to increasingly sophisticated ways of mapping disease and measuring spatial uncertainty. Many health management decisions are made for administrative areas (e.g., districts, provinces, countries) and disease maps can support these decisions by averaging their values over the regions of interest. Carrying out this aggregation in conjunction with MBG techniques has not previously been possible for very large maps, however, due largely to the computational constraints involved. This study has addressed this problem by developing a new algorithm and allows aggregation of a global MBG disease map over very large areas. It is used to estimate Plasmodium falciparum malaria prevalence and corresponding populations at risk worldwide, aggregated across regions of different sizes. These estimates are a cornerstone for disease burden estimation and are provided in full to facilitate that process.

Spatial predictions of Rhodesian Human African Trypanosomiasis (sleeping sickness) prevalence in Kaberamaido and Dokolo, two newly affected districts of Uganda

The continued northwards spread of Rhodesian sleeping sickness or Human African Trypanosomiasis (HAT) within Uganda is raising concerns of overlap with the Gambian form of the disease. Disease convergence would result in compromised diagnosis and treatment for HAT. Spatial determinants for HAT are poorly understood across small areas. This study examines the relationships between Rhodesian HAT and several environmental, climatic and social factors in two newly affected districts, Kaberamaido and Dokolo. A one-step logistic regression analysis of HAT prevalence and a two-step logistic regression method permitted separate analysis of both HAT occurrence and HAT prevalence. Both the occurrence and prevalence of HAT were negatively correlated with distance to the closest livestock market in all models. The significance of distance to the closest livestock market strongly indicates that HAT may have been introduced to this previously unaffected area via the movement of infected, untreated livestock from endemic areas. This illustrates the importance of the animal reservoir in disease transmission, and highlights the need for trypanosomiasis control in livestock and the stringent implementation of regulations requiring the treatment of cattle prior to sale at livestock markets to prevent any further spread of Rhodesian HAT within Uganda.

The risks of malaria infection in Kenya in 2009

Background

To design an effective strategy for the control of malaria requires a map of infection and disease risks to select appropriate suites of interventions. Advances in model based geo-statistics and malaria parasite prevalence data assemblies provide unique opportunities to redefine national Plasmodium falciparum risk distributions. Here we present a new map of malaria risk for Kenya in 2009.

Methods

Plasmodium falciparum parasite rate data were assembled from cross-sectional community based surveys undertaken from 1975 to 2009. Details recorded for each survey included the month and year of the survey, sample size, positivity and the age ranges of sampled population. Data were corrected to a standard age-range of two to less than 10 years (PfPR_2-10) and each survey location was geo-positioned using national and on-line digital settlement maps. Ecological and climate covariates were matched to each PfPR_2-10 survey location and examined separately and in combination for relationships to PfPR_2-10. Significant covariates were then included in a Bayesian geostatistical spatial-temporal framework to predict continuous and categorical maps of mean PfPR_2-10 at a 1 × 1 km resolution across Kenya for the year 2009. Model hold-out data were used to test the predictive accuracy of the mapped surfaces and distributions of the posterior uncertainty were mapped.

Results

A total of 2,682 estimates of PfPR_2-10 from surveys undertaken at 2,095 sites between 1975 and 2009 were selected for inclusion in the geo-statistical modeling. The covariates selected for prediction were urbanization; maximum temperature; precipitation; enhanced vegetation index; and distance to main water bodies. The final Bayesian geo-statistical model had a high predictive accuracy with mean error of -0.15% PfPR_2-10; mean absolute error of 0.38% PfPR_2-10; and linear correlation between observed and predicted PfPR_2-10 of 0.81. The majority of Kenya's 2009 population (35.2 million, 86.3%) reside in areas where predicted PfPR_2-10 is less than 5%; conversely in 2009 only 4.3 million people (10.6%) lived in areas where PfPR_2-10 was predicted to be ≥40% and were largely located around the shores of Lake Victoria.

Conclusion

Model based geo-statistical methods can be used to interpolate malaria risks in Kenya with precision and our model shows that the majority of Kenyans live in areas of very low P. falciparum risk. As malaria interventions go to scale effectively tracking epidemiological changes of risk demands a rigorous effort to document infection prevalence in time and space to remodel risks and redefine intervention priorities over the next 10-15 years.

Defining the relationship between Plasmodium falciparum parasite rate and clinical disease: statistical models for disease burden estimation

BACKGROUND

Clinical malaria has proven an elusive burden to enumerate. Many cases go undetected by routine disease recording systems. Epidemiologists have, therefore, frequently defaulted to actively measuring malaria in population cohorts through time. Measuring the clinical incidence of malaria longitudinally is labour-intensive and impossible to undertake universally. There is a need, therefore, to define a relationship between clinical incidence and the easier and more commonly measured index of infection prevalence: the "parasite rate". This relationship can help provide an informed basis to define malaria burdens in areas where health statistics are inadequate.

METHODS

Formal literature searches were conducted for Plasmodium falciparum malaria incidence surveys undertaken prospectively through active case detection at least every 14 days. The data were abstracted, standardized and geo-referenced. Incidence surveys were time-space matched with modelled estimates of infection prevalence derived from a larger database of parasite prevalence surveys and modelling procedures developed for a global malaria endemicity map. Several potential relationships between clinical incidence and infection prevalence were then specified in a non-parametric Gaussian process model with minimal, biologically informed, prior constraints. Bayesian inference was then used to choose between the candidate models.

RESULTS

The suggested relationships with credible intervals are shown for the Africa and a combined America and Central and South East Asia regions. In both regions clinical incidence increased slowly and smoothly as a function of infection prevalence. In Africa, when infection prevalence exceeded 40%, clinical incidence reached a plateau of 500 cases per thousand of the population per annum. In the combined America and Central and South East Asia regions, this plateau was reached at 250 cases per thousand of the population per annum. A temporal volatility model was also incorporated to facilitate a closer description of the variance in the observed data.

CONCLUSION

It was possible to model a relationship between clinical incidence and P. falciparum infection prevalence but the best-fit models were very noisy reflecting the large variance within the observed opportunistic data sample. This continuous quantification allows for estimates of the clinical burden of P. falciparum of known confidence from wherever an estimate of P. falciparum prevalence is available.

A world malaria map: Plasmodium falciparum endemicity in 2007

BACKGROUND

Efficient allocation of resources to intervene against malaria requires a detailed understanding of the contemporary spatial distribution of malaria risk. It is exactly 40 y since the last global map of malaria endemicity was published. This paper describes the generation of a new world map of Plasmodium falciparum malaria endemicity for the year 2007.

METHODS AND FINDINGS

A total of 8,938 P. falciparum parasite rate (PfPR) surveys were identified using a variety of exhaustive search strategies. Of these, 7,953 passed strict data fidelity tests for inclusion into a global database of PfPR data, age-standardized to 2-10 y for endemicity mapping. A model-based geostatistical procedure was used to create a continuous surface of malaria endemicity within previously defined stable spatial limits of P. falciparum transmission. These procedures were implemented within a Bayesian statistical framework so that the uncertainty of these predictions could be evaluated robustly. The uncertainty was expressed as the probability of predicting correctly one of three endemicity classes; previously stratified to be an informative guide for malaria control. Population at risk estimates, adjusted for the transmission modifying effects of urbanization in Africa, were then derived with reference to human population surfaces in 2007. Of the 1.38 billion people at risk of stable P. falciparum malaria, 0.69 billion were found in Central and South East Asia (CSE Asia), 0.66 billion in Africa, Yemen, and Saudi Arabia (Africa+), and 0.04 billion in the Americas. All those exposed to stable risk in the Americas were in the lowest endemicity class (PfPR2-10 < or = 5%). The vast majority (88%) of those living under stable risk in CSE Asia were also in this low endemicity class; a small remainder (11%) were in the intermediate endemicity class (PfPR2-10 > 5 to < 40%); and the remaining fraction (1%) in high endemicity (PfPR2-10 > or = 40%) areas. High endemicity was widespread in the Africa+ region, where 0.35 billion people are at this level of risk. Most of the rest live at intermediate risk (0.20 billion), with a smaller number (0.11 billion) at low stable risk.

CONCLUSIONS

High levels of P. falciparum malaria endemicity are common in Africa. Uniformly low endemic levels are found in the Americas. Low endemicity is also widespread in CSE Asia, but pockets of intermediate and very rarely high transmission remain. There are therefore significant opportunities for malaria control in Africa and for malaria elimination elsewhere. This 2007 global P. falciparum malaria endemicity map is the first of a series with which it will be possible to monitor and evaluate the progress of this intervention process.

A spatial national health facility database for public health sector planning in Kenya in 2008

BACKGROUND

Efforts to tackle the enormous burden of ill-health in low-income countries are hampered by weak health information infrastructures that do not support appropriate planning and resource allocation. For health information systems to function well, a reliable inventory of health service providers is critical. The spatial referencing of service providers to allow their representation in a geographic information system is vital if the full planning potential of such data is to be realized.

METHODS

A disparate series of contemporary lists of health service providers were used to update a public health facility database of Kenya last compiled in 2003. These new lists were derived primarily through the national distribution of antimalarial and antiretroviral commodities since 2006. A combination of methods, including global positioning systems, was used to map service providers. These spatially-referenced data were combined with high-resolution population maps to analyze disparity in geographic access to public health care.

FINDINGS

The updated 2008 database contained 5,334 public health facilities (67% ministry of health; 28% mission and nongovernmental organizations; 2% local authorities; and 3% employers and other ministries). This represented an overall increase of 1,862 facilities compared to 2003. Most of the additional facilities belonged to the ministry of health (79%) and the majority were dispensaries (91%). 93% of the health facilities were spatially referenced, 38% using global positioning systems compared to 21% in 2003. 89% of the population was within 5 km Euclidean distance to a public health facility in 2008 compared to 71% in 2003. Over 80% of the population outside 5 km of public health service providers was in the sparsely settled pastoralist areas of the country.

CONCLUSION

We have shown that, with concerted effort, a relatively complete inventory of mapped health services is possible with enormous potential for improving planning. Expansion in public health care in Kenya has resulted in significant increases in geographic access although several areas of the country need further improvements. This information is key to future planning and with this paper we have released the digital spatial database in the public domain to assist the Kenyan Government and its partners in the health sector.

Using remotely sensed night-time light as a proxy for poverty in Africa

BACKGROUND

Population health is linked closely to poverty. To assess the effectiveness of health interventions it is critical to monitor the spatial and temporal changes in the health indicators of populations and outcomes across varying levels of poverty. Existing measures of poverty based on income, consumption or assets are difficult to compare across geographic settings and are expensive to construct. Remotely sensed data on artificial night time lights (NTL) have been shown to correlate with gross domestic product in developed countries.

METHODS

Using national household survey data, principal component analysis was used to compute asset-based poverty indices from aggregated household asset variables at the Administrative 1 level (n = 338) in 37 countries in Africa. Using geographical information systems, mean brightness of and distance to NTL pixels and proportion of area covered by NTL were computed for each Administrative1 polygon. Correlations and agreement of asset-based indices and the three NTL metrics were then examined in both continuous and ordinal forms.

RESULTS

At the Administrative 1 level all the NTL metrics distinguished between the most poor and least poor quintiles with greater precision compared to intermediate quintiles. The mean brightness of NTL, however, had the highest correlation coefficient with the asset-based wealth index in continuous (Pearson correlation = 0.64, p < 0.01) and ordinal (Spearman correlation = 0.79, p < 0.01; Kappa = 0.64) forms.

CONCLUSION

Metrics of the brightness of NTL data offer a robust and inexpensive alternative to asset-based poverty indices derived from survey data at the Administrative 1 level in Africa. These could be used to explore economic inequity in health outcomes and access to health interventions at sub-national levels where household assets data are not available at the required resolution.

Spatial prediction of Plasmodium falciparum prevalence in Somalia

BACKGROUND

Maps of malaria distribution are vital for optimal allocation of resources for anti-malarial activities. There is a lack of reliable contemporary malaria maps in endemic countries in sub-Saharan Africa. This problem is particularly acute in low malaria transmission countries such as those located in the horn of Africa.

METHODS

Data from a national malaria cluster sample survey in 2005 and routine cluster surveys in 2007 were assembled for Somalia. Rapid diagnostic tests were used to examine the presence of Plasmodium falciparum parasites in finger-prick blood samples obtained from individuals across all age-groups. Bayesian geostatistical models, with environmental and survey covariates, were used to predict continuous maps of malaria prevalence across Somalia and to define the uncertainty associated with the predictions.

RESULTS

For analyses the country was divided into north and south. In the north, the month of survey, distance to water, precipitation and temperature had no significant association with P. falciparum prevalence when spatial correlation was taken into account. In contrast, all the covariates, except distance to water, were significantly associated with parasite prevalence in the south. The inclusion of covariates improved model fit for the south but not for the north. Model precision was highest in the south. The majority of the country had a predicted prevalence of < 5%; areas with > or = 5% prevalence were predominantly in the south.

CONCLUSION

The maps showed that malaria transmission in Somalia varied from hypo- to meso-endemic. However, even after including the selected covariates in the model, there still remained a considerable amount of unexplained spatial variation in parasite prevalence, indicating effects of other factors not captured in the study. Nonetheless the maps presented here provide the best contemporary information on malaria prevalence in Somalia.

Developing geostatistical space-time models to predict outpatient treatment burdens from incomplete national data

Basic health system data such as the number of patients utilising different health facilities and the types of illness for which they are being treated are critical for managing service provision. These data requirements are generally addressed with some form of national Health Management Information System (HMIS) which coordinates the routine collection and compilation of data from national health facilities. HMIS in most developing countries are characterised by widespread under-reporting. Here we present a method to adjust incomplete data to allow prediction of national outpatient treatment burdens. We demonstrate this method with the example of outpatient treatments for malaria within the Kenyan HMIS. Three alternative modelling frameworks were developed and tested in which space-time geostatistical prediction algorithms were used to predict the monthly tally of treatments for presumed malaria cases (MC) at facilities where such records were missing. Models were compared by a cross-validation exercise and the model found to most accurately predict MC incorporated available data on the total number of patients visiting each facility each month. A space-time stochastic simulation framework to accompany this model was developed and tested in order to provide estimates of both local and regional prediction uncertainty. The level of accuracy provided by the predictive model, and the accompanying estimates of uncertainty around the predictions, demonstrate how this tool can mitigate the uncertainties caused by missing data, substantially enhancing the utility of existing HMIS data to health-service decision-makers.

Information for decision making from imperfect national data: tracking major changes in health care use in Kenya using geostatistics

BACKGROUND

Most Ministries of Health across Africa invest substantial resources in some form of health management information system (HMIS) to coordinate the routine acquisition and compilation of monthly treatment and attendance records from health facilities nationwide. Despite the expense of these systems, poor data coverage means they are rarely, if ever, used to generate reliable evidence for decision makers. One critical weakness across Africa is the current lack of capacity to effectively monitor patterns of service use through time so that the impacts of changes in policy or service delivery can be evaluated. Here, we present a new approach that, for the first time, allows national changes in health service use during a time of major health policy change to be tracked reliably using imperfect data from a national HMIS.

METHODS

Monthly attendance records were obtained from the Kenyan HMIS for 1 271 government-run and 402 faith-based outpatient facilities nationwide between 1996 and 2004. A space-time geostatistical model was used to compensate for the large proportion of missing records caused by non-reporting health facilities, allowing robust estimation of monthly and annual use of services by outpatients during this period.

RESULTS

We were able to reconstruct robust time series of mean levels of outpatient utilisation of health facilities at the national level and for all six major provinces in Kenya. These plots revealed reliably for the first time a period of steady nationwide decline in the use of health facilities in Kenya between 1996 and 2002, followed by a dramatic increase from 2003. This pattern was consistent across different causes of attendance and was observed independently in each province.

CONCLUSION

The methodological approach presented can compensate for missing records in health information systems to provide robust estimates of national patterns of outpatient service use. This represents the first such use of HMIS data and contributes to the resurrection of these hugely expensive but underused systems as national monitoring tools. Applying this approach to Kenya has yielded output with immediate potential to enhance the capacity of decision makers in monitoring nationwide patterns of service use and assessing the impact of changes in health policy and service delivery.

Improving imperfect data from health management information systems in Africa using space-time geostatistics

BACKGROUND

Reliable and timely information on disease-specific treatment burdens within a health system is critical for the planning and monitoring of service provision. Health management information systems (HMIS) exist to address this need at national scales across Africa but are failing to deliver adequate data because of widespread underreporting by health facilities. Faced with this inadequacy, vital public health decisions often rely on crudely adjusted regional and national estimates of treatment burdens.

METHODS AND FINDINGS

This study has taken the example of presumed malaria in outpatients within the largely incomplete Kenyan HMIS database and has defined a geostatistical modelling framework that can predict values for all data that are missing through space and time. The resulting complete set can then be used to define treatment burdens for presumed malaria at any level of spatial and temporal aggregation. Validation of the model has shown that these burdens are quantified to an acceptable level of accuracy at the district, provincial, and national scale.

CONCLUSIONS

The modelling framework presented here provides, to our knowledge for the first time, reliable information from imperfect HMIS data to support evidence-based decision-making at national and sub-national levels.

Modelling distances travelled to government health services in Kenya

OBJECTIVE

To systematically evaluate descriptive measures of spatial access to medical treatment, as part of the millennium development goals to reduce the burden of HIV/AIDS, tuberculosis and malaria.

METHODS

We obtained high-resolution spatial and epidemiological data on health services, population, transport network, topography, land cover and paediatric fever treatment in four Kenyan districts to develop access and use models for government health services in Kenya. Community survey data were used to model use of government health services by febrile children. A model based on the transport network was then implemented and adjusted for actual use patterns. We compared the predictive accuracy of this refined model to that of Euclidean distance metrics. RESULTS Higher-order facilities were more attractive to patients (54%, 58% and 60% in three scenarios) than lower-order ones. The transport network model, adjusted for competition between facilities, was most accurate and selected as the best-fit model. It estimated that 63% of the population of the study districts were within the 1 h national access benchmark, against 82% estimated by the Euclidean model.

CONCLUSIONS

Extrapolating the results from the best-fit model in study districts to the national level shows that approximately six million people are currently incorrectly estimated to have access to government health services within 1 h. Simple Euclidean distance assumptions, which underpin needs assessments and against which millennium development goals are evaluated, thus require reconsideration.

Empirical modelling of government health service use by children with fevers in Kenya

An understanding of spatial patterns of health facility use allows a more informed approach to the modelling of catchment populations. In the absence of patient use data, an intuitive and commonly used approach to the delineation of facility catchment areas is Thiessen polygons. This study presents a series of methods by which the validity of these assumptions can be tested directly and hence the suitability of a Thiessen polygon catchment model explicitly assessed. These methods are applied to paediatric out-patient origin data from a sample of 81 government health facilities in four districts of Kenya. A geographical information system was used to predict the location of the catchment boundary along a transect between each pair of neighbouring facilities based on patient choice patterns. The mean location of boundaries between facilities of different type was found to be significantly displaced from the Thiessen boundary towards the lower-order facility. The affect of distance on within-catchment utilization rate was assessed by using exclusion buffers to remove the effect of neighbouring facilities. Utilization rate was found to exhibit a slight but steady decrease with distance up to 6 km from a facility. The accuracy of the future modelling of unsampled facility catchments can be increased by the incorporation of these trends.