Abstract
Background Non-adherence to antihypertensive medication treatment (AHM) is a complex health behavior with determinants that extend beyond the individual patient. The structural and social determinants of health (SDH) that predispose populations to ill health and unhealthy behaviors could be potential barriers to long-term adherence to AHM. However, the role of SDH in AHM non-adherence have been understudied.
Therefore, we aimed to define and identify the SDH factors associated with non-adherence to AHM and to quantify the variation in county-level non-adherence to AHM explained by these factors.
Methods Two cross-sectional datasets, the Centers for Disease Control and Prevention (CDC) Atlas of Heart Disease and Stroke (2014-2016 cycle) and the 2016 County Health Rankings (CHR), were linked to create an analytic dataset. Contextual SDH variables were extracted from both the CDC-CHR linked dataset. County-level prevalence of AHM non-adherence, based on Medicare fee-for-service beneficiaries’ claims data, was extracted from the CDC Atlas dataset. The CDC measured AHM non-adherence as the Proportion of Days Covered (PDC) with AHM during a 365-day period for Medicare Part D beneficiaries and aggregated these measures at the county-level. We applied Confirmatory Factor Analysis (CFA) to identify the constructs of social determinants of AHM non-adherence. AHM non-adherence variation and its social determinants were measured with structural equation models.
Results Among 3,000 counties in the US, the weighted mean prevalence of AHM non-adherence (PDC<80%) in 2015 was 25.0%, Standard Deviation (SD), 18.8%. AHM non-adherence was directly associated with poverty/food insecurity (β=0.31, P-value<0.001) and weak social supports (β=0.27, P-value<0.001), but inversely with healthy built environment (β= −0.10, P-value=0.02). These three constructs explained a third (R2=30.0%) of the variation in county-level AHM non-adherence.
Conclusion AHM non-adherence varies by geographical location, a third of which is explained by contextual SDH factors including poverty/food insecurity, weak social supports and healthy built environments.
INTRODUCTION
Nearly half of hypertension patients 1, 2 are non-adherent to AHM. Non-adherence to AHM is associated with a 27% higher risk of stroke 3 and about 56% higher risk of cardiovascular diseases (CVDs) 4. The reasons for suboptimal AHM adherence are still not well understood, despite extensive research and innovative strategies to improve this complex health behavior 5. Research on the determinants of nonadherence to AHM has focused on individual patient and provider characteristics 6-8.
The structural and social determinants of health (SDH) that predispose populations to ill health and unhealthy behaviors could be potential barriers to long-term adherence to AHM. Several studies have reported social support, food insecurity, poverty and lack of transportation as individual level SDH that are associated with medication non-adherence 9-14. However, the relationships between AHM non-adherence and contextual SDH, i.e., SDH factors measured at the community level, have not yet been well investigated 15. Emerging, but limited, data suggests that residential location and some contextual SDH are potential predictors of medication non-adherence. Nationally representative data from Medicare and commercial health plans beneficiaries in the US showed that adherence to AHM as well as antidiabetic and antilipid medications varied by geographic regions 16. A few studies, all focusing on oral antidiabetic medications, reported that medication adherence is inversely associated with neighborhood social affluence, residential stability, socioeconomic advantage and safety 17-20.
While evidence for a link between contextual SDH factors and medication non-adherence is beginning to grow, this area of research is stymied by lack of access to data and methodological limitations. Administrative healthcare claims databases are routinely used to measure non-adherence based on prescription fill data. However, these databases do not capture data on the social and structural conditions in which patients are living. There are a growing number of databases such as the Public Health Exposome 21, the CHR, and the CDC Atlas for Heart Disease and Stroke with information on contextual SDH factors. However, linking these databases with claims databases remains a challenge because of patient data privacy issues. Another challenge is the lack of a standard methodological framework defining the relationships between contextual SDH factors and medication non-adherence.
The primary objective of this analysis was to define and identify the social determinants of non-adherence to AHM and to quantify the variation in county-level non-adherence to AHM explained by these factors. We leveraged data on county-level prevalence of AHM non-adherence from the CDC Atlas of Heart Diseases and Stroke dataset and 72 unique county-level SDH variables from this dataset and the CHR datasets to implement our study objectives. A confirmatory factor analysis techniques (CFA) was applied to create constructs of social determinants of AHM non-adherence. The relationships between these constructs and AHM non-adherence were then tested with structural equation modeling.
METHODS
Data sources
The CDC Atlas of Heart Disease and Stroke (2014-2016 cycle) and the 2016 CHR datasets were linked by unique five-digit Federal Information Processing Standard (FIPS) codes of each county present in both datasets. The CDC Atlas consists of county-level estimates of all heart diseases mortality and hospitalizations based on data from the Deaths National Vital Statistics System and Centers for Medicare and Medicaid Services Medicare Provider Analysis and Review (MEDPAR) file, Part A, respectively 22.
This database also contains county-level measurements of risk factors, social and economic factors, health care delivery and insurance and health care costs data derived from multiple data sources 23. The CHR database is the most comprehensive dataset created specifically to characterize nearly all counties (3,000) by health outcomes (length of life and quality of life) and overall health factors (health behaviors, clinical care, social and economic factors, and physical environment).
Measurement of county-level non-adherence to AHM
The CDC Atlas provides medication adherence data calculated from prescription drug claims data for Medicare Advantage or Medicare fee-for-service beneficiaries aged ≥65 years with Medicare Part D coverage. AHM non-adherence was operationally defined as the proportion of days covered (PDC) with blood pressure medication for a period of 365 days. PDC<80% was considered as non-adherence to AHM. PDC is a validated measure of adherence and persistence to medications, especially among patients with repeated fills 24-26. The AHM PDC measure in the CDC Atlas data combined all individual antihypertensive agents into a single class of AHM, while still a valid measure this could potentially underestimate non-adherence since this method requires a user to stay on only a single AHM agent. All variables derived from Medicare, including PDC, were based on claims data of beneficiaries who were ≥65 years old on the extract year. Therefore, for the 2014 – 2016 CDC Atlas data, non-adherence was measured for the 2015 to 2016 period. Further details for the inclusion of Medicare data for the calculation of non-adherence are available at 22.
SDH variables
The CHR model describes communities with respect to how healthy they are (health outcomes) and existing modifiable factors that predict future health (health factors). For the purposes of ranking counties by these factors, the CHR calculated weighted composite scores for domains of both community characteristics 23. The health outcomes domains are (1) length of life and (2) quality of life; the health factors domains are: (1) health behaviors; (2) clinical care; (3) social and economic factors; and (4) physical environment. The rationale and methods for creating these domains and composite scores have already been published 27. Prior to the PCA and CFA, all variables were standardized to the means of counties within each state using methods previously applied by CHR investigators to create SDH domains 27. This standardization was necessary for the current analysis because the variables included in our analysis had different measurement scales.
Definition of social determinants of medication non-adherence (SDN) constructs
There are no standardized frameworks for identifying and studying the relationships between SDN constructs and medication non-adherence. Therefore, we used Confirmatory Factor Analysis (CFA) techniques to create a SDN measurement model to relate the observed data to four latent SDN constructs, namely: food insecurity, poverty, weak social support, and healthy built environments (Figure 1). This four-factor CFA measurement model was created on the basis of existing literature on the relationships between contextual SDH factors and medication non-adherence 15, 17, 19, as well as the availability of variables in the CHR-CDC Atlas linked dataset. We decided to define social determinants of non-adherence constructs rather than use the pre-defined SDH domains in the CHR dataset for the following reasons: 1) the CHR-based SDH factors included only 30 variables while leaving out nearly 27 other variables that may be important predictors of medication non-adherence; 2) additional potential predictors of medication non-adherence included in the CDC Atlas dataset were not included in defining SDH domains by CHR; 3) the five domains of SDH defined by CHR were based on their relative contributions to health outcomes (length and quality of life) and therefore may not have the same impact on a health behavior such as AHM non-adherence 28.
Statistical analyses
Statistical analyses were implemented with SAS (version 9.4) and IBM SPSS AMOS (version 25). The goal of these analyses was to describe and quantify the county-level variation in non-adherence to AHM that is explained by social determinants of non-adherence constructs. Structural equation models (SEMs) were used to measure the proportion of variation in county-level AHM non-adherence explained by social determinants of non-adherence constructs identified through the application of PCA and CFA techniques. Briefly, SEM is a hybridized form of analysis of variance (ANOVA), and CFA that is well suited to investigating complex relationships between dependent and independent variables. The hypothesized paths and relationships that were tested are represented in Figure 2. Because variables were standardized prior to creating social determinants of non-adherence constructs, there was no need to account for potential clustering at the state level when estimating variance and calculating p-values for the associations between constructs of social determinants of non-adherence and AHM non-adherence.
Indicators per construct. Poverty/food insecurity; Percent below poverty line (%), Uninsured (%), Food stamp/SNAP recipient (%), Food insecurity (%).
Weak social support; Children in single-parent household (%), Families with female household head (%), Residential segregation, non-white/white (%)
Healthy built environment; Population living within half a mile of a park (%), Severe housing problems (%), Access to exercise opportunities (%)
Only counties with complete data for non-adherence measurements and variables included in the creation of SDH constructs were included in the analysis. The geographic patterns of AHM non-adherence prevalence were described visually with the CDC Atlas maps tool 22. Although cross-sectional CDC and CHR datasets were used for these analyses, the primary dependent variable was measured in 2015 whereas all predictors were measured in the preceding years, 2006 – 2014.
RESULTS
The weighted mean prevalence of AHM non-adherence (PDC<80%) in 2015 among the 2067 counties (out of 3,000) included in the analysis was 25.0%, standard deviation (SD), 18.8%. The counties that were not included in the analysis did not have sufficient data on AHM non-adherence. There was a 7.1% difference in the prevalence of AHM non-adherence between the upper 90th (28.5%) and lower 10th (21.4%) percentiles of counties ranked by the prevalence of AHM non-adherence. A higher prevalence of non-adherence was concentrated within states located in the South (mean=28.4; SD, 3.2) compared to the Midwest where the prevalence was the lowest (mean=20.9; SD, 2.9) (Figure 3).
Table 1 describes the distribution of indicators of the SDN constructs and county-level demographic factors across geographic regions. The proportion of county populations aged ≥65 years were not significantly different across the four geographic regions. The rest of the demographic measures were unequally distributed across geographic regions. More than half of the counties in the South were considered to be rural, while only about a third of the Northeastern counties were rural. The county-level prevalence of all indicators of the SDN constructs differed significantly across geographical regions. The prevalence of all indicators of poverty were highest in the South. While the prevalence of households with single parents and female heads was highest in the south, the Southern (31.1%) counties were less segregated compared to the Midwestern (35.0%) and Northeastern regions (39.6%). Similarly, there were fewer parks per county (P<0.01), and access to exercise opportunities were lowest in the South (P<0.01). The prevalence of severe housing problems was highest the West (18.6%) compared to a low of 13.0% in the Midwest (P<0.01).
Confirmed social determinants of non-adherence model
The hypothesized four-factor CFA model was not confirmed by the data. Model fit was assessed based on the following indices: chi-square, Aikeke Information Criterion (AIC), Root Mean Square Error of Approximation (RMSEA), and Comparative Fit Index (CFI). e-Table 1 provides a summary of the comparative fit between the four-factor base model (Model 1) and a final three-factor model (Model 2), comprised of poverty/food insecurity, weak social support, and healthy built environments constructs. In this three-factor model the poverty and food insecurity constructs were collapsed into a single construct labelled poverty/food insecurity because of high intercorrelation between these two constructs, further explanation is provided in the next paragraph.
The factor loadings (standardized regression coefficients) and discriminant validity measures from the hypothesized CFA model were evaluated to determine how well the model captured SDN constructs. In the base model, Model 1, an indicator representing counts of pharmacies had a very low factor loading on the healthy built environment (β=-0.06). The food environment index variable also had standardized loadings >1 in the food insecurity construct (β=-1.39). Therefore, Model 1 was modified by excluding these indicators from their respective constructs to create Model 2. While all remaining indicators loaded well (from 0.23 to 0.90) on their respective constructs, the poverty and food insecurity constructs were highly inter-correlated (R2=0.98) suggesting Model 2 had a poor discriminant validity. For this reason, indicators of the food insecurity and poverty constructs were collapsed into a new construct labelled poverty/food insecurity. Indicators were sequentially loaded onto this new construct and those with the highest factor loadings were retained. The final poverty/food insecurity construct had high internal consistency (Cronbalch-alpha=0.83) and comprised of Percent below poverty line (%), Uninsured (%), Food stamp/SNAP recipient (%), Food insecurity (%) variables. This newly created construct was then added to the weak social support and healthy built environment constructs retained from Model 2 to create Model 3. This new model, consisting of three constructs – poverty/food insecurity, weak social support and healthy built environment – had a better discriminant validity based on the inter-correlations between poverty and healthy built environment (R2=-0.08), weak social support and built environment (R2=0.25) and between poverty and weak social support (R2=0.81). A test of chi-square difference (χ2=34.49 [degrees of freedom=6], P-value<0.001) between Models 2 and 3 indicated that Model 3 was a relatively better fit to the data. Therefore, Model 3 was retained for the structural path analysis. The standardized regression coefficients of indicators per each construct in Model 3 are presented in (e-Table 2).
Structural relationships between AHM non-adherence and SDN
The SEM results are summarized in Table 2 and Figure 2. The structural equation analysis confirmed that the relationships between AHM non-adherence and SDH are through associative paths involving three SDN constructs which are subdomains of SDH. All three SDN constructs were independently associated with AHM non-adherence; both poverty (β=0.31, P-value<0.001) and weak social support (β=0.27, P-value<0.001) were positively associated with AHM non-adherence. On the contrary, healthy built environment was inversely associated with AHM non-adherence (β= −0.10, P-value<0.01). Of all three, poverty/food insecurity had the strongest association with AHM non-adherence. Together, these three social determinants of medication non-adherence explained a third (R2=30.3%) of the total variation in county-level non-adherence to AHM (Figure 2).
The independence of these AHM non-adherence and SDN relationships from demographic factors were tested in sensitivity analysis. County-level distributions of race (% African-American), gender (% Female), and proportion of rural areas were included in the SEMs as exogenous variables to control for their potential confounding effects. The effects of all three SDN constructs remained statistically significant even after adjusting for these demographic factors, e-Figure 1. This suggests that the associative relationships between all three SDN constructs and AHM non-adherence are independent of these demographic factors.
DISCUSSION
Non-adherence to AHM was prevalent nationally and varied by regional and county-level locations. Using a theoretical framework and a SEM approach, poverty/food insecurity, weak social supports and healthy built environments were observed to be associated with AHM non-adherence. These findings build on a limited but growing body of literature on the relationships between residential locations and medication non-adherence.
Previous studies on this topic have been focused on testing the relationships between individual level adherence to oral antidiabetic medications and neighborhood level social and economic factors among populations with diabetes. Among 749 Mexican-Americans treated for diabetes at a University-affiliated clinic, patients who lived in neighborhoods ranked among the top quartile on neighborhood deprivation index were 60% more likely (odds ratio [OR]=1.64; 95% CI:1.12 to 2.39) not to adhere to their medications, compared to those living in the lowest quartile of neighborhoods on the deprivation index 29. In a much smaller sample (n=179), McClintock et. al. (2015) reported that the patients with diabetes who lived in neighborhoods characterized as high in social affluence, residential stability, and neighborhood advantage were about eight times (OR=8.48, 95% CI:1.71 to 42.02) more likely to adhere to oral antidiabetic medications compared to those living in neighborhoods with fewer of these environmental features 17. Among participants in the California Health Interview Survey diabetes patients living in unsafe neighborhoods were more likely to delay filling medications, OR=1.69, 95%CI:1.19 to 2.40 19. A recent study by Qato et. al. (2019) reported a significant reduction in adherence to cardiometabolic medications among patients living in low-income neighborhoods with fewer pharmacies (−7.98%; 95% CI, −8.50% to −7.47%) compared with their counterparts 20. These data, together with ours, provide evidence that the social, economic and the built environments of residential locations are potentially influential in how populations might adhere to treatment with medications.
There are several implications of our study with respect to advancing research on the topic of social determinants of AHM non-adherence. First, linking pharmaceutical claims data with aggregated databases that contain contextual SDH data is helpful in studying the potential sources of AHM non-adherence beyond individual level SDH factors. Second, owing to the complexity of large databases that capture several individual contextual SDH variables, dimension reduction methods such as CFA are highly useful in defining SDH constructs. Furthermore, the application of causal frameworks such as SEM can be a great aid in testing specific pathways and relationships between SDH constructs and AHM non-adherence. Third, our findings buttress those from previous studies to show that the social, economic and the built environments are associated with AHM non-adherence. This county-level data could be useful in formulating population-level interventions to improve AHM adherence in patient populations. Individual-level data may be liable to ecological fallacy when applied to address a population-level outcome such as county-level AHM non-adherence 30.
LIMITATIONS
This study has limitations that should be considered when interpreting findings. First is ecological fallacy - our findings may not directly translate into patient-level effects because only aggregate-level data was available for analysis. We therefore recommend future research on this topic to investigate the role of contextual SDH factors on AHM non-adherence at the individual level. Second, because only Medicare Part D beneficiaries’ data were used for calculating the prevalence of county-level non-adherence to AHM, our findings may not be generalizable to those younger than 65 years of age. Third, although we leveraged aggregate county-level data from the CHR database (one of the most comprehensive databases on SDH) and additional variables from CDC Atlas of Heart Diseases, there were limited measurements of structural SDH factors about policy (social, health and economic) and the built environment. Nearly a third of the counties did not have sufficient data to the measurement of AHM non-adherence. Future studies might well aim to augment the wealth of SDH data in both datasets with additional structural SDH factors from other databases.
STRENGTHS
Our study has some strengths over previous examinations of these issues. First, through the application of dimension reduction methods, we were able to define contextual SDH constructs that are specific to medication non-adherence. Second, the application of SEM allowed us to explore potential determinants of AHM non-adherence at the county-level. Third, our findings fill a critical gap in the lack of information on community-level determinants of AHM non-adherence. This information may be important in formulating population-level interventions to improve AHM adherence in communities. Of course, before this data can be adapted for designing interventions the interactions between contextual and individual-level SDH factors and their cumulative effects on AHM non-adherence must be investigated first.
CONCLUSIONS
AHM non-adherence varies by geographical location, a third of which is accounted for by contextual SDH factors including poverty/food insecurity, weak social supports and healthy built environments. Given that these SDH constructs accounted for about a third of the variation in AHM non-adherence at the county level, future studies would be well advised to investigate further how these contextual factors interact with individual level factors to influence AHM non-adherence.
Data Availability
All data used for this manuscript are publically available on the websites described in the manuscript.
Sources of Funding
None
Figure legends
e-Figure 1: Relationships between social determinants of non-adherence and antihypertensive medication non-adherence adjusted for demographic factors
Indicators per construct. Poverty/food insecurity; Percent below poverty line (%), Uninsured (%), Food stamp/SNAP recipient (%), Food insecurity (%).
Weak social support; Children in single-parent household (%), Families with female household head (%), Residential segregation, non-white/white (%)
Healthy built environment; Population living within half a mile of a park (%) Severe housing problems (%), Access to exercise opportunities (%)
Acknowledgments
We are greatly appreciative of Ms. McKenna Guilds, BS, College of Pharmacy, The Ohio State University, Columbus, OH, for her support in helping prepare and format tables and figures.
Footnotes
Disclosures: None of the authors have any disclosures to make.