A country bug in the city: urban infestation by the Chagas disease vector Triatoma infestans in Arequipa, Peru
- Stephen Delgado1, 2,
- Kacey C Ernst2,
- María Luz Hancco Pumahuanca3,
- Stephen R Yool1,
- Andrew C Comrie1,
- Charles R Sterling4,
- Robert H Gilman3, 5,
- César Náquira3,
- Michael Z Levy6Email author and
- the Chagas Disease Working Group in Arequipa, Peru
© Delgado et al.; licensee BioMed Central Ltd. 2013
Received: 6 August 2013
Accepted: 18 October 2013
Published: 30 October 2013
Interruption of vector-borne transmission of Trypanosoma cruzi remains an unrealized objective in many Latin American countries. The task of vector control is complicated by the emergence of vector insects in urban areas.
Utilizing data from a large-scale vector control program in Arequipa, Peru, we explored the spatial patterns of infestation by Triatoma infestans in an urban and peri-urban landscape. Multilevel logistic regression was utilized to assess the associations between household infestation and household- and locality-level socio-environmental measures.
Of 37,229 households inspected for infestation, 6,982 (18.8%; 95% CI: 18.4 – 19.2%) were infested by T. infestans. Eighty clusters of infestation were identified, ranging in area from 0.1 to 68.7 hectares and containing as few as one and as many as 1,139 infested households. Spatial dependence between infested households was significant at distances up to 2,000 meters. Household T. infestans infestation was associated with household- and locality-level factors, including housing density, elevation, land surface temperature, and locality type.
High levels of T. infestans infestation, characterized by spatial heterogeneity, were found across extensive urban and peri-urban areas prior to vector control. Several environmental and social factors, which may directly or indirectly influence the biology and behavior of T. infestans, were associated with infestation. Spatial clustering of infestation in the urban context may both challenge and inform surveillance and control of vector reemergence after insecticide intervention.
KeywordsTriatoma infestans Chagas disease Urban infestation Vector control Spatial analysis Multilevel logistic regression
Chagas disease, also known as American trypanosomiasis, is caused by the protozoan parasite Trypanosoma cruzi and is endemic in Latin America. Typically, T. cruzi is transmitted to humans via the infected excretions of various blood-feeding triatomine insect species, including Triatoma infestans. Less commonly, infection may result from congenital transmission, blood transfusion, organ transplantation, and incidental ingestion of parasite-contaminated food or drink . Chagas disease is characterized by an acute phase, which lasts 6 – 8 weeks, and a chronic phase, which persists for life. In most cases, both the acute and chronic phases of infection are asymptomatic. However, 10 – 40% of cases, depending on the geographic region, progress over a period of years to chronic disease, including potentially fatal cardiac and gastrointestinal disorders . Globally, eight million persons are infected with T. cruzi, resulting in 11,000 deaths and the loss of 430,000 disability-adjusted life years (DALYs) annually [3, 4].
Prevention and control of Chagas disease are achieved primarily via large-scale insecticide application initiatives , and T. infestans, which lives predominantly in and around human households, is a principal target for vector control . While the Southern Cone Initiative has succeeded in interrupting T. cruzi transmission by T. infestans in Brazil, Chile, and Uruguay [6, 7], T. infestans-mediated T. cruzi transmission persists in parts of Argentina, Bolivia, Paraguay, and Peru [8, 9]. Moreover, efforts to eliminate T. infestans have been complicated by the expansion of this species from sparsely populated rural regions into densely populated urban areas [10, 11].
While the ecology of T. infestans infestation in rural environments has been studied extensively, urban infestation by this insect has been investigated only preliminarily. To advance understanding of vector infestation in the urban context, this study examines infestation by T. infestans across an urban and peri-urban landscape prior to implementation of vector control. Our study focuses on Arequipa, Peru, which lies in an area with extensive T. infestans infestation and epidemic T. cruzi transmission [8, 12]. Utilizing data collected from multiple sources at multiple spatial scales, we employ spatial point pattern analysis and multilevel logistic regression to elucidate spatial patterns in infestation by T. infestans and to assess associations between several environmental and social factors and T. infestans infestation in an urban landscape. In particular, the effects of housing density, land surface temperature, and elevation were evaluated at the household level, while the effects of urban shantytowns, which have been identified as areas at higher risk for infestation by Chagas diesease vectors and vector-borne transmission of T. cruzi, were evaluated at the locality level [11, 13–15].
Results of univariate and multivariate logistic regression
Univariate logistic regression
Multivariate logistic regression
Housing density (households/hectare)
AIC = 32,993
2 – 24*
24 – 30
1.33 – 1.61
1.26 – 1.53
30 – 34
1.75 – 2.11
1.63 – 1.97
34 – 39
1.99 – 2.39
1.84 – 2.22
39 – 77
1.91 – 2.30
1.73 – 2.09
Land surface temperature (°C)
AIC = 33,300
26 – 40
1.06 – 1.10
1.08 – 1.12
Elevation (meters above sea level)
AIC = 33,134
2,120 – 2,260*
2,260 – 2,300
0.90 – 1.08
0.74 – 0.90
2,300 – 2,350
1.44 – 1.71
1.20 – 1.44
2,350 – 2,450
1.50 – 1.79
1.33 – 1.59
2,450 – 2,670
1.05 – 1.26
0 .96 – 1.16
AIC = 32,724
Results of multilevel logistic regression
Household-level fixed effects: odds ratio (95% CI)
Housing density (households/hectare)
2 – 24*
24 – 30
1.23 (1.10 – 1.37)
1.22 (1.10 – 1.37)
30 – 34
1.44 (1.29 – 1.61)
1.44 (1.29 – 1.61)
34 – 39
1.57 (1.39 – 1.76)
1.56 (1.39 – 1.75)
39 – 77
1.74 (1.54 – 1.97)
1.73 (1.53 – 1.96)
Land surface temperature (°C)
26 – 40
1.10 (1.08 – 1.13)
1.10 (1.08 – 1.13)
Elevation (meters above sea level)
2120 – 2260*
2260 – 2300
1.44 (1.18 – 1.76)
1.48 (1.21 – 1.82)
2300 – 2350
1.89 (1.49 – 2.41)
1.99 (1.57 – 2.53)
2350 – 2450
2.12 (1.63 – 2.76)
2.28 (1.75 – 2.98)
2450 – 2670
1.39 (1.07 – 1.82)
1.49 (1.14 – 1.94)
Locality-level fixed effects: odds ratio (95% CI)
1.75 (1.24 – 2.47)
IOR: (0.30 – 10.26)
Locality-level random effects
Variance (95% CI)
1.12 (0.84 – 1.49)
1.03 (0.77 – 1.37)
0.95 (0.71 – 1.27)
Change in variance
Median odds ratio
Likelihood ratio test (p-value)
v. logistic model†
v. multilevel model‡
Locality-level random effects were substantial. The median odds ratio (MOR), which is the median value of the odds ratio when comparing a higher to a lower risk locality, indicated that the median odds of infestation were two and one-half times greater in higher versus lower risk localities (Table 2).
The effect of the locality-level covariate, locality type, was also significant. Households located in shantytowns had 75% higher odds of infestation than households situated in other locality types. However, the interval odds ratio (IOR), which is the interval between the 10th and 80th percentile centered on the median value of the distribution of odds ratios for locality type, included the value one. This indicates that the effect of locality type is not as strong as the locality-level random effect (Table 2).
Household-level effects all remained statistically significant. Odds of infestation increased 10% with each 1 °C increase in land surface temperature, which was similar to the estimate from ordinary multivariate logistic regression. The highest housing density quintile had 75% higher odds of infestation relative to the lowest quintile, which is slightly diminished compared to results from the ordinary multivariate model. The fourth elevation quintile had greater than twice the odds of infestation relative to the lowest quintile, but the highest quintile showed only 50% higher odds of infestation relative to the lowest quintile. The effect of elevation was substantially increased in the multilevel versus the ordinary multivariate logistic regression model (Table 2).
Infestation by T. infestans has been found in many urban areas in Latin America, including Santiago, Chile ; Cochabamba and Sucre, Bolivia ; and Arequipa, Peru . In affected areas—urban as well as rural—prevention and control of Chagas disease relies on vector control . While infestation by and control of T. infestans has been extensively examined in the rural context, infestation in the urban milieu is less well understood. Utilizing spatial and multilevel logistic regression analysis of data collected from multiple sources at multiple spatial scales, we offer insights into the dynamics of T. infestans infestation in an urban landscape.
Prior to implementation of vector control, urban and peri-urban households in Arequipa were extensively infested by T. infestans. The intensity of infestation was spatially heterogeneous, with areas of very low and very high prevalence of infestation. Numerous clusters of infestation, small and large, were found across the six study districts, indicating that urban and peri-urban areas are conducive to the proliferation and dispersion of T. infestans. In rural landscapes, T. infestans have been shown to actively disperse by walking or flying at distances up to approximately 100 or 2,000 meters, respectively [16, 17]. In a separate study in urban Arrequipa, streets were shown to be significant barriers to the dispersion of T. infestans, and to strongly influence the spatial distribution of infestation . In contrast, flight has been observed as a main mechanism of infestation in urban San Juan, Argentina . In the present study, spatial dependence between infested households was observed at distances from 0 to 2,000 meters, suggesting that urban T. infestans may disperse by walking at shorter spans that do not cross city streets, as well as by flying at longer distances across urban blocks.
Identifying extant clusters of infestation prior to vector control may have critical consequences for implementing effective surveillance of vector reemergence subsequent to vector control. In an extensive but sparsely populated rural area in the Gran Chaco of Argentina, reinfestation by T. infestans tended to cluster in areas where infestation was aggregated prior to vector control . Infestation clusters in an extensively and densely populated urban area may be similarly problematic. The existence of numerous infestation clusters in Arequipa, many encompassing large areas and many households, should be priority areas for surveillance and control by the GRSA. Where feasible, utilization of a geographic information system to monitor T. infestans reemergence—as well as other health risks and outcomes—might be a cost effective investment for resource-constrained public health institutions in Arequipa, and elsewhere in the developing world [21, 22].
Spatial heterogeneity in urban infestation by T. infestans is likely influenced by myriad factors operating at multiple spatial scales. We evaluated only a few features, which were chosen based on ecological plausibility and data availability. In ordinary univariate and multivariate logistic regression, housing density, elevation, and land surface temperature were all positively, if not always linearly, associated with household infestation. Housing density may mediate vector dispersal. In higher density urban areas, new habitats and blood sources found in nearby houses are located at short distances from one another, thereby facilitating dispersal of refuge- or blood-seeking vectors. Also, attraction to light influences the dissemination of T. infestans, and the plentiful light sources in higher density urban areas may promote insect dispersal . Land surface temperature may affect vector biology and behavior. Both laboratory and field experiments demonstrate that T. infestans flight initiation increases at higher temperatures [17, 24], thereby promoting vector dispersal in warmer urban areas. Laboratory studies also indicate that higher temperatures increase T. infestans feeding and development rates [25, 26], and blood meal seeking is reportedly the principal cause for dispersion of triatomines . In warmer urban areas, increased feeding and development may result in increased vector dispersal. Elevation may act indirectly through socioeconomic circumstances, rather than directly through biophysical constraints. In Arequipa, lower socioeconomic status populations, often rural-to-urban migrants, typically inhabit the higher elevation hillsides, while higher socioeconomic status populations usually reside in lower elevation valleys . As such, higher infestation at higher elevation in Arequipa may be attributable to two factors: passive introduction of insects resulting from seasonal migration to and from nearby rural areas where T. infestans are prevalent, and substandard living conditions that provide habitats suitable for T. infestans infestation . The slight decrease in infestation at the highest elevations may result from the relatively recent inhabitation of these areas, leaving little time for infestation to have occurred. Elevation is unlikely to be a biophysical constraint for infestation in the currently populated areas of Arequipa, since T. infestans have been found as high as 3,682 meters above sea level in Argentina , well above the elevation of the study area.
Multilevel logistic regression revealed the importance of locality-level contextual effects and substantially diminished spatial autocorrelation present in ordinary logistic regression. The locality-level random effect, which estimates the influence of unobserved contextual effects within each locality, indicates that these unmeasured factors are associated, in median, with substantially higher risk of household infestation. In Arequipa, and elsewhere, urban shantytowns have been identified as areas with higher risk for infestation by Chagas disease vectors, and vector-borne transmission of T. cruzi[11, 13–15]. We offer further evidence that shantytowns are at higher risk for infestation by T. infestans. Controlling for locality-level effects, household-level effects for housing density, elevation, and land surface temperature all remained statistically significant and substantial.
We recognize that our study is limited in many respects. First, while we believe that household location and infestation status data are both precise and accurate, more detailed data regarding the number, life stage, and T. cruzi infection status of insects encountered during the vector control campaign were unavailable. Nor did we have in-depth data regarding households (e.g., construction materials, domestic animals) or their occupants. More detailed data would have likely improved the insights provided by our analyses. Second, we recognize that point-level household covariates are extracted from remote sensing data collected at a 30-meter scale (elevation, land surface temperature) or are spatially smoothed estimates (housing density). We also understand that land surface temperature data do not capture fine-scale temporal variability that occurs across and within days of the year, nor do they describe fine-scale spatial variation in ambient micro-climatic conditions. These issues of scale could conceivably bias the relative magnitude of observed effects. Third, at the time of data collection, portions of the six study districts were still undergoing vector control. Future analyses of areas recently reached by the vector control campaign, including districts beyond the current study area, may provide deeper and broader insights into urban and peri-urban T. infestans infestation.
The geography and ecology of T. infestans—as well as vector species for many other infectious diseases—are changing. Decreasing funding and political will and increasing insecticide resistance are endangering gains made towards interruption of vector-borne transmission of T. cruzi. For many vector-borne diseases in many parts of the world, these are not only public health concerns but also social justice issues, as economically and politically marginalized populations may suffer disproportionately. Many potentially powerful tools (e.g., Google Earth, The R Project for Statistical Computing), data sources (e.g., NASA, NOAA), and spatial and statistical methods are now freely available. Finding novel uses for these resources in conjunction with local knowledge and information—as well as increasing capacity to do so—could inspire new perspectives on and solutions to existing and emerging public health problems and their social and environmental causes and consequences.
In 2003, the Gerencia Regional de Salud de Arequipa (GRSA) initiated a vector control program to eliminate household infestation by T. infestans. The program is ongoing, is implemented on a district-by-district basis, and consists of a simple stepwise process. First, each household is assigned a unique alphanumeric code, and household codes and locations are catalogued on hand-drawn maps. Second, trained GRSA personnel go door-to-door, spray all domestic areas and peri-domestic animal enclosures in each household, and inspect for the presence of triatomines, taking advantage of the flushing effect of the insecticide. The presence or absence of T. infestans is recorded. Finally, the inspection and insecticide application process is repeated approximately six months later.
Maps of household locations, household unique identifier codes, and dates and results of household T. infestans inspections were provided by the GRSA. Within the study area 37,229 households had been inspected for infestation by T. infestans and sprayed with insecticide during the period from September 2003 through December 2008. Using GRSA maps and Google Earth imagery, we assigned geographic coordinates to these households, as well as to households within the six study districts that had not as yet been reached by the vector control campaign as of December 2008. Household geographic coordinates, household T. infestans infestation status (0/1), and household unique identifier codes were stored in a relational database management system for subsequent analysis.
Remote sensing data
Advanced Spaceborne Thermal Admission and Reflection Radiometer (ASTER) Global Digital Elevation Model Version 2 (GDEM V2) 30-meter-resolution imagery of the study area (ASTGTM2_S17W072) was acquired from the National Aeronautics and Space Administration (NASA) Earth Observing System Data and Information System (EOSDIS) . Landsat 5 Thematic Mapper (TM) 120-meter resolution thermal imagery (band 6: 10.40 – 12.50 μm) and 30-meter resolution visible (band 3: 0.63 – 0.69 μm) and near-infrared (band 4: 0.76 – 0.90 μm) imagery of the study area (WRS path 3 row 71) were obtained from the United States Geological Survey (USGS) EarthExplorer . Cloud-free images were attained for nine dates in 2008: 18 March, 19 April, 21 May, 24 July, 25 August, 26 September, 12 October, 13 November, and 15 December. A Landsat 5 TM band 3 image acquired on 25 July 1987 was obtained from NASA’s Global Orthorectifed Landsat Data Set as reference for geometric correction of 2008 images [31, 32].
In Peru, census areas are subdivided into departments, provinces, districts, and localities. A database of the locality in which each study household was located was provided by the GRSA, and a census classification of locality types was obtained from the Peru National Institute of Statistics and Informatics (INEI). Locality type categories included city (ciudad), housing development (urbanización), town (pueblo), shantytown (pueblo joven), housing association (asociación de viviendas), housing cooperative (cooperativa de viviendas), annex (anexo), hamlet (caserío), and rural community (comunidad campesina) . The GRSA and INEI databases were joined and each household was assigned a categorical variable specifying locality type for 34,275 of the 37,229 (92.1%) households with documented T. infestans inspection data.
Spatial point pattern analysis
Spatial variation in household infestation by T. infestans was evaluated by dividing the kernel-smoothed density of infested households by the kernel-smoothed density of all study households. An isotropic Gaussian smoothing kernel with a standard deviation (σ) of 45.7 meters was utilized for this analysis, where σ was selected using a likelihood cross-validation method. A 999-iteration random labeling simulation was performed to identify areas where infestation by T. infestans was significantly elevated [34, 35]. The kernel-smoothed density of all georeferenced households (n = 68,849), with σ = 27.8 meters, was estimated for evaluation as a covariate in logistic regression [35, 36].
Spatial dependence between households infested by T. infestans was assessed by computing the difference in the K-function for infested households and non-infested households at 100-meter increments from 0 to 2000 meters. A 999-iteration random labeling simulation was executed to identify distances at which spatial dependence of infested households was statistically significant .
Spatial statistical analyses were conducted using R (The R Project for Statistical Computing) . Maps of spatial variation in household infestation and household density were constructed using ArcGIS version 10 (ESRI) .
Remote sensing image analysis
Remote sensing data were utilized to derive and extract household point-level estimates of elevation and land surface temperature for evaluation as covariates in logistic regression.
The ASTER GDEM V2 image was projected to the WGS 84 UTM 19S coordinate system and resampled to a 30-meter pixel size to match the projection and spatial resolution of Landsat 5 TM imagery, and elevation data were extracted to household point locations.
Landsat 5 TM and ASTER GDEM V2 images were cropped to a 701-column by 474-row area corresponding to the rectangle bounding the six study area districts. Landsat 5 TM band 3, 4, and 6 images from 2008 were geometrically corrected utilizing a simple root mean square error minimization routine and the Landsat 5 TM band 3 image from 1985 as reference . Atmospheric correction of Landsat 5 TM band 3 and 4 images was performed using a modified dark object subtraction method , followed by topographic correction utilizing a Minnaert method . The normalized difference vegetation index (NDVI) was calculated from Landsat 5 TM band 3 and 4 images , and land surface emissivity was estimated from NDVI for each date in 2008 . Land surface temperature was derived using Landsat 5 TM band 6 thermal infrared images; land surface emissivity images; and coefficients for atmospheric transmissivity, upwelling atmospheric radiance, and downwelling atmospheric radiance [45, 46]. Land surface temperature images were overlaid, median land surface temperature was calculated for each raster pixel, and these data were extracted to household point locations.
Processing of ASTER GDEM V2 and Landsat 5 TM images was conducted using the landsat package version 1.0.8 in R version 2.15.2 . Atmospheric coefficients were obtained from the Atmospheric Correction Parameter Calculator . ArcGIS version 10 was utilized to create maps of elevation and land surface temperature, and to extract values for these variables to household point locations.
To address spatial dependence among observations and to assess locality-level effects, three multilevel logistic regression models were evaluated: a model including only locality-level random effects (Model 0); a model including household-level covariates and locality-level random effects (Model 1); and a model including household-level covariates, a locality-level covariate, and locality-level random effects (Model 2). In addition to standard odds ratios, median odds ratios (MOR) were calculated for locality-level random effects, and the interval odds ratio (IOR) were calculated for the locality-level covariate [48, 49].
Categorization scheme for locality type
Logit is the link function; Pr(y ij = 1) is the probability of household infestation; i and j indicate the ith household and the jth locality, β 1 is the vector of regression coefficients for density, where density is kernel-smoothed housing density (households/hectare) categorized by quintiles; β 2 is the vector of regression coefficients for elevation, where elevation (meters above sea level) is categorized by quintiles; β 3 is the regression coefficient for temperature, where temperature is estimated annual median land surface temperature (°C); β 4 is the regression coefficient for locality type, where locality type is a dichotomous categorization of localities into shantytown and other; β 0 is the household-level intercept; and u j is the locality-level random effect.
Regression model goodness of fit was assessed using the likelihood ratio test and the Akaike information criteria (AIC). Spatial autocorrelation in deviance residuals of the least- and best-fit models was evaluated utilizing the mark correlation function (Km) at 25-meter increments from 0 to 500 meters. A 999-iteration random labeling simulation was executed to identify distances at which regression residual spatial autocorrelation was statistically significant .
Stata/IC 12.1 was utilized for statistical analyses , and R version 2.15.2 was used for spatial analysis of regression residuals.
We gratefully acknowledge the invaluable contributions of the following institutions and individuals. The Gerencia Regional de Salud de Arequipa (GRSA) conducted vector surveillance and control and shared these data with us, and we are especially grateful to Juan G. Cornejo del Carpio and Fernando S. Málaga Chavez. The Chagas Disease Working Group in Arequipa, Peru conducted household georeferencing and assisted with database management. Dr. Denise Roe, Professor of Biostatistics at The University of Arizona, provided expert statistical consultation. The Ministerio de Salud del Perú (MINSA), Dirección General de Salud de las Personas (DGSP), Estrategia Sanitaria Nacional de Prevención y Control de Enfermedades Metaxénicas y Otras Transmitidas por Vectores (ESNPCEMOTVS), Dirección General de Salud Ambiental (DIGESA); Gobierno Regional de Arequipa; Pan American Health Organization (PAHO); and Canadian International Development Agency (CIDA) provided financial, logistical, and technical support. We wish to thank the anonymous reviewers for their insightful and invaluable comments, questions, and suggestions. ASTER GDEM V2 is a product of METI and NASA. This work was supported by the National Institutes of Health (http://nih.gov/; 5P50 AI074285-03 and 04, 3K01AI079162-02S1 to MZL, 5 T37 MD001427 to SD) and the National Institute of Allergy and Infectious Diseases (http://www.niaid.nih.gov/; 1K01AI079162-03 and R01AI101229-01A1 to MZL).
- Heymann DL (Ed): Control of communicable diseases manual. 19th edition. Washington, DC: American Public Health Association; 2008.
- World Health Organization: Control of Chagas disease: Second report of the WHO expert committee. Geneva: World Health Organization; 2002.
- World Health Organization: The global burden of disease: 2004 update. Geneva: World Health Organization; 2008.
- World Health Organization: Chagas disease: control and elimination (Report by the Secretariat). Geneva: World Health Organization; 2010.
- World Health Organization: Accelerating work to overcome the global impact of neglected tropical diseases – a roadmap for implementation. Geneva: World Health Organization; 2012.
- Pinto Dias JC: Southern Cone Initiative for the elimination of domestic populations of Triatoma infestans and the interruption of transfusional Chagas disease. Historical aspects, present situation, and perspectives. Mem Inst Oswaldo Cruz 2007, 102:11–18.View Article
- World Health Organization: Reporte del grupo de trabajo científico sobre la enfermedad de Chagas. Geneva: World Health Organization; 2007.
- Moncayo A, Silveira AC: Current epidemiologic trends for Chagas disease in Latin America and future challenges in epidemiology, surveillance and health policy. Mem Inst Oswaldo Cruz 2009, 104:17–30.PubMedView Article
- Tarleton RL, Reithinger R, Urbina JA, Kitron U, Gürtler RE: The challenges of Chagas disease – grim outlook of glimmer of hope? PLoS Med 2007, 4:e332.PubMedView Article
- Mott KE, Desjeux P, Moncayo A, Ranque P, de Raadt P: Parasitic diseases and urban development. Bull World Health Organ 1990, 68:691–696.PubMed
- Levy MZ, Bowman NM, Kawai V, Waller LA, del Carpio JG C, Cordova Benzaquen E, Gilman RH, Bern C: Periurban Trypanosoma cruzi -infected Triatoma infestans , Arequipa, Peru. Emerg Infect Dis 2006, 12:1345–1352.PubMedView Article
- Levy MZ, Small DS, Vilhena DA, Bowman NM, Kawai V, Cornejo Del Carpio JG, Cordova-Benzaquen E, Gilman RH, Caryn Bern C, Plotkin JB: Retracing micro-epidemics of Chagas disease using epicenter regression. PLoS Comput Biol 2011, 7:e1002146.PubMedView Article
- Schenone H: Enfermedad de Chagas en Chile: sectores rurales y periurbanos del area de endemo-enzootia relaciones entre condiciones de la vivienda, infestación triatomidea domiciliaria e infección por Trypanosoma cruzi del vector, del humano y de mamíferos domésticos, 1982–1985. Bol Chil Parasitol 1985, 40:58–67.PubMed
- Tibayrenc M: La maladie de Chagas en Bolivie: données préliminaires sur les cycles domestiques: méthodes simplifiées de captures des triatomes. Cah Orstom (Série Entomologie Médicale et Parasitologie) 1984, 22:51–53.
- Bowman NM, Kawai V, Levy MZ, Cornejo Del Carpio JG, Cabrera L, Delgado F, Malaga F, Cordova Benzaquen E, Pinedo VV, Steurer F, Seitz AE, Gilman RH, Bern C: Chagas disease transmission in periurban communities of Arequipa, Peru. Clin Infect Dis 2008, 46:1822–1828.PubMedView Article
- Vazquez-Prokopec GM, Ceballos LA, Kitron U, Gürtler RE: Active dispersal of natural populations of Triatoma infestans (Hemiptera: Reduviidae) in rural northwestern Argentina. J Med Entomol 2004, 41:614–621.PubMedView Article
- Schofield CJ, Lehane MJ, McEwen PK, Catalá S, Gorla DE: Dispersive flight by Triatoma infestans under natural climatic conditions in Argentina. Med Vet Entomol 1992, 6:51–56.PubMedView Article
- Barbu CM, Hong A, Manne JM, Small DS, Quintanilla Calderón JE, Sethuraman K, Quispe-Machaca V, Ancca-Juárez J, del Carpio JG C, Málaga FS, Náquira C, Levy MZ: The effects of city streets on an urban disease vector. PLOS Computational Biology 2013, 9:e1002801.PubMedView Article
- Vallvé SL, Rojo H, Wisnivesky-Colli C: Urban ecology of Triatoma infestans in San Juan, Argentina. Mem Inst Oswaldo Cruz 1996, 91:405–408.PubMedView Article
- Porcasi X, Catalá SS, Hrellac H, Scavuzzo MC, Gorla DE: Infestation of rural houses by Triatoma infestans (Hemiptera: Reduviidae) in southern area of Gran Chaco in Argentina. J Med Entomol 2006, 43:1060–1067.PubMedView Article
- Mott KE, Nuttall I, Desjeux P, Cattland P: New geographical approaches to control of some parasitic zoonoses. Bull World Health Organ 1995, 73:247–257.PubMed
- Stansfield SK, Walsh J, Prata N, Evans T: Information to improve decision making for health. In Disease Control Priorities in Developing Countries. Edited by: Jamison DT, Breman JG, Measham AR, Alleyne G, Claeson M, Evans DB, Jha P, Mills A, Musgrove P. Washington, DC: Oxford University Press and The World Bank; 2006:1017–1013.
- Zeledón R, Rabinovich JE: Chagas disease: an ecological appraisal with special emphasis on its insect vectors. Ann Rev Entomol 1981, 26:101–133.View Article
- Lehane MJ, McEwen PK, Whitaker CJ, Schofield CJ: The role of temperature and nutritional status in flight initiation by Triatoma infestans . Acta Trop 1992, 52:27–38.PubMedView Article
- Catalá S: The biting rate of Triatoma infestans in Argentina. Med Vet Entomol 1991, 5:325–333.PubMedView Article
- Gorla DE, Schofield CJ: Population dynamics of Triatoma infestans under natural climatic conditions in the Argentine Chaco. Med Vet Entomol 1989, 3:179–194.PubMedView Article
- Bayer AM, Hunter GC, Gilman RH, Cornejo del Carpio JG, Naquira C, Bern C, Levy MZ: Chagas disease, migration and community settlement patterns in Arequipa, Peru. PLoS Negl Trop Dis 2009, 3:e567.PubMedView Article
- Instituto Nacional de Estadística e Informática del Perú: XI Censo de Población y VI de Vivienda. 2007. http://www.inei.gob.pe
- National Aeronautics and Space Administration Earth Observing System Data and Information System. http://reverb.echo.nasa.gov
- United States Geological Survey EarthExplorer. http://earthexplorer.usgs.gov
- Global Orthorectifed Landsat Data Set. http://www.landsat.org/ortho/
- Tucker CJ, Grant DM, Dykstra JD: NASA’s global orthorectified Landsat data set. Photogramm Eng Remote Sensing 2004, 70:313–322.
- Instituto Nacional de Estadística e Informática del Perú: Estadística de Centros Poblados. 1993. http://www.inei.gob.pe
- Kelsall JE, Diggle PJ: Non-parametric estimation of spatial variation in relative risk. Stat Med 1995, 14:2335–2342.PubMedView Article
- Diggle PJ: Statistical analysis of spatial point patterns (Second edition). London: Arnold; 2003.
- Diggle PJ: A kernel method for smoothing point process data. Appl Stat 1985, 34:138–147.View Article
- Diggle PJ, Chetwynd AG: Second-order analysis of spatial clustering for inhomogeneous populations. Biometrics 1991, 47:1155–1163.PubMedView Article
- The R Project for Statistical Computing. http://www.r-project.org
- ESRI. http://www.esri.com
- Goslee SC: Analyzing remote sensing data in R: the landsat package. J Stat Softw 2011, 43:1–25.
- Song C, Woodcock CE, Seto KC, Lenney MP, Macomber SA: Classification and change detection using Landsat TM data: when and how to correct atmospheric effects? Remote Sens Environ 2000, 75:230–244.View Article
- Riaño D, Chuvieco E, Salas J, Aguado I: Assessment of different topographic corrections in Landsat-TM data for mapping vegetation types. IEEE Trans Geosci Remote Sens 2003, 41:1056–1061.View Article
- Jensen JR: Introductory digital image processing: a remote sensing perspective (Third edition) Upper Saddle River. NJ: Pearson Education, Inc; 2005.
- Jiménez-Muñoz JC, Cristóbal J, Sobrino JA, Sòria G, Ninyerola M, Pons X: Revision of the single-channel algorithm for land surface temperature retrieval from Landsat thermal-infrared data. IEEE Trans Geosci Remote Sens 2009, 47:339–349.View Article
- Barsi JA, Schott JR, Palluconi FD, Hood SJ: Validation of a Web-based atmospheric correction tool for single thermal band images. Proc SPIE 2005, 5882:136–142.
- Chander G, Markham BL, Helder DL: Summary of current radiometric calibration coefficients for Landsat MSS, TM, ETM+, and EO-1 ALI sensors. Remote Sens Environ 2009, 113:893–903.View Article
- Atmospheric correction parameter calculator http://atmcorr.gsfc.nasa.gov
- Larsen K, Merlo J: Appropriate assessment of neighborhood effects on individual health: integrating random and fixed effects in multilevel logistic regression. Am J Epidemiol 2005, 161:81–88.PubMedView Article
- Merlo J, Chaix B, Ohlsson H, Beckham A, Johnell K, Hjerpe P, Råstam L, Larsen K: A brief conceptual tutorial of multilevel analysis in social epidemiology: using measures of clustering in multilevel logistic regression to investigate contextual phenomena. J Epidemiol Community Health 2006, 60:290–297.PubMedView Article
- Goreaud F: Apports de l’analyse de la structure spatiale en foret tempérée à l’étude et la modélisation des peuplements complexes. PhD thesis. Paris: Ecole Nationale du Genie Rural, des Eaux et des Forêts; 2000.
- Stata http://www.stata.com/
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/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.