| Literature DB >> 22470551 |
Nicolas Rodondi1, Isabella Locatelli, Drahomir Aujesky, Javed Butler, Eric Vittinghoff, Eleanor Simonsick, Suzanne Satterfield, Anne B Newman, Peter W F Wilson, Mark J Pletcher, Douglas C Bauer.
Abstract
BACKGROUND: Guidelines for the prevention of coronary heart disease (CHD) recommend use of Framingham-based risk scores that were developed in white middle-aged populations. It remains unclear whether and how CHD risk prediction might be improved among older adults. We aimed to compare the prognostic performance of the Framingham risk score (FRS), directly and after recalibration, with refit functions derived from the present cohort, as well as to assess the utility of adding other routinely available risk parameters to FRS.Entities:
Mesh:
Substances:
Year: 2012 PMID: 22470551 PMCID: PMC3314613 DOI: 10.1371/journal.pone.0034287
Source DB: PubMed Journal: PLoS One ISSN: 1932-6203 Impact factor: 3.240
Baseline characteristics and unadjusted associations with incident CHD events (n = 2193; number of CHD events = 351).
| Variable | Mean ± SD / n (%) | HR (95% CI) | p |
| Age | 73.50±2.85 | 1.03 (1.00, 1.07) | 0.09 |
| Age (categories) | 0.03 | ||
| 70–71 | 672 (30.6) | ||
| 72–75 | 934 (42.6) | 0.99 (0.77, 1.27) | |
| 76–78 | 464 (21.2) | 1.05 (0.78, 1.41) | |
| 79 | 123 (5.6) | 1.59 (1.04, 2.43) | |
| Gender | |||
| Men | 981 (44.7) | ||
| Women | 1212 (55.3) | 0.52 (0.42, 0.65) | <0.001 |
| Race | |||
| White | 1293 (59.0) | ||
| Black | 900 (41.0) | 0.96 (0.78, 1.19) | 0.73 |
| Site | |||
| Memphis | 1125 (51.3) | ||
| Pittsburgh | 1068 (48.7) | 0.99 (0.79, 1.22) | 0.89 |
| Education | 0.29 | ||
| <high school | 532 (24.3) | ||
| High school graduate | 734 (33.6) | 0.81 (0.61, 1.07) | |
| Postsecondary | 922 (42.1) | 0.87 (0.67, 1.13) | |
| Smoking status | 0.03 | ||
| Never | 1016 (46.3) | ||
| Former | 956 (43.6) | 1.41 (1.13, 1.77) | |
| Current | 221 (10.1) | 1.49 (1.04, 2.12) | |
| Alcohol, drinks/wk | 0.41 | ||
| <1 | 1535 (70.3) | ||
| 1–7 | 482 (22.1) | 0.88 (0.68, 1.14) | |
| >7 | 166 (7.6) | 1.17 (0.81, 1.70) | |
| Physical activity, kcal/wk | 0.20 | ||
| <500 | 1148 (52.3) | ||
| 500–1500 | 598 (27.3) | 0.96 (0.74, 1.23) | |
| ≥1500 | 447 (20.4) | 1.18 (0.91, 1.54) | |
| Hypertension | 1258 (57.4) | 1.28 (1.03, 1.59) | 0.02 |
| Diabetes mellitus | 292 (13.3) | 1.63 (1.24, 2.13) | <0.001 |
| Body mass index, kg/m2 | 27.41±4.91 | 1.02 (1.00, 1.04) | 0.10 |
| Abdominal circumference | 99.43±13.54 | 1.01 (1.00, 1.02) | 0.02 |
| Systolic blood pressure, per 10 mmHg | 135.72±20.63 | 1.09 (1.04, 1.15) | <0.001 |
| Diastolic blood pressure, per 10 mmHg | 71.59±11.66 | 1.12 (1.02, 1.22) | 0.01 |
| Total cholesterol, per 10 mg/dl | 204.83±37.93 | 0.99 (0.96, 1.01) | 0.31 |
| HDL-cholesterol, per 10 mg/dl | 55.46±17.12 | 0.87 (0.81, 0.93) | <0.001 |
| Total/HDL-cholesterol | 3.98±1.22 | 1.15 (1.06, 1.25) | <0.001 |
| LDL-cholesterol, per 10 mg/dl | 122.87±34.44 | 1.01 (0.98, 1.04) | 0.50 |
| Triglycerides, mg/dl | 116 (87–160) | 1.11 (0.89, 1.39) | 0.36 |
| Glucose, per 10 mg/dl | 102.46±31.88 | 1.05 (1.03, 1.08) | <0.001 |
| Framingham risk score, % | <0.001 | ||
| <5% | 468 (21.3) | ||
| 5–9.99% | 557 (25.4) | 1.35 (0.92, 1.99) | |
| 10–19.99% | 543 (24.8) | 2.12 (1.47, 3.04) | |
| ≥20% | 625 (28.5) | 3.06 (2.17, 4.31) | |
| Creatinine, mg/dl | 1 (0.9–1.1) | 1.96 (1.33, 2.87) | 0.001 |
| GFR, ml/min/1.73 m2
| 61.15±15.05 | 1.00 (0.99, 1.01) | 0.87 |
| GFR (categories) | 0.63 | ||
| ≥80 | 525 (23.9) | 1.00 | |
| 70–79.99 | 536 (24.5) | 1.11 (0.83,1.49) | |
| 60–69.99 | 555 (25.3) | 0.87 (0.64,1.18) | |
| <60 | 576 (26.3) | 1.00 (0.75,1.35) | |
| Medication use | |||
| Lipid-lowering | 229 (10.4) | 1.04 (0.75, 1.46) | 0.79 |
| Ace inhibitors | 273 (12.4) | 1.14 (0.84, 1.54) | 0.39 |
| Hormone replacement therapy | 48 (2.2) | 0.82 (0.37, 1.85) | 0.64 |
| Aspirin | 412 (18.8) | 1.31 (1.02, 1.69) | 0.03 |
Abbreviations: SD: standard deviation; HR: hazard ratio; CI: confidence interval; HDL: high-density lipoprotein; LDL: low-density lipoprotein; GFR: glomerular filtration rate.
*p for trend.
Physical activity was assessed by questionnaire about all types of walking and exercise performed in the prior week [14].
Defined by self-report of hypertension and use of anti-hypertensive medications, or measured SBP≥140 and/or DBP≥90 mmHg.
Expressed as median (25%–75%), because of skewed distribution. The effect of the logarithm of the covariates on the CHD is measured.
Classes of CHD risk at 10 years, according to Framingham functions [2].
Glomerular filtration rate (GFR) was estimated using the MDRD equation: GFR = 175 * Creatinine−1.154 * Age−0.203 * (1.212*Iblack+Iwhite)* (0.742*Ifemale+Imen) [32].
**Quartiles were used instead of clinical cut-offs to avoid categories with few participants. In particular, categories of GFR<15 and within 15–29.99 were collapsed with the category 30–59.99 (only 0.2% in the class of GFR<15 and 0.5% in the class of GFR within 15–29.99) and a ≥80 category was replaced to the usual ≥90 (only 7% for GFR≥90).
Discrimination and calibration of Framingham functions (FRS), recalibrated FRS and Health ABC function in women (n = 1212).
| Participants with Risk Factor, N (%) | CHD events (N) | FRS | Recalibrated FRS | Refit FRS (Health ABC function ) | ||
| Coef | Coef | Coef (95%CI) | HR (95% CI) | |||
| Age, y, mean (SD) | 73.41 (2.84) | 146 | 0.33766 | 0.33766 | 0.00 (−0.06,0.06) | 1.00 (0.95,1.06) |
| Age2 | −0.00268 | −0.00268 | ||||
| TC, mg/dL | ||||||
| <160 | 78 (6%) | 11 | −0.26138 | −0.26138 | ||
| 160–199 | 364 (30%) | 46 | Referent | Referent | ||
| 200–239 | 496 (41%) | 52 | 0.20771 | 0.20771 | ||
| 240–279 | 217 (18%) | 29 | 0.24385 | 0.24385 | ||
| ≥280 | 57 (5%) | 8 | 0.53513 | 0.53513 | ||
| HDL-C, mg/dL | ||||||
| <35 | 21 (2%) | 2 | 0.84312 | 0.84312 | 0.21 (−0.30,0.71) | 1.23 (0.74 ,2.04) |
| 35–44 | 149 (12%) | 22 | 0.37796 | 0.37796 | ||
| 45–49 | 149 (12%) | 21 | 0.19785 | 0.19785 | 0.14 (−0.39,0.67) | 1.15 (0.68,1.95) |
| 50–59 | 322 (27%) | 41 | Referent | Referent | Referent | Referent |
| ≥60 | 571 (47%) | 60 | −0.42951 | −0.42951 | −0.10 (−0.51,0.30) | 0.90 (0.60,1.35) |
| Blood pressure | ||||||
| Optimal | 266 (22%) | 19 | −0.53363 | −0.53363 | −0.26 (−0.86,0.33) | 0.77 (0.42,1.40) |
| Normal | 259 (21%) | 25 | Referent | Referent | Referent | Referent |
| High normal | 254 (21%) | 38 | −0.06773 | −0.06773 | 0.41 (−0.09,0.92) | 1.51 (0.91,2.51) |
| Stage I hypertension | 296 (24) | 46 | 0.26288 | 0.26288 | 0.45 (−0.02,0.91) | 1.56 (0.98,2.49) |
| Stage II–IV hypertension | 137 (11%) | 18 | 0.46573 | 0.46573 | ||
| Diabetes | 141 (12%) | 27 | 0.59626 | 0.59626 | 0.62 (0.20,1.05) | 1.86 (1.22,2.85) |
| Smoker | ||||||
| Never | 714 (59%) | 83 | Referent | Referent | Referent | Referent |
| Former | 388 (32%) | 48 | ||||
| Current | 110 (9%) | 15 | 0.29246 | 0.29246 | 0.29 (−0.25,0.83) | 1.34 (0.78,2.29) |
| Mean survival function at t = 7.5 years, S0(t) | 0.9717 | 0.8898 | 0.8962 | |||
| C-index | 0.577 | 0.577 | 0.598 | |||
| H-L statistics | 121.43 (<0.001) | 22.73 (0.007) | 7.96 (0.539) | |||
Abbreviations: FRS: Framingham risk score; CHD: coronary heart disease; coef: coefficient; CI: confidence interval; SD: standard deviation; TC: total cholesterol; HDL-C: high-density lipoprotein cholesterol.
*Some of the Framingham risk factors categories were collapsed to avoid cells with limited numbers of events and /or unpredictive trends. Total cholesterol and age2 were omitted because they were unpredictive in these older women. The proportionality assumption was tested using the Therneau and Grambsch statistics, which is based on the Schoenfeld residuals. The assumption was accepted (p = 0.14).
Based on Wilson et al. [2].
Cholesterol categories proposed by the National Cholesterol Education Program [24].
Blood pressure categories: Optimal (Systolic<120, Diastolic>80); Normal (Systolic<130, Diastolic>85); High normal (Systolic<140, Diastolic>90); Stage I (Systolic<160, Diastolic<100); Stage II–IV (Systolic ≥160, Diastolic ≥100) [24].
Estimated from the Framingham adjusted survival rate (survival rate at the mean value of the risk factors) at 10 years: S0(10) = 0.96246 [2], as: Ŝ0(7.5) = S0(10)0.75 = 0.9717 (exponential model).
Kaplan-Meier survival function at t = 7.5 years on HABC data, similar to reference [24].
**Adjusted survival rate at t = 7.5 years obtained on the HABC cohort as the baseline survival functions of the multivariate Cox model, similar to reference [9].
After bootstrap correction for the optimism (1000 bootstrap samples from the original dataset [20]), c-index = 0.564 (p = 0.54 for comparison with Framingham function).
Adaptation to the Cox model of the Hosmer-Lemeshow test of goodness of fit [33], comparing observed and expected failures within deciles of predicted risk. Larger p values indicate better calibration [21].
Discrimination and calibration of Framingham functions (FRS), recalibrated FRS and Health ABC function in men (n = 981).
| Participants with Risk Factor, N (%) | CHD events (N) | FRS | Recalibrated FRS | Refit FRS (Health ABC function ) | ||
| Coef | Coef | Coef (95% CI) | HR (95% CI) | |||
| Age, y, mean (SD) | 73.613 (2.86) | 205 | 0.04826 | 0.04826 | 0.05 (−0.00,0.10) | 1.05 (1.00,1.10) |
| TC, mg/dL | ||||||
| <160 | 139 (14%) | 23 | −0.65945 | −0.65945 | −0.32 (−0.78,0.14) | 0.73 (0.46,1.15) |
| 160–199 | 451 (46%) | 94 | Referent | Referent | Referent | Referent |
| 200–239 | 303 (31%) | 69 | 0.17692 | 0.17692 | ||
| 240–279 | 70 (7%) | 15 | 0.50539 | 0.50539 | 0.10 (−0.20,0.39) | 1.10 (0.82,1.48) |
| ≥280 | 18 (2%) | 4 | 0.65713 | 0.65713 | ||
| HDL-C, mg/dL | ||||||
| <35 | 140 (14%) | 28 | 0.49744 | 0.49744 | ||
| 35–44 | 295 (30%) | 72 | 0.24310 | 0.24310 | Referent | Referent |
| 45–49 | 160 (16%) | 38 | Referent | Referent | ||
| 50–59 | 204 (21%) | 40 | −0.05107 | −0.05107 | −0.23 (−0.58,0.13) | 0.80 (0.56,1.13) |
| ≥60 | 182 (19%) | 27 | −0.48660 | −0.48660 | −0.60 (−1.02 ,−0.19) | 0.55 (0.36,0.83) |
| Blood pressure | ||||||
| Optimal | 214 (22%) | 27 | −0.00226 | −0.00226 | −0.47 (−0.95,0.02) | 0.63 (0.39,1.02) |
| Normal | 210 (22%) | 42 | Referent | Referent | Referent | Referent |
| High normal | 188 (19%) | 51 | 0.28320 | 0.28320 | ||
| Stage I hypertension | 258 (26%) | 59 | 0.52168 | 0.52168 | 0.18 (−0.16,0.53) | 1.20 (0.85,1.70) |
| Stage II–IV hypertension | 111 (11%) | 26 | 0.61859 | 0.61859 | ||
| Diabetes | 151 (15%) | 38 | 0.42839 | 0.42839 | 0.23 (−0.12,0.58) | 1.26 (0.88,1.79) |
| Smoker | ||||||
| Never | 302 (31%) | 56 | Referent | Referent | Referent | Referent |
| Former | 568 (58%) | 125 | ||||
| Current | 111 (11%) | 24 | 0.52337 | 0.52337 | 0.28 (−0.15,0.71) | 1.32 (0.86,2.03) |
| Mean survival function at t = 7.5 years, S0(t) | 0.9241 | 0.7929 | 0.8032 | |||
| C-index | 0.583 | 0.583 | 0.606 | |||
| H-L statistics | 16.27 (0.062) | 16.11 (0.065) | 4.89 (0.844) | |||
Abbreviations: FRS: Framingham risk score; CHD: coronary heart disease; coef: coefficient; CI: confidence interval; SD: standard deviation; TC: total cholesterol; HDL-C: high-density lipoprotein cholesterol.
*Some of the Framingham risk factors categories were collapsed to avoid cells with limited numbers of events and /or unpredictive trends. The proportionality assumption was tested using the Therneau and Grambsch statistics, which is based on the Schoenfeld residuals. The assumption was accepted (p = 0.33).
Based on Wilson et al. [2].
Cholesterol categories proposed by the National Cholesterol Education Program [24].
Blood pressure categories: Optimal (Systolic<120, Diastolic>80); Normal (Systolic<130, Diastolic>85); High normal (Systolic<140, Diastolic>90); Stage I (Systolic<160, Diastolic<100); Stage II–IV (Systolic ≥160, Diastolic ≥100) [24].
Estimated from the Framingham adjusted survival rate (survival rate at the mean value of the risk factors) at 10 years: S0(10) = 90015 [2], as: Ŝ0(7.5) = S0(10)0.75 = 0.9241 (exponential model).
Kaplan-Meier survival function at t = 7.5 years on HABC data, similar to reference [24].
**Adjusted survival rate at t = 7.5 years obtained on the HABC cohort as the baseline survival functions of the multivariate Cox model, similar to reference [9].
After bootstrap correction for the optimism (1000 bootstrap samples from the original dataset [20]), c-index = 0.580 (p = 0.90 for comparison with Framingham function).
Adaptation to the Cox model of the Hosmer-Lemeshow test of goodness of fit [33], comparing observed and expected failures within deciles of predicted risk. Larger p values indicate better calibration [21].
Figure 1Predicted risk of CHD events at 7.5 years according to original Framingham functions, recalibrated Framingham functions and Health ABC functions.