Long-term Outcomes of Peripheral Artery Disease In Veterans: Analysis of the PEripheral ARtery Disease Long-term Survival Study (PEARLS)
========================================================================================================================================

* Saket Girotra
* Qiang Li
* Mary Vaughan-Sarrazin
* Brian C. Lund
* Mohammad Al-Garadi
* Joshua A. Beckman
* Rohit Nathani
* Richard M. Hoffman
* Paul S. Chan
* Subhash Banerjee
* Shirling Tsai
* Dharam J. Kumbhani
* Nicole Minniefield-Young
* Kim G. Smolderen
* Shipra Arya
* Cathy Nguyen
* Michael E. Matheny
* Glenn T. Gobbel

## Abstract

**Background** Contemporary research in peripheral artery disease (PAD) remains limited due to lack of a national registry and low accuracy of diagnosis codes to identify PAD patients in electronic health records.

**Methods & Results** Leveraging a novel natural language processing (NLP) system that identifies PAD with high accuracy using ankle brachial index (ABI) and toe-brachial index (TBI) values, we created a registry of 103,748 patients with new onset PAD patients in the Veterans Health Administration (VHA). Study endpoints include mortality, cardiovascular (hospitalization for acute myocardial infarction or stroke) and limb events (hospitalization for critical limb ischemia or major amputation) and were identified using VA and non-VA encounters. The mean age was 70.6 years; 97.3% were males, and 18.5% self-identified as Black race. The mean ABI value was 0.78 (SD: 0.26) and the mean TBI value was 0.51 (SD: 0.19). Nearly one-third (32.4%) patients were currently smoking and 35.4% formerly smoked. Prevalence of hypertension (86.6%), heart failure (22.7%), diabetes (54.8%), renal failure (23.6%), and chronic obstructive pulmonary disease (35.4%) was high. At 1-year, 9.4% of patients had died. The 1-year incidence of cardiovascular events was 5.6 per 100 patient-years and limb events was 4.5 per 100 patient-years.

**Conclusions** We have successfully launched a registry of >100,000 patients with a new diagnosis of PAD in the VHA, the largest integrated health system in the U.S. The ncidence of death and clinical events in our cohort is high. Ongoing studies will yield important insights regarding improving care and outcomes in this high-risk group.

## Introduction

An estimated 8 to 12 million Americans suffer from lower extremity peripheral artery disease (PAD).1,2 Although PAD shares risk factors (e.g., diabetes, smoking, hypertension) with other atherosclerotic vascular diseases such as coronary artery disease (CAD), patients with PAD experience a higher rate of cardiovascular (myocardial infarction, stroke) and limb events (critical limb ischemia, amputation), which contribute to excess disability, death and medical expenditures.3–6 Importantly, racial and ethnic disparities in the incidence and outcomes of PAD have been reported, with a higher incidence of lower extremity amputations in Black, as compared with White, individuals with PAD.7,8

While clinical practice guidelines recommend a multi-pronged strategy including medications (e.g., statins), risk factor control (blood pressure [BP] and diabetes), lifestyle changes, and revascularization in select patients,9 the above therapies are unevenly applied in clinical practice.10–16 This may be due to a paucity of contemporary evidence regarding existing treatments on long-term outcomes. Unlike CAD, there are also unique challenges to reliably identifying PAD patients within electronic health records (EHR) data; existing PAD diagnosis codes historically have had low positive predictive values ranging from 27.5% to 69.4%.17–20 The limited accuracy of diagnosis codes have hampered the development of large registries to improve care and outcomes in PAD.

We recently developed and validated a natural language processing (NLP) system that can extract ankle-brachial index (ABI) and toe-brachial index (TBI) values from ABI test reports documents in the Veterans Affairs (VA) EHR and use these direct measures of ABI and TBI values to identify patients with PAD.21 In contrast to using only International Classification of Diseases (ICD) codes, our NLP-based approach had a sensitivity of 89.3%, specificity of 92.3%, and positive predictive value of 92.3% for identifying PAD.21 Leveraging this NLP system, we have launched the **Pe**ripheral **AR**tery **L**ong-term **S**urvival Study (PEARLS), a longitudinal study of PAD in the Veterans Health Administration (VHA) – the largest integrated health system in the U.S. In this paper, we describe the creation of the PEARLS registry – a contemporary cohort of newly diagnosed PAD patients identified using our novel NLP system in the VHA EHR. We also provide descriptive data on the baseline characteristics of our cohort including data on demographics, co-morbidities, medications, and longitudinal data on 1-year incidence of mortality, cardiovascular and limb events.

## Methods

This work was approved by the Institutional Review Board (IRB) and Research & Development Committees at Dallas, Iowa City and Tennessee Valley Health System VA centers as well as the IRBs at University of Texas Southwestern Medical Center and University of Iowa.

The primary objective of the study was to develop a longitudinal cohort of patients with incident PAD using our novel NLP system. In the first stage of cohort creation, we executed the NLP system on the entire corpus of ABI report documents in the VA’s EHR. Among patients identified as having PAD, we applied additional criteria to restrict our cohort to patients with a new diagnosis of PAD and those who use VA regularly for healthcare **(Figure 1)**. Next, we identified baseline demographics, co-morbidities, smoking status, medications, vital signs, and laboratory values for our cohort. We followed the cohort longitudinally for clinical events using a combination of VA and non-VA (Medicare) data sources to ensure complete follow up.

![Figure 1.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/21/2024.08.20.24312328/F1.medium.gif)

[Figure 1.](http://medrxiv.org/content/early/2024/08/21/2024.08.20.24312328/F1)

Figure 1. 
Study Flowchart

### Data Sources

The list of data sources used in the current study are summarized in **Table S1.** The VA Corporate Data Warehouse (CDW) is a national repository of clinical data for patients receiving care at a VA facility. Because the VA is a single-payer health care system, it generates data from clinical functions (e.g. clinical notes, procedural and radiology reports) and claims-based functions (e.g. clinical diagnoses and procedures using International Classification of Diseases Ninth and Tenth – Clinical Modification (ICD-9 CM and ICD-10 CM) and Current Procedural Terminology [CPT] codes). The above data are available for inpatient and outpatient encounters. In addition, data on smoking status, vital signs from each clinical encounter, laboratory test results and prescriptions filled at a VA pharmacy are also available. **Table S1** includes details of data domains within the CDW that were used for identifying study variables. Unique patient identifiers allow linkage of data across domains to create a comprehensive patient record.

Since Veterans may also receive care at a non-VA facility, relying on VA data alone may not provide complete information regarding clinical events after PAD diagnosis. Therefore, we obtained two additional sources of data to supplement clinical data available in the VHA. These include the a) Fee-basis files and Consolidated Dataset (CDS) files and b) Medicare claims. The Fee-basis (2019 and earlier) and CDS (2020 onwards) data include claims arising from services provided to Veterans at a non-VA facility (or by a non-VA provider) that are paid by the VA. Typically, the VA pays for emergent care received by a Veteran at a non-VA facility (e.g., hospitalization for acute myocardial infarction (MI) or stroke). The VA also pays for routine visits at a non-VA facility if travelling to a VA facility for the same services imposes undue burden (e.g, vascular surgery clinic at a non-VA facility near a Veteran’s home if comparable services at a VA facility require long travel). Each year the VA Information Resource Center (VIReC) obtains 100% Medicare enrollment and claims data from the Center for Medicare and Medicaid (CMS) for all Veterans. These include Part A (fee for service), Part B & Carrier (outpatient), Part C (Advantage), and Part D (medications). For our cohort, we obtained all Part A (fee-for-service), Part B & Carrier (outpatient facility & physician claims), Part C (Medicare Advantage) and Part D (medications) claims for our cohort during 2013-2022 from VIReC. Part C data were only available through 2021. Because >70% of Veterans in our cohort were older than 65 years of age (Medicare eligible) at the time of PAD diagnosis, the addition of Medicare claims ensured more comprehensive follow up from clinical encounters outside the VA especially those older than 65 years of age. Thus, the availability of these additional data sources ensures near complete follow up of our cohort of PAD patients.

### ABI Reports

Within the CDW data, we identified all reports of ABI tests performed during 2013-2020 at a VA facility. Depending on the facility, ABI tests maybe performed by vascular surgery, radiology or primary care. Therefore, ABI report documents can be found in either the Text Integration Utility (TIU) or Radiology domain within Unstructured Clinical data in CDW. In both domains, documents are organized into document types with each document type having a unique ID and an associated title (e.g., “Vascular Lab ABI report”). The document type IDs are unique to each VA facility and can be used to identify all individual documents within that facility (e.g., all ABI test results that are reported using the “Vascular Lab ABI report” title would have a common document type ID). Often, a facility uses more than one document type for reporting ABI results especially if ABI tests are performed by multiple specialties at that facility. We included all possible document type IDs used for reporting ABI results at each facility. Finally, in some instances, a facility used a ‘generic’ document type ID (e.g. “Vascular Lab Report”) for reporting both ABI studies and non-ABI vascular studies (e.g., carotid duplex, etc.). Such generic document type IDs were also included. In our validation study, the NLP system extracted no false positive ABI or TBI values from non-ABI vascular studies.21 A total of 424 document type IDs used for reporting ABI results in the VHA during 2013-2020 were identified which corresponded to 1,182,023 individual documents.

### NLP Deployment

Since the NLP system was developed using only a small random sample of individual ABI documents (800 total), certain adaptations were necessary to handle the large document corpus without losing accuracy. In our development work, we used the same text-based features for random forest model training and subsequent testing. However, this approach proved to be infeasible at scale due to temporal and spatial requirements, which requires re-training the random forest model for every document. To overcome this challenge, we first redesigned and refactored the code underlying the tool and re-tested model performance using the previous training and test corpora to ensure no change in performance. We also built a disk-based storage system that could hold millions of files and allow for rapid retrieval of documents without consuming excessive disk space. Finally, we created a database programming script that can identify PAD based on extracted ABI and TBI values (**Figure S1**). With these changes, we were able to process all of the documents in a total of 95 hours. This equated to 3.48 documents processed per second or over 300,000 documents per day.

### Creation of the PEARLS Cohort

Since our overall objective was to identify Veterans with a new diagnosis of PAD during 2015 to 2020, we intentionally added two additional years of ABI data (2013-2014) as a ‘look-back’ period. This was done to ensure that patients identified as PAD based on abnormal ABI or TBI for the first time in 2015-2020 did not also have abnormal ABI or TBI value during the prior 2 years.

Figure 1 shows the derivation of our study cohort. After processing ∼1.2 million documents and excluding non-ABI vascular studies, non-Veterans and patients missing unique IDs, a total of 336,245 patients with 617,978 ABI report documents remained. Among these patients, 161,070 (47.9%) met criteria for PAD (i.e., ABI value of <0.9 or a TBI value <0.7, **Figure S1**). For each patient, all ABI reports were arranged in chronological order, and the index date (date of PAD diagnosis) was defined as the earliest date of abnormal ABI or TBI value in any limb. From this sample, we excluded patients with evidence of PAD in the 2 years prior to 2015 or a history of lower extremity revascularization or amputation in the 2 years before the index date. Additional exclusions were non-regular use of VA services (less than 2 outpatient, or 1 inpatient and 1 outpatient visit in the prior 2 years), age less than 40 years on index date (as PAD is rare in those younger than 40 years), and death before the index date (data coding errors). Our final sample included 103,748 patients with a new diagnosis of PAD during 2015 to 2020.

### Study Outcomes

The registry collected a number of outcomes relevant to PAD. All-cause mortality was determined using the VA Death Ascertainment file. Other important outcomes that were collected included cardiovascular events (hospitalization for acute myocardial infarction [MI]) or ischemic stroke) and limb events (CLI and major amputation). Major amputation was defined as any amputation above the ankle joint. The diagnosis and procedure codes used to identify these events are listed in **Table S2.** Hospitalization for MI, stroke, and CLI were included only when they were reported as primary diagnosis on the inpatient encounter as multiple prior studies have demonstrated better accuracy only when the primary diagnosis is considered. 22–26

### Study Variables

Demographic characteristics included age, sex, race (White, Black, other), and ethnicity (Hispanic vs. non-Hispanic). Smoking status up to 12 months before the index date was determined using the health factors table in VHA CDW, which contains details on routine screenings and social factors that impact health. Co-morbidities were based on clinical diagnosis entered in inpatient and outpatient encounters from the VA, Fee-basis and CDS files during the 2 years prior to the index date, and defined using algorithms originally developed by Elixhauser, and updated by Quan.27 We used *Vital Signs* data from outpatient encounters to determine baseline systolic and diastolic blood pressure during 24 months before the index date. We used *Laboratory* data to determine baseline levels of total cholesterol, low-density lipoprotein (LDL), high-density lipoprotein (HDL), triglycerides and hemoglobin A1c levels within 12 months before the index date. We also included the ABI and TBI values extracted by the NLP system from the qualifying ABI report. We identified prescription fills using VA Pharmacy and Medicare Part D data to determine baseline use of lipid lowering drugs, antiplatelet drugs, anticoagulants, anti-hypertensive and anti-diabetic drugs. For each drug, we defined active use at baseline as filling a prescription for that drug prior to the index date within a time-period that was less than or equal to twice the days’ supply dispensed. For example, a PAD patient with an index date of August 1, 2015 who filled a 30-day supply of atorvastatin on June 15, 2015 (46 days before the index date) would be considered an active statin user as the drug was filled within 60 days (2×30 days’ supply).

### Statistical Analyses

Baseline characteristics of our study cohort were reported using means (standard deviations) and frequencies (proportions) as appropriate. Next, we examined the overall rate of each study endpoint: mortality, hospitalization for MI or stroke (cardiovascular events) and CLI and amputation (limb events) at 1-year follow-up. We also constructed survival curves as well as the cumulative index function for the secondary endpoints over a 1-year follow-up period. For non-fatal endpoints, the incidence rates were computed after accounting for the competing risk of death and expressed per 100 patient-years. Finally, for non-fatal endpoints, we also quantified the proportion of events that were identified using different data sources (VA data, Fee-basis, Medicare fee-for-service, and Medicare Advantage).All analyses were performed using SAS version 9.4 (SAS Institute, Cary, North Carolina).

## Results

**Table 1** shows baseline characteristics of our cohort. Among 103,748 eligible patients with PAD, the mean age was 70.6 years (standard deviation [SD] of 9.1 years), and the majority (97.3%) were men. Nearly 1 in 5 patients (18.5%, N=19,208) were of Black race and 3.0% (N=3,130) were of Hispanic ethnicity. The mean (SD) ABI and TBI in the worst limb was 0.78 (0.26) and 0.51 (0.19), respectively. Among 84% of patients with available data on smoking within 1 year of index data, 32.2% were currently smoking and 35.7% formerly smoked. Overall, there was a high prevalence of hypertension (86.6%), coronary artery disease (42.9%), cerebrovascular disease (23.8%), diabetes (54.8%), chronic kidney disease (23.6%) and chronic obstructive pulmonary disease (35.4%). At baseline, the mean systolic blood pressure was 133 mm Hg, and mean diastolic blood pressure was 71 mm Hg, and 58.4% (N=98,869 with available BP data) had a systolic blood pressure >130 mm Hg or diastolic blood pressure >80 mm Hg. The mean values of cholesterol, triglycerides and hemoglobin A1c are also included in **Table 1**.

View this table:
[Table 1.](http://medrxiv.org/content/early/2024/08/21/2024.08.20.24312328/T1)

Table 1. 
Baseline Characteristics of Veterans with New Diagnosis of PAD 2015-2020

**Table 2** shows baseline use of medications in our PAD cohort. Nearly two-thirds (65.4%) of PAD patients were on a statin at the time of PAD diagnosis, but only 34.2% were on a high-intensity statin. Among non-statin lipid-lowering medications, use of ezetimibe (1.1%) and Proprotein Convertase Subtilisin/Kexin-9 (PCSK-9) inhibitors (< 0.1%) at baseline was uncommon. Anti-plaletelet drugs were used in 35.0% (asprin: 24.5%, P2Y12 inhibitors: 13.3%, and phosphodiesterase inhibitors: 3.0%) while anticoagulants were used in 12.9% of patients (warfarin: 6.9%, novel oral anticoagulant 6.3%). Only 25 patients (<0.1%) were taking low dose rivaroxaban at baseline. Overall, 81.2% patients were receiving an anti-hypertensive with angiotensin converting enzyme-inhibitors, angiotensin receptor blockers and angiotensin receptor neprilysin inhibitor being the most common medications. Among the 41.3% of patients receiving anti-diabetic agents, metformin was the most common agent (24.0% patients) followed by sulfonylureas and insulin. Baseline use of glucagon-like peptide (GLP-1) agonists and sodium glucose transporter-2 (SGLT-2) inhibitors was low.

View this table:
[Table 2.](http://medrxiv.org/content/early/2024/08/21/2024.08.20.24312328/T2)

Table 2. 
Medications at Baseline in PAD Cohort

The number of clinical events and the incidence of study endpoints are reported in **Table 3** and **Figures 2 & Figure 3**, respectively. Overall, 9,723 participants (9.4%) died within 1-year of the index date. During 1-year follow-up, 2,706 (2.6%) experienced an acute MI and 2881 (2.8%) experienced a stroke, 2465 (2.4%) experienced a CLI hospitalization, and 2616 (2.5%) underwent a major amputation. After accounting for the competing risk of mortality, the 1-year incidence rate of cardiovascular events was 5.6 per 100 patient-years (Figure 3A), which included acute MI (2.8 per-100 patient-years) and stroke (3.0 per-100 patient-years). The 1-year incidence of limb events was 4.5 per-100 patient-years **(Figure 3B)** which comprised of hospitalization for CLI (2.5 per-100 patient-years) and major amputation (2.7 per-100 patient-years).

View this table:
[Table 3.](http://medrxiv.org/content/early/2024/08/21/2024.08.20.24312328/T3)

Table 3. 
The number of patients who experienced a study endpoint at 1-year

![Figure 2.](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/21/2024.08.20.24312328/F2.medium.gif)

[Figure 2.](http://medrxiv.org/content/early/2024/08/21/2024.08.20.24312328/F2)

Figure 2. Kaplan Meier Curve for Survival at 1-year.
The figure shows the Kaplan-Meier estimates of survival over 1-year follow up

![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/21/2024.08.20.24312328/F3/graphic-7.medium.gif)

[](http://medrxiv.org/content/early/2024/08/21/2024.08.20.24312328/F3/graphic-7)

![](http://medrxiv.org/https://www.medrxiv.org/content/medrxiv/early/2024/08/21/2024.08.20.24312328/F3/graphic-8.medium.gif)

[](http://medrxiv.org/content/early/2024/08/21/2024.08.20.24312328/F3/graphic-8)

Figure 3. Cumulative Incidence of Cardiovascular Events and Limb Events.
Panel A shows the cumulative incidence of cardiovascular events (hospitalization for acute myocardial infarction or stroke) through 1-year. Panel B shows the cumulative incidence of limb events (hospitalization for critical limb ischemia or amputation). The cumulative incidence function curves were estimated after accounting for competing risk of mortality.

**Table S3** shows the proportion of non-fatal events identified using different data sources. Of the 5329 patients with a cardiovascular event (acute MI or stroke), 1721 (31.9%) were identified in the VHA, 2071 (38.4%) in Fee-basis (non-VA hospital claims paid by the VHA), 1177 (21.8%) and 423 (7.8%) in Medicare fee-for-service and Advantage, respectively. Among 4307 patients with a limb event (CLI or major amputation), 2726 (63.3%) were in the VHA, 660 (15.3%) in Fee-basis, 644 (15.0%) in Medicare fee-for-service, and 277 (6.4%) in Medicare Advantage claims.

## Discussion

Using a previously validated NLP system,21 we have created a registry of >100,000 newly diagnosed patients with PAD in the VHA. In contrast to solely using administrative diagnosis and procedure codes that have modest accuracy,17–20 identification of PAD in our cohort is based on abnormal ABI and TBI values which are considered the gold standard methods for clinical diagnosis of PAD.28 Importantly, our cohort is racially and ethnically diverse with >19,000 individuals self-identified as Black, and 3,000 self-identified as Hispanic. The prevalence of smoking and other co-morbidities, including hypertension, coronary and cerebrovascular disease, diabetes, and chronic kidney disease was high. Rates of all-cause mortality approached 10% at 1-year, with high rates of cardiovascular and limb events, which underscore the enormous impact of PAD in real-world patients.

Prior studies of PAD in large health systems have been limited by the imprecision of using ICD codes for identifying PAD. These studies have reported the sensitivity of ICD codes to range from 34.7% to 76.9% and the positive predictive value to range from 27.5% to 69.4%.17–20 The low accuracy of ICD codes for PAD identification likely affects inferences drawn from such studies. To overcome this challenge, we developed an NLP system that can extract ABI and TBI values including laterality from ABI reports in the VHA.21 An algorithm based on NLP-extracted ABI and TBI values had a higher sensitivity of 83.3% and positive predictive value of 92.4% as compared ICD diagnosis and procedure codes for identifying PAD.21 Building on this prior work, we now report the findings following the implementation of our NLP system in national VA data to develop an inception cohort of newly diagnosed PAD in the VHA during 2015-2020. The PEARLS cohort includes >100,000 Veterans with a new diagnosis of PAD, making it one of the largest contemporary cohorts of PAD to our knowledge. Importantly, through linkage with an array of VA and non-VA data sources, rich information on key variables (e.g., smoking status, co-morbidities, vital signs, laboratory values and medications) as well as comprehensive follow-up is available. Thus, the PEARLS registry is well positioned to evaluate the long-term trajectory of health outcomes, evaluate the association of PAD treatments with long-term outcomes, and identify opportunities for improving care.

There are several advantages of studying PAD in the VHA. First, the VHA provides care to over 9 million Veterans across 130 VA medical centers and 1138 outpatient sites that use a standardized EHR.29 Risk factors for PAD which include smoking, hypertension, and diabetes are highly prevalent leading to a high burden of PAD in Veterans. Second, although women are under-represented, our cohort still includes ∼3,000 women representing one of the largest contemporary cohorts of women with PAD. Importantly, racial and ethnic subgroups are well-represented in our study with the inclusion of >19,000 Black individuals and ∼3000 Hispanic individuals, groups that have been under-represented in prior studies. Third, in addition to inpatient and outpatient claims, PEARLS also includes rich information on medications (VA *Pharmacy* and Part D claims), *Vital Signs* captured during clinical encounters, and *Laboratory* results that are generally unavailable in other cohorts. Most importantly, the availability of Fee-basis, CDS and Medicare files ensures comprehensive follow up data for defining clinical exposures (e.g., medications) and clinical endpoints from VA and non-VA encounters (except the small minority of Veterans insured with commercial insurance plans).30,31 This overcomes a key limitation of other EHR studies that are often limited in assessing events if they occurred in a different health system. Finally, VA data are updated nightly providing real-time access avoiding delay inherent in studies that use other insurance claims data.

The high incidence of clinical events in our cohort underscores the substantial impact of PAD on patients. Nearly 10% of patients with PAD died within 1 year of diagnosis. The incidence of cardiovascular events was 5.6% per 100 patient-years at 1-year accounting for competing risk of mortality, highlighting that PAD is a marker of malignant vascular phenotype with clinical event rates that far exceed those of patients with stable atherosclerotic cardiovascular disease.32 For instance, among 32,961 patients with stable CAD, the incidence of death or cardiovascular events was 7.6% at 5 years.32 It was also striking that the incidence of limb events -hospitalization for CLI and major amputation was high as well (4.5 per-100 patient-years) approaching that of cardiovascular events. The downstream impact of a limb event – risk of infection, limited mobility and disability – is enormous and contributes to high mortality and poor quality of life in PAD patients. Thus, our findings support guideline recommendations which consider abnormal ABI as a risk enhancer for primary prevention and the importance of developing strategies for not only evaluating therapies for reducing cardiovascular but also limb events.33

The use of an accurate method of identifying PAD patients combined with rich and comprehensive data provides a unique opportunity to enhance contemporary PAD research in important ways. Ongoing studies from our cohort will yield important insights regarding patterns of care delivery for PAD patients (e.g., medications, risk factor control and revascularization) and the extent to which these exposures are associated with long-term outcomes. Moreover, given the sizeable number of Black and Hispanic individuals in our cohort, we will also be able to evaluate whether disparities exist in care and outcomes of PAD patients by race and ethnicity. Because health insurance is not a barrier in the VHA, we will also have a unique opportunity to isolate the contribution of access to healthcare from other social determinants of health in explaining disparities. Finally, the large number of VA sites in our cohort provides an opportunity to understand variation in care patterns across sites with the goal of identifying best practices.

Our study findings should be considered in the context of the following limitations. First, ABIs can be normal in patients with non-compressible vessels. Although we used TBI values to identify PAD in patients with falsely normal ABI values, TBI measurements are not always performed. Thus, patients with non-compressible vessels who did not have TBI measured are not included. Second, data on clinical symptoms and severity (e.g., claudication, rest pain) are also not included. Since routine ABI screening of asymptomatic patients is not universally recommended, we expect that most patients in our cohort have symptomatic PAD. Third, we identified clinical endpoints using a combination of VA claims, non-VA claims paid by the VA or Medicare claims to ensure near complete follow up. However, it is possible that the incidence of clinical endpoints in younger Veterans not insured by Medicare or those with private insurance is underestimated. Fourth, medication data were based on prescriptions that were either filled at a VA pharmacy or using Medicare Part D insurance. 30,31 However, the VA provides generous pharmacy benefits for Veterans which include a) modest co-payment not to exceed $11 for all drugs including brand name medications, b) no co-payments for disabled Veterans with service connection and c) the convenience of mail-order and online pharmacy. These incentives likely minimize the potential for missing information on medication exposures in our cohort.34,35 Finally, due to its availability over the counter, use of aspirin may be underreported in our study.

## Conclusion

We have successfully developed a registry of >100,000 patients with a new diagnosis of PAD in the Veterans Health Administration using a novel NLP algorithm. Ongoing studies from our cohort will yield important insights into the association between management of PAD and clinical outcomes, with the goal of identifying opportunities for improving care.

**Funding**: The study is funded by the National Heart, Lung and Blood Institute (NHLBI, R01HL166305 and R56HL158803; PI: Girotra; Co-I: Vaughan-Sarrazin, Lund, Al-Garadi, Beckman, Hoffman, Tsai, Smolderen, Arya, Matheny, Gobbel). Dr. Lund receives funding from the VA Office of Research and Development and the VA Office of Rural Health. Dr. Chan and Dr. Girotra also receive receive funding from the NHLBI (R01HL160734) and the American Heart Association. This material is the result of work supported with resources and the use of facilities at the Iowa City, Nashville and North Texas Veterans Affairs. None of the funders had any role in the design and conduct of the study; collection, management, analysis, and interpretation of the data; preparation for publication. The views expressed in this article are those of the authors and do not necessarily reflect the position or policy of the Department of Veterans Affairs or the United States Government.

## DISCLOSURES

Dr. Smolderen receives consulting fees from Terumo, Cook Medical, Happify and grants from Johnson & Johnson, Merck, Abbott, and Shockwave

## Data Availability

Due to the proprietary nature of the data used in the study, data cannot be made available

*   Received August 20, 2024.
*   Revision received August 20, 2024.
*   Accepted August 21, 2024.


*   © 2024, Posted by Cold Spring Harbor Laboratory

This pre-print is available under a Creative Commons License (Attribution-NonCommercial 4.0 International), CC BY-NC 4.0, as described at [http://creativecommons.org/licenses/by-nc/4.0/](http://creativecommons.org/licenses/by-nc/4.0/)

## REFERENCES

1.  1.Fowkes FG, Rudan D, Rudan I, Aboyans V, Denenberg JO, McDermott MM, Norman PE, Sampson UK, Williams LJ, Mensah GA, et al. Comparison of global estimates of prevalence and risk factors for peripheral artery disease in 2000 and 2010: a systematic review and analysis. Lancet. 2013;382:1329–1340. doi: 10.1016/S0140-6736(13)61249-0
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S0140-6736(13)61249-0&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=23915883&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F21%2F2024.08.20.24312328.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000325975700023&link_type=ISI) 

2.  2.Song P, Rudan D, Zhu Y, Fowkes FJI, Rahimi K, Fowkes FGR, Rudan I. Global, regional, and national prevalence and risk factors for peripheral artery disease in 2015: an updated systematic review and analysis. Lancet Glob Health. 2019;7:e1020–e1030. doi: 10.1016/S2214-109X(19)30255-4
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/S2214-109X(19)30255-4&link_type=DOI) 

3.  3.Mahoney EM, Wang K, Cohen DJ, Hirsch AT, Alberts MJ, Eagle K, Mosse F, Jackson JD, Steg PG, Bhatt DL, et al. One-year costs in patients with a history of or at risk for atherothrombosis in the United States. Circulation Cardiovascular quality and outcomes. 2008;1:38–45. doi: 10.1161/CIRCOUTCOMES.108.775247
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiY2lyY2N2b3EiO3M6NToicmVzaWQiO3M6NjoiMS8xLzM4IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDgvMjEvMjAyNC4wOC4yMC4yNDMxMjMyOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

4.  4.Mahoney EM, Wang K, Keo HH, Duval S, Smolderen KG, Cohen DJ, Steg G, Bhatt DL, Hirsch AT, Reduction of Atherothrombosis for Continued Health Registry I. Vascular hospitalization rates and costs in patients with peripheral artery disease in the United States. Circulation Cardiovascular quality and outcomes. 2010;3:642–651. doi: 10.1161/CIRCOUTCOMES.109.930735
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6ODoiY2lyY2N2b3EiO3M6NToicmVzaWQiO3M6NzoiMy82LzY0MiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA4LzIxLzIwMjQuMDguMjAuMjQzMTIzMjguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

5.  5.Franklin H, Rajan M, Tseng CL, Pogach L, Sinha A, Mph M. Cost of lower-limb amputation in U.S. veterans with diabetes using health services data in fiscal years 2004 and 2010. J Rehabil Res Dev. 2014;51:1325–1330. doi: 10.1682/JRRD.2013.11.0249
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1682/JRRD.2013.11.0249&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25625913&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F21%2F2024.08.20.24312328.atom) 

6.  6.Hirsch AT, Hartman L, Town RJ, Virnig BA. National health care costs of peripheral arterial disease in the Medicare population. Vasc Med. 2008;13:209–215. doi: 10.1177/1358863X08089277
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1358863X08089277&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18687757&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F21%2F2024.08.20.24312328.atom) 

7.  7.Arya S, Binney Z, Khakharia A, Brewster LP, Goodney P, Patzer R, Hockenberry J, Wilson PWF. Race and Socioeconomic Status Independently Affect Risk of Major Amputation in Peripheral Artery Disease. J Am Heart Assoc. 2018;7. doi: 10.1161/JAHA.117.007425
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NToiYWhhb2EiO3M6NToicmVzaWQiO3M6MTE6IjcvMi9lMDA3NDI1IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDgvMjEvMjAyNC4wOC4yMC4yNDMxMjMyOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

8.  8.Aday AW, Duncan MS, Patterson OV, DuVall SL, Alba PR, Alcorn CW, Tindle HA, Creager MA, Bonaca MP, Damrauer SM, et al. Association of Sex and Race With Incident Peripheral Artery Disease Among Veterans With Normal Ankle-Brachial Indices. JAMA Netw Open. 2022;5:e2240188. doi: 10.1001/jamanetworkopen.2022.40188
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamanetworkopen.2022.40188&link_type=DOI) 

9.  9.Gerhard-Herman MD, Gornik HL, Barrett C, Barshes NR, Corriere MA, Drachman DE, Fleisher LA, Fowkes FG, Hamburg NM, Kinlay S, et al. 2016 AHA/ACC Guideline on the Management of Patients With Lower Extremity Peripheral Artery Disease: A Report of the American College of Cardiology/American Heart Association Task Force on Clinical Practice Guidelines. Circulation. 2017;135:e726–e779. doi: 10.1161/CIR.0000000000000471
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiRlVMTCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjExOiIxMzUvMTIvZTcyNiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA4LzIxLzIwMjQuMDguMjAuMjQzMTIzMjguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

10. 10.Willey J, Mentias A, Vaughan-Sarrazin M, McCoy K, Rosenthal G, Girotra S. Epidemiology of lower extremity peripheral artery disease in veterans. Journal of vascular surgery. 2018;68:527–535 e525. doi: 10.1016/j.jvs.2017.11.083
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jvs.2017.11.083&link_type=DOI) 

11. 11.Mentias A, Qazi A, McCoy K, Wallace R, Vaughan-Sarrazin M, Girotra S. Trends in Hospitalization, Management, and Clinical Outcomes Among Veterans With Critical Limb Ischemia. Circ Cardiovasc Interv. 2020;13:e008597. doi: 10.1161/CIRCINTERVENTIONS.119.008597
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1161/CIRCINTERVENTIONS.119.008597&link_type=DOI) 

12. 12.Arya S, Khakharia A, Binney ZO, DeMartino RR, Brewster LP, Goodney PP, Wilson PWF. Association of Statin Dose With Amputation and Survival in Patients With Peripheral Artery Disease. Circulation. 2018;137:1435–1446. doi: 10.1161/CIRCULATIONAHA.117.032361
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjExOiIxMzcvMTQvMTQzNSI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA4LzIxLzIwMjQuMDguMjAuMjQzMTIzMjguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

13. 13.Pande RL, Perlstein TS, Beckman JA, Creager MA. Secondary prevention and mortality in peripheral artery disease: National Health and Nutrition Examination Study, 1999 to 2004. Circulation. 2011;124:17–23. doi: 10.1161/CIRCULATIONAHA.110.003954
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6MTQ6ImNpcmN1bGF0aW9uYWhhIjtzOjU6InJlc2lkIjtzOjg6IjEyNC8xLzE3IjtzOjQ6ImF0b20iO3M6NTA6Ii9tZWRyeGl2L2Vhcmx5LzIwMjQvMDgvMjEvMjAyNC4wOC4yMC4yNDMxMjMyOC5hdG9tIjt9czo4OiJmcmFnbWVudCI7czowOiIiO30=) 

14. 14.Patel KK, Jones PG, Ellerbeck EF, Buchanan DM, Chan PS, Pacheco CM, Moneta G, Spertus JA, Smolderen KG. Underutilization of Evidence-Based Smoking Cessation Support Strategies Despite High Smoking Addiction Burden in Peripheral Artery Disease Specialty Care: Insights from the International PORTRAIT Registry. J Am Heart Assoc. 2018;7:e010076. doi: 10.1161/JAHA.118.010076
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1161/JAHA.118.010076&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F21%2F2024.08.20.24312328.atom) 

15. 15.Saxon JT, Safley DM, Mena-Hurtado C, Heyligers J, Fitridge R, Shishehbor M, Spertus JA, Gosch K, Patel MR, Smolderen KG. Adherence to Guideline-Recommended Therapy-Including Supervised Exercise Therapy Referral-Across Peripheral Artery Disease Specialty Clinics: Insights From the International PORTRAIT Registry. J Am Heart Assoc. 2020;9:e012541. doi: 10.1161/JAHA.119.012541
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1161/JAHA.119.012541&link_type=DOI) 

16. 16.Kumbhani DJ, Steg PG, Cannon CP, Eagle KA, Smith SC, Jr., Goto S, Ohman EM, Elbez Y, Sritara P, Baumgartner I, et al. Statin therapy and long-term adverse limb outcomes in patients with peripheral artery disease: insights from the REACH registry. Eur Heart J. 2014;35:2864–2872. doi: 10.1093/eurheartj/ehu080
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1093/eurheartj/ehu080&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24585266&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F21%2F2024.08.20.24312328.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000345367100006&link_type=ISI) 

17. 17.Fan J, Arruda-Olson AM, Leibson CL, Smith C, Liu G, Bailey KR, Kullo IJ. Billing code algorithms to identify cases of peripheral artery disease from administrative data. Journal of the American Medical Informatics Association : JAMIA. 2013;20:e349–354. doi: 10.1136/amiajnl-2013-001827
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1136/amiajnl-2013-001827&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=24166724&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F21%2F2024.08.20.24312328.atom) 

18. 18.Hong Y, Sebastianski M, Makowsky M, Tsuyuki R, McMurtry MS. Administrative data are not sensitive for the detection of peripheral artery disease in the community. Vasc Med. 2016;21:331–336. doi: 10.1177/1358863X16631041
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1358863X16631041&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=27114456&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F21%2F2024.08.20.24312328.atom) 

19. 19.Lasota AN, Overvad K, Eriksen HH, Tjonneland A, Schmidt EB, Gronholdt MM. Validity of Peripheral Arterial Disease Diagnoses in the Danish National Patient Registry. Eur J Vasc Endovasc Surg. 2017;53:679–685. doi: 10.1016/j.ejvs.2016.12.031
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.ejvs.2016.12.031&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F21%2F2024.08.20.24312328.atom) 

20. 20.Weissler EH, Zhang J, Lippmann S, Rusincovitch S, Henao R, Jones WS. Use of Natural Language Processing to Improve Identification of Patients With Peripheral Artery Disease. Circ Cardiovasc Interv. 2020;13:e009447. doi: 10.1161/CIRCINTERVENTIONS.120.009447
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1161/CIRCINTERVENTIONS.120.009447&link_type=DOI) 

21. 21.Friberg JE, Qazi AH, Boyle B, Franciscus C, Vaughan-Sarrazin M, Westerman D, Patterson OV, Parr SK, Matheny ME, Arya S, et al. Ankle- and Toe-Brachial Index for Peripheral Artery Disease Identification: Unlocking Clinical Data Through Novel Methods. Circ Cardiovasc Interv. 2022;15:e011092. doi: 10.1161/CIRCINTERVENTIONS.121.011092
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1161/CIRCINTERVENTIONS.121.011092&link_type=DOI) 

22. 22.Petersen LA, Wright S, Normand SL, Daley J. Positive predictive value of the diagnosis of acute myocardial infarction in an administrative database. J Gen Intern Med. 1999;14:555–558. doi: 10.1046/j.1525-1497.1999.10198.x
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1046/j.1525-1497.1999.10198.x&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=10491245&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F21%2F2024.08.20.24312328.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000082595600005&link_type=ISI) 

23. 23.Ammann EM, Schweizer ML, Robinson JG, Eschol JO, Kafa R, Girotra S, Winiecki SK, Fuller CC, Carnahan RM, Leonard CE, et al. Chart validation of inpatient ICD-9-CM administrative diagnosis codes for acute myocardial infarction (AMI) among intravenous immune globulin (IGIV) users in the Sentinel Distributed Database. Pharmacoepidemiol Drug Saf. 2018;27:398–404. doi: 10.1002/pds.4398
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1002/pds.4398&link_type=DOI) 

24. 24.McCormick N, Bhole V, Lacaille D, Avina-Zubieta JA. Validity of Diagnostic Codes for Acute Stroke in Administrative Databases: A Systematic Review. PLoS One. 2015;10:e0135834. doi: 10.1371/journal.pone.0135834
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1371/journal.pone.0135834&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=26292280&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F21%2F2024.08.20.24312328.atom) 

25. 25.Bekwelem W, Bengtson LG, Oldenburg NC, Winden TJ, Keo HH, Hirsch AT, Duval S. Development of administrative data algorithms to identify patients with critical limb ischemia. Vasc Med. 2014;19:483–490. doi: 10.1177/1358863X14559589
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1177/1358863X14559589&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=25447239&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F21%2F2024.08.20.24312328.atom) 

26. 26.Rayman G, Krishnan ST, Baker NR, Wareham AM, Rayman A. Are we underestimating diabetes-related lower-extremity amputation rates? Results and benefits of the first prospective study. Diabetes Care. 2004;27:1892–1896. doi: 10.2337/diacare.27.8.1892
    
    [Abstract/FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6NDoiQUJTVCI7czoxMToiam91cm5hbENvZGUiO3M6NzoiZGlhY2FyZSI7czo1OiJyZXNpZCI7czo5OiIyNy84LzE4OTIiO3M6NDoiYXRvbSI7czo1MDoiL21lZHJ4aXYvZWFybHkvMjAyNC8wOC8yMS8yMDI0LjA4LjIwLjI0MzEyMzI4LmF0b20iO31zOjg6ImZyYWdtZW50IjtzOjA6IiI7fQ==) 

27. 27.Quan H, Sundararajan V, Halfon P, Fong A, Burnand B, Luthi JC, Saunders LD, Beck CA, Feasby TE, Ghali WA. Coding algorithms for defining comorbidities in ICD-9-CM and ICD-10 administrative data. Med Care. 2005;43:1130–1139.
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1097/01.mlr.0000182534.19832.83&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=16224307&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F21%2F2024.08.20.24312328.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000233268500010&link_type=ISI) 

28. 28.Aboyans V, Ho E, Denenberg JO, Ho LA, Natarajan L, Criqui MH. The association between elevated ankle systolic pressures and peripheral occlusive arterial disease in diabetic and nondiabetic subjects. Journal of vascular surgery. 2008;48:1197–1203. doi: 10.1016/j.jvs.2008.06.005
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1016/j.jvs.2008.06.005&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=18692981&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F21%2F2024.08.20.24312328.atom) 
    
    [Web of Science](http://medrxiv.org/lookup/external-ref?access_num=000260725300018&link_type=ISI) 

29. 29.Veterans Health Administration. [https://www.vagov/health/aboutvhaasp#:~:text=The%20Veterans%20Health%20Administration%20(VHA,the%20VA%20health%20care%20program](https://www.vagov/health/aboutvhaasp#:~:text=The%20Veterans%20Health%20Administration%20(VHA,the%20VA%20health%20care%20program) *(Accessed July 11, 2024)*.
    
    

30. 30.Barnett PG, Hong JS, Carey E, Grunwald GK, Joynt Maddox K, Maddox TM. Comparison of Accessibility, Cost, and Quality of Elective Coronary Revascularization Between Veterans Affairs and Community Care Hospitals. JAMA Cardiol. 2018;3:133–141. doi: 10.1001/jamacardio.2017.4843
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1001/jamacardio.2017.4843&link_type=DOI) 

31. 31.Waldo SW, Glorioso TJ, Baron AE, Plomondon ME, Valle JA, Schofield R, Ho PM. Outcomes Among Patients Undergoing Elective Percutaneous Coronary Intervention at Veterans Affairs and Community Care Hospitals. Journal of the American College of Cardiology. 2020;76:1112–1116. doi: 10.1016/j.jacc.2020.05.086
    
    [FREE Full Text](http://medrxiv.org/lookup/ijlink/YTozOntzOjQ6InBhdGgiO3M6MTQ6Ii9sb29rdXAvaWpsaW5rIjtzOjU6InF1ZXJ5IjthOjQ6e3M6ODoibGlua1R5cGUiO3M6MzoiUERGIjtzOjExOiJqb3VybmFsQ29kZSI7czo0OiJhY2NqIjtzOjU6InJlc2lkIjtzOjk6Ijc2LzkvMTExMiI7czo0OiJhdG9tIjtzOjUwOiIvbWVkcnhpdi9lYXJseS8yMDI0LzA4LzIxLzIwMjQuMDguMjAuMjQzMTIzMjguYXRvbSI7fXM6ODoiZnJhZ21lbnQiO3M6MDoiIjt9) 

32. 32.Mesnier J, Ducrocq G, Danchin N, Ferrari R, Ford I, Tardif JC, Tendera M, Fox KM, Steg PG, Investigators C. International Observational Analysis of Evolution and Outcomes of Chronic Stable Angina: The Multinational CLARIFY Study. Circulation. 2021;144:512–523. doi: 10.1161/CIRCULATIONAHA.121.054567
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1161/CIRCULATIONAHA.121.054567&link_type=DOI) 

33. 33.Arnett DK, Blumenthal RS, Albert MA, Buroker AB, Goldberger ZD, Hahn EJ, Himmelfarb CD, Khera A, Lloyd-Jones D, McEvoy JW, et al. 2019 ACC/AHA Guideline on the Primary Prevention of Cardiovascular Disease: Executive Summary: A Report of the American College of Cardiology/American Heart Association Task Force on Clinical Practice Guidelines. Circulation. 2019;140:e563–e595. doi: 10.1161/CIR.0000000000000677
    
    [CrossRef](http://medrxiv.org/lookup/external-ref?access_num=10.1161/CIR.0000000000000677&link_type=DOI) 
    
    [PubMed](http://medrxiv.org/lookup/external-ref?access_num=http://www.n&link_type=MED&atom=%2Fmedrxiv%2Fearly%2F2024%2F08%2F21%2F2024.08.20.24312328.atom) 

34. 34.Department of Veterans Affairs Health Benefits [https://www.va.gov/health-care/refill-track-prescriptions/](https://www.va.gov/health-care/refill-track-prescriptions/) (Accessed October 2022).
    
    

35. 35.VA Prescription Refill and Tracking [https://www.va.gov/health-care/refill-track-prescriptions/](https://www.va.gov/health-care/refill-track-prescriptions/) (Accessed October 2022).