 Research article
 Open Access
 Open Peer Review
 Published:
Oblique decision trees for spatial pattern detection: optimal algorithm and application to malaria risk
BMC Medical Research Methodologyvolume 5, Article number: 22 (2005)
Abstract
Background
In order to detect potential disease clusters where a putative source cannot be specified, classical procedures scan the geographical area with circular windows through a specified grid imposed to the map. However, the choice of the windows' shapes, sizes and centers is critical and different choices may not provide exactly the same results.
The aim of our work was to use an Oblique Decision Tree model (ODT) which provides potential clusters without prespecifying shapes, sizes or centers. For this purpose, we have developed an ODTalgorithm to find an oblique partition of the space defined by the geographic coordinates.
Methods
ODT is based on the classification and regression tree (CART). As CART finds out rectangular partitions of the covariate space, ODT provides oblique partitions maximizing the interclass variance of the independent variable. Since it is a NPHard problem in R^{N}, classical ODTalgorithms use evolutionary procedures or heuristics. We have developed an optimal ODTalgorithm in R^{2}, based on the directions defined by each couple of point locations. This partition provided potential clusters which can be tested with MonteCarlo inference.
We applied the ODTmodel to a dataset in order to identify potential high risk clusters of malaria in a village in Western Africa during the dry season. The ODT results were compared with those of the Kulldorff' s SaTScan™.
Results
The ODT procedure provided four classes of risk of infection. In the first high risk class 60%, 95% confidence interval (CI95%) [52.22–67.55], of the children was infected. MonteCarlo inference showed that the spatial pattern issued from the ODTmodel was significant (p < 0.0001).
Satscan results yielded one significant cluster where the risk of disease was high with an infectious rate of 54.21%, CI95% [47.51–60.75]. Obviously, his center was located within the first high risk ODT class. Both procedures provided similar results identifying a high risk cluster in the western part of the village where a mosquito breeding point was located.
Conclusion
ODTmodels improve the classical scanning procedures by detecting potential disease clusters independently of any specification of the shapes, sizes or centers of the clusters.
Background
Since the development of warning systems and environmental hazards awareness, a wide range of statistical methods has been provided to identify disease clusters and spatial patterns. These methods have been classified into three groups [1–3]:
Tests for focused clustering where the putative source is prespecified [4, 5, 2];
Tests for global clustering with statistics using the distance between cases [6–9];
General tests for localized clusters where the putative source or potential clusters cannot be prespecified [6, 10, 11].
This paper focuses on the latter tests i.e. on general procedures for the determination of spatial patterns. These patterns allow us to localize disease clusters where the disease rate is particularly high. Since the Openshaw's Geographical Analysis Machine (GAM), numerous works have proposed extensions or modifications of this method. The GAM lays out a regular grid of points covering the region under study. Then it generates overlapping circular windows centered at each grid point with constant radii depending on the grid spacing. The procedure is repeated at different predetermined values of the radius and thus defines potential clusters. Alternative procedures use circular windows centered at the observed point locations [10] and scan the area through this irregular grid. The use of squared shaped windows has also been proposed [6]. A general review of spatial methods is provided by Waller and Gotway [12] as well as in several publications [13–15].
The Kulldorff's scan statistic is one of the most interesting and used methods for cluster analysis [16, 1, 17]. The scan statistic is a likelihood ratio based method, which Kulldorff [11] defined without any assumptions about the shape, size or collection of locations for the scanning windows. However, various algorithms are necessary to calculate the test statistics for different defined types of scanning windows. Softwares (e.g. SaTScan™ [18]) have been implemented for some of these particular windows/algorithms. The SaTScan™ imposes on the map circular windows positioned on regular (such as GAM) or irregular grid (defined by the observed point locations). For each center point, the radius varies continuously from zero to a prespecified upper bound. Each of the circular windows, moving through the different centers and with different radii, is a possible candidate for containing a cluster of cases.
It is noteworthy that the detection of potential clusters is enforced on circular shaped (or squared shaped) windows. The various algorithms applied the scan statistic method to windows centered at either grid or observed point locations. These two procedures define different sets of potential clusters and therefore may not provide exactly the same results. Furthermore, changing the windows' shape may also provide different clusters. Gangnon and Clayton introduced a bayesian approach [19] for clustering which does not require cluster's locations or shapes to be specified but which requires some prior specifications of the distribution of various cluster size and shapes (hierarchical priors). However, given the large number of potential models, the posterior distribution cannot be directly provided. Therefore, Gangnon and Clayton limit the number of models under consideration by using a randomized method to build models with high posterior distributions. They approximate the posterior distribution over the limited number of cluster models incorporating hierarchical prior. Patil and Taillie [20] proposed an adaptation of the scan statistics to detect clusters without restricted shape. It reduces the size of the potential cluster set by determining levels of the rates of cases. The potential cluster set consists on all the connected components that have rates higher than a fixed level. Each level determines a potential cluster set. But the determination of levels is datadependent. Furthermore in a practical point of view, not all of the observed rates can be used as levels in order to avoid providing a computationally impracticable number of potential cluster sets. Other procedures use stochastic optimization algorithm to reduce the number of examined potential clusters [21]. But again these methods used for the determination of potential clusters are not optimal from a classification viewpoint.
The aim of the present work is to provide an optimal partitioning procedure using Oblique Decision Trees in order to detect spatial patterns and to optimize the potential clusters determination without prior specifications. Rather than using a likelihood ratio test, this new approach, which is not a scan statistic, is based on the calculus of the interclass variance during each of many splits of the space before providing the final pattern.
Methods
CART and ODTmodels
Treebased models such as CART (Classification And Regression Trees) [22] are nonlinear and nonparametric alternatives to linear models for regression and classification problems (such as linear regression, logistic regression, linear discriminant analysis, linear proportional hazard models). CART models are fitted by binary recursive partitioning of a multidimensional covariate space, in which the dataset is successively split into increasingly homogeneous subsets until a specified criterion is satisfied. For the first partition, CART searches the best possible place to split a continuous variable into two classes and defines two subspaces which maximize overall class separation (i.e. interclass variance of the dependent variable). Each of these subspaces subsequently serves as the basis for further partitioning independently of the others and so on. At each step the variable used for each split is selected from all the explicative variables so as to provide an optimal partition given the previous actions. Partitions' sequence is summarized by a binary tree. The root node of tree corresponds to the entire data space. Partitions of the space are associated with descendants of the root node. The leaves of the tree, or terminal nodes, correspond to subspaces which are not further partitioned. The stability of the procedure can be improved using Data resampling.
While CARTmodels are widely used as exploratory techniques they are lesscommonly used for prediction. Trees generally rely on fewer assumptions than classical methods and handle a wide variety of data structures. Furthermore they are easy to use and to interpretate, and thus provide a wide range of application fields. The use of CART procedure has been considered by others in a variety of medical problems [22, 23] such as, for example, survival analysis [24–26], longitudinal analysis, diagnostic and prognostic studies or clinical trials [27–30].
One particular application is signal processing [31], in which the problem concerns the detection of multiple change points in the mean. The CART procedure can be used to estimate simultaneously the changepoints and the means by recovering an underlying piecewise constant function f(t).
If m _{ k }are the means for each piecewise k, then t _{ k }are the changepoints:
If we extend this point of view to the covariate space defined by geographic coordinates, CART estimates the "changelines" (instead of changepoints) of a piecewise constant function on R^{2}. In other words, treebased procedure can easily determine spatial patterns.
However, one limitation is that CART provides axisparallel splits i.e. rectangular spatial patterns. Oblique decision trees (ODT) deal with this problem. Those algorithms produce oblique (and then polygonal) partitioning of the covariate space. However, oblique trees are less popular than axisparallel trees because the splits are less straightforward to interpret and oblique procedures require greater computational complexity than axisparallel algorithms. Finding the best oblique tree in the covariate space is a NPHard problem [32]. Therefore, existing ODT algorithms use deterministic heuristics or evolutionary algorithms (like the OC1 system [33]) to find appropriate hyperplanes for partitioning the covariate space [22, 34, 32, 33]. Comparisons of the different procedures are provided, for example, by Murthy [33] CantuPaz [34] and Bradley [35].
Despite this difficulty in R^{N}, it is easier to find an oblique partition in the particular case of a space determined by the geographic coordinates, i.e. in R^{2}. Evolutionnary or heuristic algorithms are not robust. They provide occasionally local minima [33] and therefore are not optimal procedures in R^{2}. The ODTalgorithm we have developped is an optimal procedure to reach an optimal solution without using heuristics or evolutionary procedures.
ODT algorithm
The general purpose of the entire procedure consists on finding several partitions of the plane. We present the first step which allows finding the best oblique split of the plane. Going recursively, this algorithm will split the plane into several partitions, until reaching a specific criterion.
This subsection is organized as follow:
i. First, we will introduce how the plane is splitted into two adjacent partitions according to the interclass variance.
ii. Second, we will present how the finite set of oblique lines is determined, still within the first step of the entire procedure.
iii. Third, we will propose an optimization of this algorithm.
i. The splitting method proceeds as follows.
Consider, in the geographical space represented by the plane with an orthogonal basis {x, y} and a fixed origin O, npoints M _{ i }with coordinates {x _{ i }, y _{ i }}. These coordinates can represente the geographic coordinates of a point location provided by GPS. To each point M _{ i }a numeric random variable Z _{ i }(called explained or dependant variable) is associated with the observation z _{ i }Whereas the CART procedure partitions the plane according to a line parallel to the axis maximizing the interclass variance of z _{ i }, our procedure partitions the plane according to an oblique line maximizing in the same way the interclass variance of z _{ i }.
To find this oblique line according to the direction
we have to define the perpendicular direction uand the angle .
From a general viewpoint, for a fixed direction
the procedure has to:
Orthogonally projects the points M _{ i }on the (O, u) direction, defining the coordinate u _{ i };
Considers all the u _{ i }as potential threshold in the way to split the plane with the direction perpendicular to the direction uand going through u _{ i };
Finds the optimal split between two adjacent classes, maximizing the interclass variance of z _{ i }according to theses projections.
ii. The splitting method provides a finite set of cluster proceeding as follows.
Before detailing the algorithm, we have to study the different splitting directions
i.e. to specify wich angles θ have to be analyzed. For a global solution the algorithm can scan all the oblique directions (i.e. all the θ) between zero and π. In a heuristic way one can also discretize this interval providing a finite number of angles θ. But these two methods are not optimal.
The optimal algorithm for an optimal solution is easy to implement. Obviously, two points M _{ i }(x _{ i }, y _{ i }) and M _{ j }(x _{ j }, y _{ j }) have the same projected coordinates on the (O, u) direction if and only if M _{ i } M _{ j }is perpendicular to (O, u) (Figure 1). Then the number of critical directions, defined by the θ _{ ij }angles, exists and is finite.
For each direction
passing through two points M _{ i }(x _{ i }, y _{ i }) and M _{ j }(x _{ j }, y _{ j }), we define φ _{ ij }the angle between the line M _{ i } M _{ j }and the xaxis.
Then:
As previously defined, θ is the angle between the xaxis and the direction (O, u) perpendicular to M _{ i } M _{ j }.
Then for each couple (M _{ i }, M _{ j }), we have
Each critical angle θ _{ ij }defines an angular sector. Within each sector, the order of the coordinates projected on the (O, u) direction does not depend on this direction. For points M _{ i }and M _{ j }the difference (u _{ j } u _{ i }) of their coordinates projected on (O, u) verifies:
(u _{ j } u _{ i }) cos(φ _{ ij }) = (x _{ j } x _{ i }) sin(θ  θ _{ ij }) (1)
with (u _{ j } u _{ i }) = (y _{ j } y _{ i }) sin(θ) for x _{ i }= x _{ j }⇔ φ _{ ij }=
Thus (u _{ j } u _{ i }) depends continuously on θ. The sign of this difference cannot change within the angular sector since (u _{ j } u _{ i }) = 0 only if θ = θ _{ ij }.
It follows that the interclass variances (and then the ODT procedure) is not modified within each sector. As a direct consequence of (1) the transition from a sector to another via the critical angle θ _{ ij }(Figure 2) induces the same order except the permutation of the two adjacent elements (u _{ i }, u _{ j }).
Note that for aligned points M _{ i }, M _{ j }and M _{ k }the algorithm has to permute the adjacent element group (u _{ i }, u _{ j }, u _{ k }). Similarly for parallel directions M _{ i } M _{ j }and M _{ k } M _{ l }, the algorithm has to permute at the same time the couples of adjacent elements (u _{ i }, u _{ j }) and (u _{ k }, u _{ l }).
Note again that all these angular sectors define as much covariates. Thus the procedure comes to the usual CART procedure. But the number of different critical angles is
and using CART this way overconsumes time and space. For example, in our application the number of point locations is n= 164, hence the number of different angular sectors is N= 13270.
iii. We present now an optimization of our algorithm. A less time consuming and more efficient algorithm is a stepwise analysis of the angular sector, ordered according to the observed θ _{ ij }. At each step the algorithm uses the previous calculus.
Because only two elements between two adjacent sectors are permuted only one interclass variance has to be reloaded, related to the single different split (or some interclass variances for the group of permuted couples, related to a few different splits). The procedure inherits the calculus of the other interclass variances from the previous sector with the exception of the interclass variance related to the single permutation.
Thus, the algorithm complexity is
(n ^{2} log n) in time and (n) in space for one split. Finally, our algorithm splits the plane into two adjacent partitions as follows:
Arrange the x _{ i };
Calculate and arrange the θ _{ ij }via the a _{ ij };
Calculate
;
For each potential split of the first angular sector (corresponding to the xaxis), i.e. for each value of x _{ i }:
Calculate the ∑ z _{ i }for each class (on both sides of the threshold x _{ i }) and then the interclass variance, using the previous results;
If the calculated interclass variance is greater than the previous one, store the results;
For the next angular sector
Permute the corresponding x _{ i } x _{ j }(or the group of elements);
Calculate the ∑ z _{ i }only for the two classes generated by the split between x _{ j }and x _{ i }(or some splits for the group of permuted elements);
If the new interclass variance is greater than the previous optimum, store the results;
Until all sectors are scanned.
This algorithm goes on recursively until a specific criterion is reached and the Oblique Decision Tree is completed.
For simplicity we will not herein discuss special procedures of CART such as stopping rules, pruning algorithms or resampling methods; these are examined elsewhere [22, 36].
Dataset
Malaria is the major parasitic disease in the world affecting approximately 300–500 million individuals annually. About two percents of the individuals infected with Plasmodium falciparum die. Most of the deaths occur in children. In the last decade, the incidence of malaria has been increasing at an alarming rate in Africa representing over 90% of the reported cases in the world [37].
The study area was the whole village of Bancoumana located in the administrative circle of Kati (Mali, Western Africa). This village is located in the high Niger's valley, a Sudanese savannah area, about 60 km southwest from the capital city Bamako. The main activities are rice cultures and truck farming along the Niger river. This village is 2.5 km^{2} wide, with 8 000 inhabitants (MRTC census, 1998) and about 1 600 children under 9 years. The transmission of malaria is high during the rain season (usually from June to October, with temperatures varying between 25°C and 40°C). It decreases then, reaching a low level of transmission one or two months thereafter.
The project investigated at a villagelevel approach (using a 1–3 m resolution scale) the risk of malaria infection. The presence of P. falciparum, the main infectious agent of malaria in this area, in blood smears was investigated in 1 461 children living in 164 households during the dry season in March 2000. Among them, 474 children had a positive blood smear (32.44%, CI95% [30.09–34.89]). Localization was performed through GPS receivers. Thus, all children were geocoded at a point location (corresponding to their house). Geodatabase and cartographic displays were provided with the ArcGIS 8.3 software (ESRI, Redlands, CA).
Human subjects' research conducted in these studies was approved by the Institutional Committee on Ethics of the Mali Faculty of Medicine and Pharmacy, University of Mali. To obtain informed consents a stepwise consent process was applied as described by Doumbo [38]. First, the community informed consent was obtained before the beginning of the study. Second, the informed consent of the parents or guardians of the children were orally obtained before each clinical or biological investigation.
Data analysis
The ODTalgorithm was implemented with the Matlab Software 7.0.1 (The Mathworks Inc. 2004). We applied the ODT procedure to the dataset using the GPS coordinates of each location as independent covariates and the parasitic positivity rate (rate of positive blood smears per houses) as dependant variable. Thus ODT provided an optimal partition of the geographical area, i.e. a spatial pattern of the disease risk. We chose to use two classical stopping rules [22]. First, the ODT stopped if a class was made up of less than 15 locations. Second, we prunned a node if, after partition, one of the two resulting classes was made up of less than 3 locations.
For inference we considered the constant risk hypothesis as a model of "no pattern". Under this null hypothesis each child is at the same disease risk within the observation period regardless of his location. Thus the classes issued from the ODT displayed similar disease risk. However in keeping with many spatial health applications [12], we can not rely on asymptotic arguments to derive theoretically the associated distributions under the null hypothesis. Monte Carlo (MC) simulations were flexible tools for such assessment.
Similarly to many statistical models, we used for inference the explained variability rate R_{v}, defined as the ratio of the interclass sum of squared errors (SCE) (outcome of the ODT model) and the total SCE. We considered a Monte Carlo inference conditional on the set of all locations and on the local number of subjects. The total number of cases varied from simulation to simulation with an expected value (the total number of cases on the observed dataset). In this way, the simulations assessed spatial variations in the local proportion of cases conditional on the set of all locations. Monte Carlo simulations reflected a constant risk hypothesis similarly to the Rushton and Lolonis [39] approach. We ran 999 simulations under the constant risk hypothesis i.e. homogeneous Poisson distribution. Under this null hypothesis we applied the ODTalgorithm for each of the random dataset and calculated the empirical distribution of R_{v}. Thus the MC inference provided pvalues for testing whether or not the observed explained variability rate is a realization of the theoretical (simulated) distribution under the constant risk hypothesis. In other words, MC inference tested the ODTmodel and provided the significance of the spatial pattern issued from the oblique decision tree.
We compared the ODTmodel outputs with those of the scan statistic method. For the latter, we used the software program SaTScan™ [18] in order to test for the presence of spatial clusters of malaria infection and to estimate their locations. The identification of high risk clusters with the SaTScan™ was performed under the Poisson probability model assumption using a maximal cluster size of 50% of the total population. For statistical inference, 999 Monte Carlo replications were performed. The null hypothesis of no clustering was rejected when the simulated pvalue was lower than or equal to 0.05.
During the data analysis we calculated all confidence intervals of rates according to the Wilson method [40].
Results
Oblique Decision Tree
The ODT (Figure 3) partitioned the village into four risk classes. The explained variability rate is high, i.e. R_{v} = 83.96% of the variability is explained by the ODTmodel. The global risk of disease (Table 1) was 32.44%, CI95% [30.09–34.89]. The ODT provided two classes of high infection risk. In the first high risk class (P2), located in the western part of the village (Figure 5), the risk was 60%, CI95% [52.22–67.55]. In the second high risk class (P3), located in the southern part of the village, the risk was about 50% with a large confidence interval. Note that during the rain season about 80% of the children had a positive blood smear in the whole village. Investigations at this site pointed to a small pond located within the western high risk class, and to ricefields located in the southern part of the village, both having been identified as Anopheles (the vector of malaria) breeding places.
Monte Carlo inference provided a global test, testing the null hypothesis of a homogeneous Poisson distribution of the malaria infection cases within the study area. Under this null hypothesis we provided (999 simulated sets and one observed set) the empirical distribution of the explicated variability rate R_{v} (Figure 4). In this application the R_{v} provided by the ODTmodel significantly differed (p < 0.0001) from the one provided under the homogeneous Poisson distribution, i.e. the spatial pattern was significant.
Satscan approach
The Satscan results yielded one significant clusters (Table 2). In the first cluster (S1) the risk of disease was high with an infectious rate of 54.21% (CI95% [47.51–60.75]). Obvioulsy, his center was located within the high risk ODT partition (P2) and the risk of disease were similar in S1 and P2 (Figure 5). The second and third high clusters were not significant, totalizing only one point location each.
Discussion
For spatial cluster detection, the specification of the shape and size of the clusters is required rather than using political or administrative definitions of zones. For this purpose scanning methods provide sets of potential clusters but the problem of the choice of the shape still remains. Different scanning grid and different windows' shapes or sizes may provide different sets of potential clusters. To reduce this difficulty we introduced ODTmodels with the aim to detect spatial pattern without prespecifying windows' shape. In contrast to classical scanning procedures, neither the shape, size, nor centroid location have to be specified by the users. Thus, ODT are optimal procedures from the classification viewpoint.
Furthermore the spatial pattern obtained by the ODTmodel defines a potential clusters set which can then be tested using the classical Monte Carlo inference. Similarly to Satscan, inference analysis has to avoid multiple testing inherent to such a procedure. The Kulldorffs procedure provides first a potential cluster set. Second, this procedure performs a significance test based on the local likelihood ratio statistic for each cluster in a way that compensates for the multiple testing. In our work we provide a global inference, testing the significativity of the spatial pattern obtained by ODT. Note that, similarly to Kulldorff's inference, likelihood ratio tests can be used to test the spatial pattern.
Recently, Tango proposed a flexibly spatial scan statistic to detect noncircular clusters [41]. But the Tango's method is not practically feasible for large clusters (more than 30 point locations).
Our findings indicate that the ODTmethod is consistent with the classical Kulldorf's scan statistic. ODT procedure is thus a classification tool widely usable for spatial pattern detection. When compared to ODT, the scan statistic did not detect the second high risk cluster (P3). This is probably due to the lack of points fitting in this cluster (3 point locations and 26 childrens). The 95% confidence interval of the disease rate in this cluster is large (32.06%67.94%). Nevertheless, after investigations in the village, a putative source of disease risk has been detected at this location. The two nonsignificant high risk clusters (S2 and S3) enclose only one point location each. This might explain both the lack of significativity with the satscan method and the lack of detection by the ODT.
After detection of a significant spatial pattern, the next logical step is to test whether this pattern can be explained by known or suspected risk factors. For example in the context of malaria, environmental factors such as mosquitoes breeding sites or thatched habitations might be identified and appropriate measures can be proposed to enhance the disease's control policy.
ODTmodels allow for a flexible relationship between the variables. The relationships between covariates do not need to be linear or additive and the interactions do not need to be prespecified or to be of a particular multiplicative form. The literature about treebased models is increasing particularly for studies focusing on formal inference procedures [31, 36]. In contrast to classical ODTprocedures the algorithm herein described is optimal since it uses neither evolutionary algorithms [32, 34] nor heuristics [22]. While the problem is NPHard in R^{N}, the algorithm remains polynomial in R^{2}.
The stability of treemodels can be improved by resampling methods. It is noteworthy that scanning methods such as satscan can also benefit from resampling. Resampling methods may improve the determination of the potential cluster set when the SaTScan™ procedure uses windows centered at each point location. Among different stopping rules, criteria have to be chosen according to the tradeoff between variance and bias of prediction. The usually chosen rules are known as flexible and robust methods [31]. But as our application results indicate, less restrictive rule can be used for specific epidemiological dataset in order to improve the interpretation of the ODTmodels' output. For rare diseases it might be necessary to use less stringent stopping rules than for diseases characterized by an epidemic evolution. This is related to the definition of "cluster of cases" which depends on the epidemiological profile of the disease.
The risk of infection has a high geographic variability [42, 43] and the knowledge of this variability is essential to enhence malaria control programs' efficiency [44]. Moreover, the detection of highrisk locations is one recommendation of the 20^{th} WHO technical report [45]. In this context, the development of GIS displays data on local malaria cases and then stratification of malaria risk providing the opportunity for more focal (and then efficient) malaria control programs [42].
Conclusion
In conclusion, Oblique Decision Tree is a new approach for spatial pattern detection and has the following features:
ODT improve the classical scanning procedures by providing polygonal potential clusters;
ODT are not bound by fixed centroid locations, sizes or shapes. Thus, first they have an enhanced flexibility. Second, the results are independent from the shape size and center prespecification;
ODT provide an optimal partition in the classification viewpoint.
Thus, ODTmodels favorably compare with other cluster detection methods for spatial epidemiology.
Abbreviations
 CART:

Classification And Regression Tree
 CI95%:

95% Confidence Interval
 GIS:

Geographic Information System
 GPS:

Global Positioning System
 MC:

Monte Carlo
 MRTC:

Malaria Research and Training Center, Bamako, Mali.
 ODT:

Oblique Decision Tree
 SCE:

sum of squared errors
References
 1.
Kulldorff M, Feuer EJ, Miller BA, Freeman LS: Breast cancer in northeastern United States: a geographical analysis. Am J Epidemiol. 1997, 146: 161170.
 2.
Bithell JF: The choice of test for detecting raised disease risk near a point source. Stat Med. 1995, 14: 23092322.
 3.
Cuzick J, Edwards R: Spatial clustering for inhomogeneous populations. J R Stat Soc [Ser B]. 1990, 52: 73104.
 4.
Tango T: A class of tests for detecting 'general' and 'focused' clustering of rare diseases. Stat Med. 1995, 14: 23232334.
 5.
Diggle PJ, Morris S, Elliott P, Shaddick G: Regression modelling of disease risk in relation to point sources. J R Stat Soc [Ser A]. 1997, 160: 491505. 10.1111/1467985X.00076.
 6.
Anderson NH, Titterington DM: Some methods for investigating spatial clustering, with epidemiological applications. J R Stat Soc [Ser A]. 1997, 160: 87105. 10.1111/1467985X.00047.
 7.
Tango T: Score tests for detecting excess risks around putative sources. Stat Med. 2002, 21: 497514. 10.1002/sim.1003.
 8.
Diggle PJ, Chetwynd AG: Secondorder analysis of spatial clustering for inhomogeneous populations. Biometrics. 1991, 47: 11551163.
 9.
GomezRubio V, ferrandiz J, Lopez A: Detecting clusters of diseases with R. Proceedings of the 3rd International Workshop on Distributed Statistical Computing: March 20–22 2003; Vienna Austria. Edited by: Kurt Hornik, Friedrich Leisch, Achim Zeileis. 2003, [http://www.ci.tuwien.ac.at/Conferences/DSC2003/]
 10.
Turnbull BW, Iwano EJ, Burnett WS, Howe HL, Clark LC: Monitoring for clusters of disease: application to leukemia incidence in upstate New York. Am J Epidemiol. 1990, 132: S136143.
 11.
Kulldorff M: A spatial scan statistic. Commun Stat Theor M. 1997, 26: 14811496.
 12.
Waller LA, Gotway CA: Applied spatial statistics for public health data. 2004, Wiley: Hoboken New Jersey
 13.
Wakefield J, Elliott P: Issues in the statistical analysis of small area health data. Stat Med. 1999, 18: 23772399. 10.1002/(SICI)10970258(19990915/30)18:17/18<2377::AIDSIM263>3.3.CO;27.
 14.
Kulldorff M, Nargawalla N: Spatial disease clusters: detection and inference. Stat Med. 1995, 14: 799810.
 15.
Thomas AJ, Carlin BP: Late detection of breast and colorectal cancer in Minnesota counties: an application of spatial smoothing and clustering. Stat Med. 2003, 22: 113127. 10.1002/sim.1215.
 16.
Sheehan TJ, De Chello LM, Kulldorff M, Gregorio DI, Gershman S, Mroszczyk M: The geographic distribution of breast cancer incidence in Massachusetts 1988 to 1997, adjusted for covariates. Int J Health Geogr. 2004, 3: 1710.1186/1476072X317.
 17.
Hjalmars U, Kulldorff M, Gustafsson G, Nagarwall N: Childhood leukemia in Sweden: using GIS and spatial scan statistic for cluster detection. Stat Med. 1996, 15: 707715. 10.1002/(SICI)10970258(19960415)15:7/9<707::AIDSIM242>3.3.CO;2W.
 18.
Kulldorff M: SaTScanTM v5.lSoftware for the spatial and spacetime scan statistics. 2004, Information Management Services Inc., Silver Spring, Maryland, [http://www.satscan.org]
 19.
Gangnon RE, Clayton MK: Bayesian detection and modeling of spatial disease clustering. Biometrics. 2000, 56: 922935. 10.1111/j.0006341X.2000.00922.x.
 20.
Patil GP, Taillie C: Upper level set scan statistic for detecting arbitrarily shaped hotspots. Environ Ecol Stat. 2004, 11: 183197. 10.1023/B:EEST.0000027208.48919.7e.
 21.
Duczmal L, Assunciao RM: A simulated annealing strategy for the detection of arbitrarily shaped spatial clusters. Comput Statist Data Anal. 2004, 45: 269286. 10.1016/S01679473(02)00302X.
 22.
Breiman L, Friedman JH, Olshen RA, Stone CJ: Classification and regression trees. 1993, Chapman & Hall: New York
 23.
Segal MR, Tager IB: Trees and tracking. Stat Med. 1993, 12: 21532168.
 24.
Xu R, Adak S: Survival analysis with timevarying regression effects using a treebased approach. Biometrics. 2002, 58: 305315. 10.1111/j.0006341X.2002.00305.x.
 25.
Leblanc M, Crowley J: Relative Risk trees for censored survival data. Biometrics. 1992, 48: 411425.
 26.
Schmoor C, Ulm K, Schumacher M: Comparison of the Cox model and the regression tree procedure in analyzing a randomized clinical trial. Stat Med. 1993, 12: 23512366.
 27.
Zhang H, Holford T, Bracken MB: A treebased method of analysis for prospective studies. Stat Med. 1996, 15: 3749. 10.1002/(SICI)10970258(19960115)15:1<37::AIDSIM144>3.3.CO;2S.
 28.
Crichton NJ, Hinde JP, Marchini J: Models for diagnosing chest pain: is cart helpful?. Stat Med. 1997, 16: 717727. 10.1002/(SICI)10970258(19970415)16:7<717::AIDSIM504>3.0.CO;2E.
 29.
Fu CY: Combining loglinear model with classification and regression tree (CART): an application to birth data. Comput Statist Data Anal. 2004, 45: 865874. 10.1016/S01679473(03)000926.
 30.
McBride WJH, Mullner H, Muller R, Labrooy J, Wronski I: Determinants of dengue 2 infection among residents of charters towers, Queensland, Australia. Am J Epidemiol. 1998, 148: 11111116.
 31.
Gey S: Bornes de risque, détection de ruptures boosting: trois thèmes statistiques autour de CART en régression. PhD thesis. 2002, University of Paris XI
 32.
Heath D, Kasif M, Salzberg S: Induction of oblique decision trees. Proceedings of the 13th International Joint Conference on Artificial Intelligence: August 28September 3 1993; Chambery, France. Edited by: Ruzena Bajcsy. 1993, Morgan Kaufmann, 10021007.
 33.
Murthy SK, Kasif M, Salzberg S: A system for induction of oblique decision trees. J Artif Intell Res. 1994, 2: 132.
 34.
CantuPaz E, Kamath C: Inducing oblique decision trees with evolutionary algorithms. IEEE Trans Evol Comput. 2003, 7: 5468. 10.1109/TEVC.2002.806857.
 35.
Brodley CE, Utgoff PE: Multivatiate decision trees. COINS technical reports 9282. 1992, University of Massachusetts
 36.
Ghattas B: Agrégation d'arbres de décision binaires: application à la prévision de l'ozone dans les Bouches du Rhône. PhD thesis. 2001, University of AixMarseille II
 37.
Tanser F, Le Sueur D: The application of geographical information systems to important public health problems in Africa. Int J Health Geogr. 2002, 1: 1
 38.
Doumbo OK: It takes a village: medical research and ethics in Mali. Science. 2005, 307: 679681. 10.1126/science.1109773.
 39.
Rushton G, Lolonis P: Exploratory spatial analysis of birth defect rates in an urban population. Stat Med. 1996, 15: 717726. 10.1002/(SICI)10970258(19960415)15:7/9<717::AIDSIM243>3.0.CO;20.
 40.
Newcombe RG: Twosided confidence intervals for the single proportion: comparison of seven methods. Stat Med. 1998, 17: 857872. 10.1002/(SICI)10970258(19980430)17:8<857::AIDSIM777>3.0.CO;2E.
 41.
Tango T, Takahashi K: A flexibly shaped spatial scan statistic for detecting clusters. Int J Health Geogr. 2005, 4: 1110.1186/1476072X411.
 42.
Booman M, Durrheim DN, La Grange K, Martin C, Mabuza AM, Zitha A, Mbokazi FM, Fraser C, Sharp BL: Using a geographical information system to plan a malaria control programme in South Africa. Bull World Health Organ. 2000, 78: 14381444.
 43.
Baird JK, Agyei SO, Utz GC, Koram K, Barcus MJ, Jones TR, Fryauff DJ, Binka FN, Hoffman SL, Nkrumah FN: Seasonal malaria attack rates in infants and young children in Northern ghana. Am J Trop Med Hyg. 2002, 66: 280286.
 44.
Killeen GF, Seyoum A, Knols BGJ: Rationalizing historical successes of malaria control in africa in terms of mosquito resource availability management. Am J Trop Med Hyg. 2004, 71 (S2): 8793.
 45.
World Health Organization: Expert Committee on Malaria 20th report. World Health Organ Tech Rep. 2000, 735
Prepublication history
The prepublication history for this paper can be accessed here:http://www.biomedcentral.com/14712288/5/22/prepub
Acknowledgements
This work was supported by the French Research Ministry (PAL+ 2001 Program) and the NIH Grant MaliTulane TMRC N° AI 95002P50
We thank Dr Bernard Fichet for many valuable discussions, Dr Bernard Giusiano for program improvements, and Dr Joanny Gouvernet for his helpful comments.
Author information
Additional information
Competing interests
The author(s) declare that they have no competing interests.
Authors' contributions
J Gaudart provided the ODT procedure, implemented the algorithm, performed statistical and geographical analysis, and drafted the manuscript.
B Poudiougou, S Ranque and O Doumbo carried out the epidemiological georeferenced data and provided the parasitological analysis.
O Doumbo initiated and supervised the epidemiological study.
All authors wrote and approved the final manuscript.
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
Rights and permissions
About this article
Received
Accepted
Published
DOI
Keywords
 Malaria
 Potential Cluster
 Angular Sector
 Circular Window
 High Risk Class