SciELO - Scientific Electronic Library Online

 
vol.20 issue2Determination of simple and complex water level - discharge relationship in a river. Case of the Cauca RiverTheoretical understanding between competitiveness and productivity: firm level author indexsubject indexarticles search
Home Pagealphabetic serial listing  

Services on Demand

Journal

Article

Indicators

Related links

  • On index processCited by Google
  • Have no similar articlesSimilars in SciELO
  • On index processSimilars in Google

Share


Ingeniería y competitividad

Print version ISSN 0123-3033

Ing. compet. vol.20 no.2 Cali July/Dec. 2018

https://doi.org/10.25100/iyc.v20i1.5810 

Artículos

Analysis of spatial distribution of malaria in the department of Chocó for the year 2016

Análisis de la distribución espacial de la malaria en el departamento de chocó para el año 2016

Leonardo Hurtado1 

Manuel Aguilar2 

Miguel Ávila2 

1 Universidad Nacional de Colombia, Bogotá - Colombia. leonardohurtado1@gmail.com

2 Universidad Distrital Francisco José de Caldas, Bogotá - Colombia. josemanuel.aguilarruiz@gmail.com; maavila@udistrital.edu.co


Abstract

The following investigation aims at the estimate of a generalized linear spatial regression model of Poisson statistical distribution that allows to explain the spatial-geographic behavior of a disease transmitted by vectors such as malaria for the specific case of the municipalities of the department of Chocó in the year 2016.

In order to obtain the probability maps, spatial statistics techniques are used based on the analysis of Lattice type areas or data such as figure criteria, distance criteria and the most commonly used criteria of physical contiguity, in order to determine the existence of spatial autocorrelation To analyze the existence of malarial spatial autocorrelation among nearby municipalities, it was carried out through the statistical dynamics of sanitary variables (unsatisfied basic needs and aqueduct coverage), environmental variables (forest cover, precipitation, humidity, height and temperature) and demographic variables such as population of each municipality.

Nevertheless, it is necessary to perform a Gaussian anamorphosis process to guarantee the standardization of the data in standardized morbidity rates (SMR) and in this way identify the statistical structure of the data that for this year came from a Poisson distribution. The study allowed finding positive spatial autocorrelation confirming that in some municipalities of Chocó as Paimado, San Pablo canton, Istmina, Tado and La union panamerica malaria is more likely to occur because of its municipal closeness according to the physical relationship of certain study variables.

The result was that the variables that best explain malaria in the municipalities of Chocó are the coverage of forests and unsatisfied basic needs that give as a product the map of spatial auto correlation and the map of probability of occurrence with the highest ratio among municipalities of the region. South of the department.

Keywords: generalized linear spatial model regression; malaria; SMR; autocorrelation

Resumen

La presente investigación tiene como propósito la estimación de un modelo de regresión espacial lineal generalizada de distribución estadística Poisson que permita explicar el comportamiento espacio-geográfico de una enfermedad transmitida por vectores como la malaria para el caso puntual de los municipios del departamento de Chocó para el año 2016.

Con el fin de obtener los mapas de probabilidad se emplean técnicas de estadística espacial basada en el análisis de áreas o datos tipo Lattice tales como criterios de figuras, criterios de distancias y los más usados los criterios de contigüidad física, con el fin de determinar la existencia de autocorrelación espacial de dicho. Para analizar la existencia de autocorrelación espacial malárica entre municipios cercanos se realizó mediante la dinámica estadística de variables sanitarias (necesidades básicas insatisfechas y cobertura de acueducto), variables ambientales (cobertura boscosa, precipitación, humedad, altura y temperatura) y variables demográficas como la población de cada municipio.

Sin embargo, es necesario realizar un proceso de anamorfosis gaussiana para garantizar la normalización de los datos en tasas de morbilidad estandarizada (SMR) y de esta forma identificar la estructura estadística de los datos que para este año fueron provenientes de una distribución Poisson. El estudio permitió encontrar autocorrelación espacial positiva lo que confirma que en algunos municipios del Chocó como Paimado, El cantón de San Pablo, Istmina, Tado y La unión panamerica la malaria tiene mayor probabilidad de ocurrencia por su cercanía municipal de acuerdo a la relación física de ciertas variables del estudio.

Se obtuvo como resultado que las variables que mejor explican espacialmente la malaria para los municipios del Chocó son la cobertura de bosques y necesidades básicas insatisfechas que dan como producto el mapa de autocorrelación espacial y el mapa de probabilidad de ocurrencia con mayor relación entre los municipios del sur del departamento.

Palabras clave: Modelo de Regresión espacial lineal generalizada; malaria; SMR; autocorrelación

1. Introduction

In recent years there has been an increase in mortality rates due to the presence of malaria along the national territory. It should be highlighted that malaria is a disease caused by the protozoa of the plasmodium genus, transmitted by vectors1. According to the Ministry of Health there are four types of this parasite (P. falciparum, P. vivax, P. malariae and P. ovale) that affect approximately 12 million people 1.

The National Institute of Health (INS) in its protocol of public health security, ensures that about 85% of the national rural territory is below 1,600 meters above sea level and has the ideal geographic, climatic and epidemiological conditions for the transmission of this disease 2. According to statistics in the last decade, there is an annual average of 130 to 150 deaths, with the Pacific region of the country being the most significant 3.

Deficiency in control programs, socioeconomic conditions of the population with a high degree of vulnerability, and problems with unsatisfied basic needs affect the gradual increase of malaria transmission 4. According to the INS in its 2014 epidemiological bulletin, it reveals the quantitative and percentage distribution of cases registered throughout the national territory, in which the department of Chocó registers more than 11,000 cases which represent 41.41% of the total cases, turning into a pillar of the spread of the disease 5. The data required for this analysis were collected through the Sivigila 2.

As previously mentioned, there are several studies that involve these variables and their existing spatial relationship. Among them, the figures provided by the DANE, INS newsletters and specifically an article from the biomedical journal on malarial epidemiology in Colombia 6. Based on these studies and the data collected, it was considered of great importance to perform a spatial statistics analysis that would prove or reject the hypothesis that there is a spatial relationship of a disease such as malaria taking into account the physical and healthy conditions as variables and that validate that the closer a study area is in relation to the infection zone, the greater the probability of infection.

this way, the main objective of the research was to implement a spatial statistical model that best explained the dependence between the observations of Malaria in 2016 and the explanatory variables that possibly have an impact on the occurrence of the disease, as well as the generation of maps that allow to stratify the zones according to the number of occurrence and the risk that the spreading of the infection in the department of Chocó presents.

2. Metodology

2.1. Study zone

The department of Chocó is located on the western side of Colombia, in the Pacific region, located between the latitudes north 4°00'50" and 8°42'32" and the longitudes west 76°2'57" and 77°53'38" (Figure 1). The departmental territory is within the area of equatorial climates, characterized by high rainfall, the temperature of its valleys and coastal lowlands, is higher than 27° C, usually accompanied by high relative humidity 90% 7. In the department there are natural national parks (Utría, Katíos and Tatamá). The majority of the population is black and originated from the African slaves brought during colonial times, for mining exploitation. In the Department there are 82 safeguards, 6 of which are shared with the department of Valle del Cauca 8.

Figure 1 Departament of Chocó. Geographic Localization in Colombia. 

The territory is made up of the basins of the Atrato, San Juan and Baudó rivers. In this geographical framework, the following physiographic units are distinguished: the coastal strip, divided by cabo Corrientes; considered the most important accident of the Pacific coast. The hydrographic system is one of the most abundant and interesting in the country 9; In addition to the Atrato, San Juan and Baudó Rivers, the Andágueda, Bebará, Bebaramá, Bojayá, Docampadó, Domingodó, Munguidó and Opogodó Rivers are important 10. Chocó has 30 municipalities, 147 corregimientos and 135 police inspections. The most important cities are in their order Quibdó, Istmina, Condoto, Acandí and Bahía Solano 11.

The limits of the department are: by the north with Panama and the Caribbean Sea; on the east with the departments of Antioquia, Risaralda and Valle del Cauca; and to the west with the Pacific Ocean. (Figure 1). It has an area of 46,530 km2, a population of 3'657,821 inhabitants and its capital with a population of 500,093 inhabitants 12 for the year 2005 according to the DANE.

2.2 Data

According to the INS there are various types of physical, economic and social conditions that help this disease to spread and affect more areas than others 5. Because of this, some variables that can provide significant information in the spatial modeling of the disease are included. First, it is considered that the height above sea level allows the survival of the species of the transmitting mosquito 4 since at higher altitudes the probability of the carrier's existence is lower. The height is usually in hand with another physical variable such as temperature which in the first instance that could play a fundamental role in the reproduction of the larva of the mosquito giving the ideal conditions for its development as it has been shown in some African countries 13. Variables such as precipitation and relative humidity (the latter are average annual values) as well as the area covered by forest create an ideal environment for the incubation, proliferation and reproduction of the insects that transmit malaria 14.

Moreover, there are social and economic factors. These factors, according to their management, can reduce or accelerate the transmission of the disease, generating appropriate habitats for the transmitting mosquito. The incidence of sewage and aqueduct coverage can be highlighted since they directly affect the quality of the population's water resources, the level of Unsatisfied Basic Needs (%), which highlights that 79.19% of the total population of the population department is in precarious conditions 15.

For the elaboration of the database, information obtained from the land use planning diagram (EOT) of each of the municipalities of the department of Chocó was taken into account since in them one of the first components within the formulation study is the description of the basic variables of the municipality. Among them and the most important we have the POT of Quibdó 16, the municipality of Tadó 17, the municipality of Lloró 18, the municipality of Riosucio 19 and the municipality of Istmina 20. The data for the other municipalities were taken by the EOT of each of them.

Figure 2 Exploratory data analysis. 

Figure 2 shows the boxplot diagram in which the distribution of malaria cases in the department is observed. There is evidence of the existence of very high outliers (values above 1,100) at the upper end of the box, these values are found in municipalities south of the department while lower values are in the north. It is observed as a value within the interquartile range corresponds to the number of deaths in the municipality of Riosucio (see arrow) which shows the relationship of the graph and the municipalities in the department.

Table 1 shows the nomenclature that will be applied in the future for each of the study variables that will be taken into account in the next analyzes.

Table 1 Nomenclature of used variables. 

Variable Nomenclatura
Malaria 2016 M2016
Population 2016 Pob
Precipitation P
Wetness h
Temperature T
Forest percentage B
Toilet and sewer coverage CAA
Unsatisfied basic needs NBI

2.3 Preparation of Malaria Maps

The calculation of the standardized morbidity rate (SMR) was made, which is obtained from the quotient between the number of cases observed for each of the municipalities in the department of Chocó and the number of expected cases, the latter was obtained multiplying the population at risk presented in each municipality and the total incidence rate, which, in turn, is the quotient between the total number of cases observed and the total population present in the 26 selected municipalities of the department. This rate allowed to detect the municipalities where there was a greater or lesser number of observed and expected cases.

It is possible to represent the study of a disease through disease probability maps, spatial autocorrelation tests, case aggregations and spatial regression models 21. As for the probability maps, the function was used, which yields the probability of finding a high or low observed value compared to a comparison with the expected values, for this case of malaria for each municipality. The test used to determine the existence of self spatial correlation is the Moran Index test (equation 1), in order to identify the degree of spatial aggregation for each of the observations the LISA statistic was evaluated (Local Indicators for Spatial Autocorrelation).

I=(n/S0)i=1nj=1nwijzizj / i=1nzi2 (Equation 1)

Where, the sum of elements of the weight matrix. The observations Z are the deviations from the mean, Xi is the value of the variable in a given spatial unit and Xj is the value of the variable in another location, usually the neighbors to Xi 22.

Knowing that linear models are not the best approximation for modeling events related to the counting of cases of a certain disease, a Poisson distribution was assumed, which is part of the generalized linear models (GLMs for its acronym in English), this being one of the most appropriate for the type of variable that will be used 23 and in this way assign statistical modifications to the values so that the variables explain most of the model and statistically significant variables are obtained at 0.05 and 0.1 of level of significance. Figure 3 shows the generic formulas of the probability distributions that are part of the GLM 24.

Figure 3 Generic formulation of generalized linear models.  

3. Analysis and Discussion of Results

To select the best spatial weights matrix, self-correlation tests and map analysis, the SMR for malaria was taken into account in 2016, as well as the explanatory variables that seek to model the endemic phenomenon to environmental and socioeconomic conditions. The statistical analysis of each variable separately as well as its spatial interaction allow the selection of the most convenient spatial criteria and in this way observe with the Moran index the possible existence of self-spatial correlation that implies the abstraction of a regression model that generate the variables that best explain the malaria rate for each municipality of the department of Chocó.

3.1 AIC Criterion: Selection of the spatial weights matrix

To make the selection of the best spatial matrix, the Akaike criterion (AIC Akaike's Information Criterion) was used for each of the types of spatial weights neighborhood matrices, which can be binary or standardized, all the criteria that were evaluated are presented in Table 2, along with each of the results of AIC.

From the criteria based on graphs, the sphere of influence was selected because it had the lowest value of AIC and because it takes into account the intersection of the regions in at least two places (Figure 4), while the criteria of 1 more neighbors close and the tower contiguity criterion give the best forecasts for the generation of the regression and analysis models of the Moran scaterplot as well as the same indicator.

Table 2 AIC for weight matrices. 

MATRICES DE VECINDADES
Contigüidad Física Criterios basados en gráficas Criterios basados en distancia
Torre Reina Triangulación Delaunay Esfera de Influencia Gráfica de Gabriel Vecinos Relativos Vecino más cercano
1 2 3 4
AIC 33,274 26,575 30,973 18,397 31,579 30,864 32,67 23,31 28,02 32,04

Figure 4 Criteria to define spatial interactions. Based on graphs, distances and spatial contiguity. 

General Histogram

Figure 5 Standardized morbidity rate, comparison between the histogram and the map for values a) minor and b) higher than the rate. Histomap function 3

Figure 5a shows the set of municipalities with the lowest standardized morbidity rate, while in Figure 4b the municipality of Bagadó has a high rate, given that the number of cases observed is greater than the number of cases. Cases that could be expected for this municipality. This analysis was applied in this stage only to the variable SMR because the explanatory variables have a role in the modeling, for this point the exploratory analysis of the malaria data seeks to give a hypothetical picture of how this phenomenon explained by itself in the municipalities of the department of Chocó.

3.2 Spatial Autocorrelation

To verify the initial hypothesis of the presence or absence of spatial dependence, the significant association of similar or different values among the neighboring municipalities of the Department of Chocó was examined through the I. de Moran tests. The global tests of the I. de Moran for the SMR yielded a value of 0.4835 generated from 999 simulations, as its value is positive it is spoken of a positive self-correlation, not very strong, but it reveals the existence of the phenomenon through of the municipalities as study areas.

Map of spatial correlation auto

Taking into account the statistical results of certain variables such as the one observed in Figure 5, it is evident that there is a certain degree of spatial relationship between the incubation of this disease and the fatal cases that occur. Taking this into account, an analysis of the global Moran index is made with which the relationship of the areas in terms of malaria and the dispersion of the points between the analysis of pair of variables can be identified.

The maps generated from the local Moran index can be visualized directly in the units of measurement that for this case are the municipalities (see Figure 6) and scatter diagrams that will be explained later in figure 8. The map of the Moran index it generates the graphic representation of the dispersion of the explanatory variables vs the response variable and they are grouped in quadrants (from High to Low according to the axis of the quadrant and their combinations). High values of the phenomenon and the variables that seek to explain it indicate a positive response that indicates that there is a relationship between the municipalities and very low values indicate that the explanatory variable is not the response of malaria in a given municipality.

Figure 6 Map of Moran. Auto spatial correlation. 

In this way the municipalities with low-low values (LL) of Figure 6 are Bahía Solano, Unguia, Acandí and Riosucio. Similarly, the municipalities with high-high values (HH) of Figure 6 are Istmina, Pan-American Union, Tadó, San Pablo Canton and Paimado. This map of the malarial relationship allows us to observe how the high values of malaria are concentrated in nearby neighbors of the south of the department while the municipalities to the north and border with Panama and the ocean have a low relation of hatching. The local index shows how areas of low malarious rate are related to other municipalities, while southern municipalities share environmental variables that generate favorable circumstances for the development of the disease.

The Moran index is based on the tower criterion that best explained the spatial autocorrelation of malaria among the study municipalities. However, the spatial autocorrelation map is not the only result that yields this statistical index. A very valuable map to identify the spatial distribution of the disease is the probability map presented in Figure 7, there some municipalities of the department show a significant response to the probability of having high rates of deaths from malaria being evaluated with the level of significance less than 0.05, 0.01, 0.001 and 0.0001 respectively.

Figure 7 Moran probability map. 

The probabilistic measure indicates that municipalities such as Acandí, Riosucio, Tado and Istmina register a degree of spatial autocorrelation at levels of significance less than 0.01 and the municipality of Unguia less than 0.001. Compared to the spatial autocorrelation map of the local Moran, the municipalities that establish a spatial relationship between the variables are quite similar, and a high probability is due to the forest cover that makes other terrain variables adapt or remove the transmitting insects.

Moran Scaterplot

There are several ways to calculate the Moran index, visualize it and even interpret it, in this case the scaterplot is presented, this graph shows the position of the representative points of each municipality in a Cartesian plane where each quadrant represents (High-high, high- low, low-low and low-high) respectively (Figure 8).

Figure 8 Dispersion diagram of the Moran index. 

A high concentration of points is observed in the HH quadrant, which is equivalent to the map in Figure 6 where some municipalities are observed to the south of the department with this behavior, the distribution of these points is observed mainly in the opposite quadrants and yields a Moran index. Of 0.48 confirming the relationship of malaria among the municipalities. These municipalities in quadrant one could explain how the area of highest concentration of cases can explain how hot spots and otherwise the quadrant three that explains the northern municipalities where the cold spots are and indicates that the spatial relationship between them it exists only that it minimizes the probability of occurrence of malaria cases.

3.3 Poisson Regression Model

Once the existence of spatial autocorrelation is identified from the scatterplot and the autocorrelation map, it is proposed to generate a regression model that has the result of identifying the variables that explain malaria among the municipalities of the department of Chocó. To identify the candidate variables, individual data analyzes are performed that seek to observe the statistical behavior of the variable against the malaria rate.

Figura 9 Analysis of forest data. 

The individual analysis of the forest cover variable shown in Figure 9 shows a tendency of the data to a function of standard deviation with respect to the mean, in this way the intercept would adjust the data. Starting from this statistical concept, the standard deviation would be applied to the forest cover variable that will be entered into the regression model.

Initially, the spatial delay model shows a square R of 0.5694 with the unmodified variable, however, applying a standard deviation function, the R squared increases to 0.6305, which improves the analysis of the response variables of the model. The Poisson regression model was defined from the variables that were considered to be explanatory for the dependent variable, as previously mentioned. Next, the formula used for the model is presented:

M2016=logPob + β0+β1h+β2P+desvβ3B+β4NBI+β5CAA (2)

The regression models were selected according to some variables in the evaluation of the levels of significance, but especially taking into account the information explained by the model. In this way, Table 3 shows the results of the R squared for each of the models taking into account the totality of proposed variables to explain malaria in the department of Chocó.

Table 3 Results regression models. 

Model Initial R-squere Modify R-squere
Clasic 0.5509 0.63
Spatial lag 0.5305 0.53
Residuals 0.5759 0.58

It was found that the statistically significant variables are the unsatisfied basic needs and forest cover with a modification to the original probability of the data, with this modification the model adjusted from 0.55 to 0.63 according to the R-square. Based on these variables entered into the model, those already mentioned are significant, since they are less than α = 0.05, while the other variables, although affecting the model in statistical and spatial terms, are necessary to convert before being eliminated from modeling (Table 4).

Table 4 Summary of the spatial regression model 

Variable Coefficient St. Error Z-Value Probability
Constant -1653.91 1015.53 -1.6286 0.1033
Log_Pob -135.43 239.73 -0.5649 0.5721
High 0.169 0.1897 0.8908 0.373
Precipitation 0.01946 0.0226 0.8585 0.3905
Desv_Forest -0.1257 0.05334 -2.3566 0.0184
NBI 30.8457 10.0117 3.0809 0.002
CAA 10.7117 7.1895 1.4899 0.1362

Significant variables can express a large amount of the model, in some cases they must be transformed to adapt better to the distribution of the data, so it is not advisable to eliminate variables in the first instance only because they are not immediately significant. Some variables such as NBI can explain a lot of the model because it is related to CAA, although it is not significant, so to finalize the construction of a regression model requires an interpretation of the variables, the scope and limitations of them within the spatial context that is addressed.

4. Conclusions

The phenomenon transmitted by malaria vectors has a tendency of self-spatial correlation, where high-value groupings occur and another one at the opposite extreme with low values. The dispersion of malaria is an effect of the presence of some factors that allow the adequate incubation and proliferation of the virus.

The sanitation factor transmitted as sewer and water coverage is significant since it generates a greater incubation of the mosquito that transmits malaria, as well as the forest cover that generates its own conditions for its development. Some variables such as precipitation and height are linked to the results of forest cover, thus contributing to the model that follows a Poisson distribution affected by data in percentage rates. Some variables such as temperature were not significant under any of the regression models, so it was prudent to dispense with it since the temperature throughout the department is quite fluctuating and in the regression model could affect the final result. The quality of life of people is an agent that accelerates the proliferation of infection and insects and more in departments where the rate of NBI is considerably high.

The study allows to conclude the relationship of the phenomenon among municipalities, taking into account how the mortality of high malaria cases are related in space as can be seen in the local Moran map and is perfectly comparable with the exploratory statistics chart. Studies like this one show the current needs of health that have not been widely addressed in aspects such as the geographical one that can open a new perspective of the malaria context in Colombia. It is recommended to analyze this type of diseases over time and validate more advanced statistical models with conjugated distributions that can explain the behavior of the data in a more precise way.

5. References

1. Ministerio de Salud y Protección Social & Federación Médica Colombiana. Malaria. Memorias. 2012. [cited 2017 Oct 11]; Available from: https://www.minsalud.gov.co/sites/rid/Lists/BibliotecaDigital/RIDE/VS/TH/memorias_malaria.pdf. [ Links ]

2. Instituto Nacional de Salud. Malaria [Internet]. 2018 [cited 2018 Apr 18]. Available from: http://simposiovirologia.ins.gov.co/temas-de-interes/Paginas/malaria.aspx. [ Links ]

3. Instituto Nacional de Salud. Guía protocolo para la vigilancia en salud pública de malaria.. PRO-R02.003.0000-021. Bogota: Subdirección de Vigilancia y Control en Salud Pública, INS; 2010. Available from: https://www.minsalud.gov.co/Documents/Salud Pública/Ola invernal/Protocolo Vigilancia MALARIA.pdf. [ Links ]

4. Sarmiento MC. Enfermedades Transmisibles en Colombia: Cambios Ambivalentes. Salud Pública. 2000; 2(1): 82-93. [ Links ]

5. Rueda CNY. Boletín Alertas epidemiológicas. Bogota: Secretaría Distrital de Salud, Alcaldia de Bogota; 2016. [cited 2018 Apr 20]. Available from: http://www.saludcapital.gov.co/Boletines de alertas en salud pblica de inters int/Alertas epidemiologicas 2 2016.pdf. [ Links ]

6. Olano VA, Brochero HL, Saenz R, Quiñones ML, Molina JA. Mapas preliminares de la distribución de especies de Anopheles vectores de malaria en Colombia. Biomedica. 2001; 21(4):402-408. [ Links ]

7. Gomez AM. Planificación y estructuración del territorio desde la perspectiva de los sistemas de agua potable y saneamiento básico, las dimensiones ambiental y etnográfica. Tesis de Maestria. Maestria en Planeación Urbana y Regional, Universidad Javeriana. Bogota; 2010. 51 p. Available from: https://repository.javeriana.edu.co/bitstream/handle/10554/815/arq36.pdf?sequence=1Links ]

8. Galvis-Aponte LA, Moyano-Támara LM, Alba-Fajardo CA. La persistencia de la pobreza en el Pacífico colombiano y sus factores asociados. Documentos de Trabajo sobre economía regiona. Numero 238. Banco de la Republica, Centro de Estudios Economicos Regionales: Cartagena; 2016. Available from: http://www.banrep.gov.co/docum/Lectura_finanzas/pdf/dtser_238.pdfLinks ]

9. Arriaga RMI, Machuca CJV, Solarte ED. Perfil audiológico de niños y niñas de 7 a 12 años de una institución educativa en Quibdó-Chocó. Especialización en Audiología. Corporación Universitaria Iberoamericana: Bogotá; 2017. Available from: http://jbb-repositorio.metabiblioteca.org:8080/bitstream/001/429/1/Perfil audiológico de niños y niñas de 7 a 12 años de una institución educativa en Quibdó-Chocó.pdf [ Links ]

10. Pardo RM. Regionalizacion de indígenas Choco Datos etnohistóricos, lingüísticos y asentamientos actuales. Boletín Museo del Oro. 1987; 56(18):46-63. [ Links ]

11. Anonimo. Choco. División político-administrativa - departamentochocó. Blog [Internet]. 2016 [cited 2017 Oct 11]. Available from: https://departamentochoc.wordpress.com/2016/03/02/division-politico-administrativa/Links ]

12. todacolombia. Choco: Departamento de Choco Colombia - Informacion detallada Chocó Colombia [Internet]. 2015 [cited 2017 Oct 11]. Available from: http://www.todacolombia.com/departamentos-de-colombia/choco.htmlLinks ]

13. Aimone AM, Brown P, Owusu-Agyei S, Zlotkin SH, Cole DC. Impact of iron fortification on the geospatial patterns of malaria and non-malaria infection risk among young children: A secondary spatial analysis of clinical trial data from Ghana. BMJ Open. 2017;7(5):1-10. [ Links ]

14. Ochoa J, Osorio L. Epidemiología de malaria urbana en Quibdó, Chocó. Biomédica. 2006; 26(2):278-285. [ Links ]

15. DANE. Necesidades básicas insatisfechas (NBI) [Internet]. 2015 [cited 2017 Oct 11]. Available from: http://www.dane.gov.co/index.php/estadisticas-por-tema/pobreza-y-condiciones-de-vida/necesidades-basicas-insatisfechas-nbiLinks ]

16. Alcaldia de Quibdo. Plan de Ordenamiento Territorial -Municipio de Quibdó Diagnostico Territorial [Internet]. 2003 [cited 2017 Oct 11]. Available from: https://www.quibdo-choco.gov.co/Proyectos/NuestrosPlanes/Plan de Ordenamiento Territorial - Municipio de Quibdó.pdfLinks ]

17. Alcaldia de Tado. Esquema de ordenamiento territorial. 2000 [cited 2017 Oct 11]. Available from: http://tado-choco.gov.co/apc-aa-files/35653435636434623663326137616166/E.O.T.Municipio_de_Tad_.pdfLinks ]

18. Alcaldia de Lloro. Esquema de ordenamiento territorial [Internet]. 2000 [cited 2017 Oct 12]. Available from: http://cdim.esap.edu.co/BancoMedios/Documentos PDF/eot_proyecto_acuerdo_lloró_chocó_2000_(69_pag_190_kb).pdfLinks ]

19. Alcaldia de Riosucio. Esquema de Ordenamiento Territorial [Internet]. 2005 [cited 2017 Oct 11]. Available from: http://cdim.esap.edu.co/BancoMedios/Documentos PDF/riosuciochocoeot2005-2016.pdfLinks ]

20. Alcaldia de Itsmina. Esquema de Ordenamiento territorial [Internet]. 2005 [cited 2017 Oct 11]. Available from: http://cdim.esap.edu.co/BancoMedios/Documentos PDF/istminachocoeot2005.pdfLinks ]

21. Sanchez PDL, Santa F, Fuentes LHJ. Modelos de regresión espacial para el comportamiento de las enfermedades infecciosas dengue y malaria en Colombia para los años 2000, 2005 y 2010. Revista GEOMATICA UD.GEO. 2012; 6: 110-120. [ Links ]

22. Pfeiffer DU, Robinson TP, Stevenson M, Stevens KB, Rogers DJ, Clements ACA. Spatial Analysis in Epidemiology. OUP Oxford; 2008. 81-109 p. [ Links ]

23. Waller LA, Gotway CA. Applied Spatial Statistics for Public Health Data. Wiley Series in Probability and Statistics. John Wiley & Sons; 2004. [ Links ]

24. Brunsdon C, Fotheringham S, Charlton M. Geographically Weighted Regression. Journal of the Royal Statistical Society Series D The Statistician. 1998; 47: 431-443. [ Links ]

1 INS. The vector receives the pathogenic organism from an infected carrier, animal or human, and transmits it either to an intermediate carrier or directly to a human carrier. Fuente http://simposiovirologia.ins.gov.co/temas-de-interes/Paginas/que-son-enfermedades-transmitidas-por-vectores.aspx

2SIVIGILA: Epidemiological surveillance system.

3Function that draws a histogram given a variable. https://www.rdocumentation.org/packages/GeoXp/versions/1.6.2/topics/histomap

Received: October 12, 2017; Accepted: May 28, 2018

Creative Commons License This is an open-access article distributed under the terms of the Creative Commons Attribution License