PMCCPMCCPMCC

Search tips
Search criteria 

Advanced

 
Logo of scirepAboutEditorial BoardFor AuthorsScientific Reports
 
Sci Rep. 2017; 7: 11465.
Published online 2017 September 13. doi:  10.1038/s41598-017-11785-x
PMCID: PMC5597615

Changing smoking-mortality association over time and across social groups: National census-mortality cohort studies from 1981 to 2011

Abstract

The difference in mortality between current and never-smokers varies over time, affecting future projections of health gains from tobacco control. We examine this heterogeneity by sex, ethnicity and cause of death on absolute and relative scales using New Zealand census data. These data included smoking status, and were linked to subsequent mortality records in 1981–84, 1996–99 and 2006–11 for 25–74 year olds (16.1 million person-years of follow-up). Age-standardised mortality rates and rate differences (SRDs) were calculated comparing current to never-smokers, and Poisson regression was used to adjust for multiple socioeconomic factors and household smoking. We found that mortality declined over time in never-smokers; however, mortality trends in current-smokers varied by sex, ethnicity and cause of death. SRDs were stable over time in European/Other men, moderately widened in European/Other women and markedly increased in Māori men and women (Indigenous population). Poisson smoking-mortality rate ratios (RRs) increased from 1981–84 to 1996–99 with a moderate increase from 1996–99 to 2006–11 (RRs 1.48, 1.77, 1.79 in men and 1.51, 1.80, 1.90 in women). Socioeconomic confounding increased over time. In summary, this marked heterogeneity in smoking-mortality RRs over time has implications for estimating the future health and inequality impacts of tobacco control interventions.

Introduction

Tobacco is a causal risk factor for a substantial burden of premature mortality and it is the greatest contributor to ethnic and socioeconomic inequalities in mortality in many high-income countries1. However, the mortality burden and social patterning of smoking varies over time26. It takes decades for the peak harms of smoking to manifest due to time-lags from smoking initiation to disease-specific mortality. Once this time-lag has elapsed, the harms of smoking compared to never-smoking may reduce over time in absolute terms, due to changes in smoking behaviours (eg, number of cigarettes per day) or product toxicity. Conversely, this tobacco-harm often occurs at the same time as long-term downward trends in never-smoker mortality due to other reasons (eg, falling coronary heart disease rates). This dynamic situation means that at least one of the absolute (rate difference, RD) and relative mortality comparisons (rate ratio, RR) between current and never-smokers, will change over time as the tobacco epidemic unfolds. Figure 1 is one possible realisation of this dynamic (others are shown in Supplementary File Appendix A). The relative differences in smoking by mortality (RRs) are likely to go up both early in the tobacco epidemic (due to time-lags from smoking initiation to mortality impact) and also may increase later in the tobacco epidemic (if the fall in never smoker mortality is faster than any fall in the absolute harms of tobacco).

Figure 1
Hypothesised progression of current and never smoker mortality rates (top), and the corresponding rate ratios and rate differences (bottom), over 100 years of the tobacco epidemic.

Does this variation in the effect size from smoking, be it either or both on the absolute and relative scales, matter for tobacco control or eradication policy? Almost certainly yes. If one is making a case for tobacco reduction or eradication based on likely future health gains, then having some idea – even quantitative projections – of the future RD and RRs comparing current (and ex) to never smokers is necessary for accurate prediction of mortality reductions and (healthy) life expectancy extensions. Such information might allow policy-makers to more appropriately prioritise tobacco control interventions relative to competing ways to reduce health loss and health inequalities (eg, via obesity prevention interventions).

This theoretical expectation of at least one of the relative and absolute smoking-mortality associations changing over time is supported by recent longitudinal studies of the tobacco epidemic that report increasing relative rates of mortality in current-smokers compared to never-smokers over several decades710. For example, there were large increases in smoking RRs and RDs for all-cause mortality and lung cancer mortality in the United States (US) between the 1960s, 1980s and 2000s8. An artefactual reason for why the smoking-mortality association may vary is increased confounding over time as smoking becomes more strongly patterned by socioeconomic position (SEP), passive smoking11 and other behaviours (eg, diet) that are independently associated with mortality. Previous studies are limited in their adjustment for confounding by SEP8, 10, 12 (eg, stratifying on education only) and may be susceptible to increased residual confounding over time13. This study thoroughly adjusts for confounding by SEP and household smoking (HHS) and the large sample size enables us to examine fine-grained heterogeneity in the data by sex, ethnicity and cause of death over time. Most existing studies are limited to one cohort (with the exception of Thun et al.)8 and a subset of the population9, 10, 12.

The change over time in excess smoking mortality is expected to be differentially phased over time by social groups such as sex, SEP and ethnicity/race. For example, tobacco is classically taken up in a population by men first, then by women14. Historically the relative mortality rate in current compared to never-smokers was greater in men than in women but several studies show a convergence in recent cohorts8, 12 consistent with epidemic phasing. There are also examples where smoking-mortality rates, RDs and RRs for lung cancer are increasing among women but appear to be reaching a plateau in men8. Studies in the United States and New Zealand have also shown notable differences in the smoking-mortality RRs by ethnicity7, 15, but little difference on the absolute scale (RD) by ethnicity7. Thus while tobacco control appears to be one of the most efficient ways to reduce health inequalities16, this too will change over time with the phasing of the tobacco epidemic.

Given this background and rationale, the objectives of this study were to: 1) describe the variation over three decades in the smoking-mortality association on both the absolute and relative scale, by sex and ethnic group; 2) determine the changing contributions of specific causes of death to excess smoking mortality over time; and 3) quantify the increasing confounding by SEP and HHS (the latter as a proxy for passive smoking exposure)11, 17 over time. European/Other and Māori populations in New Zealand were selected as a case study, extending a previous study8 to examine mortality records from three national census-linked cohorts, with 16.1 million person-years of follow-up, high quality data on smoking, ethnicity and SEP and relatively thorough control of confounding including from household smoking. Māori are the Indigenous peoples of New Zealand and have higher rates of non-communicable disease for multiple reasons – not just tobacco.

Results

Variation in SRDs and SRRs over time by population group

Figure 2 shows the age-standardised all-cause mortality rates in never- and current-smokers over time. Among European/Other men, current and never-smoker mortality rates have fallen roughly in parallel with little change in SRDs over time (435, 499 and 439 per 100,000), with a corresponding increase in SRRs over time from 1.71 to 2.48 (108% increase in excess SRR (ie, SRR minus 1); Fig. 3). Among European/Other women, the SRDs were less than in men, but widened over time from 238 to 319 per 100,000 (34% increase) with more pronounced widening in SRRs from 1.65 to 2.64 (152% increase in excess RR).

Figure 2
Age-standardised mortality rates by smoking status and ethnicity for men and women aged 25–74 years in the New Zealand Census Mortality Study in 1981–84, 1996–99 and 2006–11 (95% confidence intervals are indicated ...
Figure 3
Age-standardised all-cause mortality rates per 100,000 in current-smokers and never-smokers showing the standardised rate differences (RD) and the standardised rate ratios (RR) by cohort and sex in the New Zealand Census Mortality Study (for those ...

For Māori, there was only a 13% (men) and 6% (women) higher mortality among current versus never-smokers in the early 1980s, but there was an increase in mortality in smokers from the 1980s to 1990s (in comparison to a decrease for never-smokers). This corresponded to a marked increase in both the SRDs and SRRs among Māori men and women in this time period. From 1996–99 to 2006–11, the SRDs were stable (598 and 567 per 100,000 for Māori men, and 323 and 373 for Māori women) and somewhat higher than among European/Other within each sex by time pair. By 2006–11 the SRR was 1.71 among both Māori men and women which was considerably less than among European/Other (2.48 and 2.64 respectively).

Contributions of specific tobacco-related causes of death

Figure 4 shows the age-standardised mortality rate differences disaggregated by the major tobacco-related causes of mortality (see Supplementary Table S1 for SRD and SRR figures). The IHD smoking SRDs (the height of the IHD labelled part of the total bar) narrowed over time among European/Other men and women. Conversely, the IHD SRDs increased over time among Māori. The stroke SRDs were similar to IHD for European/Other, but Māori stroke SRDs were measured with considerable imprecision (Supplementary Table S1). SRDs for lung cancer and COPD were consistently high over time in Māori and European/Other men compared to women, with the highest SRDs in Māori men. In Māori and European/Other women, SRDs doubled for lung cancer (from 44 to 90 per 100,000 in European/Other and 79 to 164 in Māori) and for COPD. The difference between men and women reflects the generally increasing lung cancer and COPD mortality in female smokers and the generally stable lung cancer and COPD mortality rates in male smokers (Fig. 2).

Figure 4
Decomposition of the age-standardised rate differences in mortality between current smokers and never smokers by mortality type in 25–74 year olds by sex and ethnicity, New Zealand Census Mortality Study.

Confounding over time

Table 1 presents Poisson regression RRs comparing mortality in current- and never-smokers by sex for each of the three time periods, with sequential adjustment for confounders. There is a moderate reduction in the RR after adjustment for SEP and HHS at any point in time, but this confounding increases over time such that by 2006–11 among men the (excess) RRs reduce by 36% from 2.24 to 1.79 when additionally adjusted for SEP and HHS; the reduction was 24% for men in 1981-84 (RRs 1.63 adjusted for age and ethnicity, and 1.48 fully-adjusted). For women, the reduction due to confounding was 31% in 2006–11 compared to 15% in 1981–84.

Table 1
Poisson regression rate ratios for the association between smoking and mortality over time in three cohorts 25–74 years old in the New Zealand Census Mortality Study (95% confidence intervals in brackets).

Changing smoking rate ratios over time by sex, cause of death, and ethnicity

Poisson smoking-mortality RR estimates fully-adjusted for SEP and HHS increased over time, mostly between 1981–84 and 1996–99 in men and women, with only modest increases from 1996–99 to 2006–11. A similar pattern was evident for lung cancer and COPD RRs (Table 1, Supplementary Figure S7) with relatively steeper increases in the RRs for IHD and stroke in women over time.

In the fully-adjusted Poisson model, there was an interaction between smoking and ethnicity with greater smoking-mortality RRs (p-value < 0.001) in European/Other than Māori at all time points for both sexes (Table 1). RRs increased in all sex-ethnicity groups.

Discussion

This study indicates an increase in smoking-mortality RRs over time with differential phasing by social group (sex and ethnicity) and cause of death. Mortality declined over time in all groups of never-smokers; however, mortality trends in current-smokers varied by sex, ethnicity and cause of death (Fig. 2). Among European/Other, all-cause mortality rates declined in male smokers but were more stable in women smokers, consistent with findings in a US study8.

Patterns of results can be aligned and contrasted with the framework presented in Fig. 1 as an example of the phases of the tobacco epidemic. All-cause mortality SRDs for Māori and European/Other men increased from the 1980s to 1990s and then decreased in the 2000s; whereas for European/Other women the absolute mortality gap between current and never smokers appeared to plateau in the 1990s and 2000s, and among Māori women they increased over all decades. These results are consistent with the delayed peak in the tobacco epidemic intensity for women (Fig. 1). There is more variability again by cause of death. COPD/lung cancer SRDs increased among women but not men. Conversely, the SRD for COPD increased in 55 + year old men in a US study8 perhaps consistent with New Zealand up to 2011 being a step further along the framework timeline. IHD and stroke SRDs declined among European/Other in New Zealand (Fig. 4), but among Māori the IHD SRD increased over time. This likely indicates that Māori are both earlier in the tobacco epidemic and have experienced slower declines in never-smoker IHD mortality rates (Fig. 2).

The relative mortality gap between current and never (RRs) increased for men and women from 1981–84 to 1996–99 – with modest increases from 1996–99 to 2006–11 (except for Māori women where the increase was greater). This increasing relative mortality gap between current and never smokers occurred even after adjusting for an increase in the amount of confounding by SEP and HHS over time. Many long-term studies report increasing RRs over time710, however, the possible plateauing recently in our New Zealand study appears to be novel.

There are three likely contributing reasons for changing RRs over time: 1) the changing intensity of smoking over time, 2) decreasing mortality rates in never-smokers and 3) the time-lag from smoking to peak mortality harm. We expand these below in general terms that should apply to most countries, and point to examples from New Zealand as a case study.

  1. Changes in smoking intensity over time
    Increased mortality in smokers is linked to smoking intensity or pack-year history8, 9, 18, including greater duration of smoking, age at which smoking was initiated, more cigarettes per day, increased inhalation per cigarette (eg, smoking down to the butt) and/or changes to the harmfulness of the product over time such as tar/carcinogen content or roll-your-own smoking. Using New Zealand as a case example, youth smoking prevalence in New Zealand was high in the 1960s and 1970s (and possibly earlier) for men and peaked in the mid-1970s for women19, consistent with a delayed phasing of the tobacco epidemic in women. Smoking prevalence in young Māori women remained at higher rates than in European/Other and decreased more slowly from the 1970s-1990s19.
  2. Decreasing mortality rates in never-smokers
    All-cause mortality has declined substantially in never-smokers. Similar absolute mortality rate declines in current and never-smokers mathematically will – necessarily – result in an increase in the smoking-mortality RR. Furthermore, there is likely to be heterogeneity between social groups with varying background (and among never smokers) mortality rates due to differential timing of epidemiological transitions. Using New Zealand as a case example, Māori never-smokers are a population with high background mortality rates, due to other risk factors for mortality such as obesity, infectious diseases, and lower access to some health services20. If absolute excess mortality from smoking is fixed then higher never-smoker mortality will result in lower smoking-mortality RRs, as seen for Māori in this study. Our results are consistent with several studies internationally that report different smoking-mortality RR associations by social group including by ethnicity15, sex8 and country21.
  3. Time-lag from smoking to peak mortality harm

Causes of death changed in their contribution to the excess smoking-mortality over time, consistent with heterogeneity by cause of death8, 12, 22. The different patterns between IHD/stroke and lung cancer/COPD are affected by the different time-lag between smoking and peak smoking-related mortality14, 23. Using New Zealand as a case example, in the early 1980s the excess lung cancer and COPD mortality in smokers (SRDs) was notably greater in men than women, but by the late 2000s SRDs in women had increased to be similar to men. This is consistent with historical differences in the phasing of the tobacco epidemic by sex (approximately 30 years is the time-lag from smoking to peak lung cancer mortality). However, throughout the study the IHD/stroke SRD trends by smoking status did not differ markedly by sex, consistent with similar effective doses of tobacco by sex in recent times and shorter time-lags to cardiovascular disease.

Confounding

Smoking is increasingly concentrated in socially disadvantaged populations in many settings (eg, as per this European study24), and these populations usually have higher rates of other risk factors for non-communicable diseases. Thus, theoretically we expect confounding of the tobacco-mortality association to increase over time (eg, by socioeconomic factors, and more proximally diet, physical activity and alcohol consumption). This New Zealand case-study is consistent with that expectation, with increasing confounding of the smoking-mortality association by SEP over time. We were ‘only’ able to adjust for six socioeconomic factors, thus meaning residual confounding by risk factors (eg, diet, obesity) is possible. However, these risk factors are largely (not necessarily fully) explained by socioeconomic factors . In a previous paper we applied quantitative bias analysis methods to evaluate the marginal effect of confounding by alcohol and obesity (over and above SEP) on smoking-cancer associations, and showed that any residual confounding over and above SEP was likely to be small25. In this study we also adjusted for HHS, but compared to other SEP variables it was found to have little additional impact as a confounder of the active-smoking-mortality relationship.

Study strengths and limitations

This New Zealand case-study included 16.1 million person-years of follow-up, three national cohort studies, high quality census data, a study timeframe representing 30 years and thorough adjustment for confounding by SEP. Sensitivity analyses suggest that selection bias (restriction to individuals with complete data in Poisson analysis) may have contributed a small degree to higher smoking RRs (Supplementary File). Mismeasurement bias and reverse causality from people quitting smoking due to illness may both have biased our results towards the null. In each cohort we had only one cross-sectional measure of an individual’s smoking status which does not allow for duration of smoking, cigarettes per day or type of tobacco consumed. Longitudinal examination of detailed smoking exposures would improve our understanding of the impact of dose on mortality. Nevertheless, trends in such biases are unlikely to change markedly over time (see Supplementary File) and thus they are unlikely to substantially impact on observed trends in RRs and RDs over time.

Possible implications

These results signal the policy and research importance of understanding the phases of the tobacco epidemic, and the variation in mortality consequences from various trajectories of smoking prevalence. The heterogeneity in the smoking-mortality association over time has implications for how researchers project and quantify future health and inequality impacts of tobacco control interventions. Such improved information might allow policy-makers to more appropriately prioritise tobacco control interventions relative to competing ways to reduce health loss and health inequalities (eg, via obesity prevention interventions or enhanced alcohol control interventions). But of course such prioritisation also requires similar quality information on the other health problems (eg, the obesity epidemic) and also the health economic aspects of all the specific interventions (ie, the relative cost-effectiveness of specific interventions such as a tobacco tax vs a sugary drinks tax). If tobacco control interventions are prioritised, there are specific interventions that appear able to reduce inequalities in smoking (including tobacco tax increases26, some smoking cessation services27 and various “tobacco endgame strategies”28).

Conclusions

This study reports an increase in the relative mortality gap between current and never smokers (RRs) over time and that this change persisted after adjusting for confounding by socioeconomic position. There was differential phasing by sex, ethnicity and cause of death. Variation in smoking-mortality RRs is likely to be related to changes in the intensity of smoking over time, underlying trends in mortality rates among never-smokers and time-lags between smoking and peak-mortality for different causes of death. Application of a framework helped provide understanding about the phasing of the tobacco epidemic. These findings have implications for estimating the future health and inequality impacts of tobacco control interventions and how these might be prioritised relative to other interventions that reduce health loss and reduce health inequalities.

Methods

Study population

Three closed cohorts of the New Zealand resident population living in a private dwelling, on census night in 1981, 1996 and 2006 were created by linking census and mortality records for three subsequent years following the 1981 and 1996 censuses and five years subsequent to the 2006 census. Probabilistic linkage was done with QualityStage TM software using an individual’s address (meshblock or census area unit), sex, date of birth, ethnicity and country of birth as matching variables. Approximately 98% of links were estimated to be true links where the mortality record was correctly linked to a previous census record29. This provided 111,000 deaths from 16.1 million person-years of follow-up (Table 2). The percentage of deaths linked to a census record ranged from 71% in 1981 to 83% in 2006. Linkage weights were used to adjust for incomplete linkage of mortality records to make the data representative of all deaths. Namely, the inverse probabilities of a mortality record being linked to the census in each age-ethnicity-deprivation-region strata were used as weights in the analysis. We present results for individuals aged 25–74 years old in the follow-up period, given more optimal measures of SEP in this group.

Table 2
Baseline characteristics of the study population by sex and years of follow-up for New Zealand 25–74 year olds in 1981, 1996, and 2006.

Smoking information from the census night was categorised as current-smoker, ex-smoker, or never-smoker with ex-smokers put aside for the majority of analyses in this paper (see census questions in Supplementary Material Appendix B). HHS was defined as living in a household with at least one other person who currently smoked, irrespective of whether the index individual was a smoker or not. The outcomes were all-cause mortality and the four largest contributors to smoking-related mortality: lung cancer, ischaemic heart disease (IHD), stroke and chronic obstructive pulmonary disease (COPD). Ethnicity was self-identified and prioritised as Māori, Pacific, Asian or European/Other. The latter group comprised individuals who did not identify as Māori, Pacific or Asian, the large majority of whom were New Zealand European.

Analysis

Age-standardised rates, standardised rate differences (SRDs) and standardised rate ratios (SRRs) were calculated to compare mortality in current-smokers with mortality in never-smokers. Direct standardisation was applied using the WHO World Standard Population to maximise international comparability.

Poisson regression was carried out in SAS30 to estimate the mortality rate ratio (RR) in current-smokers compared with never-smokers. The fully-adjusted model included age (5-year age-groups), ethnicity (prioritised by Māori, Pacific, Asian and European/Other), equivalised household income (using a New Zealand-specific index)31, neighbourhood deprivation quintile32, 33, education (none, school, post-school qualifications), labour force status (employed, job-searching, inactive), housing tenure (owner, renting, unknown), car access (0, 1 or 2+ motor vehicles) and HHS (none, lives with a current-smoker).

Sensitivity analyses were done to investigate the impact of selection bias (restricting to individuals with complete income/SEP information), reverse causation (by excluding first year of census-mortality cohort follow-up) and a negative control (smoking is expected to have a minimal association with unintentional injury). Results from sensitivity analyses did not raise major concerns about the main results presented in this paper (see Supplementary Material Appendix D).

Data availability

Supplemental information with additional methods, results and further explanation of the smoking mortality relationship and phasing of the tobacco epidemic over time is attached. Programming code and data sharing is available on request from the authors.

Ethical approval

Ethics approval was provided by the Central Regional Ethics Committee, reference number WGT/04/10/093.

Electronic supplementary material

Acknowledgements

The NZCMS is conducted in collaboration with Statistics New Zealand and within the confines of the Statistics Act 1975. Note: Access to the data used in this study was provided by Statistics New Zealand under conditions designed to give effect to the security and confidentiality provisions of the Statistics Act 1975. The results presented in this study are the work of the authors, not Statistics New Zealand. Furthermore, we wish to thank and acknowledge comments on the final draft by Professor Richard Edwards. This work was supported by the New Zealand Ministry of Health [425630/347384/00 University of Otago - New Zealand Census Mortality and Cancer Trends study 2006–2011].

Author Contributions

Author Contributions

T.B. was the principal investigator and designed the study, prepared a figure and interpreted the results. A.T. analysed the data, prepared the figures, interpreted the results and wrote the first draft of the manuscript. J.A. acquired and analysed the data. G.D. and N.W. contributed to study design and interpreted the results. All authors commented on manuscript drafts and approved the final version.

Notes

Competing Interests

The authors declare that they have no competing interests.

Footnotes

Electronic supplementary material

Supplementary information accompanies this paper at doi:10.1038/s41598-017-11785-x

Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

References

1. Gregoraci G, et al. Contribution of smoking to socioeconomic inequalities in mortality: a study of 14 European countries, 1990–2004. Tob Control. 2017;26:206–268. doi: 10.1136/tobaccocontrol-2015-052766. [PubMed] [Cross Ref]
2. Teng AM, et al. Ethnic inequalities in cancer incidence and mortality: census-linked cohort studies with 87 million years of person-time follow-up. BMC Cancer. 2016;16:755. doi: 10.1186/s12885-016-2781-4. [PMC free article] [PubMed] [Cross Ref]
3. Blakely T, Wilson N. The contribution of smoking to inequalities in mortality by education varies over time and by sex: two national cohort studies, 1981-84 and 1996-99. Int J Epidemiol. 2005;34:1054–1062. doi: 10.1093/ije/dyi172. [PubMed] [Cross Ref]
4. Mackenbach JP, et al. Inequalities in lung cancer mortality by the educational level in 10 European populations. Eur J Cancer. 2004;40:126–135. doi: 10.1016/j.ejca.2003.10.018. [PubMed] [Cross Ref]
5. Jha P, et al. Social inequalities in male mortality, and in male mortality from smoking: indirect estimation from national death rates in England and Wales, Poland, and North America. Lancet. 2006;368:367–370. doi: 10.1016/S0140-6736(06)68975-7. [PubMed] [Cross Ref]
6. Menvielle G, et al. The role of smoking and diet in explaining educational inequalities in lung cancer incidence. J Natl Cancer Inst. 2009;101:321–330. doi: 10.1093/jnci/djn513. [PMC free article] [PubMed] [Cross Ref]
7. Hunt D, Blakely T, Woodward A, Wilson N. The smoking-mortality association varies over time and by ethnicity in New Zealand. Int J Epidemiol. 2005;34:1020–1028. doi: 10.1093/ije/dyi139. [PubMed] [Cross Ref]
8. Thun MJ, et al. 50-year trends in smoking-related mortality in the United States. N Engl J Med. 2013;368:351–364. doi: 10.1056/NEJMsa1211127. [PMC free article] [PubMed] [Cross Ref]
9. Doll R, Peto R, Boreham J, Sutherland I. Mortality in relation to smoking: 50 years’ observations on male British doctors. BMJ. 2004;328:1519. doi: 10.1136/bmj.38142.554479.AE. [PMC free article] [PubMed] [Cross Ref]
10. Sakata R, et al. Impact of smoking on mortality and life expectancy in Japanese smokers: a prospective cohort study. BMJ. 2012;345:e7093. doi: 10.1136/bmj.e7093. [PubMed] [Cross Ref]
11. Hill S, Blakely T, Kawachi I, Woodward A. Mortality among “never smokers” living with smokers: two cohort studies, 1981-4 and 1996-9. BMJ. 2004;328:988–989. doi: 10.1136/bmj.38070.503009.EE. [PMC free article] [PubMed] [Cross Ref]
12. Freedman, N. D. et al. Impact of changing US cigarette smoking patterns on incident cancer: risks of 20 smoking-related cancers among the women and men of the NIH-AARP cohort. Int J Epidemiol, dyv175, doi:10.1093/ije/dyv175 (2015). [PMC free article] [PubMed]
13. Blakely T, Hunt D, Woodward A. Confounding by socioeconomic position remains after adjusting for neighbourhood deprivation: an example using smoking and mortality. J Epidemiol Commun H. 2004;58:1030–1031. doi: 10.1136/jech.2004.019737. [PMC free article] [PubMed] [Cross Ref]
14. Lopez AD, Collishaw NE, Piha T. A descriptive model of the cigarette epidemic in developed countries. Tob Control. 1994;3:242–247. doi: 10.1136/tc.3.3.242. [Cross Ref]
15. Haiman CA, et al. Ethnic and Racial Differences in the Smoking-Related Risk of Lung Cancer. N Engl J Med. 2006;354:333–342. doi: 10.1056/NEJMoa033250. [PubMed] [Cross Ref]
16. Blakely T, et al. Health, Health Inequality, and Cost Impacts of Annual Increases in Tobacco Tax: Multistate Life Table Modeling in New Zealand. PLoS Med. 2015;12:e1001856. doi: 10.1371/journal.pmed.1001856. [PMC free article] [PubMed] [Cross Ref]
17. Li, W. et al. Secondhand smoke enhances lung cancer risk in male smokers: an interaction. Nicotine Tob Res, ntw115 (2016). [PubMed]
18. Banks E, et al. Tobacco smoking and all-cause mortality in a large Australian cohort study: findings from a mature epidemic with current low smoking prevalence. BMC Med. 2015;13:1–10. doi: 10.1186/s12916-015-0281-z. [PMC free article] [PubMed] [Cross Ref]
19. McGee R, Williams S, Stanton W. Smoking among New Zealand adolescents 1960-93. NZ Med J. 1995;108:85–87. [PubMed]
20. Tobias M, Blakely T, Matheson D, Rasanathan K, Atkinson J. Changing trends in indigenous inequalities in mortality: lessons from New Zealand. Int J Epidemiol. 2009;38:1711–1722. doi: 10.1093/ije/dyp156. [PubMed] [Cross Ref]
21. Jacobs DR, et al. Cigarette Smoking and Mortality Risk: Twenty-five–Year Follow-up of the Seven Countries Study. Arch Intern Med. 1999;159:733–740. doi: 10.1001/archinte.159.7.733. [PubMed] [Cross Ref]
22. Pujades-Rodriguez M, et al. Heterogeneous associations between smoking and a wide range of initial presentations of cardiovascular disease in 1 937 360 people in England: lifetime risks and implications for risk prediction. Int J Epidemiol. 2015;44:129–141. doi: 10.1093/ije/dyu218. [PMC free article] [PubMed] [Cross Ref]
23. Lynch J, Smith GD. A life course approach to chronic disease epidemiology. Annu Rev Public Health. 2005;26:1–35. doi: 10.1146/annurev.publhealth.26.021304.144505. [PubMed] [Cross Ref]
24. Hu, Y. et al. The impact of tobacco control policies on smoking among socioeconomic groups in nine European countries, 1990-2007. Nicotine Tob Res (E-publication 16 August), doi:10.1093/ntr/ntw210 (2016). [PubMed]
25. Blakely T, et al. The association of active smoking with multiple cancers: national census-cancer registry cohorts with quantitative bias analysis. Cancer Causes Control. 2013;24:1243–1255. doi: 10.1007/s10552-013-0204-2. [PubMed] [Cross Ref]
26. Brown T, Platt S, Amos A. Equity impact of population-level interventions and policies to reduce smoking in adults: a systematic review. Drug Alcohol Depend. 2014;138:7–16. doi: 10.1016/j.drugalcdep.2014.03.001. [PubMed] [Cross Ref]
27. Nghiem N, et al. A national quitline service and its promotion in the mass media: modelling the health gain, health equity and cost-utility. Tob Control. 2017 [PubMed]
28. van der Deen, F. S. et al. Impact of five tobacco endgame strategies on future smoking prevalence, population health and health system costs: two modelling studies to inform the tobacco endgame. Tob Control (E-publication 24 June) (2017). [PubMed]
29. Blakely T, Salmond C. Probabilistic record linkage and a method to calculate the positive predictive value. Int J Epidemiol. 2002;31:1246–1252. doi: 10.1093/ije/31.6.1246. [PubMed] [Cross Ref]
30. SAS v9.4 TS Level 1M3 (Cary, North Carolina, USA).
31. Jensen, J. Income equivalences and the estimation of family expenditure on children. (Department of Social Welfare, Wellington 1988).
32. Salmond, C., Crampton, P. & Sutton, F. NZDep96 Index of Deprivation., (Health Services Research Centre, Wellington, 1998).
33. Salmond, C., Crampton, P. & Atkinson, J. NZDep2006 Index of Deprivation. (Department of Public Health, University of Otago, Wellington 2007).

Articles from Scientific Reports are provided here courtesy of Nature Publishing Group