PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
 
Sci Rep. 2017; 7: 45849.
Published online 2017 April 12. doi:  10.1038/srep45849
PMCID: PMC5388846

Impact of metric and sample size on determining malaria hotspot boundaries

Abstract

The spatial heterogeneity of malaria suggests that interventions may be targeted for maximum impact. It is unclear to what extent different metrics lead to consistent delineation of hotspot boundaries. Using data from a large community-based malaria survey in the western Kenyan highlands, we assessed the agreement between a model-based geostatistical (MBG) approach to detect hotspots using Plasmodium falciparum parasite prevalence and serological evidence for exposure. Malaria transmission was widespread and highly heterogeneous with one third of the total population living in hotspots regardless of metric tested. Moderate agreement (Kappa = 0.424) was observed between hotspots defined based on parasite prevalence by polymerase chain reaction (PCR)- and the prevalence of antibodies to two P. falciparum antigens (MSP-1, AMA-1). While numerous biologically plausible hotspots were identified, their detection strongly relied on the proportion of the population sampled. When only 3% of the population was sampled, no PCR derived hotspots were reliably detected and at least 21% of the population was needed for reliable results. Similar results were observed for hotspots of seroprevalence. Hotspot boundaries are driven by the malaria diagnostic and sample size used to inform the model. These findings warn against the simplistic use of spatial analysis on available data to target malaria interventions in areas where hotspot boundaries are uncertain.

Malaria is an important cause of global morbidity and mortality with an estimated 3.4 billion people at risk1. The past decade has seen a large reduction in the malaria burden in some areas with an estimated 47% global reduction in mortality compared to 20002. As national policies shift from control towards elimination new approaches are needed to supplement existing tools3,4. Research and programmatic activities are increasingly acknowledging the heterogeneous nature of malaria transmission at the community level.

Identifying ‘hotspots’ and targeting malaria control interventions at these, could lead to a more sustainable reduction in malaria burden5,6. Hotspots are typically defined in both public health and ecology as areas where estimates exceed those from other areas and may fuel transmission to the surrounding areas5,7,8. Malaria transmission is difficult to measure directly and several metrics are typically used to estimate malaria burden as a proxy for transmission9. However, different malaria metrics measure different facets of the transmission cycle and may lead to different conclusions on the existence, size or location of hotspots. For example, in coastal Kenya hotspots based on clinical incidence were geographically distinct and showed different temporal dynamics compared to hotspots based on the prevalence of asymptomatic infections10,11.

The detection of malaria hotspots has become increasingly prominent in the malaria literature12,13,14,15,16,17. Model-based geostatistics (MBG) are increasingly being used to identify heterogeneity in malaria transmission and can predict areas of increased disease prevalence. MBG has been effectively applied in other disease systems that exhibit both large and small-scale variation in transmission18,19. In the context of malaria, MBG has mainly been applied at the national or provincial scales, although it has yet to be widely applied for local level spatial analysis13,20,21. It allows incorporating environmental drivers of disease transmission and information on the intensity of sampling to obtain smoothed values of disease indicators to determine spatial patterns in disease occurrence. Determining the hotspot boundaries is of great public health importance if hotspot-targeted interventions are considered. Uncertainties about hotspot boundaries would complicate and potentially reduce the impact of hotspot-targeted interventions by potentially missing populations that are particularly relevant for onward transmission or misallocating resources22.

Using data collected in a large cross-sectional malaria survey carried out in the western Kenyan highlands, the aims of this study were to compare the agreement between spatial analysis based on the prevalence of molecularly detected malaria infections and serological evidence for malaria exposure and illustrate the impact of sample size on the delineation of hotspots of malaria. The results generated are not intended to provide a gold standard for hotspot detection, but to illustrate the realities of translating theoretical concepts of disease heterogeneity into actionable public health strategies.

Methods

Data sources

Epidemiological

Epidemiological data were obtained from a community cross-sectional malaria survey conducted in July 2011 in a 100 km2 rural area in the western Kenyan highlands (0°28′S, 34°51′E)23. The site is characterized by low but heterogeneous malaria transmission, with Plasmodium falciparum being the predominant species24. Factors determining local malaria transmission patterns were recently described25.

All structures in the study area were digitized using high-resolution satellite imagery (Quickbird, DigitalGlobe Services Inc, USA) and were used as a proxy for the total population size and distribution22,23. Briefly, 17,503 individuals residing in 3,213 randomly selected households (i.e. clusters of structures forming a family unit), or approximately 30% of the total population, were surveyed, with each participant providing blood spot samples on filter paper. The unit of analysis was the household with the proportion of household residents that were positive for malaria the main outcome. A random selection of 79% of the collected samples were assayed by PCR to detect the presence of a current malaria infection, corresponding to an estimated 24% of the total population26,27; all samples were tested for anti-malarial antibody response to AMA1 and/or MSP119 measured by enzyme linked immunosorbent assay to provide a measure of malaria exposure28,29. Seropositivity to each antigen was assessed using a mixture model and consisted of those individuals with optical density values greater than the mean plus three standard deviations of the distribution of those assumed negative30. An individual was considered to be seropositive if they were positive to either or both of the antigens tested. Ethical approval for collecting the epidemiological data was granted by the London School of Hygiene & Tropical Medicine (LSHTM-5721) and the Kenya Medical Research Institute (SSC-1802). All methods were performed in accordance with good research practices and written informed consent was obtained from all participants.

Environmental

Elevation for each household was derived from the ASTER [v 2.0, NASA USA] global digital elevation model (DEM). The normalized-difference vegetation index (NDVI) was calculated for the study area using the Quickbird imagery. Mean, minimum, and maximum NDVI values from a single time-point were calculated for a 500 m circular buffer around each household. Multispectral image segmentation (MIS) of the Quickbird imagery was conducted with eCognition (v 4.0, Trimble Geospatial Imaging, Germany) software and the proportion of tree cover within the 500 m buffer was determined. Fishponds were identified using a refined MIS procedure capable of detecting smaller features and manually verified against the satellite imagery. The distance from each household to the nearest fishpond was calculated in ArcGIS (ESRI, USA).

Topographic wetness index (TWI) was calculated using the DEM as previously described31. The maximum and mean TWI values for the 500 m surrounding each household were calculated. Finally, the locations of all streams in the area were determined by first locating the likely location of streams using the topographic data and then manually digitizing the more precise stream path using the satellite imagery. The distances of each household to all stream orders were calculated32.

Determining hotspots of P. falciparum infection and exposure by MBG

MBG was used to model the spatial variation in malaria parasite or antibody prevalence18,19. Two models were generated using the PrevMap package33: malaria infection was assessed using PCR positivity and exposure to malaria was assessed using seropositivity estimates for each household (see Supplementary file 1)34. Surfaces of predicted prevalence for both outcomes were generated. To guide MBG, thresholds of risks were used that resulted in 20% of the population being included in the hotspot based on the theoretical 80–20 assumption where 20% of the population constitutes 80% of the exposure and transmission events35. We acknowledge that this threshold selection is likely to be site specific and the hotspot sizes will vary based on the threshold selected: a high threshold would result in only those areas with the highest transmission being identified as a hotspot and a more granular map whereas a less stringent threshold would mean that hotspots would be more ubiquitous.

Next, the probability that any given area exceeded the threshold that encompassed 20% of the population was determined. Areas with greater than 80% probability of exceeding the threshold were considered hotspots. In an ideal scenario, the model will produce a probability surface that is polarized into areas with 100 or 0% probability of exceeding a specified threshold. The 80% probability threshold was selected to capture those areas that are almost certainly in a hotspot as well as including those most likely to be in a hotspot, a decision made to favour a higher sensitivity rather than specificity. To gauge the sensitivity of the exceedance threshold in determining hotspots, we also identified areas that had greater than 50% probability of exceeding the threshold (i.e. any likelihood of being a hotspot). This process was repeated for both outcomes to generate separate surfaces for hotspots of current infection and exposure to malaria and there were no constraints placed on hotspot size or shape. The households consistently identified, or agreement, between hotspots of infection and exposure was assessed using Cohen’s Kappa coefficient. All analyses were conducted in R v.3.0.2 (R-Project, USA).

Sample size

The impact of sample size on the ability of the model to predict hotspots and therefore confidence in delineation of hotspot boundaries, was assessed by the change in metrics of model predictive performance: the integrated mean square error (IMSE) for the predicted surface and the discrimination index (DI) for the exceedance probabilities18. To estimate the level of performance that would have been achieved had the entire (100%) population been sampled, we imputed a complete population data set using the complete set of digitized households and the predicted malaria risk surfaces to estimate household level prevalence23. Next, we selected a random sub-set of the imputed data for each of the sampling fractions 10–90% and re-fitted the geostatistical models to each sub-set. The corresponding IMSE and DI values were calculated and plotted as functions of the sampling fraction.

To determine the impact of sample size on hotspot boundaries, the geostatistical model was then re-fitted to random subsets of the collected data, with sampling fractions between 10–90%. The resulting surfaces were imported into ArcGIS, hotspot boundaries determined, and individual households were assigned as hotspot or non-hotspot accordingly. The sensitivity and specificity of the structures correctly identified, using the complete sample as the reference, were calculated and compared using the area under the receiver operator curve (AUROC)36. A meaningful change was considered to be those with non-overlapping confidence intervals for the AUROC.

Results

MBG models

The results of the geostatistical model are consistent with previous studies and suggests a positive association between parasite prevalence and maximum and mean NDVI and a negative association with mean elevation, distance from fishponds and the proportion of tree cover (see Supplementary file 2 for model validation)22,25. The optimum model fit for seroprevalence also indicated a negative association with mean elevation, distance from fishponds and tree cover. In addition, maximum TWI, minimum NDVI and distance to 2nd and 3rd order streams had negative associations with seroprevalence, while mean TWI had a positive association (Table 1). The percentage of variability explained by the covariates was 3% for parasite prevalence and 18% for seroprevalence. The spatial stochastic process, which accounts for both local and global spatial trends, accounted for 53.4% and 49.5% of the unexplained variability according to PCR and seropositivity, respectively.

Table 1
Final adjusted mixed effects logistic regression models for both outcomes.

Comparing metrics: Molecular vs. Seroprevalence

Areas with a predicted PCR prevalence (Fig. 1a) greater than 28% and predicted seroprevalence (Fig. 1b) greater than 70% encompassed 20% of the total population. These thresholds were subsequently used to determine hotspots of infection and exposure, respectively, for the MBG approach. The probability of exceeding the defined thresholds was mapped for both current infection prevalence (Fig. 2a) and previous exposure, seropositivity (Fig. 2b). The agreement between households identified as part of hotspots derived using parasite and sero-prevalence with a probability >80% of exceeding the threshold was moderate (Kappa = 0.424). Using hotspot boundaries corresponding to areas with >50% probability of exceeding the threshold resulted in only modest improvement in agreement (Kappa = 0.478).

Figure 1
Predicted malaria prevalence using model based geostatistics.
Figure 2
Probability contour maps of exceedance surfaces for malaria prevalence.

Impact of sample size

As expected, IMSE decreased proportionally to the inverse of the sample size (Fig. 3). Based on this analysis, the models using the survey data to generate predictive surfaces were generated with an estimated baseline error of a 40% relative increase in IMSE.

Figure 3
Impact of sample size on geostatistical model efficiency.

The parasite prevalence model showed a change in the number of structures correctly identified when sample size was reduced to 70% of the sampled population, or 21% of the total population (Table 2). A second significant difference in the consistency of hotspot delineation was observed with 30% of the sampled population, or 9% of the total population (Supplementary movie 1). The geostatistical model for parasite prevalence was unable to reliably detect hotspots with less than 10% of the sampled, or 3% of the total population. The impact of sample size on the models for seroprevalence showed similar trends in terms of the proportion of the population required to consistently define hotspots. (Supplementary movie 2).

Table 2
Results of the impact of sample size on the ability to consistently detect the same structures as being located inside hotspots of malaria infection (PCR prevalence) and exposure (seroprevalence). AUROC = Area Under the Receiver Operator Curve.

Discussion

We present a MBG analysis to define hotspots in an area with known highly heterogeneous transmission in the western Kenyan highlands. Different metrics and sample sizes resulted in variation in the households identified as being located in a hotspot. These results illustrate the uncertainties in determining precise boundaries at the local level that may be relevant for targeted control intervention.

This study utilized two metrics for defining malaria heterogeneity, the prevalence of current infections and the prevalence of serological markers that are indicative of previous exposure. These metrics measure different but analogous facets of malaria transmission and at an individual level are strongly associated22,25. In the current study and other studies determining spatial patterns in malaria transmission, the localization of hotspots based on the two metrics showed considerable overlap but imperfect agreement. Recent or transient malaria hotspots may be missed by serological markers of exposure if antibodies are only detectable following repeated exposure to malaria antigens10,37. However, serologically defined hotspots would be a more stable representation of areas of consistently or historically high risk but may not reflect areas with recent infections38,39. Recent advances in identifying serological markers of recent exposure may improve the agreement between serology and current parasite prevalence40 although it remains to be established whether these are as informative as the currently used antigens that show high sensitivity in low endemic settings41.

For both metrics, only a small percentage of the geographical variability that we observed in our study population was explained by environmental covariates: the spatial residuals estimated by the geostatistical model formed the main component in identifying hotspots. Obtaining a better understanding of the spatial processes driving transmission and identifying covariates accordingly would lead to more precise delineation of MBG-defined hotspots. The specific ecological processing driving transmission will be setting dependent and the necessity of delineating this spatial dynamic will depend on how well these approaches are able to identify hotspots that are meaningful for control and elimination strategies.

As expected, variations in sample size also resulted in significant changes in the boundaries of MBG-defined hotspots. In our setting, intensive sampling of the total population was conducted for accurate hotspot delineation. If less than 20.9% of the total population was sampled (70% of all samples available for our survey), a considerable loss in accuracy was experienced. An average of 37.6% (seroprevalence) and 51.4% (PCR-prevalence) of the structures in hotspots were misclassified as not belonging to hotspots while the complete sample set identified them as hotspots. This suggests that one third to half of structures may be missed by interventions designed to target hotspots or unnecessarily targeted. The sample size thresholds identified here are likely not generalizable to other settings. The purpose of this element of the current study was to illustrate that where the hotspots are drawn will be impacted by the sample size used to inform the analysis. The current findings warn against conducting opportunistic analysis on available geocoded data7,42 if these data were only available for a small fraction of the population.

Although this comparison of metrics cannot determine which is better able to accurately identify and define true hotspots of infection in the community, these results indicate that the approach and assumptions used will affect the resulting map. The MBG approach is generally used to fit a spatial residual risk surface. It allows for a greater understanding of the nature of malaria hotspots by letting the overall risk surface depend on both measured and unmeasured risk factors18. However, the underlying inferential philosophy of MBG is that it is concerned not with how likely it is that a location has an above-average prevalence, but with how likely it is that a given location has a prevalence sufficiently high to be of practical concern in a specific setting. Therefore, the thresholds defined for both the prevalence of concern as well as the probability cut-off that an area has reached or exceeded that must be identified and affects the resulting map. This feature provides useful flexibility so that this approach can easily be tailored to different settings, but also makes it difficult to identify precise hotspots if such policy thresholds do not exist19,43.

Conclusions

This research has highlighted several gaps in our ability to reliably detect hotspots of malaria. The metric and sample size used has important consequences for hotspot boundaries in this setting. The operationally most attractive approach of sampling a small fraction of the population and use the most scalable and economically attractive malaria metric, in our study serology, has limitations in terms of the precision with which hotspot boundaries can be identified. It is important to note that the study setting was characterized by low, heterogeneous but widespread malaria transmission. There may be settings where malaria is more focal, where hotspots are more readily detectable and more consistent between metrics. For our study setting with widespread heterogeneous malaria transmission, we conclude that there are too many uncertainties surrounding hotspot location, stability and boundaries to allow evidence-based targeting of malaria hotspots with the aim of reducing community-wide malaria transmission.

Additional Information

How to cite this article: Stresman, G. H. et al. Impact of metric and sample size on determining malaria hotspot boundaries. Sci. Rep. 7, 45849; doi: 10.1038/srep45849 (2017).

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Supplementary Material

Supplementary Movie 1:
Supplementary Movie 2:
Supplementary Files:

Acknowledgments

We thank all those who contributed to all of the different stages of this work. Specifically, to those at KEMRI/CDC, including Dr. Simon Kariuki, Dr. Meghna Desai and the field and laboratory staff who assisted with the data collection in the field. We would also like to acknowledge Dr. Nuno Sepulvda for providing statistical support. This manuscript has been approved by the Director of the Kenya Medical Research Institute. The data collection was supported by the Bill and Melinda Foundation under the Malaria Transmission Consortium [No: 45114] and the Grand Challenges Grant [No. OPP1024438] awarded to JC and TB, respectively. PJD acknowledges funding support from the UK Medical Research Council [No: G0902153]. This project was partly funded by Grants from the Netherlands Organization for Scientific Research (Vidi fellowship; NWO Project Number 016.158.306.

Footnotes

The authors declare no competing financial interests.

Author Contributions Community survey design T.B., C.D., J.C., J.S.; Preparation and conduct of surveys J.S., G.H.S., A.B., W.O., C.O.; Laboratory analysis A.B., W.O., S.S., E.M.; Design of analytical plan G.H.S., P.J.D., E.G., J.C., T.B., C.D.; Data Analysis G.H.S., E.G., J.C., P.K.; Contributed to manuscript preparation G.H.S., E.G., J.C., T.B., C.D., P.J.D., J.S. All authors read and approved the final manuscript.

References

  • WHO. World Malaria Report 2013. (ed Malaria Control Program) (World Health Organization, Geneva, Switzerland, 2013).
  • WHO. World Malaria Report 2014. (ed World Health Organization) (World Health Organization, Geneva, Switzerland, 2014).
  • Hay S. I., Guerra C. A., Tatem A. J., Noor A. M. & Snow R. W. The global distribution and population at risk of malaria: past, present, and future. The Lancet Infectious Diseases 4, 327–336, doi: (2004).10.1016/s1473-3099(04)01043-6 [PMC free article] [PubMed] [Cross Ref]
  • Griffin J. T. et al. . Potential for reduction of burden and local elimination of malaria by reducing Plasmodium falciparum malaria transmission: a mathematical modelling study. The Lancet Infectious Diseases, doi: (2016).10.1016/s1473-3099(15)00423-5 [PMC free article] [PubMed] [Cross Ref]
  • Bousema T. et al. . Hitting hotspots: spatial targeting of malaria for control and elimination. PLoS medicine 9, e1001165, doi: (2012).10.1371/journal.pmed.1001165 [PMC free article] [PubMed] [Cross Ref]
  • Clements A. C. A., Reid H. L., Kelly G. C. & Hay S. I. Further shrinking the malaria map: How can geospatial science help to achieve malaria elimination? The Lancet infectious diseases 13, 709–718, doi: (2013).10.1016/s1473-3099(13)70140-3 [PubMed] [Cross Ref]
  • Kangoye D. T. et al. . Malaria hotspots defined by clinical malaria, asymptomatic carriage, PCR and vector numbers in a low transmission area on the Kenyan Coast. Malaria journal 15, 213, doi: (2016).10.1186/s12936-016-1260-3 [PMC free article] [PubMed] [Cross Ref]
  • Brown V. L. et al. . Dissecting a wildlife disease hotspot: the impact of multiple host species, environmental transmission and seasonality in migration, breeding and mortality. Journal of the Royal Society, Interface/the Royal Society 10, 20120804, doi: (2013).10.1098/rsif.2012.0804 [PMC free article] [PubMed] [Cross Ref]
  • Tusting L. S., Bousema T., Smith D. L. & Drakeley C. Measuring changes in Plasmodium falciparum transmission: precision, accuracy and costs of metrics. Advances in parasitology 84, 151–208, doi: (2014).10.1016/B978-0-12-800099-1.00003-X [PMC free article] [PubMed] [Cross Ref]
  • Bejon P. et al. . Stable and unstable malaria hotspots in longitudinal cohort studies in Kenya. PLoS medicine 7, e1000304, doi: (2010).10.1371/journal.pmed.1000304 [PMC free article] [PubMed] [Cross Ref]
  • Mosha J. F. et al. . Epidemiology of subpatent Plasmodium falciparum infection: implications for detection of hotspots with imperfect diagnostics. Malaria journal 12, 221, doi: (2013).10.1186/1475-2875-12-221 [PMC free article] [PubMed] [Cross Ref]
  • Nourein A. B. et al. . Identifying residual foci of Plasmodium falciparum infections for malaria elimination: the urban context of Khartoum, Sudan. PloS one 6, e16948, doi: (2011).10.1371/journal.pone.0016948 [PMC free article] [PubMed] [Cross Ref]
  • Kleinschmidt I., Sharp B. L., Clarke G. P. Y., Curtis B. & Fraser C. Use of generalized linear mixed models in the spatial analysis of small-area malaria incidence rates in KwaZulu Natal, South Africa. American journal of epidemiology 153, 1213–1221 (2001). [PubMed]
  • Bousema T. et al. . Identification of hot spots of malaria transmission for targeted malaria control. The Journal of infectious diseases 201, 1764–1774, doi: (2010).10.1086/652456 [PubMed] [Cross Ref]
  • Pullan R. L., Sturrock H. J., Soares Magalhaes R. J., Clements A. C. & Brooker S. J. Spatial parasite ecology and epidemiology: a review of methods and applications. Parasitology 139, 1870–1887, doi: (2012).10.1017/S0031182012000698 [PMC free article] [PubMed] [Cross Ref]
  • Mosha J. F. et al. . Hot spot or not: A comparison of spatial statistical methods to predict prospective malaria infections. Malaria journal 13, 53, doi: (2014).10.1186/1475-2875-13-53 [PMC free article] [PubMed] [Cross Ref]
  • Bejon P. et al. . A micro-epidemiological analysis of febrile malaria in Coastal Kenya showing hotspots within hotspots. eLife 3, e02130, doi: (2014).10.7554/eLife.02130 [PMC free article] [PubMed] [Cross Ref]
  • Diggle P. J., Tawn J. A. & Moyeed R. A. Model-based geostatistics. Appl. Statist 47, 299–350 (1998).
  • Diggle P. J. et al. . Spatial modelling and the prediction of Loa loa risk: decision making under uncertainty. Annals of tropical medicine and parasitology 101, 499–509 (2007). [PubMed]
  • Hay S. I. et al. . A world malaria map: Plasmodium falciparum endemicity in 2007. PLoS medicine 6 (2009). [PMC free article] [PubMed]
  • Elyazar I. R. et al. . Plasmodium falciparum malaria endemicity in Indonesia in 2010. PloS one 6, e21315, doi: (2011).10.1371/journal.pone.0021315 [PMC free article] [PubMed] [Cross Ref]
  • Bousema T. et al. . The impact of hotspot targeted interventions on malaria transmission in Rachuonyo South district in the western Kenyan Highlands: A cluster-randomized controlled trial. PLoS medicine(2016). [PMC free article] [PubMed]
  • Bousema T. et al. . The impact of hotspot-targeted interventions on malaria transmission: study protocol for a cluster-randomized controlled trial. Trials 14, 36, doi: (2013).10.1186/1745-6215-14-36 [PMC free article] [PubMed] [Cross Ref]
  • Stevenson J. C. et al. . Use of different transmission metrics to describe malaria epidemiology in the highlands of western Kenya. Malaria journal 14, 418, doi: (2015).10.1186/s12936-015-0944-4 [PMC free article] [PubMed] [Cross Ref]
  • Baidjoe A. Y. et al. . Factors associated with high heterogeneity of malaria at fine spatial scale in the Western Kenyan highlands. Malaria journal 15, 307, doi: (2016).10.1186/s12936-016-1362-y [PMC free article] [PubMed] [Cross Ref]
  • Snounou G. et al. . High sensitivity of detection of human malaria parasites by the use of nested polymerase chain reaction. Mol Biochem Parasitol 61, 315–320 (1993). [PubMed]
  • Baidjoe A. et al. . Combined DNA extraction and antibody elution from filter papers for the assessment of malaria transmission intensity in epidemiological studies. Malaria journal 12, 272 (2013). [PMC free article] [PubMed]
  • Drakeley C. J. et al. . Estimating medium- and long-term trends in malaria transmission by using serological markers of malaria exposure. Proceedings of the National Academy of Sciences of the United States of America 102, 5108–5113, doi: (2005).10.1073/pnas.0408725102 [PubMed] [Cross Ref]
  • Corran P., Coleman P., Riley E. & Drakeley C. Serology: a robust indicator of malaria transmission intensity? Trends in parasitology 23, 575–582, doi: (2007).10.1016/j.pt.2007.08.023 [PubMed] [Cross Ref]
  • Sepúlveda N., Stresman G., White M. T. & Drakeley C. J. Current Mathematical Models for Analyzing Anti-Malarial Antibody Data with an Eye to Malaria Elimination and Eradication. Journal of Immunology Research 2015, 1–21, doi: (2015).10.1155/2015/738030 [PMC free article] [PubMed] [Cross Ref]
  • Garbrecht J. & Martz L. W. The assignment of drainage direction over flat surfaces in rater digital elevation models. J Hydrology 193, 204–213 (1997).
  • Gleyzer A., Denisyuk M., Rimmer A. & Salingar Y. A fast recursive GIS algorithm for computing Strahler stream order in braided and nonbraided networks. Journal of the American Water Resources Association 40, 937–946 (2004).
  • Giorgi E. & Diggle P. Package ‘PrevMap’. CRAN(2015).
  • Diggle P. J. & Giorgi E. Model-based geostatistics for prevalence mapping in low-resource settings. J Am Stat Assoc(2015).
  • Woolhouse M. E. J. et al. . Heterogeneities in the transmission of infectious agents: Implications for the design of control programs. Proceedings of the National Academy of Sciences of the United States of America 94, 338–342 (1997). [PubMed]
  • Zhang Z. & Huang Y. A Linear Regression Framework for the Receiver Operating Characteristic (ROC) Curve Analysis. Journal of biometrics & biostatistics 3, doi: (2012).10.4172/2155-6180.1000137 [PMC free article] [PubMed] [Cross Ref]
  • Griffin J. T. et al. . Gradual acquisition of immunity to severe malaria with increasing exposure. Proceedings. Biological sciences/The Royal Society 282, 20142657, doi: (2015).10.1098/rspb.2014.2657 [PMC free article] [PubMed] [Cross Ref]
  • Bejon P. et al. . Serological evidence of discrete spatial clusters of Plasmodium falciparum parasites. PloS one 6, e21711, doi: (2011).10.1371/journal.pone.0021711 [PMC free article] [PubMed] [Cross Ref]
  • Cook J. et al. . Using serological measures to monitor changes in malaria transmission in Vanuatu. Malaria journal 9, 169, doi: (2010).10.1186/1475-2875-9-169 [PMC free article] [PubMed] [Cross Ref]
  • Helb D. A. et al. . Novel serologic biomarkers provide accurate estimates of recent Plasmodium falciparum exposure for individuals and communities. Proceedings of the National Academy of Sciences of the United States of America 112, E4438–4447, doi: (2015).10.1073/pnas.1501705112 [PubMed] [Cross Ref]
  • Bousema T. et al. . Serologic markers for detecting malaria in areas of low endemicity, Somalia, 2008. Emerging infectious diseases 16, 392–399, doi: (2010).10.3201/eid1603.090732 [PMC free article] [PubMed] [Cross Ref]
  • Hardy A. et al. . Mapping hotspots of malaria transmission from pre-existing hydrology, geology and geomorphology data in the pre-elimination context of Zanzibar, United Republic of Tanzania. Parasites & vectors 8, 41, doi: (2015).10.1186/s13071-015-0652-5 [PMC free article] [PubMed] [Cross Ref]
  • Ashton R. A. et al. . Geostatistical Modeling of Malaria Endemicity Using Serological Indicators of Exposure Collected Through School Surveys. The American journal of tropical medicine and hygiene 93, 168–177, doi: (2015).10.4269/ajtmh.14-0620 [PMC free article] [PubMed] [Cross Ref]

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group