Geographical clustering of lung cancer in the province of Lecce, Italy: 1992–2001
 Massimo Bilancia^{1}Email author and
 Alessandro Fedespina^{2}
https://doi.org/10.1186/1476072X840
© Bilancia and Fedespina; licensee BioMed Central Ltd. 2009
Received: 27 November 2008
Accepted: 01 July 2009
Published: 01 July 2009
Abstract
Background
The triennial mortality rates for lung cancer in the two decades 1981–2001 in the province of Lecce, Italy, are significantly higher than those for the entire region of Apulia (to which the Province of Lecce belongs) and the national reference rates. Moreover, analyzing the rates in the threeyear periods 1993–95, 1996–98 and 1999–01, there is a dramatic increase in mortality for both males and females, which still remains essentially unexplained: to understand the extent of this phenomenon, it is worth noting that the standardized mortality rate for males in 1999–01 is equal to 13.92 per 10000 personyears, compared to a value of 6.96 for Italy in the 2000–2002 period.
These data have generated a considerable concern in the press and public opinion, which with little scientific reasoning have sometimes identified suspected culprits of the risk excess (for example, the emission caused by a number of large industrial sites located in the provinces of Brindisi and Taranto, bordering the Province of Lecce). The objective of this paper is to study on a scientifically sound basis the spatial distribution of risk for lung cancer mortality in the province of Lecce. Our goal is to demonstrate that most of the previous explanations are not supported by data: to this end, we will follow a hybrid approach that combines both frequentist and Bayesian disease mapping methods. Furthermore, we define a new sequential algorithm based on a modified version of the BesagYorkMollié (BYM) model, suitably modified to detect geographical clusters of disease.
Results
Standardized mortality ratios (SMRs) for lung cancer in the province of Lecce: For males, the relative risk (measured by means of SMR, i.e. the ratio between observed and expected cases in each area under internal standardization) was judged to be significantly greater than 1 in many municipal areas, the significance being evaluated under the null hypothesis of neutral risk on the ground of areaspecific pvalues (denoted by ρ_{ i }); in addition, it was seen that high risk areas were not randomly distributed within the province, but showed a sharp clustering. The most perceptible cluster involved a collection of municipalities around the Maglie area (Istat code: 75039), while the association among the municipalities of Otranto, Poggiardo and Santa Cesarea Terme (Istat codes: 75057, 75061, 75072) was more ambiguous. For females, it was noteworthy the significant risk excess in the city of Lecce (Istat code: 75035), where an SMR of 1.83 and ρ_{ i }< 0.01 have been registered. BYM model for the province of Lecce: For males, Bayes estimates of relative risks varied around an overall mean of 1.04 with standard deviation of 0.1, with a minimum of 0.77 and a maximum of 1.25. The posterior relative risks for females, although smoothed, showed more variation than for males, ranging form 0.74 to 1.65, around a mean of 0.90 with standard deviation 0.12. For males, 95% posterior credible intervals of relative risks included unity in every area, whereas significantly elevated risk of mortality was confirmed in the Lecce area for females (95% posterior CI: 1.33 – 2.00). BYM model for the whole Apulia: For males, internally standardized maps showed several high risk areas bordering the province of Lecce, belonging to the province of Brindisi, and the presence of a large high risk region, including the southern part of the province of Brindisi and the eastern and southern part of the Salento peninsula, in which an increasing trend in the northsouth direction was found.
Ecological correlation study with deprivation (Cadum Index): For males, posterior mean of the ecological regression coefficient β resulted to be 0.04 with 95% posterior credible interval equal to (0.01, 0.08); similarly, β was estimated as equal to 0.03 for females (95% posterior credible interval: 0.16, 0.10). Moreover, there was some indication of nonlinearly increasing relative risk with increasing deprivation for higher deprivation levels. For females, it was difficult to postulate the existence of any association between risk and deprivation.
Cluster detection: cluster detection based on a modified BYM model identified two large unexplained increased risk clusters in the centraleastern and southern part of the peninsula. Other secondary clusters, which raise several complex interpretation issues, are present.
Conclusion
Our results reduce the alleged role of the industrial facilities located around the province of Taranto: in particular, air pollution produced around the city of Taranto (which lies to the west of the province of Lecce) has been often identified as the main culprit of the mortality excess, a conclusion that was further supported by a recent study on the direction of prevailing winds on Salento. This hypothesis is contradicted by the finding that those municipalities that directly border on the province of Taranto (belonging to the socalled "JonicoSalentina" band) are those that present low mortality rates (at least for males). In the same way, the responsibilities of energy production plants located in the province of Brindisi (Brindisi province lies to the north) appear to be of little relevance. For females, given the situation observed in the city of Lecce, and given the substantial increase in mortality observed in younger age classes, further investigation is required into the role played by changes in lifestyle, including greater net propensity to smoke that women have shown since the 80s onwards (a phenomenon which could be amplified in a city traditionally cultured and modern as Lecce, as the tobacco habit is a largely cultural phenomenon). For males, the presence of high levels of deprivation throughout the eastern and southern Salento is likely to play an important role: those with lower socioeconomic status smoke more, and gender differences may be explained on the basis of the fact that in less developed areas women have less habit to tobacco smoking and alcohol drinking (and other harmful lifestyles), which are seen as purely masculine behaviour: research into the role of material deprivation and individual lifestyle differences between genders should be further developed.
Keywords
Background
This study analyses the spatial distribution of mortality from lung cancer registered in the period 1992–2001 in the province of Lecce, Italy. The motivation of this work is that the Salento peninsula (traditional name of the province of Lecce, indicating the subregion of Italy that stretches on the South of Apulia, between the Ionian Sea to the west and the Adriatic Sea to the east) represents, in Italy, a case of considerable interest: since the first report of the Epidemiological Observatory of the Region of Apulia, which examined the causes of death in a time span of five years (from 1998 to 2002), a consistent excess of lung cancer mortality among residents in the province of Lecce was found (compared to the level observed in other Apulian provinces, see [1] for data updated in 2006). We will see shortly that these results are widely confirmed by available data: it must however be pointed out that, as in previous studies, our conclusions are based on mortality rather than on incidence data (we will discuss later this key point).
The possible causes of the risk excess were the subject of intense debate (often with inadequate scientific methodology) in the local press. We can cite, for example, the alarm raised regarding emissions produced by the Enel "Federico II" powerplant of Cerano (BR), located about 30 km from the northern provincial boundary of Lecce to the north: the plant produces, according to data reported by Legambiente (data that has been given broad emphasis in local newspapers) more than one third of total national emissions of CO_{2}: 14.4 million tonnes of CO_{2} per year against a national total of about 51.6 million in 2006, a figure worsened by the presence of a second powerplant operated near Brindisi by Edipower, from which a release of nearly 3.8 million tonnes of CO_{2} into the atmosphere per year has been estimated [2]. Similar accusations were raised about megaindustrial steel factories operated by the Ilva Corporation, and located in the municipalities of the city of Taranto and Statte, about 40 km from the northwestern boundary of Lecce province: coke oven batteries create a carcinogenic risk for workers, because of exposure to benzene and asbestos, and given the vicinity to the city and the inadequacy of measures of pollution control, a risk also exists for the general population [3]. In both cases the transport mechanism of pollutants released into the atmosphere would be caused by the peculiar wind system in this area.
The environmental situation in the province of Lecce is characterized by the absence of apparent environmental risk sources. However, there are some noteworthy points that must be mentioned and that will be better discussed in the light of the results obtained. For example, the situation of environmental pollution is particularly alarming. Evidence from the State Forestry Department shows that in the period 2001–2002 there were 4800 illegal landfills in Italy, 600 in Apulia, and of these 50% (over 300) in the province of Lecce alone, compared to 50 in the province of Foggia and 15 in Brindisi [4]. The illegal disposal of toxic substances in the environment constitutes a risk also for agrofood and livestock activities: the illegal landfills surveyed by the State Forestry Department are in fact located in rural areas, where contact with the water springs and crops represents a serious health risk to consumers. To all this we can add, particularly in Salento, the serious problem of abandonment, or worse of illegal incineration of plastics used in agriculture, with all the wellknown consequences on the food chain. There is ample confirmation of this situation contained in the report of the Environmental Protection Agency on waste disposal [5]: of more than 153,000 tonnes of special hazardous waste disposed in Apulia, only 62,000 were disposed using authorized landfills (38,000), incinerators (16,000) and treatment plants (7,000); the remaining 91,000 tons and more were disposed of illegally, usually in abandoned quarries or in warehouses rented by waste traffickers, and presumably this situation still continues. With regard to industrial structures in the area, we must mention the presence of an important plant for the cement production located since the sixties in the municipality of Galatina (Istat code: 75029), a major industrial plant in Maglie (Istat code: 75039) for olive oil refining and extraction of oil from olive residues, formerly owned by the Regional government of Apulia and rented in the early 80s by a group of local entrepreneurs: finally, it is worth noting the presence of an incinerator located in the town of Surbo (Istat code: 75083), about 20 km from Lecce, that burns hospital waste and expired pharmaceuticals.
However, it is obvious that, apart from considerations based solely on intuition or hearsay, we need a thorough descriptive epidemiological analysis, which is the fundamental tool for a better understanding of the dynamics observed in incidence/mortality rates, both in time and space; consideration of spatial dimension is always required to suggest plausible aetiological hypotheses and to identify putative sources of risk. Of course, the generation of hypotheses is only a step towards further analytical studies to confirm the suspected causal relationship (see [6–8] as useful references on the role of spatial epidemiology, and on the links with geographical information systems). Most articles appearing in the press and on specialized websites, on the other hand, report data that are too general and too aggregated in order to be considered really useful.
In an aspatial setting, traditional risk factors described in the literature include primarily the habit of smoking: [9], a classic study, estimates that about 30% of the incidence of all cancers observed in the United States during the period chosen for the survey was attributable to consumption of cigarettes. As part of the alleged causal relationship smokingcancer, [10] shows that the most frequent locations are the oropharynx, the larynx, the lung, oesophagus and bladder; [11, 12] report that cigarette smoking increases up to ten times the risk of occurrence of lung cancer and up to six times that of occurrence of laryngeal cancer. The relative risk remains higher than for those who never smoked, even after a period of abstinence of more than 40 years [9]. Even passive smoking is another well studied risk factor: passive smokers inhale a complex mixture of smoke and other combustion products, a phenomenon that is commonly referred to as environmental tobacco smoke (ETS). Early studies [13, 14] have shown that the risk of contracting lung cancer is significantly higher in women married to smokers than women married to nonsmokers, while the most recent literature is definitely focused on relations between ETS and occupational exposures such as tobacco smoke inhaled at work [15].
Another important cause of lung cancer is occupational exposure to carcinogens: the risk of lung cancer in those exposed to some dangerous substances (such as benzopyrene, asbestos or metals such as hexavalent chromium, nickel and arsenic) is on average 4–8 times greater than for the general population [16–24]. It is worth noting that environmental asbestos exposures has been repeatedly reported as a main risk factor of pleura and lung cancer incidence (in Apulia for example [25], but see also [26] for a wider review). Arsenic in drinking water is another example of environmental exposure that cannot be totally avoided [26].
The real importance of air pollution as a risk factor is still unclear: some older classic papers show that a small percentage (1–2%) of lung cancer cases can be attributed to air pollution [27, 28], while the estimates of attributable risk appear decisively higher in more recent works [29]. Even the degree of urbanization and the incidence of lung cancer are sharply associated [30–32], this association could be explained by a confounding effect due to individual causes of diseases, such as smoking habits and occupational exposure, which obviously have a significantly greater importance in most densely populated areas [33]. Of course, epidemiologic evaluation has been often confounded by difficulties in defining and measuring air pollution, and evaluating the effects of lowlevel exposures in the general population. As we will soon see, our data show an excess of lung cancer in areas that can be regarded as weakly urbanized.
Even the role of the inert gas Radon as a risk factor is the subject of intense studies [34, 35]: although chemically inert, it is also radioactive and is transformed into yet other radioactive elements, usually called "children" who are electrically charged and attach onto fine particulate matter, and can then be inhaled and deposited on the surfaces of lung tissues. Moreover, radon is a ubiquitous domestic pollutant (indoor radon) as it penetrates buildings through gas found underground [11]. We will discuss further the spatial distribution of Radon in the Apulian territory, and its possible association with disease occurrence.
In this paper we have taken into account as well the fact that socioeconomic level may be a confounding variable of the arealevel spatial distribution of a given disease: this is because the socioeconomic variables tend to be associated with individual risk factors, while on a larger scale they are generally associated with zones of high pollution and massive presence of industrial plants [36]. So, even without a direct effect from environmental exposure, one can still detect a spurious association between putative sources of pollution and levels of incidence/mortality due to disease occurrence. The socioeconomic factors can be summarized by a synthetic index of material deprivation built at the area level: such an index is usually related to the prevalence of characteristics such as unemployment, low employment rates or lowquality housing and services [37]. We will see that association with deprivation is not strong for our mortality data, but it seems to be the only reasonable hypothesis that can explain the diffuse risk increase that is almost everywhere present in the province of Lecce (except for some localized 'hotspots' that cannot be explained on the basis of poverty level).
Given the discussion we have presented, the objective of this paper is to study on a scientifically sound basis the spatial distribution of risk for lung cancer mortality in the province of Lecce. Our goal is to demonstrate that most of the previous explanations are not supported by data, and that methods of descriptive epidemiology are of primary importance to generate sensible etiologic hypothesis. To this end, we will follow a hybrid approach that combines both frequentist and Bayesian disease mapping methods; furthermore, we define a new sequential algorithm based on a modified version of the BYM model, suitably modified to detect geographical clusters of disease and to confirm results obtained on the basis of posterior summaries of the "pure" BYM model.
Methods
We considered the following analyses: 1) calculation of the gross incidence rate at provincial level, as well as of standardized rates to facilitate comparison with other provinces of Apulia and the comparison with the mortality rates observed in Italy and other European nations; the reference period used for these analyses differs from that used for spatial risk estimation, for the reason that we were interested to compare patterns of temporal disease evolution in the province of Lecce with respect to other Apulian provinces as well, a task that requires a larger temporal window and that was possible due to the wide availability of data at provincial level; 2) spatial analysis to build a risk map based on the specification of an areaspecific Poisson model, where the highrisk areas are identified on the basis of pvalues associated with the null hypothesis of noincreased risk; 3) spatial analysis adjusted for the presence of spatial correlation and extraPoisson variation in areaspecific relative risks, using a BesagYorkeMollié (BYM) model estimated by Markov Chain Monte Carlo (MCMC) methods in Winbugs; 4) adjustment for material deprivation by inserting an ecological covariate in the BYM model to take account of socioeconomic score at areal level; 5) disease cluster detection using a new modelbased Bayesian paradigm based, on a suitable modification of the BYM model.
The analyses were conducted separately for both sexes: the difference in terms of incidence and/or mortality between males and females is becoming less relevant, as witnessed also by numerous references in literature. The latest data even reduce the importance of smoking as a risk factor in males, emphasizing instead the large increase in the prevalence of female smokers [38]: [39] report data about the "epidemic" of lung cancer mortality among young women in Europe. As we do not have a priori reasons to reject the hypothesis that health effects of other risk factors could be quite different between the two sexes, the adoption of separate analyses appears to be appropriate.
Data
The data at the provincial level were obtained from the Cislaghi Italian mortality atlas based on Istat data [40], considering deaths occurred in Apulia during the period 1981–2001 for the class of disease indicated by code 162 on the ICD9CM classification (i.e.: malignant neoplasm of trachea, bronchus and lung, IX^{th} Revision of the International Classification of Causes of Death, published in 1979 by the World Health Organization). To analyze the overall mortality dynamics, the reference period has been divided into seven triennials: 1981–83, 1984–86, 1987–89, 1990–92, 1993–95, 1996–98, 1999–01.
To estimate the spatial distribution of mortality, the number of deaths for lung cancer in each of the 97 municipalities in the province of Lecce was considered in the period between January 1^{st}, 1992 and December 31th, 2001 inclusive. The data were obtained once again from the Cislaghi atlas: the choice of an extensive reference period was due as much to the necessity to not have information too scattered within each stratum in which the dataset was divided, as the impossibility to obtain the data, because of possible identity disclosures if the reference period was too short. It should also be noted that the areas actually analyzed are 96 in total; due to the abovementioned privacy issue, the Cislaghi atlas considers the neighbouring municipalities Racale and Taviano as a single administrative unit. To draw the maps, the areas of these two municipalities were aggregated, creating a single area to which was assigned the name of the Taviano municipality. These operations were carried out on a Geographic Information System (GIS) using the union operation: the population of the new area was simply assumed as the sum of the populations of both the municipalities.
Calculation of mortality rates for provincial data
where D_{ t }is the number of deaths observed for the cause of death considered in the period t, t +k (expressed in years, in our case k = 2 considering three years), while R_{ t }is the sum of the population at risk in the same period. Exploiting agespecific data in the expression of M_{ t }(both the numerator as well the denominator), a specific rate M_{t, j}, for j = 1,..., J, is obtained for each agegroup considered: the division originally planned from the Cislaghi atlas includes the classes 0, 1–14, 15–34, 35–54, 55–64, 65–74, 75 +, although the first two have never been taken into account, because no deaths have ever been observed in either of the triennial considered.
where w_{ j }is the relative weight of the jth age group in the standard European population.
Spatial analysis 1: areaspecific Poisson model
where E_{ i }= ∑_{ j }N_{ ij }q_{ j }is the number of expected deaths in the ith area. The choice of the stratumspecific reference rates q_{ j }is crucial [44]: we estimated each rate with q_{ j }= ∑Y_{ ij }/∑_{ i }N_{ ij }(internal standardization, [45, 46]); this approach centres the data with respect to the map, and the areas where there is an excess of risk are those in which the number of observed cases is higher than the number of expected cases.
Pvalue based maps
that may be used to summarize the significance of the statistical hypothesis of no increased risk within the ith area H_{0}: θ_{ i }= 1, against the alternative hypothesis H_{0} : θ_{ i }> 1 (if > 1, otherwise the alternative hypothesis is H_{0} : θ _{ i }< 1 if < 1).
All of these values, for i = 1,..., N may be classified to draw a probability map, attributing to each area a colour level that denotes class membership [47]. At a significance level of α = 0.05, highrisk areas are those in which ρ_{ i }<α and > 1 occur simultaneously.
It is worth noting that probability maps may not be very informative, as pvalues alone do not give any information about the level of risk.
Spatial analysis 2: the BesagYorkMolliè (BYM) model
Apart the abovementioned issues, outcomes in spatial units are often not independent of each other. Risk estimates of areas that are close to each other will tend to be positively correlated as they share a number of spatially varying characteristics. Ignoring the overdispersion caused by spatial autocorrelation in the residual leads to incorrect inferences: in particular, an extreme value of ρ_{ i }may be more due to the lack of fit of the saturated Poisson model than to its deviation from the null hypothesis of neutral risk. This effect of overdispersion due to spatial autocorrelation is very strong only for small area (i.e. areas with very low populations), while is negligible for large municipalities [48].
As risk estimates of areas that are close to each other will tend to be positively correlated, if all such characteristics could be properly measured, then the model would be fully specified and residual spatial variation would be fully explained. However, this will never be the case and unmeasured spatial factors will introduce spatial dependence that can be described by introducing into the model suitable random effects. This is the reason why we considered a standard BYM model for a more refined analysis [49, 50]; we briefly recall the formulation of the hierarchical model:
where α is a baseline logrelative risk, and v_{ i }and ε_{ i }are random effect representing spatial clustering (autocorrelation) and unstructured extraPoisson variation respectively. Distributional forms commonly adopted at the second level for the two random effects are:
Clustering – v_{ i }is defined by the CAR (Conditional Autoregressive, [51]) specification , where , ∂_{ i }is the set of all areas neighbouring the ith area, and n_{ i }is the number of elements within the set. Hence, the CAR effect describes spatially varying risk factors, based on which neighbouring areas tend to have similar relative risks. The specification of a CAR effect forces the estimates of areaspecific logrelative risks toward a local mean (with an obvious smoothing effect and noise reduction of the map);
Heterogeneity – ε_{ i }is used instead to describe the sources of error not spatially structured: for this reason we hypothesize, as usual, the exchangeable specification .

3^{rd} level – Prior distributions of the parameters of the two random effects must be chosen [52]: let G(a, b) denote the Gamma distribution with expected value a/b and variance a/b^{2}. For each of the two precision parameters, and , we set τ_{ v }~ Gamma(a_{ v }, b_{ v }) and τ_{ ε }~ Gamma(a_{ ε }, b_{ε)}. In this paper we used a_{ v }= 0.5, b_{ v }= 0.005 for the spatially structured component [53], which corresponds to a diffuse prior that does not artificially force a spatial structure in the logrelative risk estimates when this is not actually present in the data. By simple calculations it can be proven that, with this choice, the standard deviation of spatially structured random effects is a random variable centred around 0.01, and the probability of observing values smaller than 0.01 or larger than 2.5 is equal to 0.01. For the nonspatially structured random effect we set a_{ ε }= 0.01, b_{ ε }= 0.01: these values correspond to a noninformative prior, and they were chosen on the ground of an empirical tradeoff between the need to not use an overly informative prior (given the previous absence of knowledge on spatial distribution of mortality from lung cancer), and the care taken to avoid deterioration of the convergence of the estimation algorithm, a very common issue when a flat prior is used.
Prior to model estimation we tested SMRs for the presence of overdispersion and spatial autocorrelation, as the use of BYM model needs to be motivated if the evidence for spatial autocorrelation or extraPoisson variation is not strong. We applied the following battery of tests: Pearson chisquare and PotthoffWhittinghill's statistics for assessing homogeneity of relative risks [54]; Dean's overdispersion score for testing the presence of extraPoisson variation versus the null hypothesis of Poisson distribution [55]; Moran's I and Geary's C tests to assess the presence of spatial autocorrelation, accounting for overdispersion by assuming a negativebinomial distribution as the sampling model needed to compute the null distribution by means of parametric bootstrap [56].
Parameter estimation and disease mapping
Posterior estimates of the parameters were obtained by simulating from the joint posterior by means of a Markov Chain Monte Carlo (MCMC) algorithm, using the OpenBugs 3.0.3 software together with the GeoBugs 1.3 extension [57]; 6,000 burnin iterations on two parallel chains starting from overdispersed values were simulated: the convergence was checked using the BrooksGelmanRubin diagnostic, summarized with the coefficient which tends to 1 in case the convergence is achieved [58]. Subsequently, another 3,000 iterations (for each chain) were generated: only one out of each three was considered for estimation, in order to eliminate the serial autocorrelation and to reduce the standard error estimate of the parameters. We monitored and estimated the areaspecific relative risks θ_{ i }by means of the MCMC algorithm described in the previous section. Maps were drawn by dividing the whole range of relative risk posterior estimates in five nonoverlapping subintervals delimited by quintiles, assigning a suitable colour to each interval and classifying areas accordingly.
Within the context of cluster detection, questions have been raised about the performance of the BYM model in recovering the true risk surface. For this reason, rather than insisting on the interpretation of relative risk posterior estimates, we supplemented our result by monitoring and estimating areaspecific posterior probabilities δ_{ i }= E[I(θ_{ i }> 1)Y] = Pr{θ_{ i }> 1Y} (here I(•) denotes the event indicator function). Based on those new posterior areaspecific summaries, maps were drawn dividing the interval [0,1] in ten equally spaced subintervals, and assigning a colour to each area accordingly. The resulting maps are likely to be insufficiently informative on the actual risk level (as well as pvalue based maps are), but the may be indeed useful to confirm the presence of "hotspots" of highrisk areas exploiting results given in a wide simulation experiment based on synthetic data, where a feasible benchmark for the risk calibration problem was provided [59]. The authors, defining three different loss function representing weighted tradeoffs between false positive and false negative, showed that by declaring at "high risk" those areas where > 0.8, each area with relative risk below 2 was classified as high risk with a probability of at least 75% if the expected cases in each area are between 10 and 20; this probability approximates to 1 for areas with a relative risk of 3 if the expected cases are never less than 5.
Correction for edge effects
Even when disease counts are independent, any smoothing operation applied to SMRs that borrows information from neighbouring areas will induce edge effects, that consists in biased estimates in areas located near the boundary of the investigated region because information on what happen on the other side of the boundary is missing: a larger estimation variance will be also found in boundary areas due to the low proportion of neighbouring cases [60]. Because some putative sources of pollution are located outside the study area, and the estimation of largescale patterns are likely to be affected by such statistical biases, it was necessary to adjust for edge effects. A classical methods is to employ guard areas, which are areas external to the main study window of interest and are added to the window to provide a guard area [60]: in this case, given the availability of mortality data for the whole Apulia, it was possible to estimate a global disease map in order to assess, without distortions, the presence of spatial patterns originating from the putative sources located around Brindisi and Taranto, and the largescale structure of disease risk. In order to make comparisons between the two maps feasible, the stratumspecific reference rates q_{ j }for the Apulia map were set equal to those used for the Lecce map (external standardization): in this way, expected counts for areas located inside the province of Lecce resulted to be the same for both maps. We also considered internal standardization for the Apulia map; for the province of Lecce considered in isolation, this is essentially equivalent to shift upward posterior estimates obtained by external standardization, but it may be undoubtedly useful in order to compare the average disease risk level in the province of Lecce with that occurring in the whole Apulia. Spatial risk estimation was carried out by simulating 12,000 MCMC burnin iterations, and subsequently another 6,000 iterations (for each chain) considering only one out of each three for estimation.
Spatial analysis 3: accounting for socioeconomic factors
The areaspecific deprivation measure x_{ i }considered here is the Cadum index [61]: this measure is well suited to the information flows available in Italy, and was calculated using Census data provided by Istat for the 1991 census. The amounts taken into consideration in constructing the index are as follows:

proportion of population with primary education;

proportion of rental housing;

proportion of occupied residences without bathroom;

proportion of the active workforce unemployed or looking for a first job;

proportion of singleparent families with children.
The index in question is constructed by calculating, for each area, the z score of each variable standardized with respect to the national average and to the national standard deviation, and then adding the five scores thus obtained: by design, higher scores are observed in poorer areas. A posterior summary of the importance of Cadum index in explaining the spatial pattern of disease risk can be obtained by monitoring and drawing residual spatial variation exp(v_{ i }+ ε_{ i }), an areaspecific quantity that is adjusted to eliminate the net effect due to the level of material deprivation in the ith area: the presence of a weak spatial pattern in residual spatial variation indicates a strong association with deprivation. Alternatively, the ecological coefficient β can be considered significantly positive if its posterior mean is greater then zero and its 95% percent credible interval excludes zero: this means that there is an actual risk increase in those areas where the level of poverty is higher. It is also interesting to note that when a covariate inserted in an ecological BYM model result to be appropriate to model the spatial variation of risk, random effects v_{ i }and ε_{ i }may become unidentifiable and care must be taken to avoid poor convergence and invalid inferences if an improper posterior is used [62].
Spatial analysis 4: cluster detection and inference
where a clusterspecific effect α_{ Z }enters the model as the coefficient of the dummy variable I(A_{ i }∈ Z), which assumes value 1 if area A_{ i }is in Z and 0 otherwise. In a like vein to the modelbased approach introduced in [64], where an overdispersion parameter similar to the general overdispersion parameter in Poisson regression was used to reduce the effect of extraPoisson variability on cluster location detection, our Bayesian formulation expressly addresses the problem that cluster detection nominal type I error and power of classical Scan Statistics are likely to be appreciably affected by the presence of spatial autocorrelation [65]. We propose the following sequential modelbased algorithm for cluster detection (here by "cluster" we mean any collection of neighbouring areas):

Given an initial collection Δ of candidate zones, conditionally to Z, we fit our model for all Z ∈ Δ by suitably tuned MCMC simulations. The collection of fitted models is ranked (in terms of parsimony and predictive power) on the ground of the Deviance Information Criterion (DIC, which is a hierarchical modelling generalization of Akaike Information Criterion introduced in [66]: between two competing models the one that has a lower DIC score should be preferred, see the legend of Tab. 3 for further details), and the first cluster Z* is identified accordingly, by means of the clusterspecific indicator variable entering the model that has the lowest DIC value; a noninformative Gaussian prior is usually assumed for α_{ Z };

Let be posterior estimate of the clusterspecific effect: a new model containing an additional offset term accounting for the effect due to Z* is considered, treating I(A_{ i }∈ Z*) as an explanatory variable with known coefficient equal to . Let Δ* be the collection of zones including Z* and every cluster overlapping with Z*: the previous step is iterated assuming Δ  Δ* as the collection of candidate clusters, and a second optimal cluster is identified in case;

The procedure stops when no better data explanation is possible by letting further clusterspecific terms enter the model: this is easily appreciated by means of the sequence of the DIC scores of the best model of each iteration, in the sense that the algorithm ends when such sequence becomes increasing.
It should be noted that cluster location detection and cluster inference are quite distinct: the abovedescribed algorithm reduces the number of feasible clusters by comparing a large number of models on the ground of a dataanalytic criterion. Anyway, declaring a cluster Z statistically "significant" is a different task: in our Bayesian approach this occurs when 95% posterior credible interval for the clusterspecific logrelative risk α_{ Z }excludes zero.
Results
The alarming extent of the phenomenon appears even more evident if, considering the triennial 1999–01, national and international comparisons are made with suitable data. For example, the AIRT 2006 report on cancer in Italy displays for males a standardized rate for Italy equal to 6.96 per 10,000 person years in the triennial 2000–02, while for females the corresponding rate is 1.27 per 10,000 person years [67] (rates of the AIRT study were calculated by means of standardization with respect to European standard). The same study shows, for Apulia, 6.64 deaths for males and 0.73 for females (per 10,000 personyears, triennial 2000–02): these figures are unquestionably comparable to those we calculated for the period 1999–01 for Apulia (8.35 for males and 0.85 for females). Therefore, with reference to the period 1999–01, we can conclude that death due to lung cancer in the province of Lecce is almost double for males, compared with the national and Apulian average, while for females it is much higher than the Apulian average but comparable to the national average; once again there is a disparity between the two sexes that, as we will soon see, will be confirmed in the geographical analysis.
Spatial analysis: results
SMRs for lung cancer in the province of Lecce
Maximum likelihood (SMR) and Bayesian estimates of the relative risk of mortality for lung cancer in the province of Lecce, 1992–2001, for selected areas
Istat Code  Area  Y _{ i }  E _{ i }  SMR  Significance  95% CI _{ SMR }  95% CI _{ Bayes }  

Males  
75008  Bagnolo del Salento  15  8.87  1.69  0.05  (1.02–1.81)  1.21  (0.92–1.56) 
75072  Santa Cesarea Terme  25  15.31  1.63  0.02  (1.10–2.42)  1.23  (0.96–1.55) 
75057  Otranto  34  21.24  1.60  0.01  (1.14–2.24)  1.25  (0.99–1.57) 
75025  Cursi  28  17.97  1.56  0.02  (1.08–2.26)  1.26  (0.98–1.64) 
75061  Poggiardo  41  27.69  1.48  0.02  (1.09–2.01)  1.20  (0.96–1.47) 
75018  Castrignano de'Greci  26  18.06  1.44  0.05  (0.98–2.11)  1.20  (0.95–1.50) 
75050  Morciano di Leuca  31  21.84  1.42  0.04  (1.00–2.02)  1.17  (0.92–1.51) 
75093  Vernole  49  35.39  1.38  0.02  (1.05–1.83)  1.18  (0.95–1.46) 
75051  Muro Leccese  37  26.87  1.38  0.04  (1.00–1.94)  1.19  (0.95–1.49) 
75005  Andrano  34  24.88  1.37  0.05  (0.98–1.91)  1.12  (0.90–1.41) 
75053  Neviano  43  31.73  1.36  0.04  (1.00–1.83)  1.15  (0.92–1.40) 
75064  Presicce  39  29.25  1.33  0.05  (0.97–1.82)  1.13  (0.91–1.41) 
75021  Collepasso  49  37.39  1.31  0.04  (0.99–1.73)  1.14  (0.93–1.38) 
75026  Cutrofiano  58  45.27  1.28  0.04  (0.99–1.66)  1.16  (0.95–1.39) 
75039  Maglie  86  71.28  1.21  0.05  (0.98–1.49)  1.16  (0.98–1.37) 
Females  
75095  San Cassiano  5  1.22  4.09  0.01  (1.70–9.82)  1.20  (0.69–2.20) 
75035  Lecce  100  54.62  1.83  0.01  (1.50–2.23)  1.65  (1.33–2.00) 
BYM model for the province of Lecce
Tests for assessing the presence of heterogeneity and spatial autocorrelation in the relative risks of mortality for lung cancer in the province of Lecce, 1992–2001
Test  Sampling null distribution  Number of iterations  Pvalue 

Males  
Pearson chisquared  Multinomial  9999  0.01 
PothoffWhittinghill  Multinomial  9999  0.01 
Dean's overdispersion test  0.01  
Geary's C  Negative Binomial  9999  0.03 
Moran's I  Negative Binomial  9999  0.01 
Females  
Pearson chisquared  Multinomial  9999  0.02 
Pearson chisquared *  Multinomial  9999  0.42 
PothoffWhittinghill  Multinomial  9999  0.01 
PothoffWhittinghill *  Multinomial  9999  0.42 
Dean's overdispersion test  0.01  
Dean's overdispersion test*  1.00  
Geary's C  Negative Binomial  9999  0.53 
Moran's I  Negative Binomial  9999  0.64 
BYM model for the whole Apulia
Ecological correlation study with deprivation (Cadum index)
Cluster detection
Results of the cluster detection sequential algorithm based on a suitable modification of the BYM model.
Discussion
The results of this work, reported in the previous paragraph, makes it appear that lung cancer mortality in the province of Lecce is assuming alarming dimensions. Within the province itself, however, situation differs from municipality to municipality, but it is unquestionably not homogeneous between males and females, and this raises the problem of a rather complex interpretation. For males, as previously mentioned, the presence of two large clusters is apparent, one group of the municipalities gathered around the Maglie (Istat code: 75039) and the other located in the south of peninsula; for females, however, higher mortality rates are seen in the municipality of Lecce. Furthermore, the temporal dimension must not be neglected: Figure 2 shows clearly how, for males, the problem of increased incidence of lung cancer has subsisted, according to data available to us, for at least 25 years. Data on smoking habits and other lifestyle choices, disaggregated at the municipal level, are not available: recent data on the entire region (which unfortunately are not available separately for the two sexes) show that the percentage of nonsmokers is higher than the Italian average (64% versus 56%); they show, instead, a percentage of overweight people higher than the Italian average (36.6 against 33.9), but in any case aligned with the other regions of South Italy [69].
With regard to individual risk factors, particular attention was given, in the past, to the distribution of Radon in the Apulian territory. In a major study conducted on a sample of 310 residences in nine municipalities of Apulia (Bari, Rutigliano, Foggia, Troia, Sant'Agata di Apulia, Taranto, Latiano, Lecce e Castri di Lecce) changes in the concentration of indoor Radon in these homes were assessed during the springsummer and autumnwinter periods, and considered in association with the architectural configuration of the building and to construction materials used [70]. The reported results demonstrated that in the two selected municipalities of Lecce (Lecce and Castri di Lecce) the concentration of Radon was clearly greater than that of the other municipalities studied. The causes of this excess in these homes are mainly attributable to two factors: the building type and the geological characteristics of the subsoil. The buildings made of tufa and Lecce stone are those in which the largest concentrations of Radon were measured: it is therefore easily concluded that, in these municipalities of Lecce, the concomitant presence of many old buildings constructed with the abovementioned materials and the karst subsoil which affects the process of exhalation of Radon could, combined with smoking habits, be the cause of a significant number of lung neoplasms: at present, however, scientific feedback to this statement is required, especially bearing in mind that the spatial distribution of mortality is extremely different between the sexes, as we have wellestablished.
The presence of large industrial plants in the area of Galatina (Istat code: 75029, which borders the municipality of Cutrofiano, Istat code: 75026, located around the secondary cluster near the Maglie area), raises the question about the role of occupational exposure for male workers resident in central Salento: it would be interesting to know the employment and work place of males suffering from lung cancer to check, at least for areas classified as high risk, the possibility of an occupational exposure: furthermore, we must not forget the role that these facilities could have as point sources of pollutants and carcinogens diffused throughout the territory. These studies may become possible in the near future thanks to the consolidation of data from the JonicoSalentino Cancer Register.
Despite the impressive figures concerning the environmental situation in the province of Lecce, one cannot conclude with certainty that the problem of environmental exposure caused by toxic substances released into the environment is directly responsible for the high lung cancer mortality registered in the province of Lecce: once again the differences between the sexes make it difficult to reach a similar conclusion. In addition, as already pointed out, deaths from lung cancer among males are above the Apulian average from as early as 1981, when the issue of emergency waste was probably much less relevant (if not entirely absent).
With regard to the conclusions of the local press that, as we have already mentioned in the introduction, has often attributed the causes of risk excess to air pollution produced by the factories of Taranto, it is noteworthy to observe that in support of these hypotheses it has often been cited a recent study on wind directions of Salento, which blow predominantly from the northwest, and therefore "would spread" atmospheric pollutants to west of the province of Lecce from the province of Taranto [71]. However, our conclusions disagree: the municipalities of Lecce province within the "JonicoSalentina" belt (i.e. bordering the province of Taranto) are those with the lowest overall mortality rates (at least for males). The role of facilities for energy production installed in the province of Brindisi is unclear for similar reasons.
For females, given the situation observed in the city of Lecce, the enormous increase in mortality since the 90s is attributed to their clearly greater propensity to smoking: as the tobacco habit is a predominantly cultural phenomenon of nature, it can be assumed that many women, living in a city traditionally modern and cultured as Lecce, have gradually opted to smoke cigarettes as a symbol of an emancipated lifestyle. Even these findings, though plausible, are unfortunately weak because they were not derived from statistical models of incidence and/or mortality that take into account the latency period of lung cancer, and more precise data on cigarette consumption broken by municipality and sex: new studies will be needed to verify these new hypotheses.
For males, the presence of high levels of deprivation throughout the eastern and southern Salento is likely to play an important role: the most obvious explanation for the observed differentials is that those with lower socioeconomic status smoke more [72], and that gender differences may be explained on the basis of the fact that in less developed areas women have less habit to tobacco smoking and alcohol drinking (and to other lifestyles harmful for health), which are seen as purely masculine behaviour. Of course, potential for uncontrolled confounding may be substantial, and the inquiry about the role of deprivation should be further developed and supported by data on individual lifestyles.
Study limitations
The use of mortality rather incidence data, which could be a source of bias, was necessary by the fact that a methodical collection of cancer incidence data in the province of Lecce exists only from 1999 (thanks to the "IonicoSalentino" Cancer Registry: on these issues, see the discussion reported in [73]). However, we have a substantial degree of confidence that these distortions can be regarded as negligible, given that roughly twothirds of patients are diagnosed at an advanced stage of the disease, when available treatment approaches are rarely effective: as a consequence of the fact, lung cancer, especially nonsmall cell lung cancer, still appears as one of the more incurable neoplastic diseases, and that presents among the lowest fiveyear survival rates (see for example [74] showing the relative fiveyear survival rate – namely the relationship between the proportion of patients who survived in a cohort of individuals affected by lung cancer and the proportion of survivors expected in a comparable series of individuals not affected by the disease – equal to 11% for 2006). So, the quantitative dimensions of incidence and mortality cannot be considered too dissimilar: similar conclusion may be based on the recent data reported in [75].
The use of an areal approach invites the criticism of ecological fallacy [76], that is all individuals living in an area share the same characteristics, which is clearly not the case. In addition, the aggregate nature of the study made impossible to control for potential and very important confounders, such as smoking habits (for which neither aggregate or individual data were available) or occupational exposure among others. This implies that, notwithstanding the use of advanced disease mapping methods is essential to generate new etiological hypotheses, we have no possibility of ascertaining the real causes of risk excess clusters present in the data, and we can only make conjectures. To obtain information about the contribution of individuals to disease occurrence, and include withinarea confounders, it is desirable to use both data collected on areas and on individuals [77]: this approach will be pursued in future papers as soon as individuallevel data will become available.
Conclusion
The data analyzed in this paper reveal that in Salento there are, unfortunately, all the adverse conditions for a high mortality for lung cancer, although the real reasons for this risk excess remain, until now, not fully understood: this work is a preliminary analysis that has served to estimate the spatial pattern of risk and to generate new hypotheses for study: research into the role of material deprivation and individual lifestyle differences between genders should be further developed: other findings, highlighted by the press, seem unrealistic in the light of data. The importance of geographical epidemiology as a tool for analyzing the dynamics and determinants of health and disease in the community was stressed.
Declarations
Acknowledgements
We wish to thank Dr. Giusi Graziano, from the Mario Negri Sud Institute in Santa Maria Imbaro (Chieti, Italia), for her valuable support in the use of Cislaghi mortality atlas. Comments from two anonymous reviewers contributed to improve the quality of this paper.
Authors’ Affiliations
References
 Osservatorio Epidemiologico Regione Puglia: Stato di Salute della Regione Puglia. 2006,http://www.oerpuglia.it/Stato_sez1.pdf
 Greenpeace: Italia – Emissioni di CO_{2} 2005–06: elaborazioni su dati della Commissione Europea.http://www.greenpeace.org/italy/ufficiostampa/comunicati/tabellaemissioni
 Giua R, Spartera M, Viviano G, Ziemachi G, Carbotti G: Cancer risk for cokeoven workers in the Taranto steel plant. Epidemiol Prev. 2005, 29 (5–6 Suppl): 424.PubMedGoogle Scholar
 Corpo Forestale dello Stato: Report Ambiente. 2003,http://www3.corpoforestale.it/flex/cm/pages/ServeBLOB.php/L/IT/IDPagina/393
 Agenzia per la Protezione dell'Ambiente e per i Servizi Tecnici e Osservatorio Nazionale per i Rifiuti: Rapporto Rifiuti 2004 – Vol. II: Rifiuti Speciali. 2004Google Scholar
 Beale L, Abeliano J, Hodgson S, Jarup L: Methodological issues and approaches to spatial epidemiology. Environ Health Perspect. 2008,Google Scholar
 Jarup L: Health and environment information systems for exposure and disease mapping, and risk assessment. Environ Health Perspect. 2004, 112 (9): 9957.PubMedPubMed CentralView ArticleGoogle Scholar
 Mather FJ, White LE, Langlois EC, Franklin SC, Swalm CM, Shaffer JG, Hartley WR: Statistical methods for linking health, expoures, and hazards. Environ Health Perspect. 2004, 112 (14): 14401445.PubMedPubMed CentralView ArticleGoogle Scholar
 Doll R, Peto R: Cigarette smoking and bronchial carcinoma: dose and time relationships among regular smokers and lifelong nonsmokers. J Epidemiol Community Health. 1978, 32 (4): 303313.PubMedPubMed CentralView ArticleGoogle Scholar
 Tominaga S: Major avoidable risk factors of cancer. Cancer Letters. 1999, 143: S19S23.PubMedView ArticleGoogle Scholar
 Alberg AJ, Samet JM: Epidemiology of lung cancer. Chest. 2003, 123 (90010): 21S49.PubMedView ArticleGoogle Scholar
 Alberg AJ, Brock MV, Samet JM: Epidemiology of lung cancer: looking to the future. Journal of Clinical Oncology. 2005, 23 (14): 31753185.PubMedView ArticleGoogle Scholar
 Hirayama T: Nonsmoking wives of heavy smokers have a higher risk of lung cancer: a study from Japan. Br Med J (Clin Res Ed). 1981, 282 (6259): 183185.View ArticleGoogle Scholar
 Trichopoulos D, Kalandidi A, Sparros L, MacMahon B: Lung cancer and passive smoking. Int J Cancer. 1981, 27: 14.PubMedView ArticleGoogle Scholar
 Veglia F, Vineis P, Overvad K, Boeing H, Bergmann M, Trichopoulou A, Trichopoulos D, Palli D, Krogh V, Tumino R, Linseisen J, Steindorf K, RaaschouNielsen O, Tjonneland A, Gonzalez C, Martinez C, Dorronsoro M, Barricarte A, Cirera L, Quiros J, Day N, Saracci R, Riboli E: Occupational exposures, environmental tobacco smoke, and lung cancer. Epidemiology. 2007, 18 (6): 769775.PubMedView ArticleGoogle Scholar
 Gibbs G, Sevigny M: Mortality and cancer experience of Quebec aluminum reduction plant workers, part 4: cancer incidence. J Occup Environ Med. 2007, 49 (12): 135166.PubMedView ArticleGoogle Scholar
 Marinaccio A, Scarselli A, Binazzi A, Mastrantonio M, Ferrante P, Iavicoli S: Magnitude of asbestosrelated lung cancer mortality in Italy. Br J Cancer. 2008, 99: 173175.PubMedPubMed CentralView ArticleGoogle Scholar
 Pira E, Pelucchi C, Buffoni A, Palmas A, Turbiglio M, Negri E, Piolatto P, C LV: Cancer mortality in a cohort of asbestos textile workers. British Journal of Cancer. 2005, 92: 580586.PubMedPubMed CentralView ArticleGoogle Scholar
 Doll R: Mortality from lung cancer in absestos workers 1955. Br J Ind Med. 1993, 50 (6): 48590.PubMedPubMed CentralGoogle Scholar
 Birk T, Mundt K, Dell L, Luippold R, Miksche L, SteinmannSteinerHaldenstaett W, Mundt D: Lung cancer mortality in the German chromate industry, 1958 to 1998. J Occup Environ Med. 2006, 48 (4): 42633.PubMedView ArticleGoogle Scholar
 Sorahan T, Williams SP: Mortality of workers at a nickel carbonyl refinery, 1958–2000. Occup Environ Med. 2005, 62 (2): 8085.PubMedPubMed CentralView ArticleGoogle Scholar
 Doll R: Cancer of the lung and nose in nickel workers. Br J Ind Med. 1958, 15 (4): 217223.PubMedPubMed CentralGoogle Scholar
 Lubin JH, Pottern LM, Stone BJ, Fraumeni J, Joseph F: Respiratory Cancer in a Cohort of Copper Smelter Workers: Results from More Than 50 Years of Followup. American Journal of Epidemiology. 2000, 151 (6): 554565.PubMedView ArticleGoogle Scholar
 Lee A, Fraumeni J: Arsenic and respiratory cancer in man: an occupational study. J Natl Cancer Inst. 1969, 42 (6): 104552.PubMedGoogle Scholar
 Musti M, Pollice A, Cavone D, Dragonieri S, Bilancia M: The relationship between malignant mesothelioma and an asbestos cement plant environmental risk: a spatial casecontrol study in the city of Bari (Italy). Int Arch Occup Environ Health. 2009, 82 (4): 48997.PubMedView ArticleGoogle Scholar
 BrüskeHohlfeld I: Environmental and occupational risk factors for lung cancer. Methods Mol Biol. 2009, 472: 323.PubMedView ArticleGoogle Scholar
 Doll R: Atmospheric pollution and lung cancer. Environ Health Perspect. 1978, 22: 2331.PubMedPubMed CentralView ArticleGoogle Scholar
 Doll R, Peto R: The cause of cancer: quantitative estimates of avoidable risk of cancer in the United States today. Journal of the National Cancer Institute. 1981, 66 (6): 11911308.PubMedGoogle Scholar
 Clapp RW, Jacobs MM, Loechler EL: Environmental and occupational causes of cancer: new evidence 2005–07. Rev Environ Health. 2008, 23: 137.PubMedPubMed CentralView ArticleGoogle Scholar
 Uccelli R, Mastrantonio M, Di Paola M: Distribution of causes of death in communities with different urbanization levels. Epidemiol Prev. 2000, 24: 2837.PubMedGoogle Scholar
 Schouten LJ, Meijer H, Huveneers JA, Kiemeny LA: UrbanRural Differences in Cancer Incidence in The Netherlands, 1989–1991. International Journal of Epidemiology. 1996, 25 (4): 729736.PubMedView ArticleGoogle Scholar
 Howe HL, Keller JE, Lehnherr M: Relation between Population Density and Cancer Incidence, Illinois, 1986–1990. American Journal of Epidemiology. 1993, 138: 2936.PubMedGoogle Scholar
 Chellini E, Gorini G, Martini A, Giovannetti L, Costantini AS: Lung cancer mortality patterns in women resident in different urbanization areas in central Italy from 1987–2002. Tumori. 2006, 92 (4): 2715.PubMedGoogle Scholar
 Thompson RE, Nelson DF, Popkin JH, Popkin Z: Casecontrol study of lung cancer risk from residential radon exposure in Worcester county, Massachusetts. Health Phys. 2008, 94 (3): 22841.PubMedView ArticleGoogle Scholar
 Krewski D, Lubin JH, Zielinski JM, Alavanja M, Catalan VS, Field RW, Klotz JB, Letourneau EG, Lynch CF, Lyon JL, Sandler DP, Schoenberg JB, Steck DJ, Stolwijk JA, Weinberg C, Wilcox HB: A combined analysis of North American casecontrol studies of residential radon and lung cancer. J Toxicol Environ Health A. 2006, 69 (7): 533597.PubMedView ArticleGoogle Scholar
 Wakefield J: Sensitivity analyses for ecological regression. Biometrics. 2003, 59: 917.PubMedView ArticleGoogle Scholar
 Carstairs V: Socioeconomic factors at areal level and their relationship with health. Spatial Epidemiology: Methods and Applications. Edited by: Elliot P, Wakefield J, Best NG, Briggs DJ. 2000, Oxford University Press, 5167.Google Scholar
 Bolego C, Poli A, Paoletti R: Smoking and gender. Cardiovascular Research. 2002, 53 (3): 568576.PubMedView ArticleGoogle Scholar
 Bosetti C, Levi F, Lucchini F, Negri E, Vecchia CL: Lung cancer mortality in European women: recent trends and perspectives. Annals of Oncology. 2005, 16 (10): 15971604.PubMedView ArticleGoogle Scholar
 Cislaghi C: Gis 8 – Atlante Italiano di Mortalità 1981 – 2001 Versione 8.0 betatest. 2005, ATI ESAGoogle Scholar
 Inskip H, Beral V, Fraser P: Methods for ageadjustment of rates. Stat Med. 1983, 2: 455466.PubMedView ArticleGoogle Scholar
 Pascutto C, Wakefield J, Best NG, Richardson S, Bernardinelli L, Staines A, Elliot P: Statistical issues in the analysis of disease mapping data. Statistics in Medicine. 2000, 19: 24932519.PubMedView ArticleGoogle Scholar
 Wakefield J: Disease mapping and spatial regression with count data. Biostatistics. 2007, 8 (2): 158183.PubMedView ArticleGoogle Scholar
 Bernardinelli L, Montomoli C: Empirical Bayes versus fully Bayesian analysis of geographical variation in disease risk. Statistics in Medicine. 1992, 11: 9831007.PubMedView ArticleGoogle Scholar
 Banerjee S, Carlin BP, Gelfand AE: Hierarchical Modeling and Analysis for Spatial Data. 2003, Chapman and Hall/CRCView ArticleGoogle Scholar
 Waller LA, Carlin BP, Xia H, Gelfand AE: Hierarchical SpatioTemporal Mapping of Disease Rates. Journal of the American Statistical Association. 1997, 92 (438): 607617.View ArticleGoogle Scholar
 Choynowski M: Maps based on probabilities. Journal of the American Statistical Association. 1959, 54 (286): 385388.View ArticleGoogle Scholar
 Cressie NAC: Statistics for Spatial Data (Revised Edition). 1993, WileyInterscienceGoogle Scholar
 Wakefield J, Best NG, Waller LA: Bayesian approaches to disease mapping. Spatial Epidemiology: Methods and Applications. Edited by: Elliot P, Wakefield J, Best NG, Briggs DJ. 2000, Oxford University Press, 104127.Google Scholar
 Molliè A: Bayesian mapping of disease. Markov Chain Monte Carlo in Practice. Edited by: Gilks WE, Richardson S, Spiegelhalter D. 1996, Chapman & Hall/CRC, 359379.Google Scholar
 Rue H, Held L: Gaussian Markov Random Fields. Theory and Applications. 2005, Chapman & Hall/CRCView ArticleGoogle Scholar
 Carlin BP, Louis TA: Bayes and Empirical bayes Methods for Data Analysis. 2000, Chapman & Hall/CRC, SecondView ArticleGoogle Scholar
 Kelsall JE, Wakefield J: Discussion of "Bayesian models for spatially correlated disease and exposure data". Sixth Valencia International Meetinng on Bayesian Statistics. Edited by: Bernardo JM, Berger JO, Dawid AP, Smith AFM. 1999, Oxford University PressGoogle Scholar
 Bivand RS, Pebesma EJ, GómesRubio V: Applied Spatial Data Analysis with R. 2008, SpringerGoogle Scholar
 Dean C, Lawless JF: Tests for detecting overdispersion in Poisson regression model. Journal of American Statistical Association. 1989, 84 (406): 467472.View ArticleGoogle Scholar
 GómezRubio V, FerrándizFerragud J, LópezQuilez A: Detecting clusters of disease with R. J Geograph Syst. 2005, 7: 189206.View ArticleGoogle Scholar
 Thomas A, Best NG, Lunn D, Arnold R, Spiegelhalter D: GeoBUGS User Manual – Ver. 1.3, August 2007. 2007,http://mathstat.helsinki.fi/openbugs/Google Scholar
 Brooks SP, Gelman A: Alternative methods for monitoring convergence of iterative simulations. J Comput Graph Statist. 1998, 7: 434455.Google Scholar
 Richardson S, Thomson A, Best N, Elliot P: Interpreting posterior relative risk estimates in diseasemapping studies. Environ Health Perspect. 2004, 112 (9): 10161025.PubMedPubMed CentralView ArticleGoogle Scholar
 Lawson AB, Biggeri A, Dreassi AE: Edge effects in disease mapping. Disease Mapping and Risk Assessment for Public Health. Edited by: Lawson AB, Boehning D, Lasaffree E, Biggeri A, Viel JF, Bertollini R. 1999, Wiley, ChichesterGoogle Scholar
 Cadum E, Costa G, Biggeri A, Martuzzi M: Deprivation and mortality: a deprivation index suitable for geographical analysis of inequalities. Epidemiol Prev. 1999, 23 (3): 17587.PubMedGoogle Scholar
 Besag J, Green P, Higdon D, Mengersen K: Bayesian computation and stochastic systems. Statistical Science. 1995, 10: 341.View ArticleGoogle Scholar
 Aamodt G, Samuelsen SO, Skrondal A: A simulation study of three methods for detecting disease clusters. Int J Health Geogr. 2006, 5: 15PubMedPubMed CentralView ArticleGoogle Scholar
 Zhang T, Lin G: Spatial scan statistics in loglinear models. Computational Statistics and Data Analysis. 2009, 8 (15): 28512858.View ArticleGoogle Scholar
 Loh JM, Zhu Z: Accounting for spatial correlation in the scan statistics. The Annals of Applied Statistics. 1 (2): 560584.Google Scholar
 Spiegelhalter DJ, Best NG, Carlin BP, Linde van der A: Bayesian measures of model complexity and fit (with discussion). J Roy Stat Soc B. 64 (4): 583639.Google Scholar
 AIRT Working Group: I tumori in Italia – Rapporto 2006: la mortalità per tumore in Italia. Epidemiol Prev. 2006, 30 (Supplemento 2):Google Scholar
 Ferlay J, Autier P, Boniol M, Heanue M, Colombet M, Boyle P: Estimates of the cancer incidence and mortality in Europe in 2006. Ann Oncol. 2007, 18 (3): 581592.PubMedView ArticleGoogle Scholar
 Istituto Nazionale di Statistica: Stili di Vita e Condizioni di Salute: Indagine Multiscopo sulle Famiglie "Aspetti della Vita Quotidiana". 2003, ISTATGoogle Scholar
 Lattarulo O, Martucci L, Vitucci L: Indagine Radon nelle abitazioni della Regione Puglia. Atti della Settima Conferenza Nazionale delle Agenzie Ambientali "L'Innovazione al Servizio della Conoscenza e della Prevenzione: dai Sistemi di Monitoraggio alla Diffusione della Cultura Ambientale". 2003Google Scholar
 Mangia C, Martano P, Miglietta MM, Morabito A, Tanzarella A: Modelling local winds over the Salento peninsula. Metereological Applications. 2004, 11: 231244.View ArticleGoogle Scholar
 Martikainen P, Lahelma E, Ripatti S, Albanes D, Virtamo J: Educational differences in lung cancer mortality in male smokers. Int J Epidemiol. 2001, 30 (2): 264267.PubMedView ArticleGoogle Scholar
 Pickle LW: Mapping mortality data in the United States. Spatial Epidemiology: Methods and Applications. Edited by: Elliot P, Wakefield J, Best NG, Briggs DJ. 2000, Oxford University Press, 240252.Google Scholar
 KarimKos HE, de Vries E, Soerjomataram I, Lemmens V, Siesling S, Coebergh JWW: Recent trends of cancer in Europe: a combined approach of incidence, survival and mortality for 17 cancer sites since the 1990s. European Journal of Cancer. 2008, 44 (10): 13451389.PubMedView ArticleGoogle Scholar
 Ferlay J, Autier P, Boniol M, Heanue M, Colombet M, Boyle P: Estimates of the cancer incidence and mortality in Europe in 2006. Ann Oncol. 2007, 18 (3): 581592.PubMedView ArticleGoogle Scholar
 Greenland S, Robins J: Invited commentary: ecologic studies – biases, misconceptions, and counterexamples. Am J Epidemiol. 1994, 139 (8): 747760.PubMedGoogle Scholar
 Jackson C, Best N, Richardson S: Hierarchical related regression for combining aggregate and individual data in studies of socioeconomic disease risk factors. J R Statist Soc A. 2008, 171 (1): 159178.Google Scholar
Copyright
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.