- Research article
- Open Access
- Open Peer Review
Deprivation-specific life tables using multivariable flexible modelling – trends from 2000–2002 to 2010–2012, Portugal
BMC Public Healthvolume 19, Article number: 276 (2019)
Completing mortality data by information on possible socioeconomic inequalities in mortality is crucial for policy planning. The aim of this study was to build deprivation-specific life tables using the Portuguese version of the European Deprivation Index (EDI) as a measure of area-level socioeconomic deprivation, and to evaluate mortality trends between the periods 2000–2002 and 2010–2012.
Statistics Portugal provided the counts of deaths and population by sex, age group, calendar year and area of residence (parish). A socioeconomic deprivation level was assigned to each parish according to the quintile of their national EDI distribution. Death counts were modelled within the generalised linear model framework as a function of age, deprivation level and calendar period. Mortality Rate Ratios (MRR) were estimated to evaluate variations in mortality between deprivation groups and periods.
Life expectancy at birth increased from 74.0 and 80.9 years in 2000–2002, for men and women, respectively, and to 77.6 and 83.8 years in 2010–2012. Yet, life expectancy at birth differed by deprivation, with, compared to least deprived population, a deficit of about 2 (men) and 1 (women) years among most deprived in the whole study period. The higher mortality experienced by most deprived groups at birth (in 2010–2012, mortality rate ratios of 1.74 and 1.29 in men and women, respectively) progressively disappeared with increasing age.
Persistent differences in mortality and life expectancy were observed according to ecological socioeconomic deprivation. These differences were larger among men and mostly marked at birth for both sexes.
Life tables provide information on mortality rates and probabilities of death for specific populations defined by geographical regions and/or periods of time. They are important demographic tools as they are the basis for the estimation of life expectancy at birth, an important indicator of population health and development. Many factors are known to influence overall mortality, such as age, sex, geographical region, socioeconomic deprivation or ethnicity [1,2,3,4]. While the effect of, for example, age is largely unavoidable, the gap in mortality due to socioeconomic characteristics could be reduced with policies oriented to improve population living conditions and to change the social and economic structures .
Deprivation is a combination of individual and contextual parameters. Individual measures of deprivation are rarely available at population level and area-based measures are then used as a surrogate of individual measures even though they reflect the contextual deprivation. Many studies showed the existence of socioeconomic inequalities in health outcomes including mortality, either using measures of deprivation at individual level [6,7,8] or using area-based measures [9,10,11]. The inequalities can result from individual factors such as different lifestyle behaviours, namely, smoking, alcohol, physical activity and dietary habits, different health literacy or access to health care, among other factors. There can also be ecological deprivation inequalities, defined by varying levels of availability of services or measures of the wellbeing of the population of a given area.
Previous studies have found association between health outcomes and deprivation in Portugal [12, 13]. However, no life tables have been constructed by area-based socioeconomic level in Portugal yet. In fact, deprivation-specific life tables are available for very few European countries. An ecological measure of socioeconomic deprivation, the European Deprivation Index (EDI) [14, 15] has recently become available in Portugal. The first aim of this study was thus to build deprivation-specific life tables using the Portuguese version of the EDI. The second aim was to evaluate mortality ratios between deprivation groups and trends in inequalities between 2000 and 2002 and 2010–2012 in Portugal.
The Portuguese version of the European Deprivation Index was used as deprivation indicator. This index was built using a methodology first proposed by Pornet and colleagues in 2012  and then applied to several European countries including Portugal [14, 15]. The index is based on census variables available for each country that are most associated with variables identified from the European Union Statistics on Income and Living Conditions (EU-SILC) survey . The index for Portugal based on 2001 census includes percentage of: non-owned households, households without indoor flushing, residents with low education level (≤6th grade), household with 5 rooms or less, unemployed looking for a job, female residents aged 65 years or more, households without bath/shower and percentage of residents employed in manual occupations . A score was obtained for each parish based on the census responses of its inhabitants. Most research on inequalities uses deprivation categorised according to quintiles. To build deprivation-specific life tables the EDI score was then categorized in five quintiles from the least deprived (q1) to the most deprived (q5) such that each quintile corresponded to 20% of the Portuguese population. Each deceased was assigned with the deprivation quintile corresponding to his/her parish of residence at the time of death.
Mortality and population data
Mortality rates in life tables require counts of deaths (numerator of the rates) and population (denominator) stratified by demographic variables (age, sex, others). This information is usually made available by the national statistics offices. The number of deaths by sex, age group (0, 1–4, 5–9, …, 85+), year of death and area of residence (parish) was obtained by special request to the Statistics Portugal (Instituto Nacional de Estatística). As common practice, to increase estimates’ stability, three years of data were considered centred on each census year for which the life tables were estimated (2000–2002 and 2010–2012). Population data was retrieved from the Statistics Portugal website (www.ine.pt). Number of residents by sex, age group (0, 1–4, 5–9, …, 85+) and parish was only available for census years (2001, 2011) so that the population was considered constant over the three years of each studied period. There were 4241 parishes in Portugal in 2001, with a median population of 969 inhabitants (min-max: 39–81,845), while in 2011 this number increased to 4260 (median population: 892, min-max: 31–66,250).
In this study, both the numbers of deaths and people (residing in the parishes) were summed up across the parishes for each period by sex, age group and level of deprivation.
When at subnational level, the counts of deaths and population produced by the national statistics offices are often available only by age groups (e.g. abridged) rather than by single years of age (e.g. complete) . Several methods for building complete life tables from abridged data have been in use, namely, Elandt–Johnson, Kostaki, Brass logit, and Akima spline methods . More recently, Rachet and colleagues  suggested a modelling approach to estimate smoothed mortality rates using flexible Poisson multivariable models. Death counts are modelled in the generalised linear model framework, considering a Poisson error and using splines to capture the effect of age. This method can use complete or abridged raw data allowing the estimation of complete life tables. This type of models was considered recommendable because it derives robust and unbiased estimates without making strong assumptions about age-specific mortality profiles. Also, a simulation study has shown that this method had better goodness of fit performance than other implemented methods . The age-group specific death counts were here modelled within this generalised linear model framework, considering a Poisson error with a log link function. The offset was considered the person-years at risk. Male and female death counts were modelled separately. Covariates considered in the model were age (using the mid-age of each age group), quintile of deprivation (dep), period (2000–2002 vs 2010–2012) and interactions between deprivation and age, deprivation and period and period and age. The model can be written as:
where dage,i,j denotes the number of deaths and pyrsage,i,j the number of person-years at risk for each age, deprivation group i and period j. The functions f, g1 and g2 represent restricted cubic splines. The knots positions were fixed a priori at ages 0, 1, 2 and 88 (for men) or 89 (for women). Although Rachet and colleagues considered further five knot positions selected from a set of 100 randomly simulated locations, here, we opted to consider knots at ages 10 to 50 at 10 years intervals since the other approach produced unrealistic predicted values. From these predefined knots position, the final number of knots was selected based on the Akaike Information Criterion (AIC).
Mortality rates were predicted from the fitted models by individual year of age (0–99), for each sex, period and quintile of deprivation. Life expectancies at birth were calculated from the fitted life tables. Mortality rate ratios (MRR) in terms of age were calculated from the predicted mortalities. MRR by EDI were calculated using the least deprived group as reference and the MRR by period using the period 2000–2002 as reference. The 95% confidence intervals (CI) for MRR were built assuming a normal distribution of log MRR and using the delta method. The derivations of the expressions for the CIs are presented as Additional file 1: Supplementary Material (S1).
All calculations were performed using STATA v13.1 and R v3.4.0. The STATA command mvrs was used for fitting the flexible Poisson model .
In the period 2000–2002 a total of 316,714 deaths were observed from which 219 (0.07%) were excluded due to unknown parish of residence at the time of death. In the period 2010–2012, the total number of deaths was 316,410 and only 75 deaths (0.02%) were excluded due to unknown parish or unknown age.
Both deprivation and period were found statistically significantly associated with mortality and all final fitted models included interactions between age and deprivation, period and age and period and deprivation, which were also found to be statistically significant. Higher mortality rates increasing with higher deprivation levels and decreasing with time periods were observed. The predicted mortality rates by age, sex, period and deprivation quintiles are presented in Additional file 2: Tables S1, Additional file 3: Table S2, Additional file 4: Table S3, Additional file 5: Table S4 and in Figs. 1 and 2. For all combinations period-sex-EDI analysed, the mortality rate first decreased with age reaching a minimum around 8–10 years-old and then steadily increased with age. For all deprivation quintiles q2 to q5, the mortality rates were in general significantly higher than the mortality rates of the least deprived group (q1) (Fig. 3). In men, at birth, the MRR between the most and the least deprived group was 1.62 (95%CI: 1.54–1.70) and 1.74 (95%CI: 1.65–1.83) in periods 2000–2002 and 2010–2012, respectively. In the first time period, the MRR’s decreased with age and from 81 onwards the ratio was no longer significantly different than 1. For the most recent period, this same pattern occurred from age 93 onwards (Fig. 3, top). In women, the MRRs were lower than the ones observed for men: at birth, the MRR between the two extreme deprivation groups was 1.26 (1.18–1.35) and 1.29 (1.20–1.38) in 2000–2002 and 2010–2012, respectively. The MRRs decreased with age but remained always significantly above one (Fig. 3, bottom). A reduction in the mortality rates was observed from the period 2000–2002 to the most recent period, 2010–2012. This reduction was significant for all ages (Fig. 4) and for both men and women. However, it was observed that the relative decrease in mortality was between 40 and 60% in less than 40 year olds but only 20% in older ages (over 60 year old) in men. This age gradient in the relative decrease between the two periods is less marked in women. For men, the reduction of mortality over time was less favourable for the most deprived group.
Life expectancy at birth increased from 74.0 years in 2000–2002 to 77.6 years in 2010–2012 in men and from 80.9 to 83.8 years in women. The gap in life expectancy at birth between the least and the most deprived group for men in the first period was 1.8 years. This gap slightly increased, over the ten-year period, to 2.1 years. For women, a smaller gap in life expectancy compared to men was observed. In 2000–2002 it was 1.0 year and it remained almost nearly constant over time (0.9 years in 2010–2012). The gap in life expectancy at 65 years was lower than at birth in both sexes. For men: 0.3 and 0.7 years for 2000–2002 and 2010–2012, respectively. In women, it was 0.5 years for both periods.
Persistent differences in mortality and life expectancy were observed according to ecological socioeconomic deprivation. These differences were larger among men and decreased with age for both sexes. Although mortality decreased significantly between the two periods, the deprivation gaps in mortality/life expectancy remained nearly constant from the period 2000–2002 to 2010–2012.
The smaller socioeconomic inequalities in mortality found in women have also been observed in other countries [21,22,23]. Several factors can contribute for this different pattern, including health behaviours and occupation. According to the last national health surveys made in Portugal, the prevalence of smoking was higher in men with low socioeconomic status while in women the prevalence was higher among individuals of high socioeconomic status . Ribeiro and colleagues  found no influence of deprivation on longevity after 75 for men and a weak association for women in Portugal. We observed here a decrease in the mortality rate ratios between deprivation groups with age. Similarly to Ribeiro and colleagues findings, the difference between mortality rates ceased to be significant after ages around 80 years in men, while in women a slight but still significant difference remained at all ages.
Richardson and colleagues analysed the evolution in regional gap in life expectancy at birth from 1991 to 2008 within European Union countries . No reduction in life expectancy gaps over the two decades analysed was observed, similarly to what has been observed in this study.
EUROSTAT (Statistical Office of the European Union) publishes estimates of life expectancy by age, sex and educational attainment level for several countries of the European Union including Portugal . The difference in the estimates of life expectancy at birth between the two extreme education groups, “Less than primary”, “primary and lower secondary education” and “Tertiary education”, presents a high variability between countries. While this difference, for males in 2011, was 19.3 years in Czech Republic, it was as low as 3.6 years in Turkey. In Portugal, this gap was 4.5 years, the third lowest within the 17 EU countries with published information. For women, the gaps were generally lower, ranging from 1.7 (Italy and Malta) to 8.7 years (Bulgaria). In Portugal it was 2.0 years. Although the EUROSTAT results are based on a simple individual level deprivation measure while our study is based on a composite area-based indicator, it is interesting to note that relatively low gaps were found for Portugal. Similarly, in our study the gap in life expectancy observed was low (1–2 years). Much lower than, for instance, the observed for England where deprivation-specific life tables built using area-based indices are available . Also, the smaller gap for women is in accordance with our findings.
The 2001 EDI version was used for both periods. By the time this study was developed the 2011 version was not available yet. However, preliminary results of the more recent version show that the EDI distribution by parish does not suffer major changes from 2001 to 2011. The impact of this limitation should therefore be small. For each individual, the EDI status was attributed based on the residence at the time of death. Although this can be a limitation, it is common practice in these types of studies and there was no information available that allowed matching mortality data for each individual to their residence history along their lifetime.
In this study we used a modelling approach to predict the mortality rate profiles by age. The flexible Poisson models have already been shown useful and valid to build life tables. They were used to build region-specific life tables within the CONCORD study (Cancer survival in five continents: a worldwide population-based study)  and by region, deprivation and ethnicity in England [4, 30]. The multivariable model allowed incorporating the deprivation and period effects, as well as interactions between them and with age, in a single model. Also it allowed obtaining complete life tables from abridged raw data.
Inference based on model predictions must have in consideration the correlation between the estimated parameters of the model. Ignoring this dependence and using the classical variance formulas as if the predicted values were observed ones would result in a sub estimation of the confidence intervals range. We thus derived and presented the variance estimators for the model-based mortality rate ratios taking into account this dependence.
This study is very relevant for the surveillance and monitoring of health inequalities, but it is important to highlight that these specific life tables are crucial tools to obtain reliable estimates of cancer survival within the relative survival data setting. In this setting, information on the cause of death is not available or not reliable. The disease-related survival (net survival) is then obtained indirectly by comparing the all-cause mortality of the cohort of patients with the mortality that would be experienced by individuals with the same demographic characteristics but free of the disease . The information on this expected (also called background) mortality is obtained from population life tables, assuming that the mortality due to the disease in question is negligible relatively to the overall mortality . To obtain valid net survival estimates, the population mortality should correctly reflect the expected mortality for each patient. The use of general life tables in the estimation of net survival for subgroups of the population with different overall mortality can lead to biased estimates of net survival. Estimation of net survival by deprivation is a situation where the use of general life tables can lead to overestimation of net survival in affluent groups, if these groups have a lower overall mortality than the general population, and the underestimation of net survival in the deprived groups, if these groups have a higher mortality than the general population [33, 34]. These questions arise also when stratifying by other factors that can influence overall mortality such as ethnicity .
In conclusion, this study has shown the existence of persistent socioeconomic inequalities in overall mortality in Portugal. Deprivation-specific life tables were built for Portugal. These life tables can therefore be used for monitoring inequalities and in future studies that require background mortality information in the estimation of deprivation-specific net survival from any specific disease.
Akaike Information Criterion
European Deprivation Index
Statistical Office of the European Union
European Union Statistics on Income and Living Conditions
Mortality Rate Ratio
Woods LM, Rachet B, Riga M, Stone N, Shah A, Coleman MP. Geographical variation in life expectancy at birth in England and Wales is largely explained by deprivation. J Epidemiol Community Health. 2005;59(2):115–20.
Auger N, Alix C, Zang G, Daniel M. Sex, age, deprivation and patterns in life expectancy in Quebec, Canada: a population-based study. BMC Public Health. 2010;10(1):161.
Tobias MI, Cheung J. Monitoring health inequalities: life expectancy and small area deprivation in New Zealand. Popul Health Metrics. 2003;1(1):2.
Morris M, Woods LM, Rachet B. A novel ecological methodology for constructing ethnic-majority life tables in the absence of individual ethnicity information. J Epidemiol Community Health. 2015;69(4):361–7.
Mackenbach JP, Bakker MJ. Tackling socioeconomic inequalities in health: analysis of European experiences. Lancet. 2003;362(9393):1409–14.
Pappas G, Queen S, Hadden W, Fisher G. The increasing disparity in mortality between socioeconomic groups in the United States, 1960 and 1986. N Engl J Med. 1993;329(2):103–9.
Mackenbach JP, Kunst AE. Measuring the magnitude of socio-economic inequalities in health: an overview of available measures illustrated with two examples from Europe. Soc Sci Med. 1997;44(6):757–71.
Kohler IV, Martikainen P, Smith KP, Elo IT. Educational differences in all-cause mortality by marital status - evidence from Bulgaria, Finland and the United States. Demogr Res. 2008;19(10):2011–42.
van Lenthe FJ, Borrell LN, Costa G, Diez Roux AV, Kauppinen TM, Marinacci C, et al. Neighbourhood unemployment and all cause mortality: a comparison of six countries. J Epidemiol Community Health. 2005;59(3):231–7.
Toprani A, Li W, Hadler JL. Trends in mortality disparities by area-based poverty in new York City, 1990–2010. J Urban Health. 2016;93(3):538–50.
Pickett KE, Pearl M. Multilevel analyses of neighbourhood socioeconomic context and health outcomes: a critical review. J Epidemiol Community Health. 2001;55(2):111–22.
Ribeiro AI, Krainski ET, Carvalho MS, De Fátima de Pina M. The influence of socioeconomic deprivation, access to healthcare and physical environment on old-age survival in Portugal. Geospat Health. 2017;12(2):581.
Apolinário D, Ribeiro AI, Krainski E, Sousa P, Abranches M, Duarte R. Tuberculosis inequalities and socio-economic deprivation in Portugal. Int J Tuberc Lung Dis. 2017;21(7):784–9.
Guillaume E, Pornet C, Dejardin O, Launay L, Lillini R, Vercelli M, et al. Development of a cross-cultural deprivation index in five European countries. J Epidemiol Community Health. 2015;2015:205729.
Ribeiro AI, Mayer A, Miranda A, Pina MF. The Portuguese version of the European deprivation index: an instrument to study health inequalities. Acta Medica Port. 2017;30(1):17.
Pornet C, Delpierre C, Dejardin O, Grosclaude P, Launay L, Guittet L, et al. Construction of an adaptable European transnational ecological deprivation index: the French version. J Epidemiol Community Health. 2012;66(11):982–9.
Eurostat. Access to Microdata. European Union Statistics on Income and Living Conditions (EU-SILC) [Internet]. 2015 [cited 2018 May 19]. Available from: http://ec.europa.eu/eurostat/web/microdata/european-union-statistics-on-income-and-living-conditions
Rachet B, Maringe C, Woods LM, Ellis L, Spika D, Allemani C. Multivariable flexible modelling for estimating complete, smoothed life tables for sub-national populations. BMC Public Health. 2015;15(1):1240.
Baili P, Micheli A, Montanari A, Capocaccia R. Comparison of four methods for estimating complete life tables from abridged life tables using mortality data supplied to Eurocare-3. Math Popul Stud. 2005;12(4):183–98.
Royston P, Sauerbrei W. Multivariable modeling with cubic regression splines: a principled approach. Stata J. 2007;7(1):45–70.
Stronks K, van de Mheen H, van den Bos J, Mackenbach JP. Smaller socioeconomic inequalities in health among women: the role of employment status. Int J Epidemiol. 1995;24(3):559–68.
Mackenbach JP, Kunst AE, Groenhof F, Borgan JK, Costa G, Faggiano F, et al. Socioeconomic inequalities in mortality among women and among men: an international study. Am J Public Health. 1999;89(12):1800–6.
Mustard CA, Etches J. Gender differences in socioeconomic inequality in mortality. J Epidemiol Community Health. 2003;57(12):974–80.
Alves J, Kunst AE, Perelman J. Evolution of socioeconomic inequalities in smoking: results from the Portuguese national health interview surveys. BMC Public Health. 2015;15(1):311.
Ribeiro AI, Krainski ET, Carvalho MS, Launoy G, Pornet C, de Pina MF. Does community deprivation determine longevity after the age of 75? A cross-national analysis. Int J Public Health. 2018;63(4):469–79.
Richardson EA, Pearce J, Mitchell R, Shortt NK, Tunstall H. Have regional inequalities in life expectancy widened within the European Union between 1991 and 2008? Eur J Pub Health. 2014;24(3):357–63.
EUROSTAT. Life expectancy by age, sex and educational attainment level [Internet]. 2017 [cited 2018 May 12]. Available from: https://data.europa.eu/euodp/data/dataset/94EpMjNjb6IwaNHIfxxd8w.
Public Health Outcomes Framework: Health equity report focus on ethnicity public health outcomes Framework: health equity report-focus on ethnicity. 2017.
Spika D, Bannon F, Bonaventure A, Woods LM, Harewood R, Carreira H, et al. Life tables for global surveillance of cancer survival (the CONCORD programme): data sources and methods. BMC Cancer. 2017;17(1):159.
Cancer Research UK Cancer Survival Group. UK life tables [Internet]. [cited 2018 May 19]. Available from: http://csg.lshtm.ac.uk/tools-analysis/uk-life-tables/
Perme MP, Stare J, Estève J. On estimation in relative survival. Biometrics. 2012;68(1):113–20.
Esteve J, Benhamou E, Raymond L. Techniques for survival analysis survival analysis in descriptive epidemiology. In: Statistical methods in Cancer research volume IV - descriptive epidemiology. Lyon: IARC; 1994.
Auvinen A, Karjalainen S. Possible explanations for social class differences in cancer patient survival. IARC Sci Publ. 1997;138:377–97.
Dickman PW, Auvinen A, Voutilainen ET, Hakulinen T. Measuring social class differences in cancer patient survival: is it necessary to control for social class differences in general population mortality? A Finnish population-based study. J Epidemiol Community Health. 1998;52(11):727–34.
Blakely T, Soeberg M, Carter K, Costilla R, Atkinson J, Sarfati D. Bias in relative survival methods when using incorrect life-tables: lung and bladder cancer by smoking status and ethnicity in New Zealand. Int J Cancer. 2012;131(6):974–82.
Ana Isabel Ribeiro was supported during this study by the European Commission (Horizon 2020 grant number 633666) and by FEDER through the Operational Programme Competitiveness and Internationalization and national funding from the Foundation for Science and Technology – FCT (Portuguese Ministry of Science, Technology and Higher Education) under the Unidade de Investigação em Epidemiologia - Instituto de Saúde Pública da Universidade do Porto (EPIUnit) (POCI-01-0145-FEDER-006862; Ref. UID/DTP/04750/2013). Camille Maringe was supported by Cancer Research UK grant number C7923/A18525. The funding bodies had no influence in the design of the study and collection, analysis, and interpretation of data.
Availability of data and materials
The datasets used and/or analysed during the current study are available from the corresponding author on reasonable request.
Ethics approval and consent to participate
Only aggregated data publicly available from the Statistics Portugal office was used exempting special permissions for its use. No ethics approval was needed for this study.
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary S1. Estimation of confidence intervals for model-based predicted mortality rate ratios. (PDF 111 kb)
Table S1. Life tables by deprivation quintile for men in the period 2000–2002. (PDF 411 kb)
Table S2. Life tables by deprivation quintile for women in the period 2000–2002. (PDF 411 kb)
Table S3. Life tables by deprivation quintile for men in the period 2010–2012. (PDF 411 kb)
Table S4. Life tables by deprivation quintile for women in the period 2010–2012. (PDF 411 kb)