Risk factor analysis and spatiotemporal CART model of cryptosporidiosis in Queensland, Australia
© Hu et al; licensee BioMed Central Ltd. 2010
Received: 14 September 2009
Accepted: 28 October 2010
Published: 28 October 2010
It remains unclear whether it is possible to develop a spatiotemporal epidemic prediction model for cryptosporidiosis disease. This paper examined the impact of social economic and weather factors on cryptosporidiosis and explored the possibility of developing such a model using social economic and weather data in Queensland, Australia.
Data on weather variables, notified cryptosporidiosis cases and social economic factors in Queensland were supplied by the Australian Bureau of Meteorology, Queensland Department of Health, and Australian Bureau of Statistics, respectively. Three-stage spatiotemporal classification and regression tree (CART) models were developed to examine the association between social economic and weather factors and monthly incidence of cryptosporidiosis in Queensland, Australia. The spatiotemporal CART model was used for predicting the outbreak of cryptosporidiosis in Queensland, Australia.
The results of the classification tree model (with incidence rates defined as binary presence/absence) showed that there was an 87% chance of an occurrence of cryptosporidiosis in a local government area (LGA) if the socio-economic index for the area (SEIFA) exceeded 1021, while the results of regression tree model (based on non-zero incidence rates) show when SEIFA was between 892 and 945, and temperature exceeded 32°C, the relative risk (RR) of cryptosporidiosis was 3.9 (mean morbidity: 390.6/100,000, standard deviation (SD): 310.5), compared to monthly average incidence of cryptosporidiosis. When SEIFA was less than 892 the RR of cryptosporidiosis was 4.3 (mean morbidity: 426.8/100,000, SD: 319.2). A prediction map for the cryptosporidiosis outbreak was made according to the outputs of spatiotemporal CART models.
The results of this study suggest that spatiotemporal CART models based on social economic and weather variables can be used for predicting the outbreak of cryptosporidiosis in Queensland, Australia.
Cryptosporidiosis is a diarrhoeal disease caused by microscopic parasites of the Cryptosporidium parvum . The parasite is one of the most common causes of waterborne disease in Australia and globally and is found in drinking water and recreational water . Cryptosporidiosis can also be transmitted via contaminated food, contact between people, or contact between people and animals. It is considered a drinking-waterborne disease because the largest outbreaks of cryptosporidiosis have been associated with contaminated drinking water. The individual risk factors for cryptosporidiosis in humans include drinking water from poorly treated public and private supplies, community swimming pools, day care centres and contact with farm animals. The seasonal occurrence of cryptosporidiosis has been cited as circumstantial evidence of a climatic link [3, 4].
Spatiotemporal analyses of disease have played a major role in environmental epidemiology. Logistic regression, polyclass and generalised linear models have been widely used to study relationships between diseases and various environmental risk factors. However, these statistical techniques also perform relatively poorly with high dimensionality in spatiotemporal context [5–7]. For example, incidence rates, climate and social economic factors often involve high-order interactions and collinearities, which are often difficult to cope with using these models . Moreover, infectious disease incidence data at finer spatial scales often have a preponderance of zero counts. In this case, the use of the conventionally employed standard statistical models may result in poor estimates and prediction [9–11]. An alternative approach for analysing such data with many zeros is to consider a suite of the spatiotemporal classification and regression tree (CART) models that take advantage of the natural categorisation of these types of populations into zero incidence, "normal" incidence and outbreaks. The CART models provide a non-parametric approach that can potentially better accommodancite these complex interactions since they avoid any assumption among the variables or homoscedasticity in variances . The integrated use of spatial statistics and CARTs at a variety of spatial scales has provided new insights into ecology and environmental epidemiology [13–15].
The aims of this paper are to examine the potential impact of social economic and weather factors by LGAs on the incidence of cryptosporidiosis using spatiotemporal CART models and explore their potential as a predictive model for cryptosporidiosis in Queensland, Australia.
Queensland, located in the northeast of Australia, 10-28° south latitude and 138-153° east longitude, is the second largest state (after Western Australia) with the largest habitable area in Australia. It occupies the north-eastern quarter of the continent and covers approximately 1,727,000 km2, with 7,400 km of mainland coastline (9,800 km including islands). It has typically sub-tropical climate characteristics with average temperatures of 25°C in summer and 15°C in winter. Rainfall varies regionally and seasonally, and most of the state receives over 50% of its rainfall during summer. Average rainfall varies from less than 150 mm in the southwest region to more than 4,000 mm on the far northern coast. Queensland consists of 125 statistical local government areas (LGA), with populations ranging in size from 312 to 888,449. LGA are widely used by researchers in Australia because of the availability in many datasets.
We obtained the computerised data set on the notified cryptosporidiosis cases by LGA in Queensland for the period of 1st January - 31st December 2001 from the Queensland Department of Health. Weather and socio-economic index for areas (SEIFA) data were obtained for the same period from the Australian Bureau of Meteorology and the Australian Bureau of Statistics, respectively. SEIFA is a continuum of advantage (high values) to disadvantage (low values) and takes into account variables relating to education, occupation, wealth and living conditions. Australian Bureau of Statistics does not provide SEIFA data by temporarily (month-to-month). SEIFA variables vary spatially but rarely from month-to-month during one year period (ie., 2001) and each LGA had an observation for the period of this study. Weather data comprised monthly mean maximum temperature (°C) and monthly rainfall (mm), which was an interpolated climate surface with a 0.25 × 0.25° grid resolution (about 30 × 30 km) at the equator. Climate grid cells were calibrated using the GIS software package Vertical Map to extract average pixel values of the weather variables for each LGA.
Spatial autocorrelation analysis
x i and x j denote the observed value at location i and j, is the average of the x values over the n locations, and w ij is the spatial weight measure . Moran's I can be interpreted as follows: a value close to 0 indicates randomness, while a positive (negative) value indicates positive (negative) autocorrelation.
Spatial empirical Bayes rates smoothing
We used a Bayesian random effects model to compute the spatial empirical Bayes rates smoothing and estimate the underlying distributions of incidence rates of cryptosporidiosis disease for each LGA and for Queensland overall. Spatial smoothing can be used to reduce random variation associated with small populations at each LGA. Spatial smoothing enables observations of gradients in cryptosporidiosis disease incidences that may not be apparent from the direct observation of raw incidences. Assume that the underlying rates for each LGA are drawn from a population of rates with prior distribution characterised by a mean θ and variance φ. The Bayesian posterior estimate for the underlying risk at the ith LGA is a weighted average of the raw rate p i and the prior, with weights inversely related to their variance. The equation is: , where . P i is the population at risk in area i . The smoothed incidences of cryptosporidiosis were computed from the total number of cryptosporidiosis cases in a spatial window divided by the total number of people at risk within the window. The spatial window was specified by a spatial weights file including each LGA and its neighbours. Thus the observed cryptosporidiosis rate for a LGA with a small population at risk will be adjusted considerably toward the average whereas for a larger LGA the rate will barely change.
The spatial empirical Bayes analyses were performed using GeoDa [16, 17]. The locations of cryptosporidiosis cases by LGA were geo-coded to the digital base maps of localities using MapInfo , GeoDa  and Microsoft Access software.
Spatiotemporal CART models
Since monthly cryptosporidiosis count data at the LGA level have a high incidence of zero counts (about 85% in the dataset considered here), the analysis requires special care to account for the extra variation unaccounted for with a typical Poisson assumption. Moreover, these data typically here "outbreaks" characterised as very high incidences that can be difficult to model with a single parameter Poisson distribution. The mass of extra zeros can be thought of as two subgroups: LGA which due to some characteristic could not have cryptosporidiosis disease during the study period and LGA that could have cryptosporidiosis disease but this did not appear or was not reported. We considered a suite of three spatiotemporal CART models: 1) fitting a tree to data categorised as binary: incidence/no incidence; 2) fitting a tree to just the incidences (ignoring all of the zeros). The first two models are equivalent to a zero-inflated mixture model comprising some probability of a zero response (no incidence) and some probability of a (positive) truncated Poisson model . Zero-inflated models are used when excess zeros result in a bimodal distribution which has a mixture of a mass of extra zeros and a mass that has a Poisson [10, 19]. 3) fitting a tree for outbreak of cryptosporidiosis
CART 1: fitting a tree to data categorised as binary: incidence/no incidence
Classification trees are used to predict membership of cases or objects in the classes of a categorical dependent variable based on one or more exploratory variables . A classification tree is built in accordance with a set of binary splitting rules that successively divide the data into subgroups with maximum homogeneity, where the latter is by some measure such as the change in the impurity . A classification tree model was developed to explore possible relationships between socio-economic and weather factors and monthly cryptosporidiosis incidence. A dichotomous outcome variable was defined as whether or not one or more cryptosporidiosis cases occurred in the LGA in one month. The spatial CART model was thus described as: monthly cryptosporidiosis presence/absence ~ SEIFA + maximum temperature + rainfall + factor (season). The results are reported as the probability of occurrence of cryptosporidiosis in an LGA in a month.
CART 2: fitting a tree to positive incidences
Like a classification tree, a regression tree is built by binary recursive partitioning, but the response variable is continuous . The splitting rules used in the algorithm are based on measures such as minimising the sum of the squared deviations from the mean in the two separate subgroups. The squared residuals minimisation algorithm is identical to a Gini splitting rule . Positive monthly cryptosporidiosis incidence rates (1/100,000) by LGA were used as a continuous response variable in this model. The spatial CART model was thus described as: Monthly cryptosporidiosis rate ~ SEIFA + maximum temperature + rainfall + factor (season). The results are reported as the relative risk (RR) of positive incidence of cryptosporidiosis, compared to average incidence rate of cryptosporidiosis. RR = (expected incidence - mean of incidence)/mean of incidence.
CART3: fitting a tree to outbreaks of cryptosporidiosis
We defined an outbreak if the incidence rate (without zeros) in any month exceeded the third quartile of the incidence rates per LGA in Queensland during the 12 months of 2001. The spatial CART model was thus described as: monthly cryptosporidiosis outbreak/non-outbreak ~ SEIFA + maximum temperature + rainfall + factor (season). The results are reported as the probability of occurrence of an outbreak cryptosporidiosis in a month. A map of predicted probabilities was developed using the outputs of the CART model.
Each CART analysis consisted of three basic steps. Firstly, a preliminary tree was grown by recursive data partitioning. Secondly, nested trees were formed by reducing the number of nodes in the tree (pruning). Thirdly, 10-fold cross-validation was used to address over-fitting and to identify the optimal tree with respect to its predictive ability. A minimum node deviance of 25% of the total deviance was used to prune the trees. To control for the impact of seasonality, we decomposed the cryptosporidiosis incidence into four seasonal categories (coded as Spring: September-October-November; Summer: December-January-February; Autumn: March-April-May; Winter: June-July-August) [20, 21]. Finally, we validated the model and tested the residual mean deviance and misclassification rates (predicted number of months without cryptosporidiosis occurrence/total number of months). The statistical analysis was conducted using the S-plus software package .
For each CART model, the predicted probabilities or relative risks can be presented as a map. A map of predicted probabilities of an outbreak was developed using the outputs of the CART 3 model.
Descriptive statistics of monthly cryptosporidiosis and social economic and weather variables by LGA in Queensland*
Cryptosporidiosis incidence (1/100,000)
Maximum temperature (°C)
Total rainfall (mm)
Matrix of correlation coefficients between cryptosporidiosis and social economic and weather factors*
There was a significant positive spatial autocorrelation of cryptosporidiosis incidence (including zero incidences) with a Moran's I statistic of 0.213 (p = 0.015).
Spatial distribution empirical Bayes rates smoothing
CART 1 model
CART 2 model
CART 3 model
These results of validation analysis reveal that the three CART models had reasonable accuracy, and its utility in research needs to be further explored.
The results of this study indicate that there was remarkable variation in the spatial distribution of cryptosporidiosis and suggest that temperature and SEIFA are statistically associated with the probability of occurrence and the magnitude of the incidence of the disease in Queensland, either directly or through other unmeasured variables.
Cryptosporidiosis is one of the common waterborne diseases globally. Therefore, from a public health perspective, there is a need to control and prevent this disease. There is evidence that education and behaviour might be direct risk factors for cryptosporidiosis infection . People with a lower SEIFA level often have poor health knowledge about disease transmission and personal protection from water-borne diseases. Thus, it is important to educate people in lower SEIFA areas about how to avoid contracting this infection especially in epidemic periods.
Cryptosporidium parvum thrive in warm waters of moderate salinity. Thus, with changing weather patterns, the geographic range of cryptosporidiosis may also change. For example, severe weather events appear to be correlated with outbreaks of cryptosporidiosis in hot summers in Brisbane and South Sydney, Australia [24, 25]. In the summer and autumn the cryptosporidiosis rate was positively associated with temperatures in the current and previous month in New Zealand . The autumn cryptosporidiosis peak has been linked to increased recreational water use, swimming and outdoor activities . Our research shows that high temperature (>32°C) has a significantly effect on cryptosporidiosis incidence in LGAs with high SEFIA values. This could be explained by the observation that people in these higher SEIFA classes may be more likely to undertake recreational activities such as camping and swimming in high temperature which may increase the chance of cryptosporidium . SEIFA allows ranking of regions/areas, providing a method of determining the level of social and economic well-being in each region. Although SEIFA is a useful overall index, two areas may have a similar SEIFA score for very different reasons. Thus, it may be important to require additional detailed information (eg., education and income etc) for particular analyses. Our previous research has shown that the proportion of residents with low educational attainment has a positive and significant association with cryptosporidiosis, especially in the regions with high SEIFA value .
Three-stage spatiotemporal CART models are based on simple splits of the data and thus do not require the assumptions such as the existence of linear regression among the variables . The visualized tree structure of the CART model makes the analysis results easier to understand and explain. The structure of the CART model is a set of nodes from the top to the bottom, in which the terminal nodes show the specific pattern features of the subpopulations (eg, the number of SLA, mean, risk, and probability etc) .
The suite of spatiotemporal CART models proposed have provide a comprehensive approach to modelling count data that have many zero counts and few very large counts. Because the ordinary CART do not take into account the extra variation, the test of regression parameters often fails to control Type I error rate. To overcome this problem, we developed three stage CART models. Our study shows that the three CART models can be applied in count data with extra zeros. The models can determine which predictors affect the probability of being an incidence of cryptosporidiosis, which predictors affect how positive incidence rates and which predictors affect the probability of an outbreak of cryptosporidiosis. The three models use the same social economic and weather predictors, but allow for different interactions and estimates of relative importance of the variables in describing the difference outcomes. Hence the predictors appear vastly different effects on the three models. They may not be the same predictors for the three models, or they could even have opposite effects on the three outcomes. It is thus advantageous to consider the set of models together when assessing the association between cryptosporidiosis disease and social economic and weather factors in the model.
The models proposed here can be extended in a number of lags. A lag of 0 month was considered because of biological plausibility (i.e., the incubation period for cryptosporidiosis ranges from 2 to 12 days). Longer lags (e.g., ≥ 1 month) were found to be less important in this study.
This study has four major strengths. Firstly, to our knowledge, this is the first spatiotemporal study examining the spatial variation of cryptosporidiosis disease and temperature, rainfall and SEIFA at a LGA level. Secondly, three spatiotemporal CART models were used to identify overall patterns of cryptosporidiosis transmission. Thirdly, detailed information on social economic and weather factors (temperature, rainfall and SEIFA) by LGA was incorporated in the statistical models, which may be helpful for further ecological research and development of early warning system. Finally, a map of predicted occurrence of cryptosporidiosis from this study may have important implications for public health decision-making in identifying high-risk communities to target for control and prevent cryptosporidiosis infection.
Limitations of this study should also be acknowledged. Firstly, the study suffers from the usual problems of an ecological design . Further researches need to take into account temporal change in long term period (ie., more than 5 years) during the ecological studies at fine spatial scales (ie., postal areas). Secondly, the mechanism and pathways underlying the association of cryptosporidiosis with SEIFA and temperature remain unclear . In this study, only the SEIFA, temperature and rainfall were considered. Other factors (eg., water reservoirs or large farms) may also be directly or indirectly related to cryptosporidiosis transmission. However these data were unavailable in this study. Finally, under-reporting is likely to occur when people infected by cryptosporidiosis have sub-clinical conditions and/or did not seek to see a doctor because they knew there is no effective treatment for this disease.
The identification of disease high risk areas provided the chance to explore social economic and weather factors that may be responsible for developing early warning systems of cryptosporidiosis. Early outbreak detection provides the chance to implement control measures and education campaigns. The effectiveness of public health interventions can be improved using spatiotemporal models to identify and monitor hot spots of cryptosporidiosis and identify major risk factors and the risk areas of cryptosporidiosis, and then targeting education campaigns and disease control activities at specific areas.
The findings of this study indicate that spatiotemporal CART models based on social, economic and weather variables can be used for predicting the outbreak of cryptosporidiosis in Queensland, Australia, and the research approach developed in this study may have wide applications in the surveillance and risk management of infectious diseases.
We thank the Queensland Department of Health, Australian Bureau of Meteorology, and Australian Bureau of Statistics for providing the data on notified cryptosporidiosis cases, climate, and population growth, respectively.
- Rose J, Huffman D, Gennaccaro A: Risk and control of waterborne cryptosporidiosis. FEMS Microbiol Rev. 2002, 26: 113-123. 10.1111/j.1574-6976.2002.tb00604.x.View ArticlePubMedGoogle Scholar
- Meinhardt P, Casemore D, Miller K: Epidemiologic aspects of human cryptosporidiosis and the role of waterborne transmission. Epidemiol Rev. 1996, 18: 118-136.View ArticlePubMedGoogle Scholar
- Casman E, Fischhoff B, Small M, Palmgren C, Wu F: An integrated risk model of a drinking-water-borne cryptosporidiosis outbreak. Risk Anal. 2001, 20: 493-509.Google Scholar
- Casman E, Fischhoff B, Small M, Dowlatabadi H, Rose J, Morgan G: Climate change and cryptosporidiosis: a qualitative analysis. Climate Change. 2001, 50: 219-249. 10.1023/A:1010623831501.View ArticleGoogle Scholar
- Schroder W: GIS, geostatistics, metadata banking, and tree-based models for data analysis and mapping in environmental monitoring and epidemiology. Int J Med Microbiol. 2006, 296: 23-36. 10.1016/j.ijmm.2006.02.015.View ArticlePubMedGoogle Scholar
- Munoz J, Felicisimo A: Comparison of statistical methods commonly used in predictive modelling. J Veget Sci. 2004, 15: 285-292. 10.1111/j.1654-1103.2004.tb02263.x.View ArticleGoogle Scholar
- De'ath G, Fabricius K: Classification and regression trees: a powerful yet simple technique for ecological data analysis. Ecology. 2000, 81: 3178-3192. 10.1890/0012-9658(2000)081[3178:CARTAP]2.0.CO;2.View ArticleGoogle Scholar
- Eckel S, Louis T: Identifying effect modifiers in air pollution time-series using a two-stage analysis. Dept of Biostatistics Working Papers Working Paper. 2007, 148: 1-29.Google Scholar
- Afifi A, Kotlerman J, Ettner S, Cowan M: Methods for improving regression analysis for skewed continuous or counted responses. Annu Rev Public Health. 2007, 28: 95-111. 10.1146/annurev.publhealth.28.082206.094100.View ArticlePubMedGoogle Scholar
- Cheung Y: Zero-inflated models for regression analysis of count data: a study of growth and development. Stat Med. 2002, 21: 1461-1469. 10.1002/sim.1088.View ArticlePubMedGoogle Scholar
- Hu W, Mengersen K, Fu SY, Tong S: The use of ZIP and CART to model cryptosporidiosis in relation to climatic variables. Int J Biometeorol. 2010, 54: 433-440. 10.1007/s00484-009-0294-4.View ArticlePubMedGoogle Scholar
- Breiman L, Fredman J, Olshen R, Stone C: Classification and regression trees. 1984, New York: Chapman & Hall (Wardworth, Inc)Google Scholar
- Bel L, Allard D, Laurent JM, Cheddadi R, Bar-Hen A: CART algorithm for spatial data: Application to environmental and ecological data. Computational Statistics & Data Analysis. 2009, 53: 3082-3093.View ArticleGoogle Scholar
- Cablk M, White D, Kiester A: Assessment of spatial autocorrelation in empirical models in ecology. Predicting species occurrences: issues of scale and accuracy. Edited by: Scott M, Heglund P, Morrison M, Rafael M, Wall B, Hoffer J. 2003, Washington, DC: IslandGoogle Scholar
- Hu W, Mengersen K, Tong S: Spatial analysis of notified cryptosporidiosis infections in Brisbane, Australia. Ann Epidemiol. 2009, 19: 900-907. 10.1016/j.annepidem.2009.06.004.View ArticlePubMedGoogle Scholar
- Anselin L: Exploring spatial data with GeoDa: a workbook. 2005, Urbana, USAGoogle Scholar
- Anselin L: GeoDa. 2005, Urbana, USAGoogle Scholar
- Corporation MapInfo: MapInfo professional software 7.0. 2003, New YorkGoogle Scholar
- Li Z, Knight S, Cook LJ, Hyde LK, Holubkov R, Olson LM: Modeling motor vehicle crashes for street racers using zero-inflated models. Accid Anal Prev. 2008, 40: 835-839. 10.1016/j.aap.2007.09.022.View ArticlePubMedGoogle Scholar
- Tabachnick B, Fidell L: Using multivariate statistics. 1996, New York: Harper Collins College PublishersGoogle Scholar
- Hu W, Tong S, Mengersen K, Oldenburg B: Exploratory spatial analysis of social and environmental factors associated with the incidence of Ross River virus in Brisbane, Australia. Am J Trop Med Hyg. 2007, 76: 814-819.PubMedGoogle Scholar
- S-Plus Insightful Corporation: S-Plus 6 for Windows. 2003, Seattle, WA, 6Google Scholar
- Hellard M, Hocking J, Willis J, Dore G, Fairley C: Risk factors leading to Cryptosporidium infection in men who have sex with men. Sex Transm Infect. 2003, 79: 412-414. 10.1136/sti.79.5.412.View ArticlePubMedPubMed CentralGoogle Scholar
- Stafford R, Neville G, Towner C, McCall B: A community outbreak of Cryptosporidium infection associated with a swimming pool complex. Commun Dis Intell. 2000, 24: 236-239.PubMedGoogle Scholar
- Lemmon J, McAnulty J, Bawden-Smith J: Outbreak of cryptosporidiosis linked to an indoor swimming pool. Med J Aust. 1996, 165: 613-616.PubMedGoogle Scholar
- Lake IR, Pearce J, Savill M: The seasonality of human cryptosporidiosis in New Zealand. Epidemiol Infect. 2008, 136: 1383-1387. 10.1017/S0950268807009922.View ArticlePubMedGoogle Scholar
- Naumova EN, Jagai JS, Matyas B, DeMaria A, MacNeill IB, Griffiths JK: Seasonality in six enterically transmitted diseases and ambient temperature. Epidemiol Infect. 2007, 135: 281-292. 10.1017/S0950268806006698.View ArticlePubMedGoogle Scholar
- Lake I, Harrison F, Chalmers R, Bentham G, Nichols G, Hunter P, Kovats R, Grundy C: Case-control study of environmental and social factors influencing cryptosporidiosis. Eur J Epidemiol. 2007, 22: 805-811. 10.1007/s10654-007-9179-1.View ArticlePubMedPubMed CentralGoogle Scholar
- Cressie N: Statistics for Spatial Data. 1991, New York: WileyGoogle Scholar
- The pre-publication history for this paper can be accessed here:http://www.biomedcentral.com/1471-2334/10/311/prepub
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.