The population-based mammography screening program (MSP) was implemented by the end of 2005 in Germany, and all women between 50 and 69 years are actively invited to a free biennial screening examination. However, despite the expected benefits, the overall participation rates range only between 50 and 55 %. There is also increasing evidence that belonging to a vulnerable population, such as ethnic minorities or low income groups, is associated with a decreased likelihood of participating in screening programs. This study aimed to analyze in more detail the intra-urban variation of MSP uptake at the neighborhood level (i.e. statistical districts) for the city of Dortmund in northwest Germany and to identify demographic and socioeconomic risk factors that contribute to non-response to screening invitations.
Methods
The numbers of participants by statistical district were aggregated over the three periods 2007/2008, 2009/2010, and 2011/2012. Participation rates were calculated as numbers of participants per female resident population averaged over each 2-year period. Bayesian hierarchical spatial models extended with a temporal and spatio-temporal interaction effect were used to analyze the participation rates applying integrated nested Laplace approximations (INLA). The model included explanatory covariates taken from the atlas of social structure of Dortmund.
Results
Generally, participation rates rose for all districts over the time periods. However, participation was persistently lowest in the inner city of Dortmund. Multivariable regression analysis showed that migrant status and long-term unemployment were associated with significant increases of non-attendance in the MSP.
Conclusion
Low income groups and immigrant populations are clustered in the inner city of Dortmund and the observed spatial pattern of persistently low participation in the city center is likely linked to the underlying socioeconomic gradient. This corresponds with the findings of the ecological regression analysis manifesting socioeconomically deprived neighborhoods as risk factors for low attendance in the MSP. Spatio-temporal surveillance of participation in cancer screening programs may be used to identify spatial inequalities in screening uptake and plan spatially focused interventions.
The implementation of a nation-wide, population-based mammography screening program (MSP) started in Germany by the end of the year 2005. The stepped implementation process was completed in the state of North-Rhine Westphalia in 2009. All resident women aged between 50 and 69 years are actively invited to a mammography screening examination every two years. The participation is voluntary and free of cost. Mammography screening is a procedure of secondary cancer prevention with an aim of detecting breast cancer in early stages where therapy is less invasive (e.g., breast-conserving therapies instead of mastectomy), remaining lifespans are extended and, ideally, breast cancer mortality is reduced. Despite these expected benefits and the free provision by all statutory health insurances, the overall participation rates range only between 50 and 55 % [1]. Population-based surveys demonstrate that significant gaps exist in screening mammography uptake across population subgroups [2]. These differences are believed to substantially contribute to a higher prevalence of late stage breast cancer at diagnosis among vulnerable populations, including racial and ethnic minorities [3] and low-income groups [4, 5]. More specifically, living in an economically deprived neighborhood showed a decreased likelihood of participating in cancer screening programs and an increased risk of a late-stage breast cancer diagnosis with the correspondent unfavorable prognosis [2, 6–9]. To date, few studies analyzed the intra-urban variation of participation rates in mammography screening programs [2, 10] and to our knowledge none investigated the situation in Germany. We suggest that small-area analyses may provide important insights into the processes and factors that are associated with low participation rates and that this may help to develop spatially focused approaches to improve the participation rates in disadvantaged neighborhoods.
Therefore, this study aimed to investigate the spatio-temporal distribution of the participation rates in the mammography screening program at the neighborhood-level (e.g. statistical districts) of a large city in Germany and to identify important demographic and socio-economic factors that influence the non-attendance to screening invitations.
Methods
Study region
Dortmund is a city in the federal state of North Rhine-Westphalia in northwestern Germany with a total population of 575 944 inhabitants in 2013. It is the largest city by area and population in the Ruhr district, a metropolitan area with some 5.1 million inhabitants which is the largest urban and industrial agglomeration in Germany. Dortmund is divided into 62 statistical districts with a median female population of 4086 inhabitants per statistical district. The city’s population is characterized by a high proportion of immigrants from southeast Europe and Turkey. The coal crisis in the end of 1970 led to a massive reduction of jobs in the coal and steel industry which resulted in high unemployment rates to this day. Immigrants and their descendants grew up in more socially deprived neighborhoods than many of the autochthonous population. This resulted in a strong spatial segregation of populations with a migration background and with low economic status within the city [11].
Participation rates and geo-referencing
The statistical districts were used as the geographical reference system in which the MSP participation rates were assessed. The residence addresses of all MSP participants for the years 2005 to 2013 were stored at KV.it Dortmund, the institution which administrates the MSP documentation software MaSc [12]. KV.it assigned MSP participants to one of the 62 statistical districts by linking their home addresses to a comprehensive list of street addresses for each district. A list of individual, anonymized participants who were geo-referenced to one of the statistical districts was then transferred to the Institute of Epidemiology and Social Medicine at the University Münster [13] where all subsequent analyses were carried out.
The years 2005 and 2006 were excluded from the present analyses to avoid contamination with the various organizational aspects of the stepped-up implementation of the MSP. All eligible women receive a biennial invitation to the screening program, hence, we chose to analyze three 2-year periods: 2007/2008, 2009/2010, and 2011/2012. The participation rates were calculated using the aggregated numbers of participants and the averaged female background population (age group 50–69) for each two year-period.
Spatio-temporal mapping and regression
The spatio-temporal distribution of participation rates was analyzed within a hierarchical Bayesian framework using a multivariate binomial regression model (spatio-temporal odds model): Let nit denote the number of eligible women resident in district i and period t and Yit the number of participants in breast cancer screening, with I = 1, …, 62 and t = 1, 2, 3. We assumed that the observed number of participants (Yit) had a binomial distribution with parameters nit and θit (probability of participation). At a second level, the probability of participation θit was then decomposed on the logit scale into an overall participation rate (α), main spatial effects (ui and vi) (constant in time), main temporal effects (unstructured (Φt) and structured (γt)), and a space-time interaction term (ψit).
The proposed space-time models, assuming a nonparametric time trend and a spatio-temporal interaction term, were introduced by Knorr-Held [14] and are an extension of the spatial model introduced by Besag et al. [15]. All model terms were treated as random variables: The spatially unstructured random effect (ui) was considered independent and identically distributed (iid) with zero mean and unknown precision (τu). To account for the assumption of correlated participation rates in nearby statistical districts, the spatially structured effect (vi) is modelled for each 62 districts as an intrinsic Gaussian Markov random field with unknown precision (τv). This specification is also called a conditionally autoregressive (CAR) prior and was introduced by Besag et al.[15]. In order to insure the identifiability of the intercept α (overall participation rate), a sum-to-zero constraint was imposed on the vi’s [16]. The unstructured temporal effect (Φt) was also modelled iid with zero mean and unknown precision. For the structured time effect (γt) random walks of first order were considered [17, 18]. The interaction term (ψit) can be specified in several ways [14], here it is assumed that the two unstructured effects (vi and γt) interact [17, 18]. Therefore, the interaction effect was also specified as zero mean normal with unknown precision (iid., i.e. ψit ~ N(0, τψ). The distribution of the hyperpriors was specified as follows: Minimally informative priors were specified on the log of the unstructured effect precision (log τv ~ logGamma (1, 0.001)) and on the log of the structured effect precision (log τu ~ logGamma (1,0.001)). For the unstructured time effect, a log τφ ~ logGamma (1, 0.01) hyperprior was chosen. For the structured temporal effect and the interaction term, minimally informative priors (the default priors): log τϒ, log τψ ~ logGamma (1, 0.00005) have been used. Altogether, the distribution of the hyperpriors resembles the ones used by Ugarte et. al [18].
The specified model (Equation 3) was extended to βxTi, where xTi contains the covariates with a space-time index in order to investigate potential risk factors associated with spatio-temporal variations in the participation rates. The covariates were taken from the atlas of social structure (Sozialstrukturatlas) of Dortmund which is a collection of administratively collected data reflecting social inequalities and differences in the population [19]. These are grouped into the dimensions employment status, demography, income, welfare, and housing. A full description of the explanatory variables is given in Table 1. In order to account for multicollinearity, an initial correlation matrix was examined for high correlations among the variables. Variables with a high correlation (>0.8) were excluded from further regression analyses. For the three 2-year time periods, the data of 2008, 2010 and 2012, respectively, were included in the model, and all variables were dichotomized according to their median value. Following the suggestions of Rothman et al. [20-21] each covariate was fitted separately and model fit was assessed using the changes in deviance information criterion (DIC) (smaller values of DIC indicating more explained variance and better fit). A multivariable model was fit by selectively including variables, starting with those that showed the lowest DIC in univariable analyses, until the DIC could be no further reduced. For the Bayesian inference, the integrated nested Laplace approximation (INLA) approach was used as introduced by Rue et al. [22] and implemented in the R package R-INLA [21, 23, 24]. The Bayesian inference was also used to report the resulting odd ratios (OR) as point estimate (posterior mean) and 95 % credibility intervals (CI) as a quantification of parameter uncertainty. All computations and visualizations were done in R v. 3.0.2 [25].
Table 1
Summary statistics of included variables in the 62 statistical districts and the three time periods
Dimension/Indicator
Variable
Median value (2008; 2010; 2012)
Definition
Employment/Unemployment
Employment rate [%]
48.2; 48.4; 50.1
Proportions of regular employees of the employable population (>15 to <65 years) with primary residence in Dortmund.
Employment trend [%]
3.5; 5.9; 7
Number of regular employees with primary residence in Dortmund in a five-year trend/comparison in percent.
Unemployment rate (foreigner) [%]
12.6; 12.5; 12.7
Proportion of unemployed, foreign persons of the employable, foreign population (>15 to <65 years) in percent.
Unemployment rate (<25 aged) [%]
4.2; 4.5; 4.3
Proportion of unemployed young persons (15 to < 25 years) of the population in the same age group.
Unemployment rate (long-term) [%]
44.8; 41.3; 44.7
Proportion of long-term unemployed persons (>12 months) to all unemployed persons in percent.
Demography
Female population [%]
51.7; 51.5; 51.5
Proportion of women to all inhabitants with primary residence in Dortmund.
Foreign residents [%]
7.5; 7.1; 7.7
Proportion of persons without German nationality to all inhabitants with primary residence in Dortmund.
Youth quotient
20.5; 20; 19
Number of persons aged < 15 years per 100 persons aged 15 to < 65 years.
Quotient of elderly
32.7; 32.7; 31.6
Number of persons aged > 65 years per 100 persons aged 15 to < 65 years.
Birth rate
8; 7.4; 7.5
Number of births per 1000 inhabitants with primary residence in Dortmund by December 31. of each year.
Mortality rate
10.4; 10.5; 10.5
Number of deaths per 1000 inhabitants with primary residence in Dortmund by December 31. of each year.
Social affairs
Long-term social welfare (“Hartz IV”) rate [%]
10.3; 11.9; 11.6
Proportions of unemployed (>12 months) persons (<65 years) with demand for basic financial benefits (Arbeitslosengeld II) of the residential population.
Living/habitation
Social housing rate [%]
8; 5.4; 6.2
Proportion of council-sponsored apartments on all rental appartements in percent.
Living space per person
39.6; 40.2; 40.4
The total living space in [m2] divided by the number of inhabitants with primary and secondary residence in Dortmund.
Ethics statement
KV.it [12] administrates the MSP documentation including the storage and management of the MSP participants consistent with the existing data protection legislation. For this study, KV.it aggregated participants in statistical districts so that individual women could be not identified. Data were transferred to the investigators in an anonymized form. Use of anonymized data for research purposes does not require a vote by ethics committee or an institutional review board.
Results
Mapping participation rates
The observed annual participation rates showed the overall biennial pattern of participation, i.e., one year with high number of participants and the subsequent year with lower participants (Fig. 1). Despite rising overall MSP participation rates for the three periods from 48 % (2007/08) over 50 % (2009/10) to 54 % (2011/12) (Referenzzentrum MS), a concentration of statistical districts with low participation rates persisted in the inner city of Dortmund, while the outer districts had consistently higher participation rates (Fig. 2a–c). The modeled time trends in Fig. 2d demonstrate these increasing participation rates over the three periods, where the structured time effect (γt) was more pronounced than the unstructured time effect (Φt). The spatial trends combine structured and unstructured effects (Fig. 2e) and confirm, after accounting for covariates, that lower participation rates cluster in the inner city. Finally, the interaction analyses reveals a clear space-time pattern (Fig. 2f-h) which indicates that in 2007/8, on a generally low participation level, the participation rates were particularly low in the eastern districts. This changed in 2008/9 as lower participation persisted in the western and central parts of the city, while in 2010/11 a low participation rate was found in only one inner city district.
Fig. 1
Yearly distribution of the participation rates over the period from 2007 to 2012
Fig. 2
Participation rates and random effects in the final spatio-temporal regression model. Spatial pattern of the biennial participation rates, for the periods: 2007–08 (a), 2009–10 (b), and 2011–12 (c). Black dots mark the location of the screening units in the study region. Odds ratios compared to the intercept (α) of the unstructured (Φt) and structured (γt) temporal effect (d), combined unstructured and structured spatial heterogeneity (ui + vi) (e), and spatio-temporal interaction effect (ψit) for 2007–08 (f), 2009–10 (g), and 2011–12 (h). All random effects were classified according to their quantiles
Regression analyses
Due to high correlations and close content relations, the variables: unemployment rate (total), persons with migrational background, and basic social welfare rate were excluded from the further analyses. The results of uni- and multivariable spatio-temporal regression analyses are summarized in Table 2. The odds ratios of the univariate spatio-temporal regression analyses demonstrate clearly that districts with high proportions of unemployed migrants or long-term unemployed residents showed a statistically significant impact on lower participation rates. In contrast, higher proportion of elderly population showed positive association with the participation rates. The other ecological variables clearly contained null values in their credibility intervals and were therefore considered as factors without relevant influence. In the multivariable analyses that simultaneously adjusted for spatio-temporal variation, the negative association of the proportion of unemployed migrants and long-term unemployed remained statistically significant. Proportions of migrant unemployment and long-term unemployed above the median were associated with a significant 6 % and, 3 %, respectively, increase in risk for non-attendance in the mammography screening program.
Table 2
Estimation results of univariable and multivariable analysis using the spatio-temporal model
Variable
pD
DIC
Space-time CAR without covariates
95.94
1707.09
Explanatory variables
Univariable analysis
Multivariable analysis
Odds ratio (95 % CI)
DIC (rank)
Odds ratio (95 % CI)
DIC
Unemployment rate (foreigner) [%]
0.95
(0.91; 0.99)
1706.279
0.94
(0.90; 0.98)
1706.56
Living space per person
1.03
(0.95; 1.13)
1707.15
1.01
(0.92; 1.12)
Unemployment rate (long-term) [%]
0.98
(0.94; 1.0)
1707.174
0.97
(0.94; 1.0)
Employment trend [%]
0.98
(0.95; 1.01)
1707.211
Mortality rate
1.02
(0.98; 1.06)
1707.305
Female population [%]
1.03
(0.98; 1.09)
1707.432
Social housing rate [%]
1.00
(0.96; 1.05)
1708.321
Foreign residents [%]
0.95
(0.89; 1.02)
1708.38
Quotient of elderly
1.07
(0.99; 1.15)
1708.561
Basic social welfare rate [%]
0.97
(0.91; 1.04)
1708.566
Unemployment rate (<25 aged) [%]
0.99
(0.94; 1.04)
1708.711
Youth quotient
0.99
(0.95; 1.04)
1708.801
Employment rate [%]
1.02
(0.98; 1.07)
1708.951
Birth rate
0.97
(0.94; 1.01)
1709.389
Estimates are odds ratios with their associated 95 % Credible Intervals (CI). Values below one denoting a reduced odd of participating in screening, and vice versa
Discussion
The present study analyzed small-area, intra urban variations of participation rates of the MSP in the city of Dortmund over three 2-year periods. An overall increase in the participation rates was observed over the study period, while the increase was unevenly distributed across the study area. There was a spatial concentration of statistical districts, mainly in the city center, with persistently low participation rates. Dortmund is known to have a strong gradient of socio-economical segregation [26]. Its population is characterized by a high proportion of residents which are likely to have a lower socio-economic status, because of the massive loss of jobs in mining and steel industries in this area after 1970 (also known as structural crisis) which was compounded in population segments with a migration background [11]. Muller and Berger [26] reported, in their investigation about neighborhood deprivation and prevalence of type 2 diabetes in Dortmund, that the inner city and parts of the western city are characterized by the highest level of socio-economic constraints including high proportions of immigrants, unemployed residents, residents with basic social welfare as well as a high population density and low level of incomes.
Therefore, it seems plausible that the observed spatial pattern in the participation rates is linked to the underlying socio-economic gradient. Despite the high availability of screening facilities in the inner city (Fig. 2a–c), it should be noted that these statistical districts had persistently low participation rates. In contrast, reduced participation rates in the southernmost districts may be attributable to a higher proportion of women with private health insurances who tend to abstain from public health offers. The highest rates of participation were found in the eastern parts of the city, that has more affluent statistical districts and screening facilities more nearby. The ecological regression analyses confirmed the spatio-temporal results by revealing that characteristics of disadvantage in statistical districts were related to an increased probability of non-participation in the MSP.
The identification of socio-economic risk factors at area-level as explanatory variables of non-attendance in mammography screening have been examined in previous studies [27–30] which also found that neighborhood income was an important determinant of participation [31–33]. Additionally, Peek and Han [34] reported that vulnerable groups such as the poor, the elderly, and minorities were often unaware of mammography screening programs and had a reduced awareness and a lack of information of disease prevention, diagnosis, and treatment [35]. Awareness of the program is unlikely to play a major role in Germany as all resident women were personally invited as part of a structured systematic program of early breast cancer detection; attitude towards disease prevention seems to be a more likely reason for the lower MSP attendance.
Regarding the overall increase in the participation rates over time, it should be kept in mind that the MSP is a rather recent program as compared to other European countries [36]. It became operational as part of routine care only by the end of the year 2005, and comprehensive implementation in North Rhine–Westphalia was not completed before the end of 2009 [37]. Therefore, the increasing participation rates over the study period may be mainly attributed to an increased efficiency of operational routines within the screening units which allowed for a growing numbers of screened women and it is probable that these observed spatial effects - especially within the first two study periods - are influenced by these technical and structural developments [38]. The spatio-temporal interaction effect adds to the spatial and temporal findings in that it identifies districts where the observed participation rates were reduced as compared to the entire city and throughout the study period [39]. Holding the spatial component constant confirmed the increasing overall trend of the participation rates, while holding the temporal trend constant confirmed that statistical districts in the city center had consistently reduced attendance rate in mammography screening.
This study has several strengths and limitations. An obvious strength has been the use of the Bayesian hierarchical framework in order to borrow strength from spatial and temporal neighbors to reduce the high variability inherent in the estimators, in particular, when numbers (disease counts and/or background population) are unstable [18, 39, 40]. Also, the inclusion of a space-time interaction effect is an added strength of this study, because the participation rates of mammography screening may be plausibly assumed to be dependent in space and time. Adjustment for the spatial, temporal, and space-time interaction effects depicts more clearly how the spatial pattern of the participation rates evolved over time, while the intersection of space and time is seldom considered to disentangle the complex determinants of health-related behavior and diseases [39]. Furthermore, the use of a non-parametric time trend has been a more plausible assumption than a linear time trend, because not all census tracts showed a linear increase or decrease in their participation rates. However, the analysis of effects in our study was confined to only three time periods, and hence requires caution in interpretation. The use of integrated nested Laplace approximations (INLA) reduced computing time substantially while attaining a high degree of accuracy, when fitting large, complex data sets at detailed geographic levels as used in spatio-temporal disease mapping. Given the inherently ecological nature of this study, the parameter estimates may not be used for making inferences on the individual level and therefore must not be interpreted causally. However, the results provide important hints to how social, cultural, and contextual factors may influence the attendance in mammography screening. Thus, despite the ecological nature of our study, the results may be used to provide spatially focused interventions to improve participation in disadvantaged city districts. Another limitation results from the fact that the precise number of invited women was not available due to data privacy regulations in Germany. The denominator used for calculating the participation rates contained therefore the whole female population in that age group which also included women non-eligible to screening (e.g., because of prevalent breast cancer). Thus, the participation rates may be slightly underestimated but this is not perceived as biasing the spatial associations. Another limitation results from the aggregation of the number of participants to a period of 24 month. Because the invitations to biennial screening were continuously mailed throughout the 22–26 month period, a certain amount of misclassifications is to be expected. However, as the general trend of the participation rates showed a clear biennial pattern it seems safe to assume that the main spatio-temporal process of the participation rates was captured with the temporal aggregation employed in this study.
Conclusions
This study analyzed the intra-urban participation rates of a mammography screening program within a hierarchical Bayesian framework using spatio-temporal disease models to identify regions and risk factors of low attendance. Despite a general temporal trend with increasing participation rates, spatial clustering of persistently lower participation rates was observed in the inner city districts, which are known as the socio-economically most deprived neighborhoods of Dortmund. This corresponds with the findings of the ecological regression analysis manifesting indicators of socio-economic constraint in a neighborhood as risk factors for low attendance in the MSP. The spatio-temporal interaction effect showed that the participation rates developed spatially unequally over time and that certain districts had low participation rates throughout the study period. Spatio-temporal surveillance of the participation rates and focused intervention could help identifying and reducing spatial inequalities in the uptake of mammography screening.
Declarations
Acknowledgements
We would like thank KV.it (Dortmund) for georeferencing the screening participants. Additionally, we like to thank the statistical department of Dortmund for supply of the geometrical boundary information, the background population, and the covariate data.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
Authors’ Affiliations
(1)
Institute of Epidemiology and Social Medicine, Medical Faculty, Westfälische Wilhelms-Universität Münster
(2)
Institute for Geoinformatics, Geosciences Faculty, Westfälische Wilhelms-Universität Münster
(3)
Reference Center for the Mammography Screening Program, University Hospital, Westfälische Wilhelms-Universität Münster
(4)
Epidemiological Cancer Registry North Rhine-Westphalia
References
Mammographie K. Evaluationsbericht 2011. Berlin: Zusammenfassung der Ergebnisse des Mammographie-Screening-Programms in Deutschland; 2014.Google Scholar
Lantz PM, Mujahid M, Schwartz K, Janz NK, Fagerlin A, Salem B, et al. The influence of race, ethnicity, and individual socioeconomic factors on breast cancer stage at diagnosis. Am J Public Health. 2006;96(12):2173–8. doi:10.2105/AJPH.2005.072132.View ArticlePubMedPubMed CentralGoogle Scholar
Merkin SS, Stevenson L, Powe N. Geographic socioeconomic status, race, and advanced-stage breast cancer in New York City. Am J Public Health. 2002;92(1):64–70.View ArticlePubMedPubMed CentralGoogle Scholar
Henry KA, Boscoe FP, Johnson CJ, Goldberg DW, Sherman R, Cockburn M. Breast Cancer Stage at Diagnosis: Is Travel Time Important? J Community Health. 2011;36(6):933–42. doi:10.1007/s10900-011-9392-4.View ArticlePubMedGoogle Scholar
Peipins LA, Graham S, Young R, Lewis B, Foster S, Flanagan B, et al. Time and distance barriers to mammography facilities in the Atlanta metropolitan area. J Community Health. 2011;36(4):675–83. doi:10.1007/s10900-011-9359-5.View ArticlePubMedGoogle Scholar
Dai D. Black residential segregation, disparities in spatial access to health care facilities, and late-stage breast cancer diagnosis in metropolitan Detroit. Health Place. 2010;16(5):1038–52. doi:10.1016/j.healthplace.2010.06.012.View ArticlePubMedGoogle Scholar
Stadt Dortmund. Sozialstrukturatlas 2005 - Demographische und soziale Struktur der Stadt Dortmund, ihrer Stadtbezirke und Sozialräume. In: Dezernat für Arbeit GuS, editor. Dortmund; 2005.Google Scholar
KV.IT Dortmund. IT-Gesellschaft für integrierte Services im Gesundheitswesen. 2015. http://www.kv-it-gmbh.de/. Accessed 04/21/2015.
Knorr-Held L. Bayesian modelling of inseparable space-time variation in disease risk. Stat Med. 2000;19(17–18):2555–67.View ArticlePubMedGoogle Scholar
Besag J, York J, Mollie A. Bayesian Image-Restoration, with 2 Applications in Spatial Statistics. Ann I Stat Math. 1991;43(1):1–20.View ArticleGoogle Scholar
Schrodle B, Held L. Spatio-temporal disease mapping using INLA. Environmetrics. 2010;22(6):725–34.View ArticleGoogle Scholar
Held L, Natario I, Fenton SE, Rue H, Becker N. Towards joint disease mapping. Stat Methods Med Res. 2005;14(1):61–82.View ArticlePubMedGoogle Scholar
Ugarte MD, Adin A, Goicoa T, Militino AF. On fitting spatio-temporal disease mapping models using approximate Bayesian inference. Stat Methods Med Res. 2014;23(6):507–30. doi:10.1177/0962280214527528.View ArticlePubMedGoogle Scholar
Rue H, Martino S, Chopin N. Approximate Bayesian inference for latent Gaussian models by using integrated nested Laplace approximations. J Biom Biostat. 2009;71(2):319–92. doi:10.1111/j.1467-9868.2008.00700.x.Google Scholar
Bivand RS, Gomez-Rubio V, Rue H. Spatial Data Analysis with R-INLA with Some Extensions. J Stat Softw. 2015;63(20):1–31.View ArticleGoogle Scholar
Lindgren F, Rue H. Bayesian Spatial Modelling with R-INLA. J Stat Softw. 2015;63(19):1–25.View ArticleGoogle Scholar
R Development Core Team. A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2015.Google Scholar
Dailey AB, Brumback BA, Livingston MD, Jones BA, Curbow BA, Xu X. Area-level socioeconomic position and repeat mammography screening use: results from the 2005 National Health Interview Survey. Cancer Epidemiol Biomarkers Prev. 2011;20(11):2331–44. doi:10.1158/1055-9965.EPI-11-0528.View ArticlePubMedGoogle Scholar
Dailey AB, Kasl SV, Holford TR, Calvocoressi L, Jones BA. Neighborhood-level socioeconomic predictors of nonadherence to mammography screening guidelines. Cancer Epidemiol Biomarkers Prev. 2007;16(11):2293–303. doi:10.1158/1055-9965.EPI-06-1076.View ArticlePubMedGoogle Scholar
Pornet C, Dejardin O, Morlais F, Bouvier V, Launoy G. Socioeconomic and healthcare supply statistical determinants of compliance to mammography screening programs: a multilevel analysis in Calvados. France Cancer epidemiology. 2010;34(3):309–15. doi:10.1016/j.canep.2010.03.010.View ArticlePubMedGoogle Scholar
von Euler-Chelpin M, Olsen AH, Njor S, Vejborg I, Schwartz W, Lynge E. Socio-demographic determinants of participation in mammography screening. Int J Cancer. 2008;122(2):418–23. doi:10.1002/ijc.23089.View ArticleGoogle Scholar
Kothari AR, Birch S. Individual and regional determinants of mammography uptake. Can J Public Health. 2004;95(4):290–4.PubMedGoogle Scholar
Maheswaran R, Pearson T, Jordan H, Black D. Socioeconomic deprivation, travel distance, location of service, and uptake of breast cancer screening in North Derbyshire, UK. J Epidemiol Community Health. 2006;60(3):208–12. doi:10.1136/jech.200X.038398.View ArticlePubMedPubMed CentralGoogle Scholar
Ouédraogo S, Dabakuyo-Yonli TS, Amiel P, Dancourt V, Dumas A, Arveux P. Breast cancer screening programmes: Challenging the coexistence with opportunistic mammography. Patient Educ Couns. 2014;97(3):410–7. doi:10.1016/j.pec.2014.08.016.
Peek ME, Han JH. Disparities in screening mammography. Current status, interventions and implications. J Gen Intern Med. 2004;19(2):184–94.View ArticlePubMedPubMed CentralGoogle Scholar
Chamot E, Charvet AI, Perneger TV. Who gets screened, and where: a comparison of organised and opportunistic mammography screening in Geneva, Switzerland. Eur J Cancer. 2007;43(3):576–84. doi:10.1016/j.ejca.2006.10.017.View ArticlePubMedGoogle Scholar
Heidinger O, Batzler WU, Krieg V, Weigel S, Biesheuvel C, Heindel W, et al. The incidence of interval cancers in the German mammography screening program: results from the population-based cancer registry in North Rhine-Westphalia. Dtsch Arztebl Int. 2012;109(46):781–7. doi:10.3238/arztebl.2012.0781.PubMedPubMed CentralGoogle Scholar
Bluekens AM, Karssemeijer N, Beijerinck D, Deurenberg JJ, van Engen RE, Broeders MJ, et al. Consequences of digital mammography in population-based breast cancer screening: initial changes and long-term impact on referral rates. Eur Radiol. 2010;20(9):2067–73. doi:10.1007/s00330-010-1786-7.View ArticlePubMedPubMed CentralGoogle Scholar
Papoila AL, Riebler A, Amaral-Turkman A, Sao-Joao R, Ribeiro C, Geraldes C, et al. Stomach cancer incidence in Southern Portugal 1998–2006: a spatio-temporal analysis. Biom J. 2014;56(3):403–15. doi:10.1002/bimj.201200264.View ArticlePubMedGoogle Scholar
