Age at quitting smoking as a predictor of risk of cardiovascular disease incidence independent of smoking status, time since quitting and pack-years
BMC Research Notes volume 4, Article number: 39 (2011)
Risk prediction for CVD events has been shown to vary according to current smoking status, pack-years smoked over a lifetime, time since quitting and age at quitting. The latter two are closely and inversely related. It is not known whether the age at which one quits smoking is an additional important predictor of CVD events. The aim of this study was to determine whether the risk of CVD events varied according to age at quitting after taking into account current smoking status, lifetime pack-years smoked and time since quitting.
We used the Cox proportional hazards model to evaluate the risk of developing a first CVD event for a cohort of participants in the Framingham Offspring Heart Study who attended the fourth examination between ages 30 and 74 years and were free of CVD. Those who quit before the median age of 37 years had a risk of CVD incidence similar to those who were never smokers. The incorporation of age at quitting in the smoking variable resulted in better prediction than the model which had a simple current smoker/non-smoker measure and the one that incorporated both time since quitting and pack-years. These models demonstrated good discrimination, calibration and global fit. The risk among those quitting more than 5 years prior to the baseline exam and those whose age at quitting was prior to 44 years was similar to the risk among never smokers. However, the risk among those quitting less than 5 years prior to the baseline exam and those who continued to smoke until 44 years of age (or beyond) was two and a half times higher than that of never smokers.
Age at quitting improves the prediction of risk of CVD incidence even after other smoking measures are taken into account. The clinical benefit of adding age at quitting to the model with other smoking measures may be greater than the associated costs. Thus, age at quitting should be considered in addition to smoking status, time since quitting and pack-years when counselling individuals about their cardiovascular risk.
CVD risk associated with smoking varies not only with smoking status but also with intensity and duration of smoking or smoking pack-years, time since quitting and age at quitting. Many studies have examined the lag in health benefit of smoking cessation measured by time since quitting and occurrence of a CVD event [1–8]. However age at quitting also affects the health benefits of smoking cessation [9, 10]. The risks of mortality and smoking-related disease increase with age at quitting. However, the role of age at quitting as a predictor of CVD risk in the presence of time since quitting and pack-years is unclear. As CVD is more common among elderly people it is likely that age at quitting and time since quitting, inversely correlated, influence CVD risk in opposite directions.
This study explores whether and how age at quitting influences risk of CVD incidence, using data from the Framingham Offspring Heart Study. It uses a smoking status variable with and without incorporating other smoking variables such as time since quitting among past smokers and pack-years among current smokers, while controlling for common risk factors.
Methods & Results
Study Design and Sample
The Framingham Offspring Heart Study details for design, selection criteria, examination procedures and criteria for CVD events have been described elsewhere [11–15]. Participants were eligible for the present study if at examination 4 (1988 to 1992) they were CVD-free and aged 30-74 with nonmissing data on covariates. The final study sample consisted of 3751 participants (mean age 51.61; 1937 women).
Measurement of CVD Risk Factors
The risk factors included were smoking status with various definitions (expanded below), systolic and diastolic blood pressure (SBP & DBP), total cholesterol/high-density-lipoprotein (HDL) ratio or both (depending on which provided a better prediction of outcome), age, sex, diabetes status and body-mass index (BMI). Smoking status was initially defined as a dichotomous current smoker/non-smoker variable. The other definitions of smoking status included four, six and eight categories. The four category smoking status variable was defined as: never smokers, former smokers with age at quitting below 37 years, former smokers with age at quitting of 37 years or older, and current smokers. The six category smoking status variable was defined as: never smokers, former smokers with time since quitting 5 years or less and over 5 years and current smokers with under 20, 20-39, and 40 or more pack-years where current pack-years were calculated by dividing the number of cigarettes being smoked per day by 20 to obtain an estimate of "packs" and multiplying this by the number of years a person was a smoker. The eight category smoking status variable was defined as: never smokers, past smokers into 4 groups with two levels of age at quitting (≤44, >44 years) at each of two levels of time since quitting (≤5, and >5 years) and current smokers with under 20, 20-39, and 40 or more pack-years.
Blood pressure was the average of two physician-obtained measures. Cholesterol and various smoking measures were based on standardized enzymatic methods and self-report, respectively. Diabetes was defined as a fasting glucose ≥126 mg/dL. Age at quitting and time since quitting were calculated at examination 4 by combining smoking status information at each examination with history of smoking status from examination 1 .
Development & assessment of predictive models
The Cox proportional-hazards model  was used to relate risk factors to the risk of CVD incidence during follow-up from examinations 4 to 7. The assumption of proportionality of hazards was satisfied; tested by taking interaction between a covariate and log (survival time)  and plotting Schoenfeld residuals against survival time.
To improve the interpretability of the predictive models we categorized time since quitting, age at quitting and pack-years. It was observed that the lag time for a beneficial effect of smoking cessation on risk of CVD incidence was five years after which the risk stabilized . In the literature there is no maximum age for quitting without increasing the risk of a CVD event compared to a never smoker and there was no apparent cutpoint to dichotomise this variable as the predicted time to the onset of a CVD event declined almost linearly with age at quitting (results not shown). Thus the median, shown by simulation to result in minimum loss in efficiency , was used to dichotomise age at quitting.
Four models were fitted for the outcome. Each included a composite measure of smoking status and all other risk factors found to be significantly related to the outcome. Model 1 included smoking status as a simple current smoker/non-smoker variable with current non-smoker as the reference category. To incorporate the effect of age at quitting into smoking status, quitters were separated from current non-smokers and categorized by age at quitting in Model 2, which incorporated smoking status with categories <37 and ≥37 years for age at quitting, never smokers and current smokers. To examine whether incorporating age at quitting to smoking status improved risk prediction, Model 2 was compared to Model 1. To examine whether incorporating age at quitting improves risk prediction to a model which already includes time since quitting and pack-years in smoking status, Models 3 and 4 were fitted and compared. Model 3 incorporated smoking status that included categories for never smokers, ≤5 and >5 years for time since quitting, and <20, 20-39 and 40+ for pack-years. Other categorizations for pack-years and time since quitting were found to be less effective in terms of predictive ability. Model 4 added age at quitting to Model 3 with smoking status having six categories - never smoker, current smoker, and past smokers quitting ≤5 years and whose age at quitting was ≤ 44 or >44 years, and those quitting >5 years and whose age at quitting was ≤ 44 or >44 years. Compared with Model 2 age at quitting was categorized differently because for the initial categorization of age at quitting at < 37 and ≥37 years there were inadequate numbers of cases in one of the joint categories of this variable with time since quitting resulting in inefficient estimation of its regression coefficient. The other cutpoints prior to reaching age 44 produced the same result until the cutpoint reached age 44 which did not yield inadequate number of cases in any of the joint categories. In Models 2 through 4 the reference category for smoking status was never smoker.
For assessing the discriminative ability of a model and improvement between two nested models we used Harrell's c statistic [19, 20] and a test for difference in two correlated c statistics . Large 'independent' association of the new covariate with the outcome is required to result in a meaningfully larger c statistic [22–24] for models possessing reasonably good discrimination, and the c statistic does not assist a physician in treatment decisions about an individual [25, 26] while reclassification statistics NRI  and IDI  do [25, 26]. Thus, we used the latter to supplement c-statistic analyses [27, 28]. For calculating NRI we assessed risk reclassification  by sorting the predicted risk for each model into four clinically meaningful categories (<6%, 6% to < 10%, 10% to < 20%, and ≥ 20%). The benefit and cost of using a new model compared to a baseline model can be measured by the proportions of subjects with and without subsequent events, respectively, who are classified as high risk (eg. ≥ 20%) according to the new model . There was negligible overoptimism in c and NRI estimates obtained by bootstrapping as these were less than 0.007 and 0.005 respectively.
For assessing calibration of the fitted models and improvement in global fit between two nested models we computed the Hosmer-Lemeshow statistic and its modification  and likelihood ratio test respectively. Neither Models 3 nor 4 included current age as a covariate because it had exact collinearity with time since quitting and age at quitting.
The sample risk factor characteristics at baseline examination 4 are shown in Table 1. The sample consists of 26.7% never smokers, 48.4% quitters (15.1% of whom quit within 5 years of the baseline measurement and 33.1% of whom quit before age 37 years), and 14.5% current smokers (of whom 16.6% have exposure ≥40 pack-years).
Table 2 shows that Model 2 improved predictive ability significantly compared to Model 1. Model 3 performed well in terms of model discrimination and overall fit but less well in terms of calibration. Model 4 performed well on all model performance indicators; significantly improving predictive ability compared to Model 3 (Table 3). Thus, age at quitting was an independent predictor of risk of CVD incidence regardless of including time since quitting and pack-years in the model.
Compared to never smokers, the risk of CVD incidence based on Model 2 was 7.3% higher (RR = 1.073, 95% CI 0.804 ~ 1.433) for those who quit before age 37 years and 58.1% higher (RR = 1.581, 95% CI 1.193 ~ 2.094) for those who quit at least at age 37 years (Table 4). For the former category the relative risk was not significantly different from the never smokers while for the latter category it was. Based on the final model (Model 4), the risk among those quitting more than 5 years prior to the baseline exam and whose age at quitting was 44 years or less was close to never smokers. Risk among those quitting within 5 years prior to the baseline exam and whose age at quitting was over 44 years was about three times higher than that of never smokers (Table 5).
Reclassification of subjects
This section describes how many subjects were reclassified overall and with respect to 'high risk' category of ≥ 20% when we compared the preferred full model against the reference model. Comparing Model 2 against Model 1, for participants who experienced a CVD event, the net gain in reclassification proportion was significantly different from zero (p = 0.0113) (Table 6) and significant for participants who did not experience an event (p = 0.0025) and for all participants (p = 0.0112). For those who experienced a CVD event, using Model 4 rather than Model 3 did not improve net gain in reclassification proportion significantly (p = 0.2935) (Table 7). The result was similar for participants who did not experience an event (p = 0.1545) and for all participants (p = 0.1558).
Table 8 shows that based on Model 2 instead of Model 1, 16.5% of those developing a CVD event would have moved up to the 'high risk' category of ≥20% while of those not having a CVD event 9.4% would have moved to this risk category, the difference of which is highly significant (p < 0.0001). Similarly, Table 9 shows that if we had used Model 4 rather than Model 3, 14.6% of those who develop CVD would be appropriately assessed for their cardiovascular risk while only 7.6% of those who do not develop CVD would be falsely assessed for their cardiovascular risk, the difference of which is highly significant (p < 0.0001).
Sensitivity of the results
We have adjusted for all major confounders of smoking to address confounding bias in the risk models. To address the possibility of distortion due to medical treatments affecting the risk of a CVD event we found that the regression coefficients of the models were fairly insensitive to the inclusion of cardioactive medications. To address the possibility of reverse causation, we excluded from the baseline cohort those with a cancer history and other non-CVD conditions. This did not substantially influence the results. Sub-analyses conducted by excluding from baseline cohort those smokers who quit after examination 4 and those quitters who took up smoking after examination 4, and later those current smokers from baseline cohort whose pack-years changed substantially in subsequent examinations did not influence our results.
Merits and demerits of this study
The study's key strength is that it not only evaluates improvement in predicting CVD risk when models incorporate age at quitting but also quantifies the proportions of people receiving clinical benefits and costs. However, a cost-benefit analysis of including this variable was not possible as the same number of CVD events was prevented by the full and reduced models. Also, as the Framingham cohort has an ethnically white predominance the generalizability of our models to other ethnic groups is unknown.
The incorporation of age at quitting in smoking status resulted in better prediction compared to the model which had a current smoker/non-smoker measure and to the model which incorporated both time since quitting and pack-years in smoking status. Thus, age at quitting was an independent predictor of CVD incidence even after accounting for time since quitting and pack-years.
We also showed that if we had incorporated age at quitting in smoking status instead of a current/non-smoker measure, a significantly higher proportion of those developing a CVD event would have moved up to the 'high risk' category compared to those not having a CVD event who moved up to this category. The result was similar if the model added age at quitting in smoking status which already incorporated time since quitting and pack-years. The former would be appropriately treated while the latter would be falsely treated if we included age at quitting in smoking status. Those appropriately treated can benefit from additional screening for CVD risk and would require more aggressive intervention for smoking cessation  and would thus aid in preventing more deaths. However, this benefit would be at the cost of falsely identifying people who do not develop CVD as high risk who may unnecessarily receive additional screening and may cause undue stress and burden to the smoking cessation programs. From a CVD prevention perspective the benefits associated with smoking cessation clearly outweigh the costs for CVD screening and smoking cessation programs. Age at quitting should be taken into account, as well as other smoking measures, when counselling individuals about their cardiovascular risk.
Lightwood JM, Glantz SA: Short-term economic and health benefits of smoking cessation: myocardial infarction and stroke. Circulation. 1997, 96 (4): 1089-96.
Ockene JK, Kuller LH, Svendsen K, Meilahn M: The Relationship of Smoking Cessation to Coronary Heart Disease and Lung Cancer in the Multiple Risk Factor Intervention Trial (MRFIT). AJPH. 1990, 80 (8): 10.2105/AJPH.80.8.954.
Ockene IS, Miller NH: Cigarette Smoking, Cardiovascular Disease, and Stroke: A Statement for Healthcare Professionals From the American Heart Association. Circulation. 1997, 96: 3243-3247.
Wannamethee SG, Shaper AG, Whincup PH, Walker M: Smoking cessation and the risk of stroke in middle-aged men. JAMA. 1995, 274: 155-160. 10.1001/jama.274.2.155.
Tamura U, Tanaka T, Okamura T, Kadowaki T, Yamato H, Tanaka H, Nakamura M, Okayama A, Ueshima H, Yamagata Z: Changes in weight, cardiovascular risk factors and estimated risk of coronary heart disease following smoking cessation in Japanese male workers: HIPOP-OHP study. Journal of Atherosclerosis & Thrombosis. 2010, 17 (1): 12-20.
Teo KK, Ounpuu S, Hawken S, Pandey MR, Valentin V, Hunt D, Diaz R, Rashed W, Freeman R, Jiang L, Zhang X, Yusuf S, INTERHEART Study Investigators: Tobacco use and risk of myocardial infarction in 52 countries in the INTERHEART study: a case-control study. Lancet. 2006, 368 (9536): 647-58. 10.1016/S0140-6736(06)69249-0.
Honjo K, Iso H, Tsugane S, Tamakoshi A, Satoh H, Tajima K, Suzuki T, Sobue T: The effects of smoking and smoking cessation on mortality from cardiovascular disease among Japanese: pooled analysis of three large-scale cohort. Tobacco Control. 2010, 19: 50e57-10.1136/tc.2009.029751.
Mannan H, Stevenson C, Peeters A, Walls H, McNeil J: Framingham risk prediction equations for CVD incidence using detailed measures for smoking. Heart International. 2010, 5 (2): 49-57. 10.4081/hi.2010.e11.
Doll R, Peto R, Boreham J, Sutherland I: Mortality in relation to smoking: 50 years' observations on male British doctors. BMJ. 2004, 328: 1519-10.1136/bmj.38142.554479.AE. (26 June) AE (published 22 June 2004)
Ostbye T, Taylor DH: The Effect of Smoking on Years of Healthy Life (YHL) Lost among Middle-Aged and Older Americans. Health Services Research. 2004, 39: 531-551. 10.1111/j.1475-6773.2004.00243.x.
Wood D, Wray R, Poulter N: JBS 2: Joint British Societies' guidelines on prevention of cardiovascular disease in clinical practice. Heart. 2005, 91 (Suppl 5): v1-v52.
D'Agostino S, Grundy S, Sullivan LM, Wilson P, for the CHD Risk Prediction Group: Validation of the Framingham coronary heart disease prediction scores: results of a multiple ethnic groups investigation. JAMA. 2001, 286: 180-187.
Beswick A, Brindle P: Risk scoring in the assessment of cardiovascular risk. Curr Opin Lipidol. 2006, 17: 375-386. 10.1097/01.mol.0000236362.56216.44. 1999; 159: 1197-1204
Liu J, Hong Y, D'Agostino RB, Wu Z, Wang W, Sun J, Wilson PW, Kannel WB, Zhao D: Predictive value for the Chinese population of the Framingham CHD risk assessment tool compared with the Chinese Multi-Provincial Cohort Study. JAMA. 2004, 291: 2591-2599. 10.1001/jama.291.21.2591.
D'Agostino RB, Vasan RS, Pencina MJ, Wolf PA, Cobain M, Massaro JA, Kannel WB: General Cardiovascular Risk Profile for Use in Primary Care. Circulation. 2008, 117: 743-753.
Cox DR: Regression models and life tables. J Royal Stat Soc. 1972, 34 (series B): 187-220.
Hosmer DW, Lemeshow S: Applied Survival Analysis: Regression Modeling of Time to Event Data.
Taylor JMG, Yu M: Bias and Efficiency Loss Due to Categorizing an Explanatory Variable, Journal of Multivariate Analysis. 2002, 83: 248-263.
Harrell FE, Lee KL, Mark DB: Multivariable prognostic models: issues in developing models, evaluating assumptions and adequacy, and measuring and reducing errors. Stat Med. 1996, 15: 361-387. 10.1002/(SICI)1097-0258(19960229)15:4<361::AID-SIM168>3.0.CO;2-4.
D'Agostino R, Nam BH: Evaluation of the performance of survival analysis models: discrimination and calibration measures. In: Handbook of Statistics. 2004, Amsterdam, The Netherlands: Elsevier, 1-25.
Antolini L, Nam BH, Agostino RB: Inference on correlated discrimination measures in survival analysis: a nonparametric approach. Commun Stat Theory Methods. 2004, 33: 2117-2135. 10.1081/STA-200026579.
Pepe MS, Janes H, Longton G, Leisenring W, Newcomb P: Limitations of the odds ratio in gauging the performance of a diagnostic, prognostic, or screening marker. American Journal of Epidemiology. 2004, 159: 882-890. 10.1093/aje/kwh101.
Greenland P, O'Malley PG: When is a new prediction marker useful? A consideration of lipoprotein-associated phospholipase A2 and C-reactive protein for stroke risk. Archives of Internal Medicine. 2005, 165 (21): 2454-2456. 10.1001/archinte.165.21.2454.
Ware JH: The limitations of risk factors as prognostic tools. New England Journal of Medicine. 2006, 355: 2615-2617. 10.1056/NEJMp068249.
Janes H, Pepe MS, Gu W: Assessing the value of risk predictions by using risk stratification tables. Annals of Internal Medicine. 2008, 149 (10): 751-760.
Vickers AJ, Cronin AM: Traditional statistical methods for evaluating prediction models are uninformative as to clinical value: Towards a decision analytic framework. Seminars in Oncology. 2010, 37 (1): 31-38. 10.1053/j.seminoncol.2009.12.004.
Pencina MJ, D'Agostino RB, D'Agostino RB, Vasan RS: Evaluating the added predictive ability of a new marker: from area under the ROC curve to reclassification and beyond. Stat Med. 2008, 27: 157-172. 10.1002/sim.2929.
Pencina MJ, D'Agostino RB: Overall C as a measure of discrimination in survival analysis: model specific population value and confidence interval estimation. Stat Med. 2004, 23: 2109-2123. 10.1002/sim.1802.
Executive summary of the Third Report of the National Cholesterol Education Program (NCEP) Expert Panel on Detection, Evaluation, and Treatment of High Blood Cholesterol in Adults (Adult Treatment Panel III). JAMA. 2001, 285: 2486-2497. 10.1001/jama.285.19.2486.
The Framingham Heart Study - Offspring (FHS-O) is conducted and supported by the NHLBI in collaboration with the FHS-O Study Investigators. This Manuscript was prepared using a limited access dataset obtained from the NHLBI and does not necessarily reflect the opinions or views of the FHS-O or the NHLBI. This research was supported by an NHMRC health services research grant (no. 465130), an NHMRC/NHF PhD scholarship and a Vichealth Fellowship.
The authors declare that they have no competing interests.
HRM was involved in all stages of this research as the principal author. CES, AP and HLW read the draft of the paper and provided useful suggestions. JJM is the Principal Investigator of the grant which enabled this research to be carried out. All authors read and approved the final manuscript.
About this article
Cite this article
Mannan, H.R., Stevenson, C.E., Peeters, A. et al. Age at quitting smoking as a predictor of risk of cardiovascular disease incidence independent of smoking status, time since quitting and pack-years. BMC Res Notes 4, 39 (2011). https://doi.org/10.1186/1756-0500-4-39