(Stroke. 2002;33:1776.)
© 2002 American Heart Association, Inc.
Original Contributions |
From the Diabetes Trials Unit (V.K., R.J.S., A.I.A., I.M.S., S.E.M., R.R.H.), Oxford Centre for Diabetes, Endocrinology and Metabolism, University of Oxford, Radcliffe Infirmary, and Department of Primary Health Care (H.A.N.), Institute of Health Sciences, Oxford, United Kingdom.
Correspondence to Dr Richard J. Stevens, Diabetes Trials Unit, Oxford Centre for Diabetes, Endocrinology and Metabolism, University of Oxford, Radcliffe Infirmary, Woodstock Rd, Oxford OX2 6HE, UK. E-mail richard.stevens@ dtu.ox.ac.uk
| Abstract |
|---|
|
|
|---|
Methods During 30 700 person-years of follow-up, 188 first strokes (52 fatal) occurred. Model fitting was carried out by maximum likelihood estimation using the Newton-Raphson method. Diagnostic plots were used to compare survival probabilities calculated by the model with those calculated using nonparametric methods.
Results Variables included in the final model were duration of diabetes, age, sex, smoking, systolic blood pressure, total cholesterol to high-density lipoprotein cholesterol ratio and presence of atrial fibrillation. Not included in the model were body mass index, hemoglobin A1c, ethnicity, and ex-smoking status. The use of the model is illustrated with a hypothetical study power calculation.
Conclusions This model forecasts the absolute risk of a first stroke in people with type 2 diabetes using variables readily available in routine clinical practice.
Key Words: blood pressure cholesterol diabetes mellitus statistics stroke
| Introduction |
|---|
|
|
|---|
Tools for calculation of cardiovascular risk are now frequently used to assist decisions about clinical management for both individuals and populations.15 Charts and computer programs have been developed to make calculations in the course of a patient-doctor visit, but typically use coronary heart disease risk as a surrogate for total cardiovascular risk.1619 The majority of the models currently in use for predicting vascular risk in individuals with diabetes use equations from the Framingham study, 12,2023 which included only 337 individuals with diabetes out of 5573.24
Data from the 5102 patients with type 2 diabetes in the UK Prospective Diabetes Study (UKPDS) provide an opportunity to derive a diabetes-specific model for risk of stroke.25,26 The model presented here is specifically designed for a type 2 diabetic population using variables available in routine clinical practice, including duration of diabetes and age at diagnosis.
| Materials and Methods |
|---|
|
|
|---|
For this analysis, data from 4549 patients of white, Afro-Caribbean, or Asian-Indian ethnic groups were included. Excluded were 39 patients of other ethnic groups; 37 with stroke before diagnosis of diabetes; 215 with missing data for blood pressure, lipids, or electrocardiography; and 262 with follow-up times too short (<4 years) for the model-fitting process. Characteristics of those included in the model are shown in Table 1.
|
Patients had blood pressure and biochemical measurements, including hemoglobin (Hb) A1c, lipid, and lipoprotein fractions recorded at diagnosis of diabetes, at the end of a 3- to 4-month dietary run-in period, and then at annual intervals.25 HbA1c was measured by high-performance liquid chromatography (Biorad Diamat Automated Glycosylated Hemoglobin Analyzer), nondiabetic range 4.5% to 6.2%.26,27 Cholesterol measurements were within the limits for accuracy of the lipid standardization program of the Centers for Disease Control and Prevention (CDC) in Atlanta, Ga, but when compared with the CDC reference laboratory in the United Kingdom were 2.5% lower and 1.2% lower (total cholesterol and HDL cholesterol, respectively).2729 Systolic blood pressure (SBP) was recorded each year as the mean of three measurements, with a coefficient of variation <15%, taken at the same visit.4 AF was identified clinically or by electrocardiography.
End Points
Stroke in the UKPDS is defined as a neurological deficit with symptoms or signs lasting 1 month or more.25 No distinction was made between ischemic, embolic, and hemorrhagic strokes. In patients with multiple strokes, only the first stroke is considered here. The model treats death from causes other than stroke as censored.
Statistical Model and Methods
The mathematical form of the model and its derivation have been described elsewhere.30 Table 2 shows the variables considered for inclusion in the model. Model building then proceeded stepwise using likelihood ratio tests for significance. The following well-documented risk factors for stroke were tested for inclusion at a significance level of 0.05: sex, age at diagnosis of diabetes, SBP, atrial fibrillation, smoking, and ratio of total:HDL cholesterol.7,10,31 The following additional possible risk factors were tested at a significance level of 0.005, to reduce the danger of a type I error (erroneous inclusion) occurring: HbA1c, Afro-Caribbean ethnicity, Asian-Indian ethnicity, body mass index (BMI), and duration of diagnosed diabetes. Proportional-hazards assumptions were verified with log-cumulative hazard plots32 (data not shown). For AF, there were too few cases for a reliable log-cumulative hazard plot.
|
To improve model stability, SBP, total cholesterol, HDL cholesterol, and HbA1c were each taken to be the mean of values taken 1 and 2 years after diagnosis of diabetes, as in our previous article.30 Model fitting was carried out by maximum likelihood estimation, using the Newton-Raphson method as implemented in the Numerical Algorithms Group C Library.33
In a supplementary analysis, we fitted alternative models in which the ratio T:H was replaced with each of the following: total cholesterol, HDL cholesterol, low-density lipoprotein (LDL) cholesterol, triglycerides, and the ratio LDL:HDL cholesterol.
During the early years of the study, all-cause mortality was not significantly different from that in the general population, possibly because people with life-threatening illnesses were excluded from the study.34 To avoid selection bias, years 0 to 4 of data were not used in model fitting; 262 patients with follow-up times <4 years were excluded as described above, and the fitting process for the remaining data used likelihoods conditional on surviving to 4 years without a stroke.
Model adequacy was checked with diagnostic plots, comparing modeled survival probabilities, calculated as the average across the 4549 patients, to survival probabilities for the study population calculated by nonparametric (life-table) methods.32
External validity of the model was examined by comparison with the rates of fatal stroke reported in 1370 patients of the Wisconsin Epidemiologic Study of Diabetic Retinopathy (WESDR) older-onset cohort, who were followed for an average of 8.3 years.35 The risk equation of this article was used to calculate 8.3-year risk of stroke in smoking and nonsmoking men and women with type 2 diabetes and characteristics as reported by WESDR.3638 Because the blood pressure and cholesterol values taken are the mean of a large number, correction for regression dilution was incorporated. These risk estimates were then combined in a weighted average, reflecting the male-to-female mix and the prevalence of smoking in the WESDR cohort, to give an estimate of the overall risk of stroke, for comparison with the rates of fatal stroke reported by WESDR. Rates of nonfatal stroke were not reported by WESDR.
| Results |
|---|
|
|
|---|
Table 3 shows the parameters included in the final model. Not significant were HbA1c (P=0.086), Afro-Caribbean (P>0.5) and Asian-Indian (P>0.5) ethnicity, and ex-smoking status (P>0.25). BMI achieved conventional significance (P=0.0085) but did not meet our criteria for inclusion in the model.
|
In a patient who has had no stroke in the first T years of diagnosed diabetes, the probability of a stroke in the next t years is RT(t) = exp[-qx1.145Tx(1- 1.145t)/(1-1.145)], where q=0.00186 x 1.092age-55x 0.700sex x 1.547smokingx8.554AFx 1.122(SBP-135.5)/10x 1.138T:H-5.11.
When the input values are of greater or lesser stability than those used to develop the modelfor example, blood pressure on a single occasion, rather than averaged over 2 occasionsthe effect of regression dilution must be considered, as indicated in the Appendix.
This model uses the ratio T:H for compatibility with our existing model software.30 However, it was found that models of similarly good fit to the data could be built using either total cholesterol (P=0.0044), LDL cholesterol (P=0.0037), or the ratio LDL:HDL cholesterol (P=0.0013) as the only measure of dyslipidemia. Nonnested models cannot be compared by likelihood ratio tests, but graphical comparisons found these models to be very similar to the T:H ratio model (data not shown). Triglycerides were significant (P=0.037) as the sole lipid measure, but not in models already containing LDL cholesterol, LDL:HDL cholesterol ratio, or T:H ratio. HDL cholesterol was not significant in any model (P>0.1).
Example
Consider a power calculation for a hypothetical trial of lipid-lowering therapy in primary prevention of stroke, among white men with type 2 diabetes. Based on data from WESDR, assume a cohort with mean SBP 147 mm Hg, total cholesterol 5.65 mmol/L, HDL cholesterol 1.11 mmol/L, age 55 years at diagnosis of diabetes, and with 12 years of diagnosed diabetes.35,36 For a nonsmoking man, with characteristics as above and assuming no AF, calculations are as follows: q= 0.00186x1.09255-55x 1x1x1x1.122(147- 135.5)/10x 1.138(5.65/1.11)- 5.11=0.00212, and hence R12(5)= 1-exp[-0.00212x1.14512 x(1-1.1455)/(1- 1.145)]=0.069, so that such a patient has a 6.9% probability of a stroke within the next 5 years, conditional on not dying from causes other than stroke. For a patient who smokes, the probability is 10.5%. Assuming a 30% smoking prevalence, the average 5-year risk in the population is (0.3x10.5%)+(0.7x6.9%)=8.0%. To have 80% power to detect an effect of a therapy that lowers risk by 15%, a 5-year trial would need
15 000 patients. These event rates are likely to be overestimates for the whole population, because not all patients will remain in the study for 5 years, as a result of deaths from cardiac and other causes as well as loss to follow-up. A more thorough analysis would correct also for regression dilution and account for the variation in risk factors across the population rather than use mean values.
Validation
Figure 1 shows that survival rates predicted by the model lie close to the rates observed in the UKPDS. Observed and modeled survival rates are shown in Figure 2a and 2b, split by SBP and by the ratio T:H, confirming the ability of the model to adjust for these risk factors.
|
|
For the WESDR cohort, in which 59 deaths from stroke were observed, the equation derived here predicts a 13.2% risk of fatal and nonfatal stroke in 8.3 years, corresponding to 197 strokes in 1370 patients. Previous studies have reported stroke fatality rates between 23% and 30% in diabetes.7,39,40 Applying a 26.5% fatality rate to our prediction of 197 strokes gives 52 fatal strokes. For comparison, the Framingham equations for mortality predict 33 fatal strokes in the same cohort.41
| Discussion |
|---|
|
|
|---|
Several limitations arise from the use of clinical trial rather than epidemiologic data, although we have attempted to correct for the impact of trial selection criteria, as described above. For time periods <4 years or >20 years after diagnosis of diabetes, or for ages >65 at diagnosis of diabetes, predictions from this model are extrapolations. Our definition of nonfatal stroke is more stringent than is sometimes used, requiring signs or symptoms to persist for more than a month. Our model does not distinguish between fatal and nonfatal strokes or between ischemic and hemorrhagic strokes. It would be ideal to test for any effect of antihypertensive and other therapies, but this is not practical because many patients took many different therapies over the years of UKPDS follow-up.4,43 Certain limitations are generic to all modeling exercises. Users should be aware that the risk ratio between individuals at 2 levels of a risk factor is not necessarily equal to the risk reduction obtainable for a single individual by changing levels of the risk factor. For tight blood pressure control, the risk reduction obtained in the UKPDS was greater than the risk ratio in an observational model would suggest.3 The current version of the model does not provide for estimation when data are missing. HbA1c was not found to be significant in this analysis, as it was in a previous UKPDS analysis that used more person-years of data and a more informative variable, updated mean HbA1c.6 BMI achieved conventional significance (P<0.05) but did not meet our inclusion criteria, which were designed to protect the generalizability of the risk equation. A risk equation including BMI might not generalize well, because previous studies have found that BMI is not a risk factor for stroke in type 2 diabetes or have even reported a possible protective effect of obesity.44 47
The comparison with stroke death rates in WESDR suggests that the model may generalize well to other populations. Although our derived figure, 52 fatal strokes, is lower than the 59 fatal strokes, this is to be expected because we model first stroke only. There are few other published models for the risk of stroke, and none specific to diabetes. The best-known previously published risk equations for stroke are derived from the Framingham cohort,20,41 and have been used in guidelines for clinical practice48 and in health economic modeling in the general population49 and in diabetes.2123 The 5734 individuals in Framingham were older than in the UKPDS, (mean age 66 years) and experienced a higher event rate, 472 strokes in 10 years, but fewer than 10% had diabetes. The Framingham profile adjusts for left ventricular hypertrophy. Exclusion of the latter from the present paper gives our equations greater applicability, at the cost of some discriminatory power. The British Regional Heart Study model for risk of stroke does not adjust for presence of diabetes,50 and a recent cardiovascular risk score provides only relative risk in the case of stroke.51
Although guidelines for primary prevention often use risk of coronary heart disease as a surrogate for risk of cardiovascular disease,48,52 accurate estimation of stroke as well as coronary risk is preferable, for example, in hypertension management17 and for patients of African origin.15 The model here is suitable for use in a clinical setting.
Risk estimators also have a role to play in research and in policy formulation. The example in the Results section has demonstrated the use of the model in a power calculation, indicating that in primary prevention the large numbers needed to establish a lipid-lowering effect on stroke may require large-scale meta-analyses.53,54 The literature applying the Framingham equation to cost-effectiveness analyses in diabetes demonstrates a role for a diabetes-specific risk estimator.22,23,55 We hope that the present article will be of use to researchers, planners, and clinicians in the management of cerebrovascular risk in diabetes.
| Acknowledgments |
|---|
The cooperation of the patients and staff of the UKPDS is much appreciated. The participating centers in the UKPDS were Radcliffe Infirmary, Oxford; Royal Infirmary, Aberdeen; University Hospital, Birmingham; St. Georges Hospital and Hammersmith Hospital, London; City Hospital, Belfast; North Staffordshire Royal Infirmary, Stoke-on-Trent; Royal Victoria Hospital, Belfast; St. Helier Hospital, Carshalton; Whittington Hospital, London; Norfolk & Norwich Hospital; Lister Hospital, Stevenage; Ipswich Hospital; Ninewells Hospital, Dundee; Northampton Hospital; Torbay Hospital; Peterborough General Hospital; Scarborough Hospital; Derbyshire Royal Infirmary; Manchester Royal Infirmary; Hope Hospital, Salford; Leicester General Hospital; Royal Devon & Exeter Hospital. We are grateful to Prof Peter Bennet and to Dr Andrew Farmer for helpful comments.
| Appendix 1 |
|---|
|
|
|---|
Adjustment for regression dilution is appropriate if the available data are not subject to biovariability, for example, if the blood pressure and cholesterol levels are the means over very large populations. Based on a maximum likelihood method,57 the following corrections are appropriate. For blood pressure, replace the parameter estimate 1.122 with 1.1221.25=1.155. For lipid ratio, replace the parameter estimate 1.138 with 1.1381.28=1.181. Because variables are centered around their means in the model equation (for example, blood pressure is used as SBP-135.5), no adjustment is necessary to the intercept parameter.
We use a method reported previously30 to estimate the adjustment necessary if the input data are of greater biovariability than ours. If the available blood pressure is the mean of just 3 measurements, replace 1.122 with 1.101. If blood pressure is a single measurement, replace 1.122 with 1.060. If the lipid ratio is based on single measurements of HDL and total cholesterol, replace the parameter 1.138 with 1.111.
Received November 19, 2001; revision received February 18, 2002; accepted April 3, 2002.
| References |
|---|
|
|
|---|
This article has been cited by other articles:
![]() |
P. Gaede, W. J. Valentine, A. J. Palmer, D. M.D. Tucker, M. Lammert, H.-H. Parving, and O. Pedersen Cost-Effectiveness of Intensified Versus Conventional Multifactorial Intervention in Type 2 Diabetes: Results and projections from the Steno-2 study Diabetes Care, August 1, 2008; 31(8): 1510 - 1515. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. S. Huang, Q. Zhang, N. Gandra, M. H. Chin, and D. O. Meltzer The Effect of Comorbid Illness and Functional Status on the Expected Benefits of Intensive Glucose Control in Older Patients with Type 2 Diabetes: A Decision Analysis Ann Intern Med, July 1, 2008; 149(1): 11 - 19. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. J. Hayes, P. M. Clarke, P. G. Glasziou, R. J. Simes, P. L. Drury, and A. C. Keech Can Self-Rated Health Scores Be Used for Risk Prediction in Patients With Type 2 Diabetes? Diabetes Care, April 1, 2008; 31(4): 795 - 797. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. L. Air and B. M. Kissela Diabetes, the Metabolic Syndrome, and Ischemic Stroke: Epidemiology and possible mechanisms Diabetes Care, December 1, 2007; 30(12): 3131 - 3140. [Full Text] [PDF] |
||||
![]() |
C. A. Cull, C. C. Jensen, R. Retnakaran, and R. R. Holman Impact of the Metabolic Syndrome on Macrovascular and Microvascular Outcomes in Type 2 Diabetes Mellitus: United Kingdom Prospective Diabetes Study 78 Circulation, November 6, 2007; 116(19): 2119 - 2126. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. J. Hoerger, K. A. Hicks, S. W. Sorensen, W. H. Herman, R. E. Ratner, R. T. Ackermann, P. Zhang, and M. M. Engelgau Cost-Effectiveness of Screening for Pre-Diabetes Among Overweight and Obese U.S. Adults Diabetes Care, November 1, 2007; 30(11): 2874 - 2879. [Abstract] [Full Text] [PDF] |
||||
![]() |
The Mount Hood 4 Modeling Group Computer Modeling of Diabetes and Its Complications: A report on the Fourth Mount Hood Challenge Meeting Diabetes Care, June 1, 2007; 30(6): 1638 - 1646. [Abstract] [Full Text] [PDF] |
||||
![]() |
C. B. Giorda, A. Avogaro, M. Maggini, F. Lombardo, E. Mannucci, S. Turco, S. S. Alegiani, R. Raschetti, M. Velussi, E. Ferrannini, et al. Incidence and Risk Factors for Stroke in Type 2 Diabetic Patients: The DAI Study Stroke, April 1, 2007; 38(4): 1154 - 1160. [Abstract] [Full Text] [PDF] |
||||
![]() |
T. Wada, M. Naito, H. Kenmochi, H. Tsuneki, and T. Sasaoka Chronic Nicotine Exposure Enhances Insulin-Induced Mitogenic Signaling via Up-Regulation of {alpha}7 Nicotinic Receptors in Isolated Rat Aortic Smooth Muscle Cells Endocrinology, February 1, 2007; 148(2): 790 - 799. [Abstract] [Full Text] [PDF] |
||||
![]() |
X. Yang, W.-Y. So, A. P.S. Kong, C.-S. Ho, C. W.K. Lam, R. J. Stevens, R. R. Lyu, D. D. Yin, C. S. Cockram, P. C.Y. Tong, et al. Development and Validation of Stroke Risk Equation for Hong Kong Chinese Patients With Type 2 Diabetes: The Hong Kong Diabetes Registry Diabetes Care, January 1, 2007; 30(1): 65 - 70. [Abstract] [Full Text] [PDF] |
||||
![]() |
A. Balamurugan, R. Ohsfeldt, T. Hughes, and M. Phillips Diabetes Self-management Education Program for Medicaid Recipients: A Continuous Quality Improvement Process The Diabetes Educator, November 1, 2006; 32(6): 893 - 900. [Abstract] [Full Text] [PDF] |
||||
![]() |
H. Zhou, D. J.M. Isaman, S. Messinger, M. B. Brown, R. Klein, M. Brandle, and W. H. Herman A Computer Simulation Model of Diabetes Progression, Quality of Life, and Cost Diabetes Care, December 1, 2005; 28(12): 2856 - 2863. [Abstract] [Full Text] [PDF] |
||||
![]() |
E. Selvin, J. Coresh, S. H. Golden, F. L. Brancati, A. R. Folsom, and M. W. Steffes Glycemic Control and Coronary Heart Disease Risk in Persons With and Without Diabetes: The Atherosclerosis Risk in Communities Study Arch Intern Med, September 12, 2005; 165(16): 1910 - 1916. [Abstract] [Full Text] [PDF] |
||||
![]() |
R. M. Clifford, W. A. Davis, K. T. Batty, and T. M.E. Davis Effect of a Pharmaceutical Care Program on Vascular Risk Factors in Type 2 Diabetes: The Fremantle Diabetes Study Diabetes Care, April 1, 2005; 28(4): 771 - 776. [Abstract] [Full Text] [PDF] |
||||
![]() |
W. H. Herman, T. J. Hoerger, M. Brandle, K. Hicks, S. Sorensen, P. Zhang, R. F. Hamman, R. T. Ackermann, M. M. Engelgau, R. E. Ratner, et al. The Cost-Effectiveness of Lifestyle Modification or Metformin in Preventing Type 2 Diabetes in Adults with Impaired Glucose Tolerance Ann Intern Med, March 1, 2005; 142(5): 323 - 332. [Abstract] [Full Text] [PDF] |
||||
![]() |
L. Rodriguez-Manas, M. Castro, and J. Manzarbeitia Review: Cardiovascular risk and diabetes in old people The British Journal of Diabetes & Vascular Disease, January 1, 2005; 5(1): 7 - 11. [Abstract] [PDF] |
||||
![]() |
E. Selvin, S. Marinopoulos, G. Berkenblit, T. Rami, F. L. Brancati, N. R. Powe, and S. H. Golden Meta-Analysis: Glycosylated Hemoglobin and Cardiovascular Disease in Diabetes Mellitus Ann Intern Med, September 21, 2004; 141(6): 421 - 431. [Abstract] [Full Text] [PDF] |
||||
![]() |
American Diabetes Association Consensus Panel Guidelines for Computer Modeling of Diabetes and Its Complications Diabetes Care, September 1, 2004; 27(9): 2262 - 2265. [Full Text] [PDF] |
||||
![]() |
R. J. Stevens, R. L. Coleman, A. I. Adler, I. M. Stratton, D. R. Matthews, and R. R. Holman Risk Factors for Myocardial Infarction Case Fatality and Stroke Case Fatality in Type 2 Diabetes: UKPDS 66 Diabetes Care, January 1,&nb |