Open access peer-reviewed chapter

Multiple Regression for the Schistosomiasis Positivity Index Estimates in the Minas Gerais State - Brazil at Small Communities and Cities Levels

Written By

Ricardo J.P.S. Guimarães, Corina C. Freitas, Luciano V. Dutra, Guilherme Oliveira and Omar S. Carvalho

Submitted: May 3rd, 2012 Reviewed: September 17th, 2012 Published: January 16th, 2013

DOI: 10.5772/53500

Chapter metrics overview

1,801 Chapter Downloads

View Full Metrics

1. Introduction

Schistosomiasis, caused by Schistosoma mansoni, is an endemic disease conditional on the presence of snails of aquatic habits of the genus Biomphalaria.

In Brazil, there are eleven species and one subspecies of Biomphalariagenus mollusks that have been identified: B. glabrata(Say, 1818), B. tenagophila(Orbigny, 1835), B. straminea(Dunker, 1848), B. peregrina(Orbigny, 1835), B. schrammi(Crosse, 1864), B. kuhniana(Clessin, 1883), B. intermedia(Paraense & Deslandes 1962), B. amazonica(Paraense 1966), B. oligoza(Paraense 1974), B. occidentalis(Paraense 1981), B. cousini(Paraense, 1966) and B. tenagophilaguaibensis(Paraense 1984) [1].

In Minas Gerais state, the presence of seven species: B. straminea, B. tenagophila, B. peregrina, B. schrammi, B. intermediaand B. occidentaliswas reported [1]. Among these, there are three Biomphalariaspecies (B. glabrata, B. tenagophilaand B. straminea) that have been found to be naturally infected with S. mansoni. Other three species, B. amazonica, B. peregrinaand B. cousini, were experimentally infected, being considered as potential hosts of this trematode [2-4]. B. glabratais of great importance, due to its extensive geographic distribution, high infection indices and efficiency in the schistosomiasis transmission. In endemic areas, large concentrations of these snails, together with other risk factors, favor the existence of localities with high prevalence [5-7].

The snails of the Biomphalariagenus live in a wide range of habitats, particularly in shallow and slow running waters and with floating or rooted vegetation. As these snails are distributed over large geographic areas and their populations are adapted to different environmental conditions, they can tolerate large variations in physical, chemical and biological environment in which they live [8, 9].

The intermediate hosts’ distribution of the parasite in Minas Gerais associated with favorable eco-epidemiological conditions gives the schistosomiasis expansive character not seen even in non-endemic regions [6, 10, 11].

Public health and the environment are influenced by the patterns of space occupation. Therefore, the use of geoprocessing techniques to analyze the spatial distribution of health problems allows one to determine local risks and delimit areas that concentrate the most vulnerable situations (occurrence of disease, characteristics of the environment and habitat of the intermediate host / vector). It is also possible with the use of geographic information systems to plan, schedule, control, monitor, and evaluate the diseases in groups according to their risk of transmission [12].

The use of Geographic Information Systems (GIS) and statistical tools in health has been facilitated by access to epidemiological data bases, enabling the production of thematic maps that contribute to the formulation of hypotheses about the spatial distribution of diseases and their relation to the socioeconomic variables [13].

The use of GIS and Remote Sensing (RS) are powerful tools for working complex analysis of a large number of information and viewing the results of this analysis in graphical maps. Since the seventies, RS has been applied to social sciences and health [14]. There are numerous information collected by RS data, describing some biotic and abiotic factors [15]. Application of RS and GIS techniques for mapping the risk of parasitic diseases, including schistosomiasis, has been performed over the past 15 years [16].

The estimate of schistosomiasis prevalence using GIS was first used in the Philippines and the Caribbean by [17, 18]. In Brazil, the use of GIS in schistosomiasis was first used by [19] in the state of Bahia. The authors constructed maps with environmental characteristics (total precipitation for three consecutive months, the annual maximum and minimum temperature and diurnal temperature differences), prevalence of S. mansoniand distribution of snails to study the spatial and temporal dynamics of infection and identify the environmental factors that influence the distribution of schistosomiasis. The results indicated that the snail population density and duration of annual dry season are the most important determinants for the prevalence of schistosomiasis in the study areas.

Table 1 shows a brief history of the use of GIS techniques in the study of schistosomiasis in several countries.

The main objective of the present study is to establish a relationship between schistosomiasis positivity index and the environmental and socioeconomic variables, in the Minas Gerais State, Brazil, using multiple linear regressions at small communities and cities levels.

VectorSpeciesStudy AreaSatellite-sensorTechnical-variablesReference
-Schistosoma sppPhilippines, the CaribbeanLandsat (MSS)climate[17, 18]
Oncomelania sppSchistosoma sppChinaNOAA (AVHRR), Lansdat (TM)ecological zones[20]
B. truncatus,
B. alexandrina
S. mansoni,
S. haematobium
EgyptNOAA (AVHRR)temperature, NDVI[21]
Phlebotomus papatasiSchistosoma sppSoutheast AsiaNOAA (AVHRR)NDVI[22]
B. stramineaS.mansoniBrazil-spatial distribution[23]
B. pfeifferiS.mansoniKenya-linear regression, mapping techniques, cluster analysis[24]
B. alexandrinaS. mansoni,
S. haematobium
EgyptNOAA (AVHRR), Lansdat (TM)dT, NDVI, MDE[25]
B. glabrata, B. straminea, B. tenagophilaS.mansoniBrazil-temperature, precipitation, DEM, soil type, vegetation type[19]
Oncomelania sppSchistosoma sppChinaLandsat (TM)classification, GIS[26]
S. mansoni,
S. haematobium
Tanzania-GIS, logistic regression[27]
B. alexandrinaS. mansoniEgyptNOAA (AVHRR), Lansdat (TM)dT, BED, NDVI[28]
B. glabrata, B. stramineaS. mansoniBrazilNOAA (AVHRR)NDVI, dT[29]
Bulinus sppS. haematobiumTanzaniaNOAA (AVHRR)LST, NDVI, DEM, precipitation, logistic regression[30]
B. pfeifferiS. mansoniEthiopiaNOAA (AVHRR)LST, NDVI[31]
B. pfeifferiS. mansoniEthiopiaNOAA (AVHRR)NDVI, temperature, logistic regression[32]
Oncomelania sppS. japonicumChinaNOAA (AVHRR), Lansdat (TM)TNDVI[33]
-Schistosoma sppChad, CameroonNOAA (AVHRR)ecology[34]
B. pfeifferi, B. senegalensisS. mansoni,
S. haematobium
Africa (sub-Saharan Africa)NOAA (AVHRR)SIG[35]
-S. haematobiumChad-environmental data[36]
-S. mansoni,
S. haematobium
CameroonNOAA (AVHRR), EROSlogistic regression[37]
Oncomelania sppSchistosoma sppChinaLandsat (TM)RS[38]
Oncomelania hupensisS. japonicumChinaLandsat (TM)LU[39]
S. japonicum,
S. mansoni,
S. haematobium
ChinaLandsat (TM)SIG[40]
Oncomelania hupensisS. japonicumChinaLansdat (TM)TNDVI[41]
Bulinus spp.S. haematobiumKenyaNOAA (AVHRR)Tmax[42]
B. glabrataS. mansoniBrazil-GPS[43]
-S. mansoni,
S. haematobium
UgandaLandsat (TM)ecological zones[44]
-Schistosoma sppUgandaNOAA (AVHRR)LST[45]
Oncomelania hupensisS. japonicumJapan-PDA[46]
Oncomelania hupensisS. japonicumChinaIkonos, ASTERMDE[47]
-S. mansoniCôte d'IvoireLandsat (ETM), NOAA (AVHRR), EROS, MODISenvironmental and socioeconomic data[16, 48]
O. hupensisS. japonicumChinaLandsat (TM)NDVI[49]
Oncomelania hupensisS. japonicumChinaNOAA (AVHRR)LST[50]
Oncomelania hupensisS. japonicumChinaLandsat (TM)SAVI[51]
-S. mansoniBrazil-logistic regression models and Bayesian spatial models[52]
Biomphalaria spS.mansoniBrazilMODIS, SRTMregression, elevation, mixture model, NDVI[53]
B. glabrataS.mansoniBrazil-spatial analysis, GPS, immunological data[54]
Biomphalaria spS.mansoniBrazil-social and environmental data, regression[55]
Biomphalaria sppS. mansoniBrazil-GPS and GIS[56]
B. glabrataS. mansoniBrazil-kernel[57]
-Schistosoma sppAfrica-ecology, GIS, RS, geostatistics[58]
B. pfeifferiS. mansoniCôte d'Ivoire-socioeconomic data, logistic regression, Bayesian model[59]
B. sudanica, B. stanleyiS. mansoniUganda-spatial analysis[60]
-S. haematobiumTanzania-social and ecological data, Bayesian models, logistic regression, NDVI, elevation, cluster analysis[61]
Biomphalaria spS.mansoniBrazilMODIS, SRTMsocial and environmental data, RS, NDVI, temperature, regression[62]
Biomphalaria spS.mansoniBrazilMODIS, SRTMlinear regression, imprecise classification, regionalization and pattern recognition[63]
Oncomelania hupensisS. japonicumChinaSPOTecological data, land use, land cover, classification, Bayesian model, RS, NDVI, slope, LST[64]
-Schistosoma sppChina-GIS, spatial analysis and clustering, Bayesian model,[65]
-S. mansoni,
S. haematobium
MaliNOAA (AVHRR)Bayesian models, NDVI, LST, GIS, logistic regression[66]
Biomphalaria spS.mansoniBrazil-kriging, spatial distribution[67, 68]
Biomphalaria spS.mansoniBrazil-Fuzzy logic[69]
-S. japonicumChina-GIS, spatial analysis, clustering, kernel[70]
Biomphalaria spS. mansoniBrazilMODISmeteorological data, socioeconomic, sanitation, RS, regression[71]
B. stramineaS.mansoniBrazil-kernel, GPS, spatial distribution[72]
B. glabrataS.mansoniBrazilMODISmixture model[73]
Biomphalaria spS.mansoniBrazilMODIS, SRTMsocial and environmental data, sanitation, biological, RS, NDVI, temperature, regression, kriging[1, 74]
Biomphalaria spS.mansoniBrazilMODISdecision tree, environmental data, RS[75]
-S. mansoni,
S. haematobium
East Africa-Bayesian geostatistics, logistic regression, Markov chain Monte Carlo simulation,[76]
Biomphalaria spS.mansoniAfricaMODISClimate change, spatial distribution, temperature, precipitation, MaxEnt, soil[77]
-S.mansoniEthiopia, KenyaNOAA (AVHRR)geostatistics, LST, NDVI, elevation, environmental data, LQAS, LpCP[78]
Biomphalaria sppS. mansoniBrazil-GPS, GIS, spatial distribution[79]

Table 1.

Use of geoprocessing techniques in the study of schistosomiasis.


2. Material and methods

2.1. Materials

The study area includes 4,846 small communities (called localities) in the entire State of Minas Gerais, Brazil. The dependent variable is the schistosomiasis positivity index (Ip). Ipwere obtained from the Brazilian Schistosomiasis Control Program (PCE) through the annual reports of the Secretary of Public Health Surveillance (SVS) and the Secretary of Health in the State of Minas Gerais (SESMG). From the 4,846 locations mentioned above, only 1,590 of them have information on the positivity of the disease. Since schistosomiasis is a disease characterized by environmental and social factors, environmental and socioeconomic variables were used as explanatory variables, as well as a variable containing information about presence of intermediate hosts. A brief description of these variables is given below.

2.2. Schistosomiasis positivity index

Schistosomiasis positivity index (Ip) values were obtained in 1,590 localities from the Brazilian Schistosomiasis Control Program (PCE) through the Annual Reports of the Secretary of Public Health Surveillance (SVS) and the Secretary of Health in the State of Minas Gerais (SESMG). The Ipdata were obtained from the database SISPCE (Information System of the Brazilian Schistosomiasis Control Program) from 1996 to 2009. The Kato-Katz technique is the methodology used to determine positivity index, examining one slide per person.

These Ipwere determined for each locality iby:


where: riis the number of infected people and niis the total population in locality i.

2.3. Intermediate hosts

Information about the existence of Biomphalariasnails were provided at a municipality basis by the Laboratory of Helminthiasis and Medical Malacology of the Rene Rachou Research Center (CPqRR/Fiocruz-MG).

The distribution of Biomphalaria snails used for this study was defined as: B. glabrata, B. tenagophila, B. straminea, B. glabrata + B. tenagophila, B. glabrata + B. straminea, B. tenagophila + B. straminea, B. glabrata + B. tenagophila + B. stramineaand No Biomphalaria. The class “No Biomphalaria” includes information about the non-occurrence of Biomphalariaspecies or information about non-transmitter species in Brazil, such as B. peregrina, B. schrammi, B. intermedia, B. occidentalis, etc.

The spatial distribution of the schistosomiasis Ipand the Biomphalariaspecies data are presented in Fig. 1.

2.4. Environmental data

Twenty eight environmental variables were obtained from remote sensing and meteorological sources.

The remote sensing variables were derived from Moderate Resolution Imaging Spectroradiometer (MODIS) and from the Shuttle Radar Topography Mission (SRTM) sensor.

The variables of MODIS sensor used were collected in two seasons, summer (from 17/Jan/2002 to 01/Feb/2002 period) and winter (from 28/Jul/2002 to 12/Aug/2002 period). MODIS data were composed by the blue, red, near and middle infrared bands and also the vegetation indices (NDVI and EVI) [73].

Figure 1.

Spatial distribution of the (a) schistosomiasis and (b)Biomphalariaspecies in Minas Gerais State, Brazil.

The Linear Spectral Mixture Model (LSMM) is an image processing algorithm that generates fraction images with the proportion of each component (vegetation, soil, and shade) inside the pixel, which is estimated by minimizing the sum of square of the errors. In this work, the so called vegetation, soil, and shade fraction images were generated using the MODIS data, and the estimated values for the spectral reflectance components were also used as an input to the regression models [73].

Others variables obtained from SRTM were also used in this study: the digital elevation model (DEM) and slope (derived from DEM). Based on the SRTM data, a drainage map of Minas Gerais was generated, and the variables: water percentage in municipality (QTA) and water accumulation (WA) were derived. Six meteorological variables consisting of total precipitation (Prec), minimum (Tmin) and maximum (Tmax) temperature average for summer and winter seasons were obtained from the Center for Weather Forecast and Climate Studies (CPTEC), in the same date of MODIS images.

2.5. Socioeconomic data

Socioeconomic variables obtained by The Brazilian Institute of Geography and Statistics (IBGE) census for the year 2000 were also used as explanatory variables. The variables used in this work are those related to the water quality (percentage of domiciles with access to the general net of water supply, access to the water through wells or springheads, and with other access forms to the water), and to the sanitary conditions (the percentage of domiciles with bathroom connected to rivers or lakes, connected to a ditch, to rudimentary sewage, to septic sewage, to a general net, to other sewerage type, with bathroom or sanitarium and without bathroom or sanitarium).

2.6. Methods

Indicator kriging and multiple linear regressions were employed to estimate the presence of the intermediate host and the schistosomiasis disease, respectively.

2.7. Indicator Kriging

Since information about existence of Biomphalariais only available on municipality basis, indicator kriging was used in this study to make inferences, in a grid basis, about the presence of the Biomphalariaspecies (B. glabrata, B. tenagophilaand/or B. straminea), intermediate hosts of S. mansoni. The method allows spatialization of the data conditioned to the sample set of categorical attributes, aiming at the spatial distribution and production of maps.

The categorical attributes (classes) used for this study were defined as: B. glabrata, B. tenagophila, B. straminea, B. glabrata + B. tenagophila, B. glabrata + B. straminea, B. tenagophila + B. straminea, B. glabrata + B. tenagophila + B. stramineaand No Biomphalariatotalizing eight probable classes.

The snail attributes (class of species and localization) were distributed along the drainage network of 15 River Basins (Buranhém, Doce, Grande, Itabapoana, Itanhém, Itapemirim, Jequitinhonha, Jucuruçu, Mucuri, Paraíba do Sul, Paranaíba, Pardo, Piracicaba/Jaguari, São Francisco and São Mateus), according to the methodology used by [67].

In [1], however the indicator kriging was used only at municipalities’ level, but in this study it was used for localities level.

Indicator kriging procedures were applied to obtain an approximation of the conditional distribution function of the random variables. Based on the estimated function, maps of snail spatial distributions along with the corresponding uncertainties for the entire state and also map of estimated prevalence of schistosomiasis were built.

The indicator kriging result was used as a variable in multiple regression models.

2.8. Multiple linear regressions

Multiple linear regressions are a form of regression analysis in which data are modeled by a least squares function which is a linear combination of the model parameters and depends on one or more independent variables.

The regression analysis was applied with the schistosomiasis Ipas dependent variable, in addition to 93 quantitative variables (28 environmental variables and 65 socioeconomic variables), and one qualitative variable resulting from the kriging (presence or absence of B. glabrata) as explanatory variables.

The dependent variable was randomly divided into two sets: one with 852 cases (localities) for variables selection and model definition, and another with 738 cases for model validation.

Due to the high number of independent variables, some procedures were performed for variables selection. The relations among the dependent and the independent variables were analyzed in terms of correlation, multi co-linearity, and possible transformations that better explain the dependent variable.

A logarithmic transformation for the dependent variable (denoted by lnIp) was made as it improved the correlation with independent variables.

The analysis of the correlation matrix showed that some variables had non-significative correlations with lnIpat 95% confidence level, and also some variables were highly correlated among themselves, indicating that those variables could be excluded from future analysis.

Since multi co-linearity effects among the remaining independent variables were detected, variables selection techniques were used in order to choose a set of variables that better explain the dependent variable. Variable selection was performed by the R2 criterion using all possible regressions [80].

This selection technique consists in the identification of a best subset with few variables and a coefficient of determination R2 sufficiently close to that when all variables are used in the model.

Interaction effects were also analyzed to be included in the model. After performing the residual analysis, the chosen regression model was then validated. The final estimated regression function was computed using the entire data set (definition and validation), and it was applied to all localities to build a risk map for schistosomiasis positivity index.

The multiple regressions were developed based on two approaches: a global model (throughout the state) and a regional model (regionalization).

Regionalization is a classification procedure using the SKATER algorithm (Spatial ‘K’luster Analysis by Tree Edge Removal) applied to spatial objects with an areal representation (municipalities), which groups them into homogeneous contiguous regions [81].

Regionalization was applied in Minas Gerais to divide the state into four homogeneous regions. The choice of the number of regions was based on the spatial distribution of localities (Figure 1b) in order to achieve an adequate number of localities in each region.

The regional model was developed by doing a regression model separately in each of the four regions formed by first applying the SKATER algorithm using environmental variables [74].

The models validation was performed using the Root Mean Square Error (RMSE) and the Mean Squared Prediction Error (MSPR), given by [80].


where Ipi and I^pirepresent, respectively, the observed and predicted positivity index in the i-th observation and nis the number of observations of the data model definition (i =1,..., n).

The RMSE measures the variation of the observed values around the estimated values. Ideally, the values of RMSE are close to zero. The MSPR is computed the same way as the RMSE, but using validation samples.

The final models were applied in all 4,846 localities to estimate the positivity index.

2.9. Simple average interpolator

The simple average interpolator (SAI) algorithm of the software SPRING [82] was used to estimate the value of Ipat each point (x,y) of the grid. This estimative is based on the simple average of the variable values in the eight nearest neighbors of this point, according to equation (3).


where I^piis the estimated positivity index of the 8 neighbors of the point (x,y) and f(x,y) is the interpolation function.

The file generated by interpolation was a grid with spatial resolution of 1 km. The purpose of using this tool was to determine which of the mesoregions presented estimated values above 15% (class with high positivity index).


3. Results and discussion

The GeoSchisto Database ( was created containing all variables used in this study.

The indicator kriging result was a regular grid of 250 x 250 meters with the estimate of Biomphalariaspecies class for the entire Minas Gerais State. The indicator kriging result is presented in Fig. 2a. The variable B. glabrataused in regression models is presented in Fig. 2b.

Figure 2.

a) Kriging and (b) estimated presence ofB. glabrataby indicator kriging.

3.1. Global model

The five variables selected were: presence or absence of the B. glabrata, summer precipitation (PCs), summer minimum temperature (TNs), winter Enhanced Vegetation Index (EVIw) and households with a bathroom or toilet and sewage from septic tank type (V31).

The final model, with R2 = 0.18, was:


Fig. 3a shows the estimated Ipfor all 4,842 localities in Minas Gerais using the estimated regression equation (4). Figure 3b shows the plot of the residuals, resulting from the difference between observed and estimated Ipfrom 1,590 locations. In Figure 3b, dark colors (red and blue) represent overestimated values, light colors (red and blue) underestimated ones, and in white are the municipalities where the estimated prevalence differs very little from the true values.

The precipitation, minimum temperature, EVI and sanitation were positively correlated with Ip. This is consistent with the adequate environmental conditions for the transmission of schistosomiasis. Also, the transmission depends on the presence of B. glabrata.

Figure 3.

Global model: (a) estimatedIpand (b) residuals.

The result of this model has the same variables (BG, TNs, Eviwand sanitation) obtained by [74] when estimatives were done on a municipality basis, indicating a great similarity between the two global models. The difference is in the sanitation variable where the variable obtained by [74] was related to the type of water (well or spring) and this study to the type of sewage system (septic tank).

3.2. Regional model

The Minas Gerais State was divided into four regions using the SKATER algorithm. Table 2 presents the number of localities in each region used for model generation and for model validation. The regionalization can be seen in Figure 4 a.

Model GenerationModel ValidationTotal
Region 1 (R1)10466170
Region 2 (R2)428262690
Region 3 (R3)220338558
Region 4 (R4)10072172

Table 2.

Number of localities in each region used for model generation and for model validation.

Regression models were developed for each of the four regions with the same 94 variables used in the global model, and the same selection procedure. Different numbers of variables were selected in each region to determine the best regression model.

The final models generated for each region (Fig. 4c) and their R2 were:


where: PCW(winter precipitation), V25(percentage of households with another form of access to water), ΔTS(difference of summer maximum and minimum temperature), BG(presence or not of the B. glabrata), TNS(summer minimum temperature), EVIW(winter Enhanced Vegetation Index), V261(percentage of residents in households with another form of water supply), NDVIS(summer Normalized Difference Vegetation Index), V33(percentage of housing with bathroom or toiled connected to a ditch), V283(percentage of households without bathrooms), EVIS(summer Enhanced Vegetation Index), PCS(summer precipitation), QTA(water percentage in municipality), V254(percentage of households with water supply network general) and V272(percentage of households without toilet or sanitation).

Figure 4.

Regional model: (a) estimatedIp, (b) residuals. (c) regionalization.

Figure 4a shows the estimated values of Ipfor all 4,842 localities in the Minas Gerais State using equations (3, 4, 5 and 6). Also, Figure 4b shows the residues from 1,590 locations. In this figure, red and blue represent overestimates, cyan and magenta represent the underestimated values and in the white localities with good estimate.

The regional model for Region 1 (R1) reflects the effect of sanitation (households with other forms of water than tap water, wells or springs) and the influence of weather (precipitation and temperature of summer). Region 1 achieved a R2 value of 0.35. The model obtained by [74] for the same Region 1 also has the same sanitation variable (percentage of homes with another type of access to water). The relationship between temperature and disease was also obtained by [29] and [55, 74].

The models for Regions 2 and 3 (R2 and R3) show the presence of B. glabrataassociated with the effect of vegetation (Eviw) and sanitation. Among the regional models, Region 2 had the lowest R2 (0.21) and Region 3 had the highest R2 (0.38).

The model for Region 4 (R4) shows that Ipwas associated with vegetation (Evis), weather (precipitation and temperature) and sanitation (type of water and sewage). The R2 found for this model was 0.22.

In all models the presence of B. glabrata, sanitation, vegetation index and temperature were the most important variables. These characteristics are the same as environmental conditions for the presence and development of snails (infection of the intermediate host) and sanitation (water contamination - presence of S. mansonicercariae) obtained by [74] which were obtained at municipalities level.

[29] also showed that the distribution of schistosomiasis in Bahia, at municipalities level, is related to the vegetation index (NDVI) and temperature (ΔTs) using sensor data from low spatial resolution (AVHRR/NOAA).

3.3. Simple Averages Interpolator (SAI)

Table 3 presents the mean square error (RMSE) and Mean Squared Error of Prediction (MSPR) for the global and regional models, for each region. From this table we can observe that the mean square decreased from 10.739 to 9.979 when we used separate models for each region. It was also noted that the RMSE of the Regional Model was smaller than the RMSE of the Global Model for all four regions, highlighting the importance of using different equations and different variables for each region. Since the Regional Model can be considered a better model the Simple Averages Interpolator (SAI), was applied using the known positivity index of the 1,590 localities (Fig. 5a), and using the regression estimated positivity index of all 4,842 localities (Fig. 5b). The objective of applying SAI to all estimated Ipvalues is to indicate current and potential local transmission of schistosomiasis.


Table 3.

Residual analysis of the dependent variable (Ip) for the models.

n (number of localities)

Figure 5a shows clusters presence in six mesoregions (Norte de Minas, Jequitinhonha, Vale do Mucuri, Vale do Rio Doce, Metropolitana de Belo Horizonte and Zona da Mata) with the highest Ipvalues. In Figure 5b the same six mesoregions can be noticed; however two news clusters in Sul/Sudoeste de Minas and Triângulo Mineiro/Alto Parnaíba mesoregions presented, respectively high and middle Ipvalues.

Figure 5.

The averages interpolator: (a)Ipand (b) estimatedIpby Regional Model. (c) Mesoregions of Minas Gerais State.

Thus, the Norte de Minas, Jequitinhonha, Vale do Mucuri, Vale do Rio Doce, Metropolitana de Belo Horizonte and Zona da Mata mesoregions are endemics areas.

Sul/Sudoeste de Minas and Triângulo Mineiro/Alto Parnaíba mesoregions are not endemic areas, but have a schistosomiasis focus (Itajubá municipality in Sul mesoregion). The Sul/Sudoeste de Minas mesoregion has 146 municipalities representing about 20% of municipalities in Minas Gerais State and is a non-endemic area for schistosomiasis. Due to the high concentration of cities in an area of 49,523.893 km2 (which represents less than 10% of the area of Minas Gerais State) and a high agricultural economy, it is a region with high risk of schistosomiasis transmission. Therefore, it would be interesting to do a detailed study in the Sul mesoregion to determine the schistosomiasis Ip.

Also, it would be interesting to keep surveillance in the municipalities of the Triângulo Mineiro/Alto Parnaíba mesoregion that presented B. glabratapresence.


4. Conclusions and future work

This study shows the importance of a joint use of GIS and RS to study the risk of schistosomiasis. Moreover, it can be concluded that the combined use of GIS and statistical techniques allowed the estimation of schistosomiasis Ip. Results of the regression models confirmed the importance of the use of environmental variables to characterize the snail habitat in the endemic area of the state of Minas Gerais.

Results of the regression models show that regionalization improves the estimation of the disease in Minas Gerais. Based on this model, a schistosomiasis risk map was built for Minas Gerais. [74] and [75] also obtained a better model with the use of regionalization when estimating schistosomiasis at a municipality level.

The Simple Averages Interpolator is a technique that may indicate possible local to transmission and surveillance of schistosomiasis.

It is recommended the use of GPS for field surveys together and the application of this methodology with images of better spatial resolution (10-30m) in other states for validation. Also, we recommend using a smaller area (municipality or mesoregion) estimate for the schistosomiasis.

The methodology used in this study can be utilized to control schistosomiasis in the areas with occurrence of the disease and also it can be used to take preventive measures to prevent the disease transmission.

Next step will be to utilize data from the PCE by localities to study other diseases such as ascariasis, hookworm, trichuriasis, etc, using data from CBERS and/or Landsat and new methodologies (Geographically Weighted Regression, Generalized Additive Model, etc).



The authors woud like to acknowledge the support of Sandra da Costa Drummond (Fundação Nacional de Saúde) and the support of CNPq (Conselho Nacional de Desenvolvimento Científico e Tecnológico) (grants # 300679/2011-4, 384571/2010-7, 302966/2009-9, 308253/2008-6).


  1. 1. Guimarães RJPS, Fonseca FR, Dutra LV, Freitas CC, Oliveira GC, Carvalho OS. A study of schistosomiasis prevalence and risk of snails presence spatial distributions using geo-statistical tools. In: MB R, (ed.). Schistosomiasis. Rijeka: InTech, 2012. p255-80.
  2. 2. Correa LR, Paraense WL. Susceptibility ofBiomphalaria amazonicato infection with two strains ofSchistosoma mansoni. Rev Inst Med Trop São Paulo. 1971;13 387-90.
  3. 3. Paraense WL, Correa LR. Susceptibility ofBiomphalaria peregrinafrom Brazil and Ecuador to two strains ofSchistosoma mansoni. Rev Inst Med Trop São Paulo. 1973;15 127-30.
  4. 4. Caldeira RL, Teodoro TM, Gomes MFB, Carvalho OS. Preliminary studies investigating the occurrence ofBiomphalaria cousiniin Brazil. Mem Inst Oswaldo Cruz. 2010;105(4) 485-7.
  5. 5. Carvalho OS, Rocha RS, Massara CL, Katz N. Expansão da esquistossomose mansoni em Minas Gerais. Mem Inst Oswaldo Cruz. 1987;82(Suppl. IV) 295-8.
  6. 6. Carvalho OS, Massara CL, Rocha RS, Katz N. Esquistossomose Mansoni no Sudoeste do Estado de Minas Gerais (Brasil). Rev Saúde Pública. 1989;23 341-4.
  7. 7. Carvalho OS, Massara CL, Silveira Neto HV, Guerra HL, Caldeira RL, Mendonça CLF, et al. Re-evaluation of Schistosomiasis Mansoni in Minas Gerais, Brazil. II. Alto Paranaiba Mesoregion. Mem Inst Oswaldo Cruz. 1997;92 141-2.
  8. 8. WHO. Study Group on the Ecology of Intermediate Snail Hosts of Bilharziasis. Geneva: World Health Organization; 1957.
  9. 9. Juberg P, Schall VT, Barbosa JV, Gatti MJ, Soares MS. Behavior ofBiomphalaria glabrata, the intermediate host snail ofSchistosoma mansoni, at different depths in water in laboratory conditions. Mem Inst Oswaldo Cruz. 1987;82 179-208.
  10. 10. Katz N, Carvalho OS. Introdução recente da esquistossomose mansoni no sul do estado de Minas Gerais, Brasil. Mem Inst Oswaldo Cruz. 1983;78 281-4.
  11. 11. Carvalho OS, Rocha RS, Massara CL, Katz N. Primeiros casos autóctones de esquistossomose mansonica em região do noroeste do Estado de Minas Gerais (Brasil). Rev Saúde Pública São Paulo. 1988;22 237-9.
  12. 12. Carvalho MS, Pina MF, Santos SM. Conceitos Básicos de Sistemas de Informação Geográfica e Cartografia Aplicados à Saúde. Brasília: OPS/Ministério da Saúde, 2000.
  13. 13. Rojas LI, Barcellos C, Petter P. Utilização de mapas no campo da epidemiologia no Brasil. Informe Epidemiológico do SUS. 1999;8 25-35.
  14. 14. Cline BL. New eyes for epidemiologists: aerial photography and other remote sensing techniques. Am J Epidemiol. 1970;92 85-9.
  15. 15. Beck LR, Lobitz BM, Wood BL. Remote Sensing and Human Health: New Sensors and New Opportunities. Emerg Infect Dis. 2000;6(3) 217-27.
  16. 16. Raso G, Matthys B, N'Goran EK, Tanner M, Vounatsou P, Utzinger J. Spatial risk prediction and mapping ofSchistosoma mansoniinfections among schoolchildren living in western Côte d'Ivoire. Parasitology. 2005;131(01) 97-108.
  17. 17. Cross ER, Bailey RC. Prediction of areas endemic for schistosomiasis through use of discriminant analysis of environmental data. Mil Med. 1984;149 28-30.
  18. 18. Cross ER, Sheffield C, Perrine R, Pazzaglia G. Predicting areas endemic for schistosomiasis using weather variables and a Landsat data base. Mil Med. 1984 Oct;149(10) 542-4.
  19. 19. Bavia ME, Hale L, Malone JB, Braud DH, Shane SM. Geographic information systems and the enviromental risk of Schistosomiasis in Bahia, Brazil. Am J Trop Med Hyg. 1999;60(4) 566-72.
  20. 20. Chen S, Hu J. Geo-ecological zones and endemic diseases in China -- a sample study by remote sensing. Preventive Veterinary Medicine. 1991;11(3-4) 335-44.
  21. 21. Malone JB, Huh OK, Fehler DP, Wilson PA, Wilensky DE, Holmes RA, et al. Temperature data from satellite imagery and distribution of schistosomiasis in Egypt. Am J Trop Med Hyg. 1994;51(3) 714-22.
  22. 22. Cross ER, Newcomb WW, Tucker CJ. Use of weather data and remote sensing to predict the geographic and seasonal distribution ofPhlebotomus papatasiin southwest Asia. Am J Trop Med Hyg. 1996;54(5) 530-6.
  23. 23. Teles HMS. Distribuição deBiomphalaria stramineaao Sul da Região Neotropical, Brasil. Rev Saude Publica. 1996;30(4) 341-9.
  24. 24. Kloos H, Fulford AJC, Butterworth AE, Sturrock RF, Ouma JH, Kariuki HC, et al. Spatial patterns of human water contact andSchistosoma mansonitransmission and infection in four rural areas in Machakos District, Kenya. Social Science & Medicine. 1997;44(7) 949-68.
  25. 25. Malone JB, Abdel-Rahman MS, El Bahy MM, Huh OK, Shafik M, Bavia M. Geographic information systems and the distribution ofSchistosoma mansoniin the Nile delta. Parasitol Today. 1997;13(3) 112-9.
  26. 26. Gong P, Spear R, Seto E, Zhou Y, Xu B, Maxzle D, et al. Remote sensing and GIS for schistosomiasis control in Sichuan, China, an overview. Proceedings of Geoinformatics'99; 1999 19-21 JUNE; Ann Arbor, MI; 1999. p. 1-9.
  27. 27. Lwambo NJS, Siza JE, Brooker S, Bundy DAP, Guyatt H. Patterns of concurrent hookworm infection and schistosomiasis in schoolchildren in Tanzania. Transactions of the Royal Society of Tropical Medicine and Hygiene. 1999;93(5) 497-502.
  28. 28. Abdel-Rahman MS, El-Bahy MM, Malone JB, Thompson RA, El Bahy NM. Geographic information systems as a tool for control program management for schistosomiasis in Egypt. Acta Trop. 2001;79(1) 49-57.
  29. 29. Bavia ME, Malone JB, Hale L, Dantas A, Marroni L, Reis R. Use of thermal and vegetation index data from earth observing satellites to evaluate the risk of schistosomiasis in Bahia, Brazil. Acta Trop. 2001;79(1) 79-85.
  30. 30. Brooker S, Hay SI, Issae W, Hall A, Kihamia CM, Lwambo NJ, et al. Predicting the distribution of urinary schistosomiasis in Tanzania using satellite sensor data. Trop Med Int Health. 2001 Dec;6(12) 998-1007.
  31. 31. Kristensen TK, Malone JB, McCarroll JC. Use of satellite remote sensing and geographic information systems to model the distribution and abundance of snail intermediate hosts in Africa: a preliminary model forBiomphalaria pfeifferiin Ethiopia. Acta Trop. 2001;79(1) 73-8.
  32. 32. Malone JB, Yilma JM, McCarroll JC, Erko B, Mukaratirwa S, Zhou X. Satellite climatology and the environmental risk ofSchistosoma mansoniin Ethiopia and east Africa. Acta Trop. 2001;79(1) 59-72.
  33. 33. Zhou XN, Malone JB, Kristensen TK, Bergquist NR. Application of geographic information systems and remote sensing to schistosomiasis control in China. Acta Trop. 2001;79(1) 97-106.
  34. 34. Beasley M, Brooker S, Ndinaromtan M, Madjiouroum EM, Baboguel M, Djenguinabe E, et al. First nationwide survey of th health of schoolchildren in Chad. Trop Med Int Health. 2002;7(7) 625-30.
  35. 35. Brooker S. Schistosomes, snails and satellites. Acta Tropica. 2002;82(2) 207-14.
  36. 36. Brooker S, Beasley M, Ndinaromtan M, Madjiouroum EM, Baboguel M, Djenguinabe E, et al. Use of remote sensing and a geographic information system in a national helminth control programme in Chad. Bull World Health Organ. 2002;80 783-9.
  37. 37. Brooker S, Hay SI, Tchuem Tchuenté LA, Ratard R. Using NOAA AVHRR data to model human helminth distribuitions in planning disease control in Cameroon. Photogrametric Engineering & Remote Sensing. 2002;68(2) 175-9.
  38. 38. Davis GM, Wu WP, Chen HG, Liu HY, Guo J-G, Lin D-D, et al. A baseline study of importance of bovines for humanSchistosoma japonicuminfections around Poyang Lake, China: villages studied and snail sampling strategy. Am J Trop Med Hyg. 2002;66 359-71.
  39. 39. Seto E, Xu B, Liang S, Gong P, Wu W, Davis GM, et al. The use of remote sensing for predictive modeling of schistosomiasis in China. Photogramm Eng Rem Sens. 2002;68(2) 167-74.
  40. 40. Zhou X, Acosta L, Willingham AL, Leonardo LR, Minggang C, Aligui G, et al. Regional Network for Research, Surveillance and Control of Asian Schistosomiasis (RNAS). Acta Trop. 2002;82(2) 305-11.
  41. 41. Zhou X, Dandan L, Huiming Y, Honggen C, Leping S, Guojing Y, et al. Use of landsat TM satellite surveillance data to measure the impact of the 1998 flood on snail intermediate host dispersal in the lower Yangtze River Basin. Acta Trop. 2002;82(2) 199-205.
  42. 42. McNally K. Developing Risk Assessment Maps ForSchistosoma HaematobiumBased On Climate Grids And Remotely Sensed Data: Louisiana State University; 2003.
  43. 43. Barbosa CS, Araújo KC, Antunes L, Favre T, Pieri OS. Spatial distribution of schistosomiasis foci on Itamaracá Island, Pernambuco, Brazil. Mem Inst Oswaldo Cruz. 2004;99(Suppl. I) 79-83.
  44. 44. Brooker S, Kabatereine NB, Clements ACA, Stothard JR. Schistosomiasis control. The Lancet. 2004;363 658-9.
  45. 45. Kabatereine NB, Brooker S, Tukahebwa EM, Kazibwe F, Onapa AW. Epidemiology and geography ofSchistosoma mansoniin Uganda: implications for planning control. Trop Med Int Health. 2004;9(3) 372-80.
  46. 46. Nihei N, Kajihara N, Kirinoki M, Chigusa Y, Saitoh Y, Shimamura R, et al. Fixed-point observation ofOncomelania nosophorain Kofu Basin--establishment of monitoring system of schistosomiasis japonica in Japan. Parasitology International. 2004;53 199-205.
  47. 47. Xu B, Gong P, Biging G, Liang S, Seto E, Spear R. Snail Density Prediction for Schistosomiasis Control Using Ikonos and ASTER Images. Photogram Eng Rem S. 2004;70(11) 1285-94.
  48. 48. Raso G, Vounatsou P, Singer BH, N'Goran EK, Tanner M, Utzinger J. An integrated approach for risk profiling and spatial prediction ofSchistosomamansoni– hookworm coinfection. PNAS. 2006;103(18) 6934-9.
  49. 49. Guo J-G, Penelope V, Cao C-L, Jurg U, Zhu H-Q, Daniel A, et al. A geographic information and remote sensing based model for prediction ofOncomelania hupensishabitats in the Poyang Lake area, China. Acta Tropica. 2005;96(2-3) 213-22.
  50. 50. Yang G-J, Vounatsou P, Zhou X-N, Tanner M, Utzinger J. A Bayesian-based approach for spatio-temporal modeling of county level prevalence ofSchistosoma japonicuminfection in Jiangsu province, China. International Journal for Parasitology. 2005;35(2) 155-62.
  51. 51. Zhang Z-Y, Xu D-Z, Zhou X-N, Zhou Y, Liu S-J. Remote sensing and spatial statistical analysis to predict the distribution ofOncomelania hupensisin the marshlands of China. Acta Trop. 2005;96(2-3) 205-12.
  52. 52. Brooker S, Alexander N, Geiger S, Moyeed RA, Stander J, Fleming F, et al. Contrasting patterns in the small-scale heterogeneity of human helminth infections in urban and rural environments in Brazil. International Journal for Parasitology. 2006;36(10-11) 1143-51.
  53. 53. Freitas CC, Guimarães RJPS, Dutra LV, Martins FT, Gouvea EJC, Santos RAT, et al. Remote Sensing and Geographic Information Systems for the Study of Schistosomiasis in the State of Minas Gerais, Brazil. Geoscience and Remote Sensing Symposium, 2006 IGARSS 2006 IEEE International Conference on; 2006 July 31 2006-Aug. 4 2006; 2006. p. 2436-9.
  54. 54. Gazzinelli A, LoVerde PT, Haddad JPA, Pereira WR, Bethony J, Correa-Oliveira R, et al. The spatial distribution ofSchistosoma mansoniinfection before and after chemotherapy in the Jequitinhonha Valley in Brazil. Mem Inst Oswaldo Cruz. 2006;101(Suppl. I) 63-71.
  55. 55. Guimarães RJPS, Freitas CC, Dutra LV, Moura ACM, Amaral RS, Drummond SC, et al. Analysis and estimative of schistosomiasis prevalence for Minas Gerais state, Brazil, using multiple regression with social and environmental spatial data. Mem Inst Oswaldo Cruz. 2006;101(Suppl. I) 91-6.
  56. 56. Tibiriçá SHC. Identificação morfológica e molecular, biometria, abundância e distribuição geográfica deBiomphalariaspp (Preston, 1910) (Mollusca, Planorbidae) no município de Juiz de Fora, Minas Gerais. MSc thesis.Universidade Federal de Juiz de Fora; 2006.
  57. 57. Araújo KCGM, Resendes APC, Souza-Santos R, Silveira Júnior JC, Barbosa CS. Análise espacial dos focos deBiomphalaria glabratae de casos humanos de esquistossomose mansônica em Porto de Galinhas, Pernambuco, Brasil, no ano 2000. Cadernos de Saúde Pública. 2007;23 409-17.
  58. 58. Brooker S. Spatial epidemiology of human schistosomiasis in Africa: risk models, transmission dynamics and control. Trans R Soc Trop Med Hyg. 2007;101(1) 1-8.
  59. 59. Matthys B, Tschannen AB, Tian-Bi NT, Comoé H, Diabaté S, Traoré M, et al. Risk factors forSchistosoma mansoniand hookworm in urban farming communities in western Côte d'Ivoire. Tropical Medicine & International Health. 2007;12(6) 709-23.
  60. 60. Moira AP, Fulford AJC, Kabatereine NB, Kazibwe F, Ouma JH, Dunne DW, et al. Microgeographical and tribal variations in water contact andSchistosoma mansoniexposure within a Ugandan fishing community. Tropical Medicine & International Health. 2007;12(6) 724-35.
  61. 61. Clements ACA, Brooker S, Nyandindi U, Fenwick A, Blair L. Bayesian spatial analysis of a national urinary schistosomiasis questionnaire to assist geographic targeting of schistosomiasis control in Tanzania, East Africa. International Journal for Parasitology. 2008;38(3-4) 401-15.
  62. 62. Guimarães RJPS, Freitas CC, Dutra LV, Moura ACM, Amaral RS, Drummond SC, et al. Schistosomiasis Risk Estimation in Minas Gerais State, Brazil, using Environmental Data and GIS techniques. Acta Trop. 2008;108 234-41.
  63. 63. Martins FT. Mapeamento do risco da esquistossomose no estado de Minas Gerais, usando dados ambientais e sociais. MSc thesis. INPE; 2008.
  64. 64. Yang K, Wang X-H, Yang G-J, Wu X-H, Qi Y-L, Li H-J, et al. An integrated approach to identify distribution ofOncomelania hupensis, the intermediate host ofSchistosoma japonicum, in a mountainous region in China. International Journal for Parasitology. 2008;38(8-9) 1007-16.
  65. 65. Zhang Z, Carpenter TE, Chen Y, Clark AB, Lynn HS, Peng W, et al. Identifying high-risk regions for schistosomiasis in Guichi, China: A spatial analysis. Acta Trop. 2008;107 217-23.
  66. 66. Clements ACA, Bosqué-Oliva E, Sacko M, Landouré A, Dembélé R, Traoré M, et al. A Comparative Study of the Spatial Distribution of Schistosomiasis in Mali in 1984-1989 and 2004-2006. PLoS Negl Trop Dis. 2009;3(5) e431.
  67. 67. Guimarães RJPS, Freitas CC, Dutra LV, Felgueiras CA, Moura ACM, Amaral RS, et al. Spatial distribution ofBiomphalariamollusks at São Francisco River Basin, Minas Gerais, Brazil, using geostatistical procedures. Acta Trop. 2009;109 181-6.
  68. 68. Guimarães RJPS, Freitas CC, Dutra LV, Felgueiras CA, Drummond SC, Oliveira G, et al. Use of Indicator Kriging to Investigate Schistosomiasis in Minas Gerais State, Brazil. J Trop Med. 2012;2012(Article ID 837428) 10.
  69. 69. Martins-Bedé FT, Godo L, Sandri S, Dutra LV, Freitas CC, Carvalho OS, et al. Classification of Schistosomiasis Prevalence Using Fuzzy Case-Based Reasoning. Bio-Inspired Systems: Computational and Ambient Intelligence: Berlin/Heidelberg, 2009. p8.
  70. 70. Zhang Z, Clark AB, Bivand R, Chen Y, Carpenter TE, Peng W, et al. Nonparametric spatial analysis to detect high-risk regions for schistosomiasis in Guichi, China. Transactions of the Royal Society of Tropical Medicine and Hygiene. 2009;103(10) 1045-52.
  71. 71. Carvalho OS, Scholte RGC, Guimarães RJPS, Freitas CC, Drummond SC, Amaral RS, et al. The Estrada Real project and endemic diseases: the case of schistosomiasis, geoprocessing and tourism. Mem Inst Oswaldo Cruz. 2010;105(4) 532-6.
  72. 72. Galvão AF, Favre TC, Guimarães RJPS, Pereira APB, Zani LC, Felipe KT, et al. Spatial distribution ofSchistosoma mansoniinfection before and after chemotherapy with two praziquantel doses in a community of Pernambuco, Brazil. Mem Inst Oswaldo Cruz. 2010;105(4) 555-62.
  73. 73. Guimarães RJPS, Freitas CC, Dutra LV, Scholte RGC, Amaral RS, Drummond SC, et al. Evaluation of a linear spectral mixture model and vegetation indices (NDVI and EVI) in a study of schistosomiasis mansoni andBiomphalaria glabratadistribution in the state of Minas Gerais, Brazil. Mem Inst Oswaldo Cruz. 2010;105(4) 512-8.
  74. 74. Guimarães RJPS, Freitas CC, Dutra LV, Scholte RGC, Martins-Bedé FT, Fonseca FR, et al. A geoprocessing approach for studying and controlling schistosomiasis in the state of Minas Gerais, Brazil. Mem Inst Oswaldo Cruz. 2010;105(4) 524-31.
  75. 75. Martins-Bedé FT, Dutra LV, Freitas CC, Guimarães RJPS, Amaral RS, Drummond SC, et al. Schistosomiasis risk mapping in the state of Minas Gerais, Brazil, using a decision tree approach, remote sensing data and sociological indicators. Mem Inst Oswaldo Cruz. 2010;105(4) 541-8.
  76. 76. Schur N, HürSchur N, Stensgaard A, Chimfwembe K, Mushinge G, Simoonga C, et al. Spatially explicit Schistosoma infection risk in eastern Africa using Bayesian geostatistical modeling. Acta Trop. 2011. In Press.
  77. 77. Stensgaard A-S, Utzinger J, Vounatsou P, Hürlimann E, Schur N, Saarnak CFL, et al. Large-scale determinants of intestinal schistosomiasis and intermediate host snail distribution across Africa: Does climate matter? Acta Trop. 2011. In Press.
  78. 78. Sturrock HJW, Gething PW, Ashton RA, Kolaczinski JH, Kabatereine NB, Brooker S. Planning schistosomiasis control: investigation of alternative sampling strategies forSchistosomamansonito target mass drug administration of praziquantel in East Africa. International Health.3(3) 165-75.
  79. 79. Tibiriça SHC, Mittherofhe A, Castro MF, Lima AC, Gonçalves M, Pinheiro IO, et al. Malacological survey ofBiomphalariasnails in municipalities along theEstrada Realin the southeast of the State of Minas Gerais, Brazil. Rev Soc Bras Med Trop. 2011;44(2) 163-7.
  80. 80. Neter J, Kutner MH, Nachtssheim CJ, Wasserman W. Applied linear statistical models. Boston: WCB/McGraw-Hill, 1996.
  81. 81. Assunção RM, Neves MC, Camara G, Freitas CC. Efficient regionalization techniques for socio-economic geographical units using minimum spanning trees. Int J Geogr Inf Sci. 2006;20 797-811.
  82. 82. INPE. SPRING V.4.1 Sistema de Processamento de Informações Georeferenciadas 2005.

Written By

Ricardo J.P.S. Guimarães, Corina C. Freitas, Luciano V. Dutra, Guilherme Oliveira and Omar S. Carvalho

Submitted: May 3rd, 2012 Reviewed: September 17th, 2012 Published: January 16th, 2013