Impact of HbA1c testing on cardiovascular morbidity and mortality: A Regression Discontinuity Design study
PROJECT GROUP Irene Petersen (Department of Primary Care & Population Health, University College London and Dept. of Clinical Epidemiology, Aarhus University, Denmark) Federico Ricciardi, Gianluca Baio (Department of Statistical Science, University College London) Sia Kromann Nicolaisen, Reimar W Thomsen and Lars Pedersen (Department of Clinical Epidemiology, Aarhus University, Denmark)
ABOUT Young, modern university from 1928 44,500 students, 10% international students A graduate university. More than 50% of students study at master or PhD level
REGRESSION DISCONTINUITY DESIGN (RDD) In epidemiology, there is a great need for designs that can support generation of causal evidence To date, RDD has not been implemented much in epidemiological research using real world data Therefore, we aimed to explore the potential for RDD to generate causal evidence using existing healthcare data in Denmark
AIM In particular, our aims were: To review the necessary RDD design assumptions in relation to a real study on HbA1c testing and cardiovascular morbidity and mortality in Denmark. To run simple intention-to-treat (ITT) analyses and explore the robustness of the effect estimates for different bandwidths. If possible, incorporate treatment in the analyses and derive causal complier estimates (CCHR), and compare the findings against the existing literature.
HBA1C The term HbA1c refers to glycated haemoglobin HbA1c reflects the average blood glucose level over the last 8-12 weeks In Denmark, HbA1c is used to diagnose patients with type2 diabetes
DISTRIBUTION OF HBA1C MEASUREMENTS IN DENMARK AMONG PEOPLE WITH NO PREVIOUS DIABETES DIAGNOSIS 12.5 10.0 Percent 7.5 5.0 2.5 0.0 20 25 30 35 40 45 50 55 60PROFESSOR, 65 MSC, 70 PHD 75 HbA1c, mmol/mol
EPIDEMIOLOGICAL QUESTION We want to evaluate the effect (ITT) of treatment eligibility (HbA1c 48 mmol/mol) on allcause mortality and cardiovascular events Subsequently, we want to estimate the effect (CCHR) of early metformin initiation versus deferred or no metformin initiation Note: According to Danish guidelines, metformin treatment initiation is recommended for patients newly diagnosed with diabetes, i.e., patients should be treated with metformin if their HbA1c level exceeds 48 mmol/mol
THE RDD TIME TO EVENT MODEL Y = Mortality and cardiovascular events Z = Value of the first HbA1c measurement (LABKA) c = 48 mmol/mol Threshold indicator (0 if Z < c, 1 if Z c) T = Treatment indicator (T = 0 if the patient does not receive metformin, T = 1 otherwise) Z c = Centered assignment variable
DATA AND POPULATION Data from North and Central Denmark Regions Total population ~1.8M Complete coverage Free, tax-supported healthcare for everyone in Denmark Record-linkage at an individual level
LABKA The clinical laboratory information system (LABKA) Holds test results from every blood sample analyzed in clinical chemistry department laboratories in the North and Central Denmark regions (pop of ~1.8M) Includes results from any patient visiting a hospital or a general practitioner Few exceptions are some results from small, rapid point-of-care devices used by medical staff or patients themselves Data available since 1992 (only selected departments). Since late 1990s, full geographical coverage Data are updated annually
STUDY POPULATION We included individuals aged 40-80 years at the time of their first HbA1c measurement between 2006 and 2014 (N = 525,266) We excluded individuals with a previous diagnosis of type I and type II diabetes, individuals treated with metformin or other glucose lowering drugs prior to the initial HbA1c measurement, and individuals with a cardiovascular event before the initial HbA1c measurement After exclusions, we had 290,333 patients
HISTOGRAM OF THE ASSIGNMENT VARIABLE 12.5 10.0 Percent 7.5 5.0 2.5 0.0 20 25 30 35 40 45 50 55 60PROFESSOR, 65 MSC, 70 PHD 75 HbA1c, mmol/mol
PROBABILITY OF METFORMIN INITIATION 1.0 Cumulative incidence 0.8 0.6 0.4 0.2 0.0 25 30 35 40 45 50 55 60 65 70 75 80 85 90 HbA1c, mmol/mol 1 month 3 months
BASELINE CHARACTERISTICS Range 1 Range 2 Range 3 40-47 48-56 Std diff 42-47 48-53 Std diff 44-47 48-51 Std diff Total 76,031 8,011 36,360 6,710 14,216 5,416 Male 47.6% 52.5% -0.10 48% 52% -0.08 48.4% 51.6% -0.06 Median age 62.0 62.9-0.21 62.9 63.0-0.02 63.4 63.3 0.07 Heart failure 861 (1.1%) 202 (2.5%) -0.10 531 (1.5%) 179 (2.7%) -0.08 253 (1.8%) 134 (2.5%) -0.05 Peripheral vascular disease 1,418 (1.9%) 197 (2.5%) -0.04 787 (2.2%) 173 (2.6%) -0.03 338 (2.4%) 145 (2.7%) -0.02 Obesity 2,274 (3.0%) 405 (5.1%) -0.11 1,257 (3.5%) 337 (5.0%) -0.08 558 (3.9%) 257 (4.7%) -0.04 Cancer 2,211 (2.9%) 350 (4.4%) -0.08 1,201 (3.3%) 285 (4.2%) -0.05 513 (3.6%) 214 (4.0%) -0.02
BASELINE CHARACTERISTICS Range 1 Range 2 Range 3 40-47 48-56 Std diff 42-47 48-53 Std diff 44-47 48-51 Std diff Total 76,031 8,011 36,360 6,710 14,216 5,416 Male 47.6% 52.5% -0.10 48% 52% -0.08 48.4% 51.6% -0.06 Median age 62.0 62.9-0.21 62.9 63.0-0.02 63.4 63.3 0.07 Heart failure 861 (1.1%) 202 (2.5%) -0.10 531 (1.5%) 179 (2.7%) -0.08 253 (1.8%) 134 (2.5%) -0.05 Peripheral vascular disease 1,418 (1.9%) 197 (2.5%) -0.04 787 (2.2%) 173 (2.6%) -0.03 338 (2.4%) 145 (2.7%) -0.02 Obesity 2,274 (3.0%) 405 (5.1%) -0.11 1,257 (3.5%) 337 (5.0%) -0.08 558 (3.9%) 257 (4.7%) -0.04 Cancer 2,211 (2.9%) 350 (4.4%) -0.08 1,201 (3.3%) 285 (4.2%) -0.05 513 (3.6%) 214 (4.0%) -0.02
BASELINE CHARACTERISTICS Range 1 Range 2 Range 3 40-47 48-56 Std diff 42-47 48-53 Std diff 44-47 48-51 Std diff Total 76,031 8,011 36,360 6,710 14,216 5,416 Male 47.6% 52.5% -0.10 48% 52% -0.08 48.4% 51.6% -0.06 Median age 62.0 62.9-0.21 62.9 63.0-0.02 63.4 63.3 0.07 Heart failure 861 (1.1%) 202 (2.5%) -0.10 531 (1.5%) 179 (2.7%) -0.08 253 (1.8%) 134 (2.5%) -0.05 Peripheral vascular disease 1,418 (1.9%) 197 (2.5%) -0.04 787 (2.2%) 173 (2.6%) -0.03 338 (2.4%) 145 (2.7%) -0.02 Obesity 2,274 (3.0%) 405 (5.1%) -0.11 1,257 (3.5%) 337 (5.0%) -0.08 558 (3.9%) 257 (4.7%) -0.04 Cancer 2,211 (2.9%) 350 (4.4%) -0.08 1,201 (3.3%) 285 (4.2%) -0.05 513 (3.6%) 214 (4.0%) -0.02
BASELINE CHARACTERISTICS (MEDICATION) Range 1 Range 2 Range 3 40-47 48-56 Std diff 42-47 48-53 Std diff 44-47 48-51 Std diff Total 76,031 8,011 36,360 6,710 14,216 5,416 Statins 17,039 (22.4%) 1,969 (24.6%) -0.05 8,924 (24.5%) 1,683 (25.1%) -0.01 3,647 (25.7%) 1,389 (25.6%) 0.00 NSAIDS 19,726 (25.9%) 2,363 (29.5%) -0.08 9,803 (27.0%) 1,945 (29.0%) -0.05 3,948 (27.8%) 1,556 (28.7%) -0.02 Glucocorticoids 5,461 (7.2%) 859 (10.7%) -0.12 3,003 (8.3%) 721 (10.7%) -0.08 1,354 (9.5%) 582 (10.7%) -0.04 Diuretics 16,045 (21.1%) 2,610 (32.6%) -0.26 8,965 (24.7%) 2,186 (32.6%) -0.18 3,959 (27.8%) 1,760 (32.5%) -0.10 Antidepressants 10,438 (13.7%) 1,277 (15.9%) -0.06 5,199 (14.3%) 1,082 (16.1%) -0.05 2,108 (14.8%) 871 (16.1%) -0.03 Antibiotics 25,094 (33.0%) 2,951 (36.8%) -0.08 12,654 (34.8%) 2,496 (37.2%) -0.05 5,185 (36.5%) 1,991 (36.8%) -0.01
BASELINE CHARACTERISTICS (MEDICATION) Range 1 Range 2 Range 3 40-47 48-56 Std diff 42-47 48-53 Std diff 44-47 48-51 Std diff Total 76,031 8,011 36,360 6,710 14,216 5,416 Statins 17,039 (22.4%) 1,969 (24.6%) -0.05 8,924 (24.5%) 1,683 (25.1%) -0.01 3,647 (25.7%) 1,389 (25.6%) 0.00 NSAIDS 19,726 (25.9%) 2,363 (29.5%) -0.08 9,803 (27.0%) 1,945 (29.0%) -0.05 3,948 (27.8%) 1,556 (28.7%) -0.02 Glucocorticoids 5,461 (7.2%) 859 (10.7%) -0.12 3,003 (8.3%) 721 (10.7%) -0.08 1,354 (9.5%) 582 (10.7%) -0.04 Diuretics 16,045 (21.1%) 2,610 (32.6%) -0.26 8,965 (24.7%) 2,186 (32.6%) -0.18 3,959 (27.8%) 1,760 (32.5%) -0.10 Antidepressants 10,438 (13.7%) 1,277 (15.9%) -0.06 5,199 (14.3%) 1,082 (16.1%) -0.05 2,108 (14.8%) 871 (16.1%) -0.03 Antibiotics 25,094 (33.0%) 2,951 (36.8%) -0.08 12,654 (34.8%) 2,496 (37.2%) -0.05 5,185 (36.5%) 1,991 (36.8%) -0.01
FIRST HBA1C MEASUREMENT AND THE EVENT RATE (PER PY) 0.08 0.06 Rate 0.04 0.02 0.00 25 30 35 40 45 50 55 60 65 70 75 80 85 90 HbA1c, mmol/mol
FIRST HBA1C MEASUREMENT AND THE EVENT RATE (PER PY) 0.08 0.06 N~7,000 N~3,000 Rate 0.04 N~ 290,000 0.02 N~36,000 0.00 25 30 35 40 45 50 55 60 65 70 75 80 85 90 HbA1c, mmol/mol
RDD MODEL FOR THE EVENT RATE where log(hazard(y i Z i )) = β 0 + β 1 Z i c + β 2 1 Z i < c + β 3 Z i c 1(Z i c) β 1 is the slope and β 0 +β 2 the intercept below the threshold. β 1 +β 3 is the slope and β 0 the intercept above the threshold. β 2 is the difference in intercepts in the threshold.
MODEL FIT RANGE 42-47 MMOL/MOL VERSUS 48-53 MMOL/MOL 0.08 0.06 Rate 0.04 0.02 0.00 25 30 35 40 45 50 55 60 65 70 75 80 85 90 HbA1c, mmol/mol
INTENTION-TO-TREAT ESTIMATES EFFECT OF TREATMENT ELIGIBILITY ON CARDIOVASCULAR OUTCOMES AND MORTALITY Strata Bandwidth Mortality Stroke and AMI Combined All 40-56 0.7690 (0.6817 ; 0.8675) 0.9773 (0.8111 ; 1.1775) 0.8236 (0.7410 ; 0.9153) 42-53 0.7282 (0.6264 ; 0.8467) 0.9623 (0.7625 ; 1.2145) 0.7852 (0.6880 ; 0.8960) 44-51 0.7246 (0.5789 ; 0.9070) 0.8182 (0.5758 ; 1.1626) 0.7412 (0.6085 ; 0.9029) 2006-2011 40-56 0.7860 (0.6916 ; 0.8934) 0.9708 (0.7946 ; 1.1861) 0.8303 (0.7416 ; 0.9296) 42-53 0.7578 (0.6450 ; 0.8903) 0.9758 (0.7588 ; 1.2550) 0.8067 (0.6998 ; 0.9298) 44-51 0.7661 (0.5994 ; 0.9790) 0.8258 (0.5607 ; 1.2161) 0.7736 (0.6228 ; 0.9608) 2012-2014 40-56 0.7876 (0.5465 ; 1.1351) 0.9759 (0.5673 ; 1.6790) 0.8820 (0.6473 ; 1.2016) 42-53 0.6584 (0.4254 ; 1.0191) 0.9282 (0.4857 ; 1.7741) 0.7751 (0.5361 ; 1.1207) 44-51 0.6560 (0.3709 ; 1.1603) 0.8540 (0.3555 ; 2.0517) 0.7232 (0.4438 ; 1.1785)
INTENTION-TO-TREAT ESTIMATES STABLE ESTIMATES EVEN IN SMALLER BANDWIDTHS Strata Bandwidth Mortality Stroke and AMI Combined All 40-56 0.7690 (0.6817 ; 0.8675) 0.9773 (0.8111 ; 1.1775) 0.8236 (0.7410 ; 0.9153) 42-53 0.7282 (0.6264 ; 0.8467) 0.9623 (0.7625 ; 1.2145) 0.7852 (0.6880 ; 0.8960) 44-51 0.7246 (0.5789 ; 0.9070) 0.8182 (0.5758 ; 1.1626) 0.7412 (0.6085 ; 0.9029) 2006-2011 40-56 0.7860 (0.6916 ; 0.8934) 0.9708 (0.7946 ; 1.1861) 0.8303 (0.7416 ; 0.9296) 42-53 0.7578 (0.6450 ; 0.8903) 0.9758 (0.7588 ; 1.2550) 0.8067 (0.6998 ; 0.9298) 44-51 0.7661 (0.5994 ; 0.9790) 0.8258 (0.5607 ; 1.2161) 0.7736 (0.6228 ; 0.9608) 2012-2014 40-56 0.7876 (0.5465 ; 1.1351) 0.9759 (0.5673 ; 1.6790) 0.8820 (0.6473 ; 1.2016) 42-53 0.6584 (0.4254 ; 1.0191) 0.9282 (0.4857 ; 1.7741) 0.7751 (0.5361 ; 1.1207) 44-51 0.6560 (0.3709 ; 1.1603) 0.8540 (0.3555 ; 2.0517) 0.7232 (0.4438 ; 1.1785)
INTENTION-TO-TREAT ESTIMATES ALL ESTIMATES POINT IN THE SAME DIRECTION Strata Bandwidth Mortality Stroke and AMI Combined All 40-56 0.7690 (0.6817 ; 0.8675) 0.9773 (0.8111 ; 1.1775) 0.8236 (0.7410 ; 0.9153) 42-53 0.7282 (0.6264 ; 0.8467) 0.9623 (0.7625 ; 1.2145) 0.7852 (0.6880 ; 0.8960) 44-51 0.7246 (0.5789 ; 0.9070) 0.8182 (0.5758 ; 1.1626) 0.7412 (0.6085 ; 0.9029) 2006-2011 40-56 0.7860 (0.6916 ; 0.8934) 0.9708 (0.7946 ; 1.1861) 0.8303 (0.7416 ; 0.9296) 42-53 0.7578 (0.6450 ; 0.8903) 0.9758 (0.7588 ; 1.2550) 0.8067 (0.6998 ; 0.9298) 44-51 0.7661 (0.5994 ; 0.9790) 0.8258 (0.5607 ; 1.2161) 0.7736 (0.6228 ; 0.9608) 2012-2014 40-56 0.7876 (0.5465 ; 1.1351) 0.9759 (0.5673 ; 1.6790) 0.8820 (0.6473 ; 1.2016) 42-53 0.6584 (0.4254 ; 1.0191) 0.9282 (0.4857 ; 1.7741) 0.7751 (0.5361 ; 1.1207) 44-51 0.6560 (0.3709 ; 1.1603) 0.8540 (0.3555 ; 2.0517) 0.7232 (0.4438 ; 1.1785)
INTENTION-TO-TREAT ESTIMATES (STRATIFIED) EFFECT OF TREATMENT ELIGIBILITY ON CARDIOVASCULAR OUTCOMES AND MORTALITY Strata Bandwidth Mortality Stroke and AMI Combined Women 40-56 0.7200 (0.6019 ; 0.8614) 0.8397 (0.6196 ; 1.1380) 0.7687 (0.6546 ; 0.9026) 42-53 0.6817 (0.5453 ; 0.8522) 0.8038 (0.5549 ; 1.1643) 0.7327 (0.6003 ; 0.8944) 44-51 0.7109 (0.5109 ; 0.9892) 0.5773 (0.3347 ; 0.9958) 0.6916 (0.5149 ; 0.9290) Men 40-56 0.8090 (0.6874 ; 0.9522) 1.0628 (0.8392 ; 1.3461) 0.8638 (0.7507 ; 0.9939) 42-53 0.7660 (0.6243 ; 0.9398) 1.0817 (0.8017 ; 1.4594) 0.8265 (0.6928 ; 0.9859) 44-51 0.7343 (0.5407 ; 0.9972) 1.0462 (0.6602 ; 1.6578) 0.7855 (0.6025 ; 1.0240) Age <60y 40-56 0.7695 (0.5640 ; 1.0497) 1.3664 (0.9491 ; 1.9672) 0.9797 (0.7690 ; 1.2482) 42-53 0.7157 (0.4842 ; 1.0578) 1.5215 (0.9489 ; 2.4395) 0.9912 (0.7289 ; 1.3478) 44-51 0.6522 (0.3615 ; 1.1766) 0.9299 (0.4504 ; 1.9197) 0.7853 (0.4923 ; 1.2527) Age >=60y 40-56 0.8303 (0.7282 ; 0.9468) 0.9197 (0.7395 ; 1.1438) 0.8471 (0.7529 ; 0.9530) 42-53 0.7735 (0.6566 ; 0.9112) 0.8644 (0.6604 PROFESSOR, ; 1.1315) MSC, PHD 0.7840 (0.6770 26 JUNE ; 0.9078) 2017 44-51 0.7788 (0.6108 ; 0.9932) 0.8016 (0.5356 ; 1.1997) 0.7679 (0.6175 ; 0.9550)
INTENTION-TO-TREAT ESTIMATES (STRATIFIED) EFFECT OF TREATMENT ELIGIBILITY ON CARDIOVASCULAR OUTCOMES AND MORTALITY Strata Bandwidth Mortality Stroke and AMI Combined Women 40-56 0.7200 (0.6019 ; 0.8614) 0.8397 (0.6196 ; 1.1380) 0.7687 (0.6546 ; 0.9026) 42-53 0.6817 (0.5453 ; 0.8522) 0.8038 (0.5549 ; 1.1643) 0.7327 (0.6003 ; 0.8944) 44-51 0.7109 (0.5109 ; 0.9892) 0.5773 (0.3347 ; 0.9958) 0.6916 (0.5149 ; 0.9290) Men 40-56 0.8090 (0.6874 ; 0.9522) 1.0628 (0.8392 ; 1.3461) 0.8638 (0.7507 ; 0.9939) 42-53 0.7660 (0.6243 ; 0.9398) 1.0817 (0.8017 ; 1.4594) 0.8265 (0.6928 ; 0.9859) 44-51 0.7343 (0.5407 ; 0.9972) 1.0462 (0.6602 ; 1.6578) 0.7855 (0.6025 ; 1.0240) Age <60y 40-56 0.7695 (0.5640 ; 1.0497) 1.3664 (0.9491 ; 1.9672) 0.9797 (0.7690 ; 1.2482) 42-53 0.7157 (0.4842 ; 1.0578) 1.5215 (0.9489 ; 2.4395) 0.9912 (0.7289 ; 1.3478) 44-51 0.6522 (0.3615 ; 1.1766) 0.9299 (0.4504 ; 1.9197) 0.7853 (0.4923 ; 1.2527) Age >=60y 40-56 0.8303 (0.7282 ; 0.9468) 0.9197 (0.7395 ; 1.1438) 0.8471 (0.7529 ; 0.9530) 42-53 0.7735 (0.6566 ; 0.9112) 0.8644 (0.6604 PROFESSOR, ; 1.1315) MSC, PHD 0.7840 (0.6770 26 JUNE ; 0.9078) 2017 44-51 0.7788 (0.6108 ; 0.9932) 0.8016 (0.5356 ; 1.1997) 0.7679 (0.6175 ; 0.9550)
CSAL TREATMENT EFFECT ESTIMATES (BOR ET AL., EPIDEMIOLOGY, 25(5), PP. 729-737) RDD analysis on the combined endpoint Bandwidth 42-53 Paramter Interpretation Estimate h_above Hazard among eligibles 0.04 h_below Hazard among non-eligibles 0.07 h_above_treat Hazard in always-takers and treated compliers 0.04 h_below_treat Hazard in always-takers 0.00 TreatProb_above Treatment probability among eligibles 0.19 TreatProb_below Treatment probability among non-eligibles 0.15 Prob_complier Difference in treatment probability in the threshold 0.03 ITT_frd Difference in hazards in the threshold (ITT) -0.02 CACE_frd Complier causal difference in the threshold (CACE) -0.73 h_above_complier Treated complier hazard 1.24 h_below_complier Control complier hazard 1.97 CCHR_frd Complier causal hazard ratio (CCHR) 0.63
INTENTION-TO-TREAT (ITT) AND COMPLIER CSAL HAZARD RATIOS (CCHR) Strata Bandwidth Combined event ITT Combined event CCHR All 40-56 0.8236 (0.7410 ; 0.9153) 0.83 42-53 0.7852 (0.6880 ; 0.8960) 0.63 44-51 0.7412 (0.6085 ; 0.9029) 0.44
SUMMARY The effect (ITT) of treatment eligibility showed a 20-30% reduced risk of cardiovascular events and mortality. The reduced risk was mainly driven by mortality. The reduced risk was also observed across age, gender and calendar year. The effect (CCHR) of early metformin treatment showed a more than 20% reduced risk of cardiovascular events and mortality. The estimated difference in treatment probabilities at the threshold is small, i.e., the complier effects (CCHR) are unstable and should be interpreted with caution.
STRENGTH OF THE STUDY Data collection independent of research question A large study population and complete follow-up on everyone High validity and completeness of endpoints
LIMITATIONS Fuzzy treatment assignment No data on primary care diagnoses We did not use information on treatment or HbA1c measurements during follow-up
EVIDENCE FROM CLINICAL TRIALS In a meta-analysis of twelve RCTs among metformin initiators (n = 5,455 ) and comparators (n=8,996), metformin was associated with a significant reduction in cardiovascular events compared to placebo or no therapy (OR 0.79, 95% CI 0.64-0.98) In a small subgroup of overweight patients in the UKPDS trial, metformin-treated patients had a 30% reduced risk of macrovascular complications and a 36% reduced all-cause mortality
CONCLUSION We were able to evaluate the RDD assumption on routinely recorded healthcare data The data showed acceptable balance in observed covariates around the threshold There was no evidence of bunching around the threshold Below the threshold, around 95% were untreated within 3 months Above the threshold, less than 20% were treated within 3 months
CONCLUSION (CON T) Simple ITT analyses showed acceptable robustness in estimates across all bandwidths. The analyses showed that patients eligible for treatment seem to have a reduced risk of cardiovascular events and mortality Complier estimates (CCHR) showed a beneficial effect of early metformin treatment. The CCHR estimates should be interpreted with caution. However, estimates were comparable with the existing literature