Lung Cancer Mortality in Tuscany from 1971 to 2010 and Its Connections with Silicosis: A Space-Cohort Analysis Based on Shared Models

Dreassi, Emanuela

doi:https://doi.org/10.1155/2018/4964569

Computational and Mathematical Methods in Medicine

On this page

Abstract Introduction Data Results Discussion and Conclusions Conflicts of Interest Acknowledgments References Copyright Related Articles

Research Article | Open Access

Volume 2018 | Article ID 4964569 | https://doi.org/10.1155/2018/4964569

Lung Cancer Mortality in Tuscany from 1971 to 2010 and Its Connections with Silicosis: A Space-Cohort Analysis Based on Shared Models

Emanuela Dreassi¹

Academic Editor: Konstantin G. Arbeev

Received22 Aug 2017

Revised13 Dec 2017

Accepted28 Dec 2017

Published28 Jan 2018

Abstract

Lung cancer mortality in Tuscany (Italy) for males, from 1971 and 2010, is investigated. A hierarchical Bayesian model for space-time disease mapping is introduced. Such a model belongs to the class of shared random effect models and exploits the birth-cohort as the relevant time dimension. It allows for highlighting common and specific patterns of risk for each birth-cohort. The results show that different birth-cohorts exhibit quite different spatial patterns, even if the socioeconomic status is taken into account. In fact, there were different occupational exposures before and after the Second World War. The birth-cohort 1930–35 exhibits high relative risks related to particular areas. This fact could be connected with occupational exposure to risk factors for silicosis, perhaps a prognostic status for lung cancer.

1. Introduction

Descriptive epidemiology focuses on the variation of disease occurrence among populations. Studies of time-space variations in mortality are widely used in etiologic research in order to evaluate environmental exposure. Models that jointly describe space and time variation of disease risk have been extensively proposed. All of them are essentially extensions of the hierarchical Bayesian model by Besag et al. [1]. See, for example, Bernardinelli et al. [2], Waller et al. [3], Assunção et al. [4], Sun et al. [5], and Knorr-Held and Besag [6].

In particular, Knorr-Held [7] suggests a generalized linear mixed model where space, time, and different specifications of interaction terms between space and time effects are involved. Following the latter, Lagazio et al. [8] proposed the same model where the birth-cohort is the main time dimension.

Taking birth-cohort as time dimension is motivated by several biological reasons, which make it preferable to calendar period when considering the consequences of a change in the prevalence of human exposure to risk factors in the population, even if it requires the mortality or incidence data for a long time span. As a rule, age effects are related to the natural history of the disease, while birth-cohort and/or calendar period effects depend on different temporal patterns of prevalence of human exposure to risk factors. Therefore, an appropriate selection of the time scales yields important clues about the biology of the disease, and it is fundamental to correctly understand the evolution of the epidemics.

Lung cancer is the most common type of cancer in the world since 1985 and its spatiotemporal variation is examined by an extensive epidemiological literature. See, for example, Boyle and Ferlay [9] and Bray and Weiderpass [10]. The adopted time scale is the calendar period (e.g., [11]) or birth-cohort [12, 13].

As regards mortality for lung cancer, socioeconomic factors are assumed to be associated with individual exposure to risk factors (e.g., tobacco smoking or occupational exposure). In this case, the biology of the process of carcinogenesis suggests that more than 10–15 years (the latency time) should run between exposure and mortality (see [14]). See the foundational papers by Glick [15] and Mayer [16] for a discussion and a motivation of spatial analysis of cancer mortality as well as for the theoretical models that take into account the process of carcinogenesis. A strong connection between smoking habits and material deprivation has been stressed by several studies (see, e.g., [17]). Material deprivation indicators usually refer to the occurrence of subject states such as unemployment, low education, living in a very small dwelling, overcrowding, and not having a car (e.g., see [18–22]). So far, material deprivation indicators have been used as aggregate-level covariates to adjust ecological regression coefficients in small area studies [23]. In fact, a strong association with area-based deprivation and mortality was repeatedly found on one side and area-based deprivation and exposure to environmental/individual hazards on the other [24–26]. Relevant risk factors could be considered in the model with the inclusion of a material deprivation index (see [27–30]).

In this paper, a space-cohort analysis for disease mapping is developed. A model, more parsimonious than Lagazio et al. [8], is introduced in Section 2. Such a model, called Shared Model for Birth Cohort (SMBC) in the sequel, is actually a simple version of the usual intrinsic conditional autoregressive (CAR) models, but it shows reasonably good performance when comparing a disease among the strata of a given population. So far, shared models have been used for various reasons, such as to compare two or more diseases [31–33], to compare the same disease for two genders in space and time [34, 35], and for ecological regression analysis [36]. Using shared models to compare the same disease between two different birth-cohorts, as done by SMBC, seems to be new.

To assess its behavior, SMBC is then applied to a case study. The data, described in Section 3, are the lung cancer death certificates collected for males residing in the 287 municipalities of the Tuscany Region (Italy) from 1971 to 2010. The information on material deprivation is based on census data collected from the Italian Statistical Institute (ISTAT) in the years 1961, 1971, 1981, and 1991. Four variables are used as indicator of socioeconomic inequalities, namely, unemployment, low education (less than 6 years of schooling), being a tenant, and the absence of a bathroom in the house.

Finally, the results obtained by applying SMBC to the above data are discussed in Section 4. Overall, the performances of SMBC look encouraging and endorse some etiologic hypotheses of mortality for lung cancer highlighted by the medical literature. We refer to Section 4 for a detailed discussion. Here, we just note that occupational exposure to risk factors for silicosis, perhaps a prognostic status for lung cancer, is strongly supported by our application of SMBC.

2. Space-Cohort Models

Let denote the observed number of deaths for the th area and the th cohort, with and . In the sequel, following Besag et al. [1], each is assumed to have a Poisson distribution with mean value , where is the number of expected cases and the relative risk. Also, the random variables are conditionally independent given the matrix .

2.1. Space-Time Interaction Models

In Lagazio et al. [8], is modelled aswhere is a time-dependent covariate, and are, respectively, the unstructured and structured spatial variability terms [1], is an unstructured effect for the th cohort, and is an area specific cohort effect (the interaction term spatially structured). The prior distributions for , , and are multivariate normal with mean zero and precision matrices , , and , respectively. Proper Gamma priors with very high dispersion have been assumed for the hyperparameters , , .

The structure matrices [37] , , and are specified following the different nature of effects. In particular, and are taken to be identity matrices, while the structured spatial component is modelled as an intrinsic Gaussian autoregression [1, 2, 37]. Precisely, the entries of are given bywhere is the number of areas contiguous to the th one.

With this definition of , the mean of the conditional distribution of given all the other terms is , where indicates that areas and are spatially contiguous. The conditional precision is given by .

Interaction terms are assumed to be structured in space. The prior distribution is multivariate normal with mean zero and precision matrix . Again, a highly dispersed proper Gamma prior was used for . The structure matrix was defined as the Kroneker product between the matrices and an identity matrix [7, 37].

The precision is for or and for .

To complete the model specification, highly dispersed normal prior has been assumed for the regression coefficients and .

Posterior distributions of the parameters of interest have been approximated using Gibbs sampling. After a burn-in of 100,000 iterations, we retained 1,000 samples taken from the last 100,000 iterations. The posterior distributions have been summarised using the posterior mean. Because of the high number of terms in the model, convergence has been assessed only on a subset of the identifiable parameters. Gelman and Rubin [38] test and partial autocorrelation plots have been used to check for achieved convergence of relative risks and of the hyperparameters.

The interaction terms represent the differences of the spatial pattern between two or more birth-cohorts.

2.2. Shared Models for Birth-Cohort Analysis

As noted in Section 1, in this paper, shared models are used for a space-time (birth-cohort) analysis of relative risk and the resulting model is denoted by SMBC.

Shared models allow highlighting the differences between the spatial patterns of two or more birth-cohorts and represent an effective alternative to space-time interaction models. In fact, usually, shared models provide reliable information and lead to more parsimonious models.

In SMBC, following Besag et al. [1] and Knorr-Held and Best [31], it is assumed thatwhere is a birth-cohort specific intercept (overall risk level) and a spatial term. In turn, is decomposed into a shared and birth-cohort specific spatial effect; namely,where (i)the shared component is the convolution of an unstructured spatial effect and a spatially structured term ;(ii)the specific component is the convolution of an unstructured spatial effect and a spatially structured term ;(iii)the scale parameter allows the shared component to vary per birth-cohort by a constant factor.

The prior distributions are as follows. The intercept has a flat “noninformative” prior (a centered normal distribution with “large” variance). The terms , constrained to , are jointly normal with zero mean and covariance matrix The heterogeneity terms and are independent centered normal with precision parameters and . In order to cope with the spatial structure, the random effects , are modelled as intrinsic CAR models, as discussed, for example, in Lee [39], with precision parameters and .

Hyperprior distributions for the precision parameters , , , and are assumed to be Gamma .

Posterior distributions of the parameters of interest have been sampled by Gibbs sampling.

Once SMBC is implemented, an inspection of the estimated values of and provides valuable information about common and uncommon risk factors for different birth-cohorts. This procedure is exemplified in Section 4, where SMBC is applied to a real situation.

3. Data

Lung cancer death certificates were collected for males residing in the 287 municipalities of the Tuscany Region (Italy) from 1971 to 2010. For the 40 years analyzed, amounting to a total of 34,666,951 person-years, the number of recorded death certificates is 66856.

The data were made available by the Tuscany Regional Government under the research project Tuscany Atlas of Mortality 1971–1994 (see [40]) and by the Regional Mortality Register for the period 1995–2010. Deaths and corresponding populations for each municipality were cross-classified in 10 age-classes and over; the first eight age-classes being omitted) and 8 calendar periods .

The 1971–2010 database of lung cancer mortality, divided by gender at a municipality level, allows us to investigate several birth-cohorts.

The expected number of cases for each municipality was evaluated by the predicted age-specific reference rates calculated by the age-cohort model (see [41]) for the whole Tuscany Region. Figure 1 shows the epidemic curve for the considered age-classes (a) and birth-cohorts (b). The reference birth-cohort is 1920–30, when the epidemic reached its maximum. Presumably, the exposure to smoking for such a cohort began in the 1940s, after the Second World War. This model allows disentangling age effects from birth-cohort components present in longitudinal mortality data. In this way, the expected number of cases is adjusted by age and is not affected by cohort effects.

(a)

(b)

Observed and expected cases were subsequently aggregated along the diagonals of the Lexis diagram representing birth-cohorts, thus collapsing on the age dimension. For the space-time analysis, we focused on the six birth-cohorts . Table 1 illustrates incidence rates and observed cases for the whole Tuscany Region, as well as the correspondence among birth-cohorts, age-classes, and calendar periods.

Data on material deprivation were taken from census data collected from Italian Statistical Institute (ISTAT) in the years 1961, 1971, 1981, and 1991. Four variables were used as indicator of socioeconomic inequalities: unemployment, low education (less than 6 years of schooling), being a tenant, and the absence of a bathroom in the house. A composite index was then computed as the sum of the -score of the proportion of people with the conditions listed above for each municipality (see [28], for details).

The space-time dynamics of the association between mortality and deprivation level were studied using a deprivation score which results in a meaningful absolute level across the years. We considered socioeconomic factors observed in 1961, 1971, 1981, and 1991 and imputed, using a first-order random walk process, in 1956, 1966, 1976, 1986, and 1996. Further, we considered a temporal lag of 10 years between these and mortality.

4. Results

To compare birth-cohorts 1905–15 and 1930–40, we describe estimates for space-time interaction terms from the Lagazio et al. [8] model and the specific terms from the proposed model (SMBC). We focused on these birth-cohorts since their comparison seemed to be the most meaningful. In fact, the comparison highlights differences before and after Second World War. In addition, 1905–15 and 1930–40 are the youngest and oldest cohorts for which age-classes between 60 and 75 years were observed. Consequently, such cohorts have a substantial observed number of events.

Figure 2 shows the estimates for the space-time interaction terms. The areas with the highest mortality (north-western part of the region) exhibit a decrease in relative risk, while the low-risk areas in the south-eastern part of the region behave in the opposite way, still showing an increase in relative risk.

(a)

(b)

The birth-cohorts 1905–15 and 1930–40 are also compared through the more parsimonious SMBC. Figure 3 shows relative risks from SMBC for the 1905–15 and 1930–40 birth-cohorts. Figure 4(a) shows the estimates of the exponentialized parameters , that is, the common spatial distribution of relative risk for lung cancer for males in Tuscany (Italy) for 1905–15 and 1930-40. There is a high-risk area in the north-western part of the region and a low-risk area near the south-eastern border. This pattern has a strong connection with the spatial distribution of socioeconomic characteristics over the region. Indeed, the heavily industrialized and developed areas within Tuscany are located in the river Arno valley, from the capital, Florence, to Pisa and the main port of Livorno. At the beginning of 1900, the first industrial settlements were mainly close to the mountain areas. In the map, this pattern corresponds to the north-western part of the region.

(a)

(b)

(a)

(b)

(c)

(d)

Figure 4(b) also shows the estimates of the parameters , namely, the relative importance of the common pattern for the two birth-cohorts 1905–15 and 1930–40. The box-plots suggest that the epidemic of lung cancer clearly declines.

Figures 4(c) and 4(d) illustrate the estimates of the exponentialized parameters , that is, the specific spatial distribution of relative risk for birth-cohorts 1905–15 and 1930–40. For the birth-cohort 1905–15, there are some high values for the cities of Florence, Siena, Pisa, Livorno, Grosseto and for the industrial municipalities of Massa, Carrara, and Piombino. On the contrary, there is a prominent pattern for the 1930–35 birth-cohort, with higher relative risks for the municipalities of Mount Amiata (southern Tuscany). Furthermore, a decrease in the risk is observed for the cities (with the exception of Arezzo that moves in the opposite direction), mainly for Florence and its metropolitan municipalities (Scandicci, Signa, and Lastra a Signa). High risks are also observed in the municipality of Colle Val D’Elsa.

Consequently, an opposite trend can be observed in Arezzo, Colle Val D’Elsa, and the Mount Amiata areas. Results from the two different models are comparable. The proposed model seems to highlight more directly and parsimoniously the differences between two and more cohorts. The maps for the specific components are not oversmoothed, since there is minor information from the data but also a less structured model.

After the Second World War, the mines of Mount Amiata were supplying 50% of the world production of Cinnabar (a mineral consisting of mercury) and they remained the top supplier for decades. Cinnabar, with its vermilion red color, was used mainly as a primary colorant of walls, fabrics, and crockery. Before the advent of mining on a large scale, the economy of almost all the municipalities in the territory of Mount Amiata was based on agro-forestry-pastoral activities. The arrival of the mines leads to a striking transformation of this type of economy. The mining centers of Abbadia San Salvatore, Castellazzara, Piancastagnaio, and Santa Fiora reported a marked improvement in living standards. But the two occupational diseases typical of mining and metallurgical activities linked to mercury were silicosis and mercury-poisoning. Silicosis is a pneumoconiosis caused by prolonged inhaling (20–30 years) of dust containing small crystalline free silica particles. In the seventies, the mercury crisis hit the mines of Mount Amiata. In any case, the 1930–40 cohort of workers had been exposed for a long period from the postwar era when the mine was “in the fullest activity” to the closure of the mines in 1970, that is, 20–30 years of exposure. Moreover, exposure to crystalline free silica also involved the workers in the gold casting sector of Valtiberina and Arezzo and the glassware sector of Colle Val D’Elsa.

The units most at risk include those working in the metal mines, ceramics and glass factories, foundries, and the sandstone mining and granite industries. Workers with high dust exposure who develop silicosis are at increased risk of lung cancer; however, high levels of exposure to silica dust on its own are significant in the pathogenesis of lung cancer and silicosis is coincidental.

5. Discussion and Conclusions

In this paper, the space-time epidemic of lung cancer in Tuscany (Italy) is investigated. The analysis is performed via a hierarchical Bayesian model with shared and cohort-specific risk components.

It should be mentioned that this work is mainly a descriptive epidemiological investigation. In addition, three possible criticisms are the ecological bias, the misleading specification of the deprivation measure, and considering the lag of association between deprivation and mortality of 10 years. Despite such (and possibly other) limitations, our analysis yields some information on the space-time behavior of lung cancer in Tuscany and gives a few hints on the etiology of this disease.

It should be also noted that socioeconomic factors, when used at aggregate level, could give rise to the problem of “ecological fallacy.” However, they allow for identifying homogeneous groups, to study large populations at low costs, and to address questions of environmental health that might be difficult or impossible to study using other approaches (see [42], for a review of ecological analysis).

Lung cancer mortality for males in Tuscany has increased over the last 40-year period (1971–2010), showing a strong epidemic curve with the highest rates for the birth-cohorts born between 1920 and 1930. The geographical distribution of relative risk is very spatially structured and the north-western areas have a higher risk (see also [43]) due to the location of industrial plants [44]. The resulting pattern has an interesting interpretation as it highlights small towns with industrial plants in the north-west of Tuscany, the Massa-Carrara, Piombino, and the Island of Elba (with coke and iron industries since 1905). The spatial patterns suggest a role played by occupational and environmental factors. For males, the specific high-risk municipalities correspond to early industrialized areas. Lung cancer epidemics are now decreasing in all the municipalities of Tuscany. The overall time trend is moving toward greater homogeneity among areas. However, some areas (the municipalities of Mount Amiata, Arezzo, and Colle Val D’Elsa) show the opposite behavior for the younger cohorts: a new occupational risk factor, namely, silica, could be the cause of lung cancer in workers in mines and the gold and glass sectors.

The role of silicosis in the development of lung cancer, associated with silica exposure, remains controversial. Some studies indicate an increased risk of lung cancer from exposure to silica, and others limit such association to individuals with silicosis, while others show no association at all. The association has been studied for many decades. See Lynge et al. [45], Hnizdo and Sluis-Cremer 1991 [46], Amandus and Costello [47], NIOS [48], EC 2006 [49], Maciejewska [50], Vida et al. [51], MH [52], Chen et al. [53], and Liu et al. 2013 [54]. Although controversial, the decision taken in 1997 by the International Agency for Research on Cancer (IARC) [55] classified the respirable fraction of crystalline silica as a carcinogen of the first group for humans. This was followed by several epidemiological studies in different industrial fields, meta-analyses and reviews. However, unequivocal conclusions were not always reached. In 2002, the Italian Silica Network (consisting of the regions, provinces, INAIL, ISPESL, National Institute of Health, local health authorities, and scientific research centers) was set up in the aim of intervening in the field of prevention which lacked the required completeness and clarity of legislation at a European level. An international estimation conducted in 2003 by CAREX (CARcinogen EXposure, Canada) suggested that in Italy more than 254 thousand people were exposed to word-related crystalline silica risks, 12–15 thousand of whom in Tuscany. In 2004, the Region of Tuscany launched a project aimed at the prevention, measurement, and intervention in critical sectors such as quarries, construction sites, glass factories, and cement plants. Silicosis, the terrible plague which in the past, and only until a few decades ago, accompanied the industrial and economic development of our country, is not almost completely eradicated. However, still far from being won, in Italy as in the rest of Europe, is the battle against another big risk that those who work with crystalline free silica are exposed to lung cancer.

Conflicts of Interest

The author declares that they have no conflicts of interest.

Acknowledgments

The author would like to express their gratitude to Mariangela Vigotti (University of Pisa, Italy) and to Elisabetta Chellini (Istituto per lo Studio e la Prevenzione Oncologica, Firenze, Italy) for kindly making available the datasets 1971–1994 and 1995–2010, respectively.

References

J. Besag, J. York, and A. Mollié, “Bayesian image restoration, with two applications in spatial statistics,” Annals of the Institute of Statistical Mathematics, vol. 43, no. 1, pp. 1–59, 1991.
View at: Publisher Site | Google Scholar | MathSciNet
L. Bernardinelli, D. Clayton, C. Pascutto, C. Montomoli, M. Ghislandi, and M. Songini, “Bayesian analysis of space—time variation in disease risk,” Statistics in Medicine, vol. 14, no. 21-22, pp. 2433–2443, 1995.
View at: Publisher Site | Google Scholar
L. A. Waller, B. P. Carlin, H. Xia, and A. E. Gelfand, “Hierarchical Spatio-Temporal Mapping of Disease Rates,” Journal of the American Statistical Association, vol. 92, no. 438, pp. 607–617, 1997.
View at: Publisher Site | Google Scholar
R. M. Assunção, I. A. Reis, and C. D. L. Oliveira, “Diffusion and prediction of Leishmaniasis in a large metropolitan area in Brazil with a Bayesian space-time model,” Statistics in Medicine, vol. 20, no. 15, pp. 2319–2335, 2001.
View at: Publisher Site | Google Scholar
D. Sun, R. K. Tsutakawa, H. Kim, and Z. He, “Spatio-temporal interaction with disease mapping,” Statistics in Medicine, vol. 19, no. 15, pp. 2015–2035, 2000.
View at: Publisher Site | Google Scholar
L. Knorr-Held and J. Besag, “Modelling risk from a disease in time and space,” Statistics in Medicine, vol. 17, no. 18, pp. 2045–2060, 1998.
View at: Publisher Site | Google Scholar
L. Knorr-Held, “Bayesian modelling of inseparable space-time variation in disease risk,” Statistics in Medicine, vol. 19, no. 17-18, pp. 2555–2567, 2000.
View at: Publisher Site | Google Scholar
C. Lagazio, E. Dreassi, and A. Biggeri, “A hierarchical Bayesian model for space-time variation of disease risk,” Statistical Modelling, vol. 1, no. 1, pp. 17–29, 2001.
View at: Publisher Site | Google Scholar
P. Boyle and J. Ferlay, “Cancer incidence and mortality in Europe, 2004,” Annals of Oncology, vol. 16, no. 3, pp. 481–488, 2005.
View at: Publisher Site | Google Scholar
F. I. Bray and E. Weiderpass, “Lung cancer mortality trends in 36 European countries: Secular trends and birth cohort patterns by sex and region 1970-2007,” International Journal of Cancer, vol. 126, no. 6, pp. 1454–1466, 2010.
View at: Publisher Site | Google Scholar
D. R. Youlden, S. M. Cramb, and P. D. Baade, “The international epidemiology of lung cancer: geographical distribution and secular trends,” Journal of Thoracic Oncology, vol. 3, no. 8, pp. 819–831, 2008.
View at: Publisher Site | Google Scholar
F. Bray, J. E. Tyczynski, and D. M. Parkin, “Going up or coming down? The changing phases of the lung cancer epidemic from 1967 to 1999 in the 15 European Union countries,” European Journal of Cancer, vol. 40, no. 1, pp. 96–125, 2004.
View at: Publisher Site | Google Scholar
K. Shibuya, M. Inoue, and A. D. Lopez, “Statistical modeling and projections of lung cancer mortality in 4 industrialized countries,” International Journal of Cancer, vol. 117, no. 3, pp. 476–485, 2005.
View at: Publisher Site | Google Scholar
P. Armitage and R. Doll, “The age distribution of cancer and a multi-stage theory of carcinogenesis,” British Journal of Cancer, vol. 8, no. 1, pp. 1–12, 1954.
View at: Publisher Site | Google Scholar
B. J. Glick, “The Spatial Organization of Cancer Mortality,” Annals of the Association of American Geographers, vol. 72, no. 4, pp. 471–481, 1982.
View at: Publisher Site | Google Scholar
J. D. Mayer, “The role of spatial analysis and geographic data in the detection of disease causation,” Social Science & Medicine, vol. 17, no. 16, pp. 1213–1221, 1983.
View at: Publisher Site | Google Scholar
A. E. J. M. Cavelaars, A. E. Kunst, J. P. Mackenbach et al., “Educational differences in smoking: international comparison,” British Medical Journal, vol. 320, no. 7242, pp. 1102–1107, 2000.
View at: Publisher Site | Google Scholar
P. Townsend, “Deprivation,” Journal of Social Policy, vol. 16, no. 2, pp. 125–146, 1987.
View at: Publisher Site | Google Scholar
P. Townsend, P. Phillimore, and A. Beattie, Health and Deprivation: Inequalities and the North, Croom Helm, London, UK, 1988.
B. Jarman, “Identification of underprivileged areas,” British Medical Journal, vol. 286, no. 6379, pp. 1705–1709, 1983.
View at: Publisher Site | Google Scholar
R. Morris and V. Carstairs, “Which deprivation? a comparison of selected deprivation indexes,” Journal of Public Health, vol. 13, no. 4, pp. 318–326, 1991.
View at: Google Scholar
V. Carstairs, “Deprivation indices: their interpretation and use in relation to health,” Journal of Epidemiology and Community Health, vol. 49, 2, pp. S3–S8, 1995.
View at: Publisher Site | Google Scholar
S. St Leger, “Use of deprivation indices in small area studies of environment and health,” Journal of Epidemiology & Community Health, vol. 49, no. 2, pp. 1–88, 1995.
View at: Google Scholar
J. Eachus, M. Williams, P. Chan et al., “Deprivation and cause specific morbidity: Evidence from the Somerset and Avon survey of health,” British Medical Journal, vol. 312, no. 7026, pp. 287–292, 1996.
View at: Publisher Site | Google Scholar
G. Davey Smith, C. Hart, G. Watt, D. Hole, and V. Hawthorne, “Individual social class, area-based deprivation, cardiovascular disease risk factors, and mortality: the Renfrew and Paisley study,” Journal of Epidemiology and Community Health, vol. 52, no. 6, pp. 399–405, 1998.
View at: Publisher Site | Google Scholar
J. P. Pell, A. C. H. Pell, J. Norrie, I. Ford, and S. M. Cobbe, “Effect of socioeconomic deprivation on waiting time for cardiac surgery: Retrospective cohort study,” British Medical Journal, vol. 320, no. 7226, pp. 15–18, 2000.
View at: Publisher Site | Google Scholar
E. Dreassi, “A space-time analysis of the relationship between material deprivation and mortality for lung cancer,” Environmetrics, vol. 14, no. 5, pp. 511–521, 2003.
View at: Publisher Site | Google Scholar
A. Biggeri, E. Dreassi, and M. Marchi, “A multilevel Bayesian model for contextual effect of material deprivation,” Statistical Methods & Applications. Journal of the Italian Statistical Society, vol. 13, no. 1, pp. 89–103, 2004.
View at: Google Scholar | MathSciNet
D. Catelan, A. Biggeri, E. Dreassi, and C. Lagazio, “Space-cohort Bayesian models in ecological studies,” Statistical Modelling. An International Journal, vol. 6, no. 2, pp. 159–173, 2006.
View at: Publisher Site | Google Scholar | MathSciNet
E. Dreassi, A. Biggeri, and D. Catelan, “Space-time models with time-dependent covariates for the analysis of the temporal lag between socioeconomic factors and lung cancer mortality,” Statistics in Medicine, vol. 24, no. 12, pp. 1919–1932, 2005.
View at: Publisher Site | Google Scholar | MathSciNet
L. Knorr-Held and N. G. Best, “A shared component model for detecting joint and selective clustering of two diseases,” Journal of the Royal Statistical Society: Series A (Statistics in Society), vol. 164, no. 1, pp. 73–85, 2001.
View at: Publisher Site | Google Scholar | MathSciNet
L. Held, I. Natário, S. E. Fenton, H. v. Rue, and N. Becker, “Towards joint disease mapping,” Statistical Methods in Medical Research, vol. 14, no. 1, pp. 61–82, 2005.
View at: Publisher Site | Google Scholar | MathSciNet
E. Dreassi, “Polytomous disease mapping to detect uncommon risk factors for related diseases,” Biometrical Journal, vol. 49, no. 4, pp. 520–529, 2007.
View at: Publisher Site | Google Scholar | MathSciNet
S. Richardson, J. J. Abellan, and N. Best, “Bayesian spatio-temporal analysis of joint patterns of male and female lung cancer risks in Yorkshire (UK),” Statistical Methods in Medical Research, vol. 15, no. 4, pp. 385–407, 2006.
View at: Publisher Site | Google Scholar | MathSciNet
A. Biggeri, D. Catelan, and E. Dreassi, “The epidemic of lung cancer in Tuscany (Italy): A joint analysis of male and female mortality by birth cohort,” Spatial and Spatio-Temporal Epidemiology, vol. 1, no. 1, pp. 31–40, 2009.
View at: Publisher Site | Google Scholar
Y. C. MacNab, “On Bayesian shared component disease mapping and ecological regression with errors in covariates,” Statistics in Medicine, vol. 29, no. 11, pp. 1239–1249, 2010.
View at: Publisher Site | Google Scholar | MathSciNet
D. G. Clayton, “Generalized linear mixed models,” in Markov chain Monte Carlo in practice, Interdiscip. Statist., pp. 275–301, Chapman & Hall, London, UK, 1996.
View at: Google Scholar | MathSciNet
A. G. Gelman and D. B. Rubin, “Inference from iterative simulation using multiple sequences,” Statistical Science, vol. 7, no. 4, pp. 457–472, 1992.
View at: Publisher Site | Google Scholar
D. Lee, “A comparison of conditional autoregressive models used in Bayesian disease mapping,” Spatial and Spatio-Temporal Epidemiology, vol. 2, no. 2, pp. 79–89, 2011.
View at: Publisher Site | Google Scholar
M. A. Vigotti, A. Biggeri, E. Dreassi, M. A. Protti, and C. Cislaghi, Atlas of Mortality in Tuscany 1971-1994, Università degli studi di Pisa: Edizioni Plus, 2001.
D. Clayton and E. Schifflers, “Models for temporal variation in cancer rates. I: Age–period and age–cohort models,” Statistics in Medicine, vol. 6, no. 4, pp. 449–467, 1987.
View at: Publisher Site | Google Scholar
H. Morgenstern, “Ecologic studies,” in Modern Epidemiology, K. J. Rothman and S. Greenland, Eds., pp. 459–480, Philadelphia, Penn, USA, Lipincott-Raven, 2nd edition, 1998.
View at: Google Scholar
A. Biggeri, M. Bonannini, D. Catelan, F. Divino, E. Dreassi, and C. Lagazio, “Bayesian ecological regression with latent factors: atmospheric pollutants emissions and mortality for lung cancer,” Environmental and Ecological Statistics, vol. 12, no. 4, pp. 397–409, 2005.
View at: Publisher Site | Google Scholar | MathSciNet
G. Mori, “Dall’unità alla guerra: aggregazione e disgregazione di un’area regionale,” in Storia d’Italia. Le Regioni dall’Unità a Oggi, G. Mori, Ed., Einaudi, La Toscana, Italy, 1986.
View at: Google Scholar
E. Lynge, K. Kurppa, L. Kristofersen, H. Malker, and H. Sauli, “Silica dust and lung cancer: Results from the Nordic occupational mortality and cancer incidence registers,” Journal of the National Cancer Institute, vol. 77, no. 4, pp. 883–889, 1986.
View at: Google Scholar
E. Hnizdo and G. K. Sluis-Cremer, “Silica exposure, silicosis, and lung cancer: A mortality study of South African gold miners,” British Journal of Industrial Medicine, vol. 48, no. 1, pp. 53–60, 1991.
View at: Google Scholar
H. Amandus and J. Costello, “Silicosis and lung cancer in u.S. Metal miners,” Archives of Environmental & Occupational Health, vol. 46, no. 2, pp. 82–89, 1991.
View at: Publisher Site | Google Scholar
NIOS: National Institute for Occupational Safety and Health, Health Effects of Occupational Exposure to Respirable Crystalline Silica, National Institute for Occupational Safety and Health, Washington, DC, USA, 2002.
“EC: European Commission. Brussels, Belgium: European Commission; 2006”.
View at: Google Scholar
A. Maciejewska, “Occupational exposure assessment for crystalline silica dust: Approach in poland and worldwide,” International Journal of Occupational Medicine and Environmental Health, vol. 21, no. 1, pp. 1–23, 2008.
View at: Publisher Site | Google Scholar
S. Vida, J. Pintos, M.-É. Parent, J. Lavoué, and J. Siemiatycki, “Occupational exposure to silica and lung cancer: Pooled analysis of two case-control studies in Montreal, Canada,” Cancer Epidemiology, Biomarkers & Prevention, vol. 19, no. 6, pp. 1602–1611, 2010.
View at: Publisher Site | Google Scholar
MH: Ministry of Health, China’s Health Statistics Yearbook 2011, Peking Union Medical College Press, Beijing, China, 2011.
W. Chen, Y. Liu, H. Wang et al., “Long-term exposure to silica dust and risk of total and cause-specific mortality in Chinese workers: a cohort study,” PLoS Medicine, vol. 9, no. 4, article e1001206, 2012.
View at: Publisher Site | Google Scholar
Y. Liu, K. Steenland, Y. Rong et al., “Exposure-response analysis and risk assessment for lung cancer in relationship to silica exposure: A 44-year cohort study of 34,018 workers,” American Journal of Epidemiology, vol. 178, no. 9, pp. 1424–1433, 2013.
View at: Publisher Site | Google Scholar
IARC: International Agency for Research on Cancer (1997) Monographs on the Evaluation of Carcinogenic Risks to Humans, Vol. 68: Silica, Some Silicates, Coal Dust and Para-aramid Fibrils. Lyon, France:International Agency for Research on Cancer.

Copyright

Copyright © 2018 Emanuela Dreassi. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

PDF Download Citation

Download other formats

Order printed copies

Views

1563

Downloads

1191

Citations