Dynamic effects of smoking cessation on disease incidence, mortality and quality of life: The role of time since cessation
 Rudolf T Hoogenveen^{1}Email author,
 Pieter HM van Baal^{1},
 Hendriek C Boshuizen^{1} and
 Talitha L Feenstra^{1}
https://doi.org/10.1186/1478754761
© Hoogenveen et al; licensee BioMed Central Ltd. 2008
Received: 08 March 2007
Accepted: 11 January 2008
Published: 11 January 2008
Abstract
Background
To support health policy makers in setting priorities, quantifying the potential effects of tobacco control on the burden of disease is useful. However, smoking is related to a variety of diseases and the dynamic effects of smoking cessation on the incidence of these diseases differ. Furthermore, many people who quit smoking relapse, most of them within a relatively short period.
Methods
In this paper, a method is presented for calculating the effects of smoking cessation interventions on disease incidence that allows to deal with relapse and the effect of time since quitting. A simulation model is described that links smoking to the incidence of 14 smoking related diseases. To demonstrate the model, health effects are estimated of two interventions in which part of current smokers in the Netherlands quits smoking.
To illustrate the advantages of the model its results are compared with those of two simpler versions of the model. In one version we assumed no relapse after quitting and equal incidence rates for all former smokers. In the second version, incidence rates depend on time since cessation, but we assumed still no relapse after quitting.
Results
Not taking into account time since smoking cessation on disease incidence rates results in biased estimates of the effects of interventions. The immediate public health effects are overestimated, since the health risk of quitters immediately drops to the mean level of all former smokers. However, the longterm public health effects are underestimated since after longer periods of time the effects of past smoking disappear and so surviving quitters start to resemble never smokers. On balance, total health gains of smoking cessation are underestimated if one does not account for the effect of time since cessation on disease incidence rates. Not taking into account relapse of quitters overestimates health gains substantially.
Conclusion
The results show that simulation models are sensitive to assumptions made in specifying the model. The model should be specified carefully in accordance with the questions it is supposed to answer. If the aim of the model is to estimate effects of smoking cessation interventions on mortality and morbidity, one should include relapse of quitters and dependency on time since cessation of incidence rates of smokingrelated chronic diseases. A drawback of such models is that data requirements are extensive.
Background
Smoking is a risk factor for many major chronic diseases and reduces both length and quality of life [1]. In the Netherlands, about 28% of the Dutch population aged 15 or above smokes [2]. About 13% of the Dutch burden of disease and 3.7% of Dutch health care costs in 2003 could be attributed to smoking [3]. An effective tobacco control policy could thus substantially reduce the burden of disease. However, most of the effects of smoking on morbidity and mortality become manifest only after many years and not in all smokers. Since randomized trials to investigate the effect of smoking cessation interventions on disease occurrence and mortality are practically impossible, models based on current epidemiological knowledge and synthesizing data from many different sources are needed to estimate the effects of smoking cessation over time [4].
It was our aim to quantify the effects of smoking interventions on public health, taking into account the time since cessation of quitters. The latter is important, since many quitters relapse [5] and for most smoking related diseases the increased risks of former smokers only decrease gradually over time since cessation [6–8] Several models have been presented in the literature that deal with relapse and that include mortality risks of former smokers declining with time since quitting [9–13]. However, the potential use of these models is limited because they do not include all characteristics that are necessary to calculate the impact of smoking cessation on summary measures of population health combining morbidity and mortality, such as quality or disability adjusted life years. The SimSmoke model, for example, describes only all cause mortality and mortality from lung cancer [11, 12]. The effects of changing smoking class prevalence rates on other chronic diseases is calculated afterwards at each time step.
In this paper we present a model that is capable of describing the effects of smoking cessation on morbidity and mortality over time, and that overcomes the limitations of the models mentioned above. Our model describes the life course of quitters after smoking cessation taking into account relapse. It distinguishes the most important smoking related chronic diseases, and the incidence rates of former smokers depend on time since cessation. Moreover, it explicitly describes morbidity from these diseases by modeling the change of the disease prevalence rates over time. Morbidity, in turn, determines mortality in the model. By modeling diseases, effects of smoking cessation on quality of life and health care costs can simply be estimated by coupling health care costs and quality of life figures to diseases [14] instead of coupling them directly to smoking status.
The model presented in this paper is part of the RIVM Chronic Disease Model (CDM) [15]. The CDM is a statetransition Markovtype model that was designed to describe the effects of epidemiological risk factors on morbidity and mortality from several chronic diseases in the Dutch population. It includes 28 chronic diseases and several risk factors amongst which smoking, Body Mass Index, and physical inactivity. In modeling diseases explicitly, the structure of the model is similar to the Prevent model [16] and the recently presented Quit Benefits model [10, 17]. An important difference with the Prevent model is that also different risk factor classes are modeled. In comparison with the Quit Benefits Model our model includes more diseases, allows for comorbidity and has the ability to track health effects over a longer period. In the current study we describe how time since cessation is taken into account as an additional model parameter. In the next section, the model structure with respect to smoking is explained in detail. To illustrate the model, health effects are estimated of two interventions in which part of Dutch smokers quits smoking. To illustrate the strengths of the model its results are compared with the results of two simpler model versions. In the most simple version we assumed no relapse after quitting and equal incidence rates for all former smokers. To this were added incidence rates that depend on time since cessation. The final model was arrived at by also including relapse of quitters. In the discussion section, we will discuss the strengths and weaknesses of the model and elaborate on possible applications. Details on the mathematical structure of the model and its input data can be found in Appendices.
Methodology
Basic structure of the RIVM Chronic Disease Model (CDM)
The CDM is a model that describes the effects of risk factors, including smoking and overweight, on the incidence and mortality of chronic diseases in the population. It describes the effects for the total Dutch population taking into account birth and migration [15]. The CDM has been used for future projections of risk factor and disease prevalence numbers [18, 19, 19–21], cost effectiveness analyses [14] and estimates of healthy life expectancy [22]. The model describes the life course of cohorts in terms of transitions between risk factor classes and transitions between disease states over time. Risk factors and diseases are linked through relative risks of disease incidence. The CDM was formulated mathematically as a set of timecontinuous differential equations [23]. The model equations describe the 1year changes of the probability values for all risk factor classes and disease states in cohorts, specified by gender and age. The main model outcome variables are numbers of incident and prevalent cases and numbers of deaths, specified by disease, age and gender. To keep the number of model states manageable, the model describes the changes of the risk factor distributions and disease probabilities separately, i.e. as marginal distributions, but not the joint probability distribution function of all risk factor classes and diseases simultaneously.
The set of model equations consists of three components (see Appendix 1). In the initialization component the parameter values and the initial distribution of the population over all model states are calculated. In the simulation component the 1year changes of the model state prevalence numbers are calculated. These changes are the result of transitions between the risk factor classes and disease states. The transition numbers are computed as the 1year transition probabilities times the model state prevalence numbers at the start of the 1year timeinterval. Finally, in the postprocessing component the values of the output variables are calculated from the results of the simulation component.
Demographic data such as all cause mortality rates and initial population numbers were available from Statistics Netherlands [24]. To estimate incidence, prevalence and mortality rates in the general population, three types of data sources were used: general practitioner registrations for noncancer diseases, national cancer registries, and cohort studies for diabetes [25–27]. Noncancer mortality rates were estimated using a three state transition model [28]. To compute health effects in terms of qualityadjusted life years (QALYs), the CDM couples disability weights from the Dutch Burden of Disease Study to disease prevalences [22]. Disability weights reflect the severity and impact of a disease relative to death and optimal health, defined as absence of disability, and ranges from 0 (no disability) to 1 (death). The Dutch Burden of Disease Study estimated disability weights for 48 different disease categories, using the person tradeoff method. For the QALY calculations it was assumed that comorbidity reduces quality of life but that the effects are less than the sum of disability from the individual diseases[29, 30]. Health care costs were calculated by coupling estimated disease prevalence numbers to costs per patient per year per disease available from the Dutch Costs of Illness study [14, 31].
Model structure with respect to smoking
The CDM relates smoking to increased incidence rates of 14 smokingrelated chronic diseases, i.e. coronary heart disease (acute myocardial infarction (AMI) and other coronary heart disease), congestive heart failure, stroke, chronic obstructive pulmonary diseases (COPD), diabetes, and cancer of the lung, stomach, larynx, oral cavity, esophagus, pancreas, bladder and kidney. The incidence rates of smoking related diseases are increased in current smokers as well as in former smokers, with the relative risks of former smokers declining from the risk of a smoker immediately after stopping smoking to that of a never smoker as a function of time since cessation. Smoking status was defined as 'never smoker', 'current smoker', and 'former smoker'. The latter category was further subdivided into several classes based on the number of years since smoking cessation: less than 1 year ago, quitted between 1 and 2 years ago, and so forth, up to 20 years or more. The transitions involved were initiation (from 'never' to 'current' smoker), cessation (from 'current' to 'former, and quitted less than 1 year ago'), relapse (from any 'former' smoking class back to 'current'), and continuation as former smoker (e.g. from 'quitted between 2 and 3 years ago' to 'quitted between 3 and 4 years ago'). For all smokingrelated diseases we distinguished two states, i.e. without and with the disease. For each disease the transition involved was 'disease incidence', i.e. from 'diseasefree' to 'with the disease'. Finally, the model distinguishes outflow due to mortality from any cause. Mortality rates depend on the disease, i.e. persons with a disease have higher mortality rates than diseasefree persons. The parameters used in the model are the 1year probabilities of each transition between model states described.
In the next subsection, we describe in detail how the model parameters were estimated that depend on the time since cessation: relapse probabilities for former smokers and relative risks of smoking related diseases for former smokers.
Input data and estimation of model parameters dependent on time since cessation
The initial distribution over all smoking classes and the 1year smoking class transition probabilities were estimated from survey data from the Dutch Foundation on Smoking and Health [2]. This is a national survey on smoking which includes questions on current and past smoking status. Initiation and quit probabilities for neversmokers and current smokers were estimated using the retrospective data and are described elsewhere [32]. The initial smoking class probability values were calculated from the current smoking status. Since we had no input data on the initial distribution of all former smokers over time since smoking cessation, we generated this distribution by running the model initially for a birth cohort of nonsmokers (see Appendix 1). The generated distribution was used as input data.
It was assumed that the relapse rates of quitters decrease over the time since smoking cessation, and that the decrease is according to a negativeexponential curve: λ(s) = α β exp(β s) (1)
with:
s time since smoking cessation
λ smoking relapse rate
α, β regression coefficients.
The 1year probability of relapse was calculated by applying this hazard rate to a 1year time interval. The parameters α and β were estimated by fitting the probability of relapse to the data on the past history of former smokers [2, 33]. Details and results of the estimation method are given in Appendix 2.
For the risk factor class 'current smoking' relative risks, specified by gender and age, were used (see the online appendix of [22]). The relative risks of former smokers gradually decrease over time since smoking cessation in the following way. Immediately after cessation the relative risks are similar to those of current smokers. In the long run and conditional on survival the negative health effects of smoking disappear, and the relative risks are similar to those of never smokers, i.e. they attain value 1. Therefore, it was assumed that for each disease the relative risk decreases according to a negativeexponential curve starting from the value for current smokers and converging to the value 1. Furthermore, it was assumed that the rate of convergence decreases with age. This was done to account for the cumulative nature of the effects of smoking. We estimated these relative risks using the following regression model: RR_{former}(a,s) = 1 + (RR_{current}(a)  1) exp( γ(a) s) γ(a) = γ _{0} exp(η a) (2)
with:
a age
s time since smoking cessation
γ regression coefficient of time dependency
γ _{0}, η intercept and regression coefficient respectively of age dependency
RR_{former}(a) relative risks of disease incidence for former smokers
RR_{current}(a) relative risks of disease incidence for current smokers
The parameters γ _{0} and η were estimated in a twostep procedure, using relative risk data from major cohort studies reported in literature. In the first step we estimated the regression coefficient γ(a) for each study separately. Next, all regression coefficients γ(a) calculated were plotted in one graph against the mean age of the cohort. In the second step the regression coefficient γ _{0} and η were estimated from these γ(a) values. The details and results of the estimation method are given in Appendix 3.
Results

'young quitters' scenario: 10% of all current smokers aged 20–44 quit smoking in year 1;

'old quitters' scenario: 10% of all current smokers aged 45–70 quit smoking in year 1.

'simple' model: no relapse of quitters and relative risks of smoking related diseases do not depend on time since cessation.

'time dependent' model: relative risks of smoking related diseases depend on time since smoking cessation, but no relapse of quitters.

'time dependent & relapse' model: relapse of quitters and relative risks of smoking related diseases depend on time since smoking cessation.
Cumulative differences in (quality adjusted) life years between intervention and current practice scenarios over a period of 100 years
'young quitters' scenario  'old quitters' scenario  

'simple' model  'time dependent' model  'time dependent & relapse' model  'simple' model  'time dependent' model  'time dependent & relapse' model  
Life years gained (* 1,000)  525  746  439  277  317  187 
QALYs gained (* 1,000)  425  588  348  200  221  131 
Conclusion and Discussion
In this paper, we showed that not taking into account duration of smoking cessation time results in biased effect measures of interventions. The immediate public health effects are overestimated, since in reality quitters start with having the health risks of current smokers. The longterm public health effects are underestimated since after longer time periods the effects of past smoking have disappeared and so surviving quitters resemble more never smokers. On balance these two counteracting forces cause that health gains of smoking cessation are underestimated when one does not allow disease risks to depend on time since smoking cessation. This effect is more pronounced in younger populations.
These results show that simulation models are sensitive to assumptions regarding the model specification, in our case transition rates for former smokers. Therefore sensitivity analyses of simulation models should not only involve parameter values, but also model specifications should be taken into consideration. The model should be specified carefully to coincide with the questions it is supposed to answer. If the aim of the model is solely to project future health status in a business as usual scenario, a model without duration dependency might suffice. In the total population of former smokers, age correlates strongly with time since cessation. Therefore, using age specific relative risks without taking into account time since cessation would not introduce major errors if the aim is to make projections for the population of former smokers. Including time since smoking cessation becomes crucial when the aim of the modelling is to estimate effects of smoking cessation interventions instead of making projections.
A major drawback of the methodology employed is that more complex models require more detailed data, which are harder to get. For example, in the 'time dependent' models we used relative risks that were specified by both age and time since smoking cessation. The number of studies available to estimate these relative risks was much smaller than the number of studies available to estimate a single agedependent relative risk value for all former smokers combined. The reduction of absolute health risks after smoking cessation depends both on the absolute and relative disease risks for smoking. Since the absolute risks are largest for elderly ages, the reductions are most sensitive to the relative risk values for these ages. This was illustrated, for example, by the discussion on the attribution of US mortality numbers to obesity, where the attributed number of deaths were initially overestimated by assuming ageconstant relative risk values [34]. This was reason for us to analyze the effect of age on the reduction of health risks over time since smoking cessation. Indeed, empirical data (see Appendix 3) showed that the rate of decline decreased with age. As a result, the relative effects of smoking cessation are smaller for higher ages. Analogously, we analyzed the effect of age on the reduction of smoking relapse rates. However, the data available did not enable us to regress these rates on both time since cessation and age simultaneously. Moreover, these data were retrospectively obtained measurements, with the first measurement being the smoking status for the period 'more than 2 years ago'. As a result, we could do no better than assuming that the proportional reduction of the relapse rates was constant over the time since cessation and over age.
Including specific chronic diseases and relapse rates and relative risks that depend on time since smoking cessation makes the model presented well suited to calculate the impact of smoking interventions on morbidity and mortality results, when it is appropriate to take into account several effects simultaneously. Examples of these effects are dependency on age, dependency on time since smoking cessation, and competing morbidity and mortality risks. It also enabled us to calculate generic health measures such as disabilityadjusted life years. The model calculates the difference in quality adjusted life years between smokers and never smokers, and hence shows the health gain that can be gained by policy measures. By using a disease specific approach, the causal effects on quality of life and health care costs can be estimated and used in cost effectiveness analysis.
Appendix 1: Mathematical model equations
The model consists of three components. In the initialization component the parameter values and the initial distribution of the population over all model states are calculated from the input data. In the simulation component the 1year changes of the model state prevalence numbers are calculated for each cohort separately. Each cohort is defined by its initial age value, i.e. a(0). These changes in the number of persons for each state are the result of transitions between the risk factor classes and between the disease states. The numbers of transitions are computed as the 1year transition probabilities times the state prevalence numbers at the start of the 1year timeinterval. Finally, in the postprocessing component the model output variables are calculated from the results of the simulation component. All model parameters and variables are specified by gender and age, but we omit the index for gender below for reasons of readability.
Initialization component
with:
a age
c index over all smoking classes, c = never, current, former
i_{d,0}(a) data disease d incidence rate at age a
i_{d,base}(a) disease d incidence rate for never smokers at age a
n_{c,0}(a) initial smoking class probability values at age a
RR_{d}(c,a) relative risk of incidence of disease d for smoking class c
with:
d index over diseases
μ _{tot}(a) all cause mortality rates at age a
μ _{other}(a) mortality rates for other causes of death at age a
p_{d,0}(a) initial disease d prevalence rate values at age a
am_{d}(a) disease d related attributable mortality rates at age a
with:
a(0) age of cohort at initial time point t = 0
N_{c}(t) number of persons in smoking class c at time t,
for c = never, current
N_{pop,0}(a) initial total population numbers
p_{d}(t) disease d prevalence rates at time t
n_{c,0}(a) input proportions of population in smoking classes c
p_{d,0}(a) input disease d prevalence rates at age a
The initial values of the numbers of former smokers, stratified by time since cessation, were generated in a preprocessing step by running the model for a birth cohort without any disease included. In this way we calculated the distribution of former smokers n_{former}(s,a) over all cessation classes (s), specified by gender and age a. Doing so, we implicitly assumed that all smoking class transition probabilities are constant over time. Thus: N _{ former }(0,s) = n _{ former,0}(a(0))n _{ former }(s,a(0)))N _{ pop,0}(a(0)) (6)
with:
N_{former}(t,s) number of former smokers at time t in former smoking class s
s index over classes for time since smoking cessation, s = 1,...,S
E.g., s = 2 means former smoker stopped 1–2 years ago.
n_{former,0}(a) data proportion of former smokers in the population at age a
Simulation component
The simulation component describes the changes of the prevalence numbers in all smoking classes distinguished, as well as the changes of the prevalence rates for all chronic diseases included in the model. These changes are formulated as differential equations with 1year time steps.
1year changes of smoking class prevalence numbers
with:
t time parameter, with 1year steps
RR_{d}(s,a) disease d relative risk of former smoking class s
at age a, see Appendix 3
E(RR_{d}(former),t) mean disease d relative risk value at time t for former
smokers, at time t
E(RR_{d},t) mean disease d relative risk value in entire population at time t
with μ _{c}(t) all cause mortality rates for smoking class c at time t
with
λ _{start}(a), λ _{stop}(a) 1year start and stop smoking probabilities respectively,
at age a
λ _{relapse}(s) smoking relapse probabilities that depend on time since cessation class s, see appendix 2
The number of former smokers in the last class S at the end of the year are the sum of the numbers in the last and second last class at the start of the year that do not relapse or die.
1year changes in disease prevalence rates
We describe the 1year change in prevalence rates instead of numbers for each disease included. Since the mortality rates for other causes of death are assumed equal for persons with and without the disease, the change in the rate values depends only on the disease incidence and disease related excess mortality rates. The current disease incidence rates are the baseline disease incidence rates times the current mean relative risk value. i _{ d }(t) = E(RR _{ d },t)i _{ d,base }(a(t)) (11)
with:
i_{d}(t) disease d incidence rate at time t
i_{d,base}(a) baseline disease d incidence rate at age a
1year event incidence probabilities were calculated from the rate values assuming constant values over the year. These probabilities are still denoted as i_{d}(t). The prevalence rates change as a result of incidence and mortality. This equation is known as the DisModequation [28]. p _{ d }(t + 1) = p _{ d }(t) + i _{ d }(t)(1  p _{ d }(t))  em _{ d }(a(t))p _{ d }(t)(1  p _{ d }(t)) (12)
with: em_{d}(a) disease d related excess mortality probability at age a.
The parameter em_{d}(a) is the excess mortality related to disease d. It describes the additional mortality in the population with disease d as compared to the population without disease d (see the online appendix of [22]). The parameter am_{d}(a), that was defined before, describes the additional mortality rate, conditional on all risk factors and other disease states. The parameter am_{d}(a) can be interpreted as the mortality that uniquely can be attributed to disease d. It adjusts the excess mortality for mortality due to comorbid diseases. The part (1  p_{d}(t)) in the last term of the equation comes from describing changes of prevalence rates instead of numbers.
Model postprocessing component
The model output variables are computed from the results in the simulation component. These are the following:
with:
I_{d}(t) disease d incidence numbers during 1year period [t,t+1)
QALY qualityadjusted life years
w_{d}(a) disease d weight coefficients at age a
The weight coefficients w_{d}(a) describe the relative loss of quality of life value due to the disease. A value 0 means there is no loss of quality of life due to the disease; a value 1 indicates there is complete loss of quality of life, and that having the disease is no better than being dead.
Appendix 2: Relapse rates that depend on time since smoking cessation
First we present the statistical regression model that describes how relapse rates depend on the time since smoking cessation. Second, we describe how the regression parameters were estimated based on data available from a retrospective study on smoking behavior. Third, the resulting parameter estimates are presented.
Statistical model
The relapse rates are the rates at which former smokers restart smoking. They depend on the time since smoking cessation. Our formal model is defined in continuous time. Subsequently, the 1year probabilities of relapse used in our simulation model were calculated by integrating the rate values over time. Following the literature [1] we assumed that the proportional decrease of these relapse rates is constant over time, meaning that the relapse rates are highest shortly after cessation, and diminish in the long run. As a result of this assumption, the relapse rates follow a negativeexponential curve. λ _{ relapse }(s) = αβ exp( βs) (15)
with:
s time since smoking cessation [months]
λ _{relapse} smoking relapse rate
α, β regression coefficients.
Estimation procedure
The parameters α and β were estimated from a series of crosssectional population surveys on smoking behavior. The samples were representative for the Dutch population aged 15 years and above [2, 33]. Data were collected in the years 2000–2003. The complete questionnaire used is reproduced in Cappanesi et al.[32], and includes questions regarding current smoking status as well as past smoking behavior of the responders. These data were used to estimate the model parameters α and β in the following way.
Firstly, we calculated the numbers of quitters. The 1year quit probabilities were calculated from the survey by dividing the number of current nonsmokers that smoked 1 year ago by the total number of smokers 1 year ago, based on answers in the retrospective survey. These probabilities were multiplied with the proportion of current smokers available from the survey and the total population numbers available from Statistics Netherlands [24] to get the numbers of quitters for each age class. N _{ stop }(a) = λ _{ stop }(a)p _{ current }(a)N(a) (18)
with:
a age
N_{stop}(a) 1year number of quitters at age a
p_{current}(a) proportion current smokers at age a
N(a) total population numbers at age a
λ _{stop}(a) 1year quit probability of current smokers
Given trial values α and β, we calculated for each group of quitters specified by age at cessation how many quitters were still abstinent over time using the probability of relapse function. N _{ former }(a,dα,β) = N _{ stop }(a  d)S(dα,β) (19)
with: N_{former}(a,dα,β) number of persisting former smokers at age a that
quitted d months ago
with:
p_{calc}(a,dα,β) calculated proportion of former smokers at age a who
quitted d months ago
with:
p_{emp}(da) the empirical distribution of all former smokers at age a over
time since cessation classes d
SS(α,β) weighted sum of squares
At first, we estimated the parameters α and β for several age classes separately to check for different declines of smoking relapse rates over age. However, since the number of former smokers was too small for lower and middle ages, we concluded that we could not identify any effect of age, and thus assumed parameter values that were different between both sexes but constant over age.
Results
Estimated parameter values for $\widehat{\alpha}$ and $\widehat{\beta}$ for men are 1.177 and 0.150. Estimated parameter values for $\widehat{\alpha}$ and $\widehat{\beta}$ for women are 1.197 and 0.113.
Appendix 3: Relative risks that depend on time since smoking cessation
In this appendix the methods used to calculate relative risks of former smokers that depend on time since smoking cessation are described. Two different methods were used to estimate the regression parameters, depending on the data available. The first method was used for diseases with sufficient data on relative risks of former smokers specified by time since cessation. The second method was used when these data were not available. In the latter case we estimated the regression parameter from the distribution of all former smokers over time since cessation (see Appendix 1) and mean relative risks of all former smokers.
Statistical model

The relative risk of quitters equals the relative risk of current smokers.

The relative risk of former smokers approaches the relative risk of never smokers, i.e. value 1.

Relative risks of former smokers show a timeconstant proportional decrease.

The proportionality coefficients that describe the rate of decrease over time of the relative risks decrease proportionally over age
with:
a age
a*(a) transformation of a, a*(a) = (a50)^{+}: the nonnegative value of a50
s time since smoking cessation
γ regression coefficient of time dependency
η regression coefficient of age dependency
RR_{current}(a) relative risks of current smokers at age a
RR_{former}(a) mean relative risks of all former smokers at age a
The parameter γ _{0} is the reference value of parameter γ, i.e. the value of γ for age 50 years. We chose age 50 years as the origin, since we could not identify any age gradient for lower ages (see below). As a consequence, we assumed ageconstant rates of decrease over time for these ages. Thus, in the statistical model applied in our simulation model the term (a50) was replaced by the term (a50)^{+}.
Estimation of regression parameters from cohort studies stratified by time since cessation
To estimate the regression coefficients γ and η for each smokingrelated relative risk, we used data from major cohort studies presented in literature [6–8, 35–42]. To be included, studies had to present relative risks that were adjusted for at least gender and age, and the time since smoking cessation had to be subdivided in at least three classes with reported cutoff points. We estimated the regression coefficient γ by the method of least squares. For each study we checked whether the published relative risk of current smokers matched the extrapolated one based on the curve of relative risks of former smokers. If so, we added this relative risk value (with time since smoking cessation 0) to the reported ones for former smokers, and reestimated the rate of proportional decrease. We included these relative risks of current smokers to improve the precision of our estimation. Only for some forms of cancer among which lung cancer we found differences between reported and extrapolated relative risks of current smokers. These differences point at a reversal of causality: getting lung cancer is a reason to quit smoking. We plotted the calculated relative risk curves and data points regressed on the time since smoking cessation to check the model fit for each study separately.
Next, all regression coefficients γ calculated were plotted regressed on age to check the assumption of ageconstant proportional decrease. Our statetransition model uses 1year timesteps, whilte most cohort studies have much longer followup time periods. Thus, in the cohort studies the age at baseline is not the age at event. Therefore we used the estimated age at event for age a in the formula on γ(a) instead of the reported age at baseline. The age at event was calculated assuming a Gompertztype event rate. I.e., we assumed that the event rate λ increases exponentially with age a, i.e. λ(a) ∝ exp(a c). We assumed that the rate of increase of all event rates with age was equal to the one for all cause mortality. We estimated this rate value using Dutch mortality data that were available from Statistics Netherlands.
If we found an agedependent proportional decrease of the rate of change of relative risks with time since cessation, we estimated the regression coefficient η by the method of weighted least squares. For ages lower than 50 years we had not enough data to identify any relation with age, and thus assumed constant values of parameter γ. Only for all cause mortality we had sufficient data to identify the model parameters specified by age. Based on the data we could find no differences between men and women.
Estimation of regression parameters from one relative risk value for all former smokers
If no data were available from cohort studies we estimated the regression parameters from reported relative risk values of all former smokers together, and used the calculated distribution of all former smokers over time since cessation (see appendix 1). The parameters were estimated by the method of weighted least squares, where the weights were the number of former smokers for each age class.
with:
η _{0} transformed parameter value; η _{0} = η γ _{0}
RR_{former}(a) the data mean relative risk value of all former smokers
The latter equation was solved within the framework of weighted linear regression: (γ _{0}, η _{0}) = (XWX^{ T })^{1} X^{ T } Wy (26)
with: X design matrix, i.e. with as columns the independent variables E_{former}(sa)
and a*(a) E_{former}(sa) for (rows) all age classes
W diagonal weight matrix, with values n_{former,0}(a) N_{pop,0}(a), the number of
former smokers for all age classes
y response vector, i.e. with as elements the empirical values of the right hand side of equation (25).
We assumed that the decline of the relative risks of former smokers decreased with age. As a result, if equation (24) resulted in a negative value of parameter η _{0} we reestimated parameter γ _{0} with fixed value η _{0} = 0.
Results
Results of parameter estimates per smoking related disease
γ _{0}  η  

AMI  0.24228  0.05822 
Other CHD  0.24228  0.05822 
CHF*  0.0421371  0 
CVA  0.31947  0.01648 
COPD  0.20333  0.03087 
Diabetes*  0.024811  0 
Lung cancer  0.15637  0.02065 
Stomach cancer*  0.0264112  0 
Esophagus cancer*  0.0537424  0 
Larynx cancer*  0.0279918  0 
Uriny/bladder cancer  0.05417  0 
Kidney cancer*  0.0385957  0 
Pancreas cancer  0.09279  0 
Oral cavity cancer*  0.0493028  0 
Declarations
Acknowledgements
We would like to thank Ardine de Wit, Geert Jan Kommer, Peter Engelfriet and both reviewers for their critical readings of the manuscript.
Authors’ Affiliations
References
 US Department of Health and human services PHSCfDCCfCDPaHP: Reducing Tobacco Use: A Report of the Surgeon General. 2000.Google Scholar
 STIVORO: Roken, de harde feiten: Volwassenen. Den Haag: STIVORO; 2003.Google Scholar
 de Hollander A, Hoeymans N, Melse J, van Oers J, Polder Je: Zorg voor gezondheid – Volksgezondheid Toekomst Verkenning 2006. RIVM Rapport 270061003. 2006.Google Scholar
 Buxton MJ, Drummond MF, Van Hout BA, Prince RL, Sheldon TA, Szucs T, Vray M: Modelling in economic evaluation: an unavoidable fact of life. Health Econ 1997, 6: 217–27. 10.1002/(SICI)10991050(199705)6:3<217::AIDHEC267>3.0.CO;2WPubMedView ArticleGoogle Scholar
 Hughes JR, Keely J, Naud S: Shape of the relapse curve and longterm abstinence among untreated smokers. Addiction 2004, 99: 29–38. 10.1111/j.13600443.2004.00540.xPubMedView ArticleGoogle Scholar
 Enstrom JE: Smoking cessation and mortality trends among two United States populations. J Clin Epidemiol 1999, 52: 813–25. 10.1016/S08954356(99)000402PubMedView ArticleGoogle Scholar
 Kawachi I, Colditz GA, Stampfer MJ, Willett WC, Manson JE, Rosner B, Speizer FE, Hennekens CH: Smoking cessation and time course of decreased risks of coronary heart disease in middleaged women. Arch Intern Med 1994, 154: 169–75. 10.1001/archinte.154.2.169PubMedView ArticleGoogle Scholar
 Kawachi I, Colditz GA, Stampfer MJ, Willett WC, Manson JE, Rosner B, Speizer FE, Hennekens CH: Smoking cessation and decreased risk of stroke in women. JAMA 1993, 269: 232–6. 10.1001/jama.269.2.232PubMedView ArticleGoogle Scholar
 Ahmad S: Increasing excise taxes on cigarettes in California: a dynamic simulation of health and economic impacts. Prev Med 2005, 41: 276–83. 10.1016/j.ypmed.2004.10.024PubMedView ArticleGoogle Scholar
 Hurley SF, Matthews JP: The Quit Benefits Model: a Markov model for assessing the health benefits and health care cost savings of quitting smoking. Cost Eff Resour Alloc 2007, 5: 2. 10.1186/1478754752PubMed CentralPubMedView ArticleGoogle Scholar
 Levy DT, Friend KB: A simulation model of tobacco youth access policies. J Health Polit Policy Law 2000, 25: 1023–50. 10.1215/036168782561023PubMedView ArticleGoogle Scholar
 Levy DT, Cummings KM, Hyland A: Increasing taxes as a strategy to reduce cigarette use and deaths: results of a simulation model. Prev Med 2000, 31: 279–86. 10.1006/pmed.2000.0696PubMedView ArticleGoogle Scholar
 Tengs TO, Osgood ND, Lin TH: Public health impact of changes in smoking behavior: results from the Tobacco Policy Model. Med Care 2001, 39: 1131–41. 10.1097/0000565020011000000010PubMedView ArticleGoogle Scholar
 van Baal PH, Feenstra TL, Hoogenveen RT, Ardine de Wit G, Brouwer WB: Unrelated medical care in life years gained and the cost utility of primary prevention: in search of a 'perfect' costutility ratio. Health Econ 2007.Google Scholar
 Hoogenveen RT, de Hollander AEM, van Genugten MLL: The chronic disease modelling approach. RIVM Rapport 266750001. 1998.Google Scholar
 Barendregt JJ: PREVENT: the technical background. In Public Health Models: Tools for Health Policy making at National and European level. Nijkerk: Callenbach; 1999.Google Scholar
 Barendregt JJ, Bonneux L, van der Maas PJ: The health care costs of smoking. N Engl J Med 1997, 337: 1052–7. 10.1056/NEJM199710093371506PubMedView ArticleGoogle Scholar
 Struijs JN, van Genugten ML, Evers SM, Ament AJ, Baan CA, van den Bos GA: Modeling the future burden of stroke in The Netherlands: impact of aging, smoking, and hypertension. Stroke 2005, 36: 1648–55. 10.1161/01.STR.0000173221.37568.d2PubMedView ArticleGoogle Scholar
 Mulder I, Genugten MLLv, Hoogenveen RT, Hollander AEMd, BuenodeMesquita HB: The impact of smoking on future pancreatic cancer: a computer simulation. Ann Oncol 1999,10(Suppl 4):74–8. 10.1023/A:1008397005267PubMedView ArticleGoogle Scholar
 Struijs JN, van Genugten ML, Evers SM, Ament AJ, Baan CA, van den Bos GA: Modeling the future burden of stroke in The Netherlands: impact of aging, smoking, and hypertension. Stroke 2005, 36: 1648–55. 10.1161/01.STR.0000173221.37568.d2PubMedView ArticleGoogle Scholar
 Feenstra TL, van Genugten ML, Hoogenveen RT, Wouters EF, Ruttenvan Molken MP: The impact of aging and smoking on the future burden of chronic obstructive pulmonary disease: a model analysis in the Netherlands. Am J Respir Crit Care Med 2001, 164: 590–6.PubMedView ArticleGoogle Scholar
 van Baal PH, Hoogenveen RT, de Wit GA, Boshuizen HC: Estimating healthadjusted life expectancy conditional on risk factors: results for smoking and obesity. Popul Health Metr 2006, 4: 14. 10.1186/14787954414PubMed CentralPubMedView ArticleGoogle Scholar
 Lambert JD: Numerical methods for ordinary differential systems: the initial value problem. Chichester, Wiley; 1991.Google Scholar
 Statistics Netherlands: Statline. [http://www.cbs.nl]
 Hoogenveen RT, Gijsen R: Dutch DisMod for several types of cancer. RIVM Rapport 260751004. 2000.Google Scholar
 Hoogenveen RT, Gijsen R, van Genugten MLL, Kommer GJ, Schouten JSAG, de Hollander AEM: Dutch DisMod. Constructing a set of consistent data for chronic disease modelling. RIVM Rapport 260751001. 2000.Google Scholar
 Hoogenveen RT, Gijsen R, van Genugten MLL, Kommer GJ, Schouten JSAG, de Hollander AEM: Dutch DisMod. Constructing a set of consistent data for chronic disease modelling. Bilthoven: National Institute of Public Health and the Environment (RIVM) 2000.Google Scholar
 Barendregt JJ, Van Oortmarssen GJ, Vos T, Murray CJ: A generic model for the assessment of disease epidemiology: the computational basis of DisMod II. Popul Health Metr 2003, 1: 4. 10.1186/1478795414PubMed CentralPubMedView ArticleGoogle Scholar
 Melse JM, EssinkBot ML, Kramers PG, Hoeymans N: A national burden of disease calculation: Dutch disabilityadjusted lifeyears. Dutch Burden of Disease Group. Am J Public Health 2000, 90: 1241–7.PubMed CentralPubMedView ArticleGoogle Scholar
 Stouthard MEA, EssinkBot ML, Bonsel GJ, Barendregt JJ, Kramer PGN, van de Water HPA, GunningScheepers LJ, van der Maas PJ: Disability Weights for Diseases in The Netherlands. Department of Public Health. Erasmus University Rotterdam 1997.Google Scholar
 Slobbe LCJ, Kommer GJ, Smit JM, Groen J, Meerding WJ, Polder JJ: Kosten van ziekten in Nederland 2003. Zorg voor euro's – 1. RIVM Rapport 270751010. 2006.Google Scholar
 Capannesi M, Boshuizen HC, Willemsen MCvHJC: How to obtain long term projections for smoking behavior: a case study in the Dutch population. , in press.Google Scholar
 STIVORO: Roken, de harde feiten: Volwassenen. Den Haag: STIVORO; 2002.Google Scholar
 Flegal KM, Graubard BI, Williamson DF: Methods of calculating deaths attributable to obesity. Am J Epidemiol 2004, 160: 331–8. 10.1093/aje/kwh222PubMedView ArticleGoogle Scholar
 Ando M, Wakai K, Seki N, Tamakoshi A, Suzuki K, Ito Y, Nishino Y, Kondo T, Watanabe Y, Ozasa K, Ohno Y, JACC Study Group: Attributable and absolute risk of lung cancer death by smoking status: findings from the Japan Collaborative Cohort Study. Int J Cancer 2003, 105: 249–54. 10.1002/ijc.11043PubMedView ArticleGoogle Scholar
 BenShlomo Y, Smith GD, Shipley MJ, Marmot MG: What determines mortality risk in male former cigarette smokers? Am J Public Health 1994,84(8):1235–1242.PubMed CentralPubMedView ArticleGoogle Scholar
 Bosetti C, Franceschi S, Levi F, Negri E, Talamini R, la VecchiaC: Smoking and drinking cessation and the risk of oesophageal cancer. Br J Cancer 2000, 83: 689–91. 10.1054/bjoc.2000.1274PubMed CentralPubMedView ArticleGoogle Scholar
 Coughlin SS, Neaton JD, Sengupta A: Cigarette smoking as a predictor of death from prostate cancer in 348,874 men screened for the Multiple Risk Factor Intervention Trial. Am J Epidemiol 1996,143(10):1002–1006.PubMedView ArticleGoogle Scholar
 Kawachi I, Colditz GA, Stampfer MJ, Willett WC, Manson JE, Rosner B, Speizer FE, Hennekens CH: Smoking cessation and decreased risk of stroke in women. JAMA 1993, 269: 232–6. 10.1001/jama.269.2.232PubMedView ArticleGoogle Scholar
 Wakai K, Seki N, Tamakoshi A, Kondo T, Nishino Y, Ito Y, Suzuki K, Ozasa K, Watanabe Y, Ohno Y, JACC Study Group. Japan Collaborative Cohort Study: Decrease in risk of lung cancer death in males after smoking cessation by age at quitting: findings from the JACC study. Jpn J Cancer Res 2001, 92: 821–8.PubMedView ArticleGoogle Scholar
 Wannamethee SG, Shaper AG, Whincup PH, Walker M: Smoking Cessation and the Risk of Stroke in Middleaged Men. JAMA 1995, 274: 155–60. 10.1001/jama.274.2.155PubMedView ArticleGoogle Scholar
 Zeegers MPA, Goldbohm RA, van den BrandtPA: A prospective study on active and environmental tobacco smoking and bladder cancer risk (The Netherlands). Canc Caus Cont 2002, 13: 83–90. 10.1023/A:1013954932343View ArticleGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.