Background The healthy user bias is usually overlooked as an explanation in studies in which a strong association is found between poor patient medication adherence and worse disease outcomes. Such studies are increasing in frequency across disease states and influence clinical practice. Adherence to antihypertensive medications was studied to illustrate confounding in such studies.
Methods Using data from veterans with hypertension starting antihypertensive treatment, causal models were developed that predicted the risks of hospitalisation, myocardial infarction (MI) and death associated with poor adherence (<80%) while adjusting for patient demographics, baseline disease severity and disease comorbidity. In a second set of otherwise identical models, adjustment was made for time-varying blood pressure (BP), thus controlling for adherence effects that were mediated through the main pharmacological effects of the drugs. It was hypothesised that the second set of models would reveal a positive association between poor adherence and adverse disease outcomes that is largely explained by unmeasured confounders, including health-related behaviours.
Results The models that did not adjust for time-varying BP levels showed that patients with poor adherence had statistically significantly increased risks of 3.7% for hospitalisation, 28.1% for MI and 23.3% for death. These estimates exceed the benefits of these drugs demonstrated by clinical trials. When controlling for time-varying BP, the increased risks were similar (3.4% for hospitalisation, 27.7% for MI and 23.4% for death). The findings were consistent across a range of adherence thresholds (50–90%) and when allowing disease status variables to vary.
Conclusions The associations between poor adherence and outcomes are largely independent of the pharmacological effects of the drugs on BP control as well as commonly measured patient covariates. This finding suggests that even carefully designed observational adherence studies using rich clinical data are impossibly confounded and probably overestimate the true magnitude of the effect. Clinical practice guidelines based on reported adherence effects should be reconsidered.
- angiotensin receptor blockers
- β blockers
- calcium channel blockers, epidemiology
Statistics from Altmetric.com
Observational studies attempting to show increased risks of adverse outcomes due to poor adherence are becoming increasingly popular across a wide range of disease states.1–10 These types of studies are often used to support decisions in clinical care. For example, three observational studies that showed increased risks of death, myocardial infarction (MI) and stent thrombosis in patients who discontinued thienopyridine therapy following stent placement11–13 prompted a new clinical guideline recommending that the duration of antiplatelet therapy be extended following stent placement.14 In hypertension, reported effect sizes of low versus high adherence include a 15–194% increased risk of any-cause hospitalisation,15–18 43% increased risk of cardiovascular hospitalisation,18 15–61% increased risk of acute cardiovascular disease events3 19 and 50–74% increased risk of death.17 20 This is in contrast to intention-to-treat results of clinical trials where patients who took antihypertensive medications versus no exposure experienced no statistically significant effect on hospitalisation, only 11–30% decreased risk of cardiovascular events and only 11–17% decreased risk of death.21 When we examine the inference that imperfect adherence to antihypertensive agents could put patients at greater risk than no exposure at all, it suggests that the causal effect of adherence on outcomes is overestimated in these observational studies. Consequently, decisions to change care recommendations may be misguided.
One potential explanation for the overestimation in such studies is the ‘healthy user’ hypothesis,22 which suggests that adherence to pharmacotherapy may be a surrogate marker for other important health-related behaviours, as illustrated in figure 1. This bias has been proposed as an explanation for the now outmoded belief that hormone replacement therapy could protect postmenopausal women from cardiovascular disease events,23 the supposed effect of lipid-modifying therapy on preventing hip fractures24 and an apparent protective benefit of statins in community-acquired pneumonia.25 Perhaps the most compelling evidence of the healthy user effect comes from a meta-analysis showing decreased mortality in clinical trial patients who were adherent to placebo compared with those who were non-adherent to placebo.26 In most studies, including adherence studies, healthy user bias is overlooked as a potential explanation for findings.
A study was undertaken to explore the role of healthy user effects in these types of studies. We used the best available methods and secondary observational data to create causal models for the associations between poor patient adherence to antihypertensive medications and three outcomes including MI, any-cause hospitalisation and death. We then contrasted two versions of these models—those that adjusted for time-varying blood pressure (BP) levels and those that did not. We theorised that adjusting out the indirect effects of adherence on outcomes that were mediated through effects on BP—while also adjusting for patient demographics, disease severity and comorbid conditions—would show that unmeasured confounding was present since there is no plausible explanation for direct effects of adherence on these outcomes.
Study design and patients
A cohort design was used to study the risks of any-cause hospitalisation, MI and any-cause death in a nationally representative sample of veterans from six regions of the USA (Northwest, Southwest, West, Midwest, Northeast and Southeast). From these data we randomly sampled 100 000 patients who (1) had received one of five antihypertensive classes considered by the JNC-7 to be cardioprotective and first-line27 (ie, thiazide diuretics, ACE inhibitors, angiotensin receptor blockers, β-blockers and calcium channel blockers) between 1 October 2002 and 31 December 2004; (2) were at least 18 years old on the index date (the date of the first prescription for an antihypertensive agent); and (3) had no prescription for an antihypertensive agent in the year prior to the index date. Patients were then excluded if they did not have at least one BP reading on or before the index date and at least one during the follow-up period.
Sample size was determined by calculating the number of observations needed to avoid over-fitting predictive multivariable models using the heuristic that 10 events were required for each covariate. From a similar analysis conducted using data from a subset of veterans in our western region,28 we knew that we might have up to 50 levels of covariates in our final models, meaning that we would need at least 500 events. From the same data, we also calculated that the rarest of our three endpoints of interest (MI) occurred in about two patients per 1000 person-years with an average follow-up time of 4.1 years, suggesting we would need approximately 65 000 patients to avoid over-fitting. We also estimated that we would exclude about 25% of our cohort because of a lack of baseline and follow-up BP measures and determined that randomly sampling 100 000 patients would yield enough patients to avoid over-fitting.
The primary independent variable was adherence to antihypertensive medications, which was calculated using a cumulative measure of medication availability (CMA).29 The CMA was based on pharmacy refill data from the VA's Decision Support System datasets and was calculated as the sum of days on which medications were available in the numerator divided by the number of days between the first fill and the last day of supply dispensed with the last fill of an antihypertensive therapy in the denominator. The CMA was dichotomised at an adherence threshold of 80%, with patients whose ratio was ≥0.8 considered adherent and patients whose ratio was <0.8 considered non-adherent. The adherence measure was calculated for all five classes of antihypertensive medications together; if patients switched between classes they were considered to be adherent as long as their overall adherence ratio for antihypertensive medications was ≥0.8. Adherence was treated as cumulative and time-varying; ratios were calculated for each quarter following the initiation of treatment and updated at each subsequent quarter. Adherence was lagged so that cumulative adherence in each quarter predicted events in the subsequent quarter until patient censorship.
Confounding was handled by direct adjustment for other independent variables. These were selected based on their known or theoretical associations with the outcomes of interest and their availability in the dataset. Covariates included sociodemographic characteristics (sex, marital status and socioeconomic status (SES)); disease severity (baseline BP and levels of ischaemic disease); and comorbid conditions (chronic obstructive pulmonary disorder, diabetes, cancer, psychological disorders, physical disabilities, smoking disorder, sleep apnoea and thyroid disorder). Covariates were obtained using vital signs from the VA's Corporate Data Warehouse; labs, diagnosis codes and procedure codes from the VA's Medical SAS datasets (inpatient and outpatient records so named because they are kept in SAS format); and pharmacy refill data from VA's Decision Support System datasets. The primary models included only baseline adjustment for each of these covariates.
The outcomes of interest were any-cause hospitalisation, MI and any-cause death. Any-cause hospitalisation information was captured from the inpatient files of the Medical SAS datasets. MI was identified from inpatient and outpatient files from Medical SAS datasets using the ICD-9 codes 410.01 (acute MI of the anterolateral wall, initial episode of care), 410.11 (acute MI of other anterolateral wall, initial episode of care) and 410.21 (acute MI of inferolateral wall, initial episode of care). The occurrence of death was identified using the Beneficiary Identification Records Locator Subsystem death file kept by the VA Benefits Administration which contains information on veterans known to be deceased. The file is updated weekly by matching with the Social Security Administration Death Master File. The file is also updated when a hospital, cemetery, relative or acquaintance of a veteran reports a veteran's death.
Descriptive statistics were used to characterise the baseline demographic data of the cohort. Univariate and multivariable Cox proportional hazards models with time-varying exposures were used to measure the association between independent variables and each outcome. Multivariable models were constructed by reviewing the literature and considering the causal mechanisms for variables that can potentially confound the relationships of interest, regardless of whether the associations were found to be statistically significant in this dataset.
We developed models that incorporated all of the available variables that would theoretically be expected to be important for unmasking the relationship between patient adherence to antihypertensive medications and hospitalisation, MI or death. For each endpoint we contrasted two models that were identical in all respects except that one did not adjust for time-varying BP level and one did. By adjusting for time-varying BP, we were essentially blocking the pathway between adherence and outcomes that is mediated through the pharmacological effects of the drug. We theorised that any association that remained after blocking that pathway must be mediated through other pathways. Because we had also adjusted for patient demographics, disease severity and comorbidity in both sets of models, the remaining pathways must be limited and must include unmeasured healthy behaviours.
All models included baseline BP so the analyses controlled for the follow-up level of BP in addition to the changes. In the models that adjusted for time-varying BP as an explanatory variable, follow-up BPs were coded as five categorical variables: low (diastolic ≤70), normal (systolic <120 and diastolic and >70 but <80); pre-hypertension (systolic 120–139 or diastolic 80–89); stage 1 hypertension (systolic 140–159 or diastolic 90–99); and stage 2 hypertension (systolic ≥160 or diastolic ≥100). Patients were censored at the occurrence of an event, in the quarter following discontinuation of the drug, or 31 December 2007, whichever occurred first. All statistical tests were performed using Stata SE V.10 and SAS V.9.
In accordance with guidelines from the International Society for Pharmacoeconomics and Outcomes Research,30 multiple alternative adherence thresholds were evaluated in sensitivity analyses including 50%, 60%, 70% and 90%. Sensitivity analyses were also conducted that (1) did not allow BP control to vary over time (ie, BP was coded as the hypertension level at the end of the entire follow-up period); (2) allowed other explanatory variables to vary over time, with patients contributing person-time to different levels of each covariate on a quarterly basis if their conditions changed over time; and (3) excluded the days supplied with the last fill from both the numerator and the denominator for the calculation of the adherence ratio, since the algorithm we used assumes 100% adherence in that interval. Time-varying independent variables were lagged by one quarter so that the statistical analysis captured the effect of patient characteristics on outcomes in the subsequent quarter. By lagging these variables we avoided the situation in which a change in disease status that occurred at the end of a quarter predicted an event that occurred earlier in the same quarter.
Patients and blood pressure readings in the follow-up period
Of 100 000 randomly sampled veterans who met our inclusion criteria, a total of 80 359 had at least one BP recorded at baseline and in the follow-up period. Of these, 4475 were excluded because of missing zip codes (from which the SES index was calculated), leaving 75 884 patients in the final analysis. Baseline demographic characteristics for the cohort are summarised in table 1. The mean age of the veterans was 65.7 years and they were predominantly male (96.7%). By definition, patients were approximately evenly distributed among quartiles of SES index categories, an index based on several demographic variables obtained by linking the zip code with 2000 census data.31 The mean (SD) number of post-index BPs per patient was 30.2 (62.2); in each quarter each patient had a mean of 2.1–5.0 BP observations. The mean (SD) change from baseline for systolic BPs in each quarter following the index date ranged from −3.05 (25.2) mm Hg in the first quarter to −15.7 (30.0) mm Hg in the penultimate quarter for systolic BP and from −1.8 (14.8) mm Hg in the first quarter to −8.9 (16.0) mm Hg in the last quarter for diastolic BP.
A total of 14 315 patients in our study cohort had a hospitalisation, 1063 had an MI and 1132 died between the index date and discontinuation of the index regimen over an average follow-up time of 3.2–3.5 years. This corresponds to incidence rates of 61.4 hospitalisations, 4.0 MIs and 4.2 deaths per 1000 person-years. Table 2 shows a summary of the risks for each outcome associated with poor adherence to antihypertensive medications for the two models in the base case as well as sensitivity analyses. In the first set of models (those that did not control for time-varying BP levels), a cumulative adherence ratio of <80% was associated with a 3.7% increased risk of hospitalisation (95% CI 1.01 to 1.07), a 28.1% increased risk of MI (95% CI 1.08 to 1.52) and a 23.3% increased risk of death (95% CI 1.06 to 1.43) compared with patients with good adherence, while adjusting for baseline demographics, disease severity and comorbidity. When also controlling for time-varying BP, the magnitude of the association between poor adherence and each outcome was very similar. In this second set of models, poor adherence was associated with a 3.4% increased risk of hospitalisation (95% CI 1.00 to 1.07), a 27.7% increased risk of MI (95% CI 1.08 to 1.52) and a 23.4% increased risk of death (95% CI 1.06 to 1.44). Sensitivity analyses consistently showed similar results.
Our findings support the conclusion that much of the association between poor adherence and worse outcomes is mediated through a pathway other than the pharmacological effects of antihypertensive medications on BP. Similar to the findings of other researchers,3 15–19 we showed a substantial increase in the risk of MI (28%) and death (23%) associated with poor adherence to antihypertensive agents over an average follow-up of more than 3 years. Counter to intuition, this difference was present whether or not we controlled for time-varying BP, a finding that was consistently reproduced in all sensitivity analyses. If a considerable amount of the association between poor adherence to antihypertensive medications and patient outcomes was mediated through the pharmacological action pathway illustrated in the unshaded area of figure 1, then we would have observed a considerable difference between the models that did and did not adjust for time-varying BP. Instead, our results are consistent with the hypothesis that most of this association is independent of the pharmacological actions of the drug on BP and must have other causes, such as the ‘healthy user’ pathway illustrated in the shaded area of figure 1.
The problem we illustrate here is conceptually similar to the problem of reporting ‘per-protocol’ or ‘as-treated’ analyses versus ‘intention-to-treat’ (ITT) analyses in the clinical trial literature. Essentially, it has been repeatedly shown that non-adherent patients in randomised controlled trials (RCTs) have worse outcomes and that it is both biologically implausible and contradicted by the results of ITT analyses to suppose that this usually strong association is a treatment effect.32–34 This conclusion has resulted in a CONSORT (Consolidated Standards of Reporting Trials) guideline that discourages exclusive reporting of per-protocol analyses in RCT publications.35 Adherence studies are analogous to per-protocol analyses in which patient characteristics determine both the extent to which patients will experience the exposure and their degree of risk for adverse outcomes, resulting in highly confounded effect estimates.
We considered a variety of alternative observational methods that might avoid the confounding that we demonstrated: First, similar causal models developed using logistic regression have the same structural confounding shown in figure 1. Second, the case crossover design and other within-subject designs are able automatically to address time-invariant confounding. These methods have been developed to evaluate the acute effects of intermittent exposure, but their use is not as well understood under time-varying treatment and cumulative effects of drug exposures.36 37 Moreover, they are not suited to the adherence problem because a strong and probably untenable assumption must be made that other important healthy user factors do not covary with adherence. Third, instrumental variables may be useful to minimise unmeasured confounders, but there are no known instruments that are useful for adherence studies in the observational setting. Finally, newer techniques such as marginal structure models and inverse probability weighting may be better at capturing both the exposure and the measured confounders by treating the exposure and covariates as processes that evolve over time, but they do not adequately control for unmeasured patient confounders. It is most likely that the only way to assess the causal effect of adherence on outcomes is in the RCT setting and with the use of instrumental variables and nested structural models,38–40 methods that do not directly apply to the observational setting.
We were surprised that the association between poor adherence and risk of hospitalisation did not change discernibly between analyses with and without time-varying BP. This finding was consistently reproduced throughout all of the sensitivity analyses. Upon examination of our data, we theorised that this lack of difference could be explained by a possible floor effect. The mean decrease in BP from baseline was greater among non-adherent patients by about 3 mm Hg in every quarter. The adherent patients, whose BPs were lower at baseline, could not experience reductions of the same magnitude as the non-adherent patients because their lower baseline BP afforded less room to improve. Nonetheless, the fact that we saw a consistent decrease in BP in both groups suggests that the antihypertensive medications did have a clinical effect on the patients.
A central question for our analysis is whether or not it is valid to assume that BP reduction mediates the effects of antihypertensive agents on outcomes. There is strong clinical trial evidence that antihypertensives both reduce BP and prevent the outcomes of interest,21 but it is difficult to show a strong causal link between BP and outcomes. Vittenhoff recommends that three conditions should be present for concluding a variable is a mediator: (1) the predictor of interest (adherence) predicts the mediator (BP); (2) the mediator predicts the outcome (MI, hospitalisation or death) in a model controlling for the predictor of interest; and (3) adding the mediator to a multivariable model for the outcome attenuates the association for the predictor of interest.41 The first condition is difficult to show in the observational setting because of the unmeasured confounding that we highlight here. However, the matter is really a question of dose-response: poorly adherent patients are exposed to lower doses on average than adherent patients. We can infer from phase 2 clinical trials that higher levels of exposure are associated with greater reductions in BP than lower levels,42–46 thus meeting the first condition. When we examined our data for the presence of the second condition we found that, indeed, patients whose treated BP level was in the normal range had a lower risk of events than patients with high treated BP levels, when adjusting for adherence. The third condition—that the addition of the mediator to a multivariable analysis attenuates the estimated coefficient for the predictor of interest—was also met, but the magnitude of the difference was very small. We contend that it was small because most of the effect appears to be mediated through other pathways, including the healthy user pathway illustrated in figure 1.
There could be other explanations for our findings. One could be that all of the pharmacological effects on outcomes are mediated through other direct effects of the drugs that are independent of changes in BP. For example, ACE inhibitors have renal protective effects in patients with diabetes,47 48 and β blockers after MI improve outcomes in ways that are not mediated through BP.49 The magnitudes of the associations between these direct effects and outcomes of interest are not well quantified, but it would be highly speculative to conclude that these alternative effects would overshadow the effects mediated through BP. Second, adherent patients may differ systematically from non-adherent patients in their access and utilisation of care outside the VA. In that case, VA-based outcomes such as hospitalisation and MI would be captured preferentially in one group, which could potentially overwhelm the effects of BP drugs in our models. However, this effect would not be present for the death outcome, which is reported by the Social Security Administration and is not dependent on a patient's care source. The fact that our results for the death endpoint are consistent with the other two endpoints suggests that differential identification of the outcomes is not driving our findings.
Limitations of the study
This analysis has some limitations. By introducing follow-up BP into the models we may have introduced a new type of confounding—uncontrolled common causes of both BP level and outcome risk—that is not present in analyses directly relating outcomes to adherence. Some of these causes—those that relate to behavioural factors—are included in figure 1; others may be independent of healthy user characteristics. For example, we could postulate a whole host of genetic or environmental factors that lead both to high BP and worse outcomes as well as low BP and better outcomes. Since this might result in many patients with good adherence having high BPs and worse outcomes as well as many others with bad adherence having low BPs and good outcomes, then the overall analysis when pooling these patient types within the two exposure categories would have the effect of attenuating the treatment effect in the models where we adjust for BP. Thus, our results are robust to this issue. Moreover, there may be known confounders that are not adequately measured in our data. For example, although we are able to adjust for diagnoses such as angina (a marker of disease severity), the data afford us little insight into symptom frequency, which would have improved discrimination of baseline severity. However, this latter is a limitation of all such secondary database studies.
We found that the associations between poor adherence to antihypertensive medications and outcomes such as hospitalisation, MI and death are largely independent of the pharmacological effects of these drugs on BP control while adjusting for patient demographics, disease severity and comorbidity. This suggests that other unmeasured confounders must explain most of the association between patient medication adherence and outcomes. Using the current array of epidemiological methods, observational studies that attempt to assess the effect of adherence on outcomes are impossibly confounded by healthy user bias and are likely to overestimate or misidentify any effect.
The authors thank Dr Tom Greene for his thoughtful comments on drafts of this paper and Kristin Knippenberg for her careful editorial assistance. The authors would like to acknowledge the Agency for Healthcare Research and Quality (AHRQ) for a career development award supporting Joanne LaFleur (5K08HS018582), the National Cancer Institute for it's support of Richard E. Nelson (1KM1CA156723), and the VA's Health Services Research and Development (HSR&D) office for a career development award supporting Brian C Sauer.
Funding This project was supported by the Veterans Affairs (VA) Geriatric Research, Education, and Clinical Center (GRECC), the University of Utah Primary Care Research Center Scholarship and Novartis Pharmaceuticals.
Competing interests JL has support from the United States (US) Veterans Affairs (VA) Geriatric Research, Education, and Clinical Center (GRECC), the University of Utah Primary Care Research Center Scholarship and Novartis Pharmaceuticals for the submitted work. JL, REN, BCS and JRN have no relationships with any other companies that might have an interest in the submitted work in the previous 3 years or any non-financial interests that may be relevant to the submitted work. Their spouses, partners, or children have no financial relationships that might be relevant to the submitted work.
Ethics approval This study was conducted with the approval of the University of Utah Institutional Review Board and Salt Lake City Veterans Affairs Office of Research and Development.
Provenance and peer review Not commissioned; externally peer reviewed.
If you wish to reuse any or all of this article please use the link below which will take you to the Copyright Clearance Center’s RightsLink service. You will be able to get a quick price and instant permission to reuse the content in many different ways.