Large-scale drivers of malaria and priority areas for prevention and control in the Brazilian Amazon region using a novel multi-pathogen geospatial model
- Denis Valle^{1}Email author and
- Joanna M Tucker Lima^{1}
https://doi.org/10.1186/1475-2875-13-443
© Valle and Tucker Lima; licensee BioMed Central Ltd. 2014
Received: 26 August 2014
Accepted: 15 November 2014
Published: 20 November 2014
Abstract
Background
Most of the malaria burden in the Americas is concentrated in the Brazilian Amazon but a detailed spatial characterization of malaria risk has yet to be undertaken.
Methods
Utilizing 2004-2008 malaria incidence data collected from six Brazilian Amazon states, large-scale spatial patterns of malaria risk were characterized with a novel Bayesian multi-pathogen geospatial model. Data included 2.4 million malaria cases spread across 3.6 million sq km. Remotely sensed variables (deforestation rate, forest cover, rainfall, dry season length, and proximity to large water bodies), socio-economic variables (rural population size, income, and literacy rate, mortality rate for children age under five, and migration patterns), and GIS variables (proximity to roads, hydro-electric dams and gold mining operations) were incorporated as covariates.
Results
Borrowing information across pathogens allowed for better spatial predictions of malaria caused by Plasmodium falciparum, as evidenced by a ten-fold cross-validation. Malaria incidence for both Plasmodium vivax and P. falciparum tended to be higher in areas with greater forest cover. Proximity to gold mining operations was another important risk factor, corroborated by a positive association between migration rates and malaria incidence. Finally, areas with a longer dry season and areas with higher average rural income tended to have higher malaria risk. Risk maps reveal striking spatial heterogeneity in malaria risk across the region, yet these mean disease risk surface maps can be misleading if uncertainty is ignored. By combining mean spatial predictions with their associated uncertainty, several sites were consistently classified as hotspots, suggesting their importance as priority areas for malaria prevention and control.
Conclusion
This article provides several contributions. From a methodological perspective, the benefits of jointly modelling multiple pathogens for spatial predictions were illustrated. In addition, maps of mean disease risk were contrasted with that of statistically significant disease clusters, highlighting the critical importance of uncertainty in determining disease hotspots. From an epidemiological perspective, forest cover and proximity to gold mining operations were important large-scale drivers of disease risk in the region. Finally, the hotspot in Western Acre was identified as the area that should receive highest priority from the Brazilian national malaria prevention and control programme.
Keywords
Malaria Brazilian amazon Bayesian model Geospatial model Hotspot Plasmodium falciparum Malaria risk maps Forest cover Gold mining MigrationBackground
The Brazilian Amazon region plays a critical role in the Americas, both in terms of the number of malaria cases and fatalities. Approximately 40% of deaths due to malaria in the Americas occurred in Brazil alone [1] and nearly all malaria cases in Brazil originate from within the Amazon region [1, 2]. Despite the existence of extensive data regularly collected by the Brazilian malaria surveillance system, analyses of large-scale patterns of malaria incidence and its drivers in the Brazilian Amazon are rare. For instance, Olson et al.[3] found that the relationship between precipitation and malaria incidence tends to be negative in areas with extensive wetlands, whereas no discernible pattern was evident in counties with few or no wetlands. Hahn et al.[4] revealed that roads, forest fires and selective logging are important risk factors for malaria, but they found no association between deforestation rate and malaria. On the other hand, Achcar et al. [5] found a positive relationship between deforestation rate and malaria incidence as well as a negative relationship between malaria incidence and human development index and population density. Finally, Valle and Clark [6] found that higher malaria incidence is frequently observed in areas with greater forest cover. Determining the key drivers of malaria risk in this region is particularly important given current and future large-scale environmental transformations due to ongoing expansion of the road network and the construction of multiple hydro-electric dams [7, 8].
Currently, the spatial characterization of malaria risk in the Brazilian Amazon is based on the calculation of the incidence per year per thousand inhabitants (i.e., annual parasite index - API) for each county, which is categorized into low, medium and high API groups and displayed in choropleth maps [2, 9, 10]. These results are used by the Brazilian Health Ministry to prioritize counties for malaria prevention and control measures. Unfortunately, it is widely known that this type of data might have substantial noise when associated with small population sizes (as is the case for most Amazonian counties) and small case counts (as is the case for Plasmodium falciparum) [11]. The accurate mapping of P. falciparum is particularly critical from a public health perspective because, despite its relative rarity, P. falciparum infection tends to result in more severe outcomes than infections by Plasmodium vivax, the dominant type of malaria in this region [12]. Geospatial models have been widely used to ameliorate some of these problems, allowing for estimates of disease risk at a particular site to be influenced by neighbouring sites (i.e., borrowing of information across sites through spatial correlation).
In this article, the concept of borrowing information is taken one step further by allowing for estimates of P. falciparum disease risk to be influenced by P. vivax disease risk. Given that both Plasmodia are typically transmitted by the same set of vectors, sharing information across pathogens (i.e., modelling between-pathogen correlation) should improve predictions of P. falciparum disease risk. A novel multi-pathogen geospatial Bayesian model is introduced that accommodates for both between-site and between-pathogen correlation. This model aims to circumvent the problems associated with determining the API of counties with small population sizes and for relatively rare pathogens (e.g., malaria caused by P. falciparum). Model results uncover large-scale drivers of malaria incidence, spatially characterize malaria risk across the Brazilian Amazon and help identify key priority areas for malaria prevention and control interventions in the region.
Methods
Data
Several remotely sensed environmental covariates were gathered from multiple sources, including annual deforested and forested area, rainfall and length of dry season, and proximity to large water bodies. Annual deforested and forested area are based on a semi-automated annual analysis of Landsat imagery, provided by the Brazilian Institute of Space Research (INPE) [16]. Monthly rainfall data for 0.25^{o} × 0.25^{o} pixels come from the Tropical Rainfall Measuring Mission (TRMM) [17]. Precipitation data was used to calculate average monthly rainfall and the length of dry season (number of months with less than 100 mm of rainfall) [18, 19] for each year and county. To map large water bodies, a 250 m resolution global raster water mask, derived using MODIS and Shuttle Radar Topography Mission (SRTM) data was used [20]. In addition to these covariates, proximity to gold mining operations [21] and proximity to other large-scale infrastructure, such as roads and hydro-electric dams in operation, were also considered [22]. All covariates were assessed within a 50-km radius of each main city, encompassing the area where most of the population likely resides (see Additional file 1). Proximity-based covariates were transformed into dummy variables, taking on the value of 1 if the item (e.g., large water body, gold mining operation, road, or hydro-electric dam) was present within the 50-km buffer and zero otherwise.
Several socio-economic covariates based on the 2000 Brazilian Census, including income (i.e., average income of rural households), health (i.e., mortality rate for children under five years old), education (i.e., proportion of the rural population that is literate), and rural population size were incorporated in the model to understand the underlying large-scale drivers of malaria risk [23]. Finally, 2010 Census data were used to assess migration patterns, given by the proportion of people that arrived in the county over the previous three to five years. The population size for each county in each month was estimated by linearly interpolating the population counts from the 2000, 2007 and 2010 Brazilian Censuses. These population estimates were then used as offsets in the statistical model.
Statistical model
where ${\mathbf{x}}_{iy}^{T}$ is the design vector, containing covariate information, β_{ j } is the vector of regression parameters for pathogen j, and e_{ iyj } is the residual. Different assumptions can be made regarding the nature of these residuals. For instance, if one assumes that ${e}_{iyj}~N\left(0,{\sigma}_{j}^{2}\right)$ and that these residuals are independent, the resulting model would be equivalent to two independent regressions, one for each pathogen j.
Similar to most geospatial models, the correlation of the residuals for sites i and k was assumed to be given by Cor[e_{ iyj }, e_{ kyj }] = exp(-ϕD_{ ik }), where D_{ ik } is the distance between these sites. This expression describes the correlation among residuals within the same year y and pathogen j as a decay function of the distance between sites i and k. The innovation in the proposed model consists of further assuming that the correlation for pathogens j and l is given by Cor[e_{ iyj }, e_{ iyl }] = ρ. This assumption allows for information on pathogen j to influence, and potentially improve, inference on pathogen l. Such a feature might be useful if incidence of these diseases is likely to be correlated (e.g., because they are transmitted by the same vector).
where ⊗ denotes the Kronecker product. The original design matrix X_{ y } contains information regarding the covariates for all sites in year y, where each continuous covariate was standardized to have zero mean and standard deviation of one and enters the model with a linear and quadratic term. Based on the design matrix, ${\mathbf{X}}_{y}^{*}$ is defined to be equal to $\left[\begin{array}{cc}\hfill {\mathbf{X}}_{y}\hfill & \hfill \mathbf{0}\hfill \\ \hfill \mathbf{0}\hfill & \hfill {\mathbf{X}}_{y}\hfill \end{array}\right]$ and, correspondingly, the associated vector of slope parameters β* is defined to be equal to $\left[\begin{array}{c}\hfill {\mathbf{\beta}}_{1}\hfill \\ \hfill {\mathbf{\beta}}_{2}\hfill \end{array}\right]$, where β_{1} and β_{2} are the vectors for pathogens 1 and 2, respectively, This specification for ${\mathbf{X}}_{y}^{*}$ and β* allows for the different pathogens to have different regression parameters. The covariance matrix of this multivariate normal, given by σ^{2}P(ρ) ⊗ S(ϕ), is separable into an overall variance σ^{2} parameter, a pathogen correlation matrix P(ρ), and a spatial correlation matrix S(ϕ). Similar separable correlation matrices are often used for spatial-temporal models [24].
where I is the identity matrix. Priors on the variance parameters v, σ were chosen based on the recommendation from [25]. Finally, similar to [26], a discrete support for ϕ was assumed to aid mixing of a Monte Carlo Markov Chain (MCMC) algorithm. This was implemented by assuming a uniform prior over a set of values for φ, where these values were chosen based on a preliminary analysis.
To determine the large-scale drivers of malaria risk in the region, an adaptation of the methodology described in [27] was used to first determine which socio-economic and environmental covariates were significantly associated with malaria incidence. This helped avoid over-fitting the model given the large number of covariates being evaluated. The model described above was then fit using the predictors identified to be important by that analysis. Spatial predictions, on the other hand, required values for all covariates at all prediction locations. Thus, only spatially explicit environmental covariates were used for spatial predictions, ignoring the socio-economic covariates. Finally, a cross-validation exercise (described below) was performed to determine the out-of-sample predictive skill of the model when applied to create interpolated surfaces of P. falciparum disease risk. Again, only spatially explicit environmental covariates were used for this validation exercise. For both the spatial predictions and the cross-validation exercise, linear and quadratic terms of these spatially explicit environmental variables were included, whether or not they were significant.
Cross-validation
Description of the models compared in the ten-fold cross-validation exercise
Models | Assumptions | Parameters |
---|---|---|
1 – No correlation | No spatial and no pathogen correlation | ρ = 0; Φ > > 0 |
2 – Spatial | Accounts for spatial correlation but not for pathogen correlation | ρ = 0; Φ estimated |
3 – Pathogen | Accounts for pathogen correlation but not for spatial correlation | ρ is estimated; Φ > > 0 |
4 – Pathogen + Spatial | Accounts for pathogen and spatial correlation | ρ and Φ are both estimated |
To determine the out-of-sample predictive skill of these models for P. falciparum, a ten-fold cross-validation was performed. This procedure entails the random exclusion of 10% of P. falciparum observations and the use of the rest of the data to fit the four models described in Table 1. Then, left-out observations were predicted with the fitted models. This was repeated ten times, each time leaving out a distinct subset of the data. Predictive skill was assessed using two criteria: (1) the proportion of left-out observations that were within the 95% predictive credible interval (95% CI coverage). For good models, the proportion of left-out observations within this interval should be close to 95%; and, (2) mean-squared-error (MSE) of the predicted malaria incidence per person. This criterion takes into account both bias and uncertainty in predictions and smaller MSE values indicate better models. All models were fit using a customized Gibbs sampler, with Metropolis-within-Gibbs steps for non-conjugate conditional distributions. All analyses were conducted in R [28] and the R code to fit the statistical model is available upon request.
Hotspot and coldspot definitions
To complement the geospatially modelled malaria risks maps, hotspots and coldspots were mapped as a way to take into account both mean disease risk and its associated uncertainty. Hotspots are defined as sites that were predicted to have significantly higher disease risk than average. In other words, site i was deemed to be a hotspot for pathogen j if the 95% credible interval of α_{ iyj } did not straddle ${\overline{\alpha}}_{j}$ (i.e., the average disease risk across all years and counties for pathogen j). Furthermore, persistent hotspots are those areas classified as a hotspot in all five years, intermittent hotspots are areas classified as a hotspot in at least one but not all years, and potential hotspots are areas predicted to have very high malaria risk (mean risk greater than 95% of the risk in all locations) but whose 95% CI encompassed the average risk. Similar definitions were applied for coldspots (i.e., areas with significantly lower disease risk).
Results
Cross-validation results
Validation results suggest that between-pathogen correlation is more important for out-of-sample predictive skill for this dataset than spatial correlation (i.e., compare the drop in out-of-sample predictive ability of the ‘Spatial’ model that neglects pathogen correlation relative to that of the ‘Pathogen’ model that neglects spatial correlation). These findings indicate that data on P. vivax incidence can substantially improve predictions for P. falciparum incidence.
Malaria findings
Average income of rural households and length of dry season were positively associated with malaria risk. Regions with higher average rural income may have greater access to health posts, potentially increasing malaria reporting rates relative to poorer counties. Findings regarding dry season length suggest that the formation of small and/or isolated water bodies when rivers recede during the dry season likely increases mosquito breeding habitat. Deforestation rate, annual rainfall, rural population size, rural literacy rate, child mortality rate (< five years), and proximity to roads or large water bodies, had relatively little effect on malaria risk. Surprisingly, no significant effect of proximity to hydro-electric dams was found, which could be due to the relatively coarse scale of the analysis or because of the enhanced malaria prevention and control programmes in these regions. Parameter estimates with associated 95% CIs for the final model are given in Table A1 in the Additional file 1. The model’s between-pathogen correlation parameter was high (E[ρ] = 0.86 with 95% CI 0.85-0.88) supporting the conclusion from the validation exercise that models accounting for between-pathogen correlation improve predictive performance. On the other hand, spatial correlation was small (< 0.16) for sites 100 km apart and negligible (< 0.05) for sites more than 160 km apart.
Approximate location of persistent malaria hotspots by state
State (acronym) | Both | P. vivax only | P. falciparum only |
---|---|---|---|
Acre (AC) | Rodrigues Alves, Mancio Lima, Cruzeiro do Sul | ||
Amapa (AP) | Calcoene, Oiapoque, Tartarugalzinho | Ferreira Gomes, Mazagao, Porto Grande, Santana | |
Amazonas (AM) | Careiro, Guajara, Itapiranga | Iranduba, Presidente Figueiredo, Silves | Sao Sebastiao do Uatuma, Urucara |
Para (PA) | Jacareacanga | Pacaja | Anajas, Chaves, Itaituba, Novo Progresso |
Rondonia (RO) | Porto Velho, Rio Crespo, Alto Paraiso, Campo Novo de Rondonia, Cujubim | Machadinho d’Oeste, Candeias do Jamari, Jamari | |
Roraima (RR) | Bonfim, Canta, Mucajai | Amajari, Normandia, Rorainopolis |
Discussion
Multiple large-scale infrastructure projects have been planned for future years to fully integrate the Amazon region (e.g., Planos de Crescimento Acelerado and Iniciativa para a Integracao de Infra-estrutura Regional da America do Sul – IIRSA). These projects include the expansion of the road network and the construction of multiple new hydro-electric dams, both in the Brazilian Amazon and in neighbouring countries [7, 8]. Interestingly, despite the widespread concern regarding increased malaria incidence due to the construction of dams (e.g., the Brazilian Government requires malaria studies prior to dam construction in the region), there was no association between proximity to these dams and malaria incidence. This lack of association might be due to the coarse scale of analysis and/or because of enhanced malaria surveillance and control in the vicinity of dams, typically initiated and sponsored by companies that build and operate these dams as a component of their impact mitigation plan. Not only did proximity to dams show no effect, forest cover results suggest that large-scale development projects likely induce a long-term decline of malaria incidence as a result of substantial deforestation associated with infrastructure development, similar to the conclusion in [6]. Nevertheless, these results should be viewed with caution due to the inherent limitations of observational studies and absence of long-term data. Furthermore, increased deforestation may have several other detrimental public-health effects. Finally, considering only the proximity to a dam reservoir or road fails to account for possible increases in malaria risk related to dam construction and/or road creation.
Proximity to gold mining operations exhibited a strong effect on malaria risk for both P. vivax and P. falciparum. Although this association has been described repeatedly in the literature e.g., [29–31], the detection of large-scale association between malaria incidence and gold mining is relatively rare [32]. It is important to note that the data consist of officially registered gold mining operations only. These data are used as a proxy for the location of gold mining activities in general because no comprehensive data on the location of small-scale artisanal gold mining operations (i.e., garimpos) at the scale of this analysis currently exists. Findings linking migration patterns with increased malaria incidence further support the importance of gold mining in increasing malaria risk, given that these activities typically spur considerable migration to the region.
Mapping disease risk is fundamental to public health, enabling strategic decision-making, such as spatially prioritizing resource allocation and deciding where to conduct more intensive interventions. Model predictions revealed substantial spatial heterogeneity across the Brazilian Amazon, highlighting regions that consistently have very high malaria risk. While these spatial predictions can be useful, relying on mean predictions alone to identify priority areas might be unreliable if their associated uncertainty is ignored. By identifying areas that are significantly and persistently above average, a handful of priority regions were highlighted (see Table 2). The importance of accounting for uncertainty in defining hotspots is particularly evident when comparing hotspots in Figure 5 with the traditional choropleth maps used by the Ministry of Health to identify areas of high malaria risk e.g., [33]. These discrepancies emphasize the critical role of uncertainty when creating interpolated surfaces of disease risk and the importance of carefully choosing the outcome to display in maps to effectively inform and guide public decision-making (see more discussion on this topic in [34]). The identified priority sites merit further scrutiny from the scientific community to better determine additional driving factors behind malaria risk in these persistent hotspots.
Results suggest that malaria incidence, both for P. vivax and P. falciparum, generally declined between 2004 and 2008 throughout most of the Brazilian Amazon region, except for Amazonas state. Sites with significant increases in the estimated malaria risk from 2004 to 2008 tended to be sites where transportation is typically limited to boats, suggesting that logistical factors may be hindering malaria prevention and control. Together with the hotspot analysis, this temporal analysis suggests that the Brazilian Government should prioritize the region in western Acre (i.e., Mancio Lima, Rodrigues Alves and Cruzeiro do Sul) for malaria prevention and control.
In this work, each county is represented by a single point at its urban centre (see Additional file 1). An alternative geospatial modelling approach that avoids the problem of determining the exact location of people and malaria infections relies on a conditionally autoregressive (CAR) prior for the spatial random effects. In such a model, spatial correlation is modelled based on an adjacency matrix, where two counties are deemed neighbours if they share a boundary. This approach is problematic in the Brazilian Amazon since human populations are highly clustered around county seats. Because counties are large, assuming that two counties are neighbours regardless of the distance between county seats seems unrealistic (e.g., [35]). Another problem with a CAR model is that the implied correlation matrix is often very counterintuitive [36]. Finally, CAR models do not allow for fine-scale geographical prediction, implicitly assuming that risk is the same throughout the county, and the resulting display of results using choropleth maps tends to over-emphasize larger counties. On the other hand, an important limitation of models like the one presented here is that implementation via Markov Chain Monte Carlo (MCMC) algorithms is computationally expensive due to the size of the covariance matrix [37]. While this problem is alleviated by using a separable correlation matrix, computational burden remains a critical issue, particularly if the model is expanded to include more than two pathogens and/or correlation between years. In future research, this problem might be tackled through the use of low-rank representations for spatial processes [38, 39].
An important limitation of the data used in this paper is that it lacks information on the total number of malaria examinations performed [40]. As a consequence, areas with a more comprehensive malaria surveillance network may apparently have higher malaria risk because of the detection of a greater number of cases. Thus, the high correlation between P. falciparum and P. vivax incidence identified by the model may be partly attributable to site-specific differences in surveillance intensity (i.e., enhanced surveillance will increase reported incidence of malaria cases from both pathogens). Regardless of the source of correlation between pathogens, joint modelling of data from multiple pathogens is still likely to be beneficial, particularly for diseases that are closely related (e.g., transmitted by the same vector, such as dengue, yellow fever and chikungunya). For instance, cross-validation results suggest that P. vivax incidence provides considerable information on P. falciparum incidence, potentially improving the ability to infer disease risk for this less common but more severe type of malaria. Only a few statistical models exploit information from data on multiple pathogens [41, 42], suggesting that the idea of integrating data from multiple pathogens to improve model predictions has been underexplored.
The gross underestimation of malaria symptomatic cases is not an issue in this study because malaria drugs are only available after obtaining positive microscopy detection in a government health facility. However, asymptomatic individuals and subpatent infections are likely to be missed by the data [43–45], and thus a map of infection risk might be substantially different from a disease risk map. Furthermore, P. vivax and Plasmodium malariae are commonly confounded in microscopic examination of thick blood smears [46–48], and as a result, findings regarding P. vivax are expected to be generally accurate but may falter in locations with high P. malariae incidence. Finally, in the absence of more systematically collected data on this topic, the role of concurrent local malaria interventions has been ignored.
Conclusion
This article provides both methodological and epidemiological contributions. From a methodological perspective, the benefits of jointly modelling multiple pathogens for spatial predictions were illustrated, where abundant data from one pathogen can improve predictions of incidence of a related, rarer pathogen. Maps of mean disease risk (the usual end product of geospatial analyses) were contrasted with that of statistically significant disease clusters, highlighting the critical importance of uncertainty in determining disease hotspots. From an epidemiological perspective, large-scale drivers and spatial patterns of malaria incidence were revealed for the Brazilian Amazon. For instance, areas with higher forest cover and close to gold mining operations tended to have higher malaria incidence both for P. vivax and P. falciparum, suggesting that these are important large-scale drivers of disease risk in the region. Finally, key priority areas for malaria prevention and control were identified in this article. In particular, using both spatial and temporal results, the hotspot comprised of Mancio Lima, Rodrigues Alves, and Cruzeiro do Sul (Acre) should receive highest priority from the Brazilian national malaria prevention and control programme.
Declarations
Acknowledgements
We thank Phil Lounibos and Paul Psychas for providing numerous comments on an earlier version of this manuscript. This research was partially funded by a University of Florida Research Opportunity Seed Fund.
Authors’ Affiliations
References
- Partners RBM: The global malaria action plan for a malaria-free world. 2008, Geneva: World Health OrganizationGoogle Scholar
- Oliveira-Ferreira J, Lacerda MVG, Brasil P, Ladislau JLB, Tauil PL, Daniel-Ribeiro CT: Malaria in Brazil: an overview. Malar J. 2010, 9: 115-10.1186/1475-2875-9-115.PubMed CentralView ArticlePubMedGoogle Scholar
- Olson SH, Gangnon RE, Elguero E, Durieux L, Guegan J-F, Foley JA, Patz JA: Links between climate, malaria, and wetlands in the Amazon Basin. Emerging Infect Dis. 2009, 15: 659-662. 10.3201/eid1504.080822.PubMed CentralView ArticlePubMedGoogle Scholar
- Hahn MB, Gangnon RE, Barcellos C, Asner GP, Patz JA: Influence of deforestation, logging, and fire on malaria in the Brazilian Amazon. PLOS One. 2014, 9: e85725-10.1371/journal.pone.0085725.PubMed CentralView ArticlePubMedGoogle Scholar
- Achcar JA, Martinez EZ, de Souza ADP, Tachibana VM, Flores EF: Use of Poisson spatiotemporal regression models for the Brazilian Amazon forest: malaria count data. Rev Soc Bras Med Trop. 2011, 44: 749-754. 10.1590/S0037-86822011000600019.View ArticlePubMedGoogle Scholar
- Valle D, Clark J: Conservation efforts may increase malaria burden in the Brazilian Amazon. PLOS One. 2013, 8: e57519-10.1371/journal.pone.0057519.PubMed CentralView ArticlePubMedGoogle Scholar
- Finer M, Jenkins CN: Proliferation of hydroelectric dams in the Andean Amazon and implications for Andes-Amazon connectivity. PLOS One. 2012, 7: e35126-10.1371/journal.pone.0035126.PubMed CentralView ArticlePubMedGoogle Scholar
- Ministerio de Minas e Energia: Plano decenal de expansao de energia 2020. 2011, Brasilia: Distrito FederalGoogle Scholar
- Barreto ML, Teixeira MG, Bastos FI, Ximenes RAA, Barata RB, Rodrigues LC: Successes and failures in the control of infectious diseases in Brazil: social and environmental context, policies, interventions, and research needs. Lancet. 2011, 377: 1877-1889. 10.1016/S0140-6736(11)60202-X.View ArticlePubMedGoogle Scholar
- WHO: World Malaria Report 2011. 2011, Geneva: World Health OrganizationGoogle Scholar
- Waller LA, Gotway CA: Applied spatial statistics for public health data. 2004, Hoboken, NJ: John Wiley & SonsView ArticleGoogle Scholar
- Carlton J, Sina B, Adams J: Why is Plasmodium vivax a neglected tropical disease?. PLOS Negl Trop Dis. 2011, 5: e1160-10.1371/journal.pntd.0001160.PubMed CentralView ArticlePubMedGoogle Scholar
- Laboratorio de Informacao em Saude ICICT-Fiocruz: Observatorio Climate e Saude.http://www.climasaude.icict.fiocruz.br/,
- Parry L, Day B, Amaral S, Peres C: Drivers of rural exodus from Amazonian headwaters. Popul Environ. 2010, 32: 137-176. 10.1007/s11111-010-0127-8.View ArticleGoogle Scholar
- Parry L, Peres C, Day B, Amaral S: Rural-urban migration brings conservation threats and opportunities to Amazonian watersheds. Conserv Lett. 2010, 3: 251-259. 10.1111/j.1755-263X.2010.00106.x.View ArticleGoogle Scholar
- INPE, Instituto Nacional de Pesquisas Espaciais: Projeto Prodes: Monitoramento da floresta Amazonica Brasileira por satelite.http://www.obt.inpe.br/prodes/index.php,
- NASA: 3B43: Monthly 0.25 × 0.25 degree merged TRMM and other sources estimates.http://mirador.gsfc.nasa.gov/cgi-bin/mirador/presentNavigation.pl?tree=project&project=TRMM&dataGroup=Gridded&dataset=3B43:%20Monthly%200.25%20×%200.25%20degree%20merged%20TRMM%20and%20other%20sources%20estimates&version=006,
- Baker IT, Prihodko L, Denning AS, Goulden M, Miller S, da Rocha HR: Seasonal drought strees in the Amazon: reconciling models and observations. J Geophys Res. 2008, 113: G00B01-Google Scholar
- Saleska SR, Miller SD, Matross DM, Goulden ML, Wofsy SC, da Rocha HR, de Camargo PB, Crill P, Daube BC, de Freitas HC, Hutyra L, Keller M, Kirchhoff V, Menton M, Munger JW, Pyle EH, Rice AH, Silva H: Carbon in Amazon forests: unexpected seasonal fluxes and disturbance-induced losses. Science. 2003, 302: 1554-1557. 10.1126/science.1091165.View ArticlePubMedGoogle Scholar
- Carroll M, Townshend J, DiMiceli C, Noojipady P, Sohlberg R: A new global raster water mask at 250 m resolution. Int J Digital Earth. 2009, 2: 291-308. 10.1080/17538940902951401.View ArticleGoogle Scholar
- DNPM: Departamento Nacional de Producao Mineral.http://sigmine.dnpm.gov.br/webmap/,
- ANEEL: Brazilian Electricity Regulatory Agency.http://sigel.aneel.gov.br/kmz.html,
- SIDRA - IBGE: Censo Demografico e Contagem da Populacao.http://www.sidra.ibge.gov.br,
- Banerjee S, Carlin BP, Gelfand AE: Hierarchical modeling and analysis of spatial data. 2004, Boca Raton, FL: Chapman & Hall/CRCGoogle Scholar
- Gelman A: Prior distributions for variance parameters in hierarchical models. Bayesian Anal. 2006, 1: 515-533.View ArticleGoogle Scholar
- Wikle CK: Hierarchical modeling with spatial data. Handbook of Spatial Statistics. Edited by: Gelfand AE, Diggle PJ, Fuentes M, Guttorp P. 2010, Boca Raton, FL: CRC PressGoogle Scholar
- Zhao K, Valle D, Popescu S, Zhang X, Mallick B: Hyperspectral remote sensing of plant biochemistry using Bayesian model averaging with variable and band selection. Remote Sens Environ. 2013, 132: 102-119.View ArticleGoogle Scholar
- R Core Team: R: A language and environment for statistical computing. 2013, Vienna, Austria: R Foundation for Statistical ComputingGoogle Scholar
- Crompton P, Ventura AM, de Souza JM, Santos E, Strickland GT, Silbergeld E: Assessment of mercury exposure and malaria in a Brazilian Amazon riverine community. Environ Res. 2002, 90: 69-75. 10.1006/enrs.2002.4358.View ArticlePubMedGoogle Scholar
- Barbieri AF, Sawyer DO, Soares-Filho BS: Population and land use effects on malaria prevalence in the Southern Brazilian Amazon. Hum Ecol. 2005, 33: 847-874. 10.1007/s10745-005-8213-8.View ArticleGoogle Scholar
- de Oliveira EC, dos Santos ES, Zeilhofer P, Souza-Santos R, Atanaka-Santos M: Spatial patterns of malaria in a land reform colonization project, Juruena municipality, Mato Grosso, Brazil. Malar J. 2011, 10: 177-10.1186/1475-2875-10-177.PubMed CentralView ArticlePubMedGoogle Scholar
- Braz RM, Duarte EC, Tauil PL: Caracterizacao das epidemias de malaria nos municipios da Amazonia Brasileira em 2010. Cad Saude Publica. 2013, 29: 935-944. 10.1590/S0102-311X2013000500011.View ArticlePubMedGoogle Scholar
- da Silva-Nunes M: Impacto de alteracoes ambientais na transmissao da malaria e perspectivas para o controle da doenca em areas de assentamento rural da Amazonia Brasileira. Oecologia Australis. 2010, 14: 603-622. 10.4257/oeco.2010.1403.02.View ArticleGoogle Scholar
- Gelman A, Price P: All maps of parameter estimates are misleading. Stat Med. 1999, 18: 3221-3234. 10.1002/(SICI)1097-0258(19991215)18:23<3221::AID-SIM312>3.0.CO;2-M.View ArticlePubMedGoogle Scholar
- Waller L, Carlin B: Disease Mapping. Handbook of spatial statistics. Edited by: Gelfand AE, Diggle PJ, Fuentes M, Guttorp P. 2010, Boca Raton, FL: Taylor & Francis/CRC, 217-244.View ArticleGoogle Scholar
- Wall MM: A close look at the spatial structure implied by the CAR and SAR models. J Stat Plan Inference. 2004, 121: 311-324. 10.1016/S0378-3758(03)00111-3.View ArticleGoogle Scholar
- Best N, Richardson S, Thomson A: A comparison of Bayesian spatial models for disease mapping. Stat Methods Med Res. 2005, 14: 35-59. 10.1191/0962280205sm388oa.View ArticlePubMedGoogle Scholar
- Wikle CK: Low-rank representations for spatial processes. Handbook of Spatial Statistics. Edited by: Gelfand AE, Diggle PJ, Fuentes M, Guttorp P. 2010, Boca Raton, FL: Taylor & Francis/CRC, 107-118.View ArticleGoogle Scholar
- Banerjee S, Gelfand AE, Finley AO, Sang H: Gaussian predictive process models for large spatial data sets. J R Stat Soc Series B Stat Methodol. 2008, 70: 825-848. 10.1111/j.1467-9868.2008.00663.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Valle D, Clark J: Improving the modeling of disease data from the government surveillance system: a case study on malaria in the Brazilian Amazon. PLOS Comput Biol. 2013, 9: e1003312-10.1371/journal.pcbi.1003312.PubMed CentralView ArticlePubMedGoogle Scholar
- Gelfand AE, Vounatsou P: Proper multivariate conditional autoregressive models for spatial data analysis. Biostatistics. 2003, 4: 11-25. 10.1093/biostatistics/4.1.11.View ArticlePubMedGoogle Scholar
- Jin X, Carlin BP, Banerjee S: Generalized hierarchical multivariate CAR models for areal data. Biometrics. 2005, 61: 950-961. 10.1111/j.1541-0420.2005.00359.x.View ArticlePubMedGoogle Scholar
- Alves FP, Durlacher RR, Menezes MJ, Krieger H, da Silva LHP, Camargo EP: High prevalence of asymptomatic Plasmodium vivax and Plasmodium falciparum infections in native Amazonian populations. Am J Trop Med Hyg. 2002, 66: 641-648.PubMedGoogle Scholar
- Ladeia-Andrade S, Ferreira MU, de Carvalho ME, Curado I, Coura JR: Age-dependent acquisition of protective immunity to malaria in riverine populations of the Amazon Basin of Brazil. Am J Trop Med Hyg. 2009, 80: 452-459.PubMedGoogle Scholar
- Valle D, Clark JS, Zhao K: Enhanced understanding of infectious diseases by fusing multiple datasets: a case study on malaria in the Western Brazilian Amazon region. PLOS One. 2011, 6: e27462-10.1371/journal.pone.0027462.PubMed CentralView ArticlePubMedGoogle Scholar
- Scopel KKG, Fontes CJF, Nunes AC, Horta MF, Braga EM: High prevalence of Plasmodium malariae infections in a Brazilian Amazon endemic area (Apiacas - Mato Grosso State) as detected by polymerase chain reaction. Acta Trop. 2004, 90: 61-64. 10.1016/j.actatropica.2003.11.002.View ArticlePubMedGoogle Scholar
- Cavasini MTV, Ribeiro WL, Kawamoto F, Ferreira MU: How prevalent is Plasmodium malariae in Rondonia, Western Brazilian Amazon. Rev Soc Bras Med Trop. 2000, 33: 489-492.View ArticlePubMedGoogle Scholar
- Galardo AKR, Arruda M, Couto AARD, Wirtz R, Lounibos LP, Zimmerman RH: Malaria vector incrimination in three rural riverine villages in the Brazilian Amazon. Am J Trop Med Hyg. 2007, 76: 461-469.PubMedGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Comments
By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.