Advertisement
AJM online Clinical research study| Volume 129, ISSUE 9, P1001.e9-1001.e18, September 2016

Validation of Risk Assessment Models of Venous Thromboembolism in Hospitalized Medical Patients

      Abstract

      Background

      Patients hospitalized for acute medical illness are at increased risk for venous thromboembolism. Although risk assessment is recommended and several at-admission risk assessment models have been developed, these have not been adequately derived or externally validated. Therefore, an optimal approach to evaluate venous thromboembolism risk in medical patients is not known.

      Methods

      We conducted an external validation study of existing venous thromboembolism risk assessment models using data collected on 63,548 hospitalized medical patients as part of the Michigan Hospital Medicine Safety (HMS) Consortium. For each patient, cumulative venous thromboembolism risk scores and risk categories were calculated. Cox regression models were used to quantify the association between venous thromboembolism events and assigned risk categories. Model discrimination was assessed using Harrell's C-index.

      Results

      Venous thromboembolism incidence in hospitalized medical patients is low (1%). Although existing risk assessment models demonstrate good calibration (hazard ratios for “at-risk” range 2.97-3.59), model discrimination is generally poor for all risk assessment models (C-index range 0.58-0.64).

      Conclusions

      The performance of several existing risk assessment models for predicting venous thromboembolism among acutely ill, hospitalized medical patients at admission is limited. Given the low venous thromboembolism incidence in this nonsurgical patient population, careful consideration of how best to utilize existing venous thromboembolism risk assessment models is necessary, and further development and validation of novel venous thromboembolism risk assessment models for this patient population may be warranted.

      Keywords

      Clinical Significance
      • Several venous thromboembolism (VTE) risk assessment models (RAMs) have been developed, but lack external validation in hospitalized medical patients.
      • VTE incidence in medical patients is low (1%).
      • Results from an external validation study in a cohort of over 60,000 medical patients indicated poor model discrimination for VTE RAMs assessed.
      • Existing VTE RAMs have limited utility in identifying the highest-risk subset of medical patients for whom pharmacologic prophylaxis is warranted.
      Identifying medical patients at increased risk for venous thromboembolism using an individualized approach is of increasing importance and is emphasized in available guidelines.
      • Kahn S.R.
      • Lim W.
      • Dunn A.S.
      • et al.
      Prevention of VTE in nonsurgical patients: Antithrombotic Therapy and Prevention of Thrombosis, 9th ed: American College of Chest Physicians Evidence-Based Clinical Practice Guidelines.
      While several risk assessment models exist for venous thromboembolism in medical patients,
      • Barbar S.
      • Noventa F.
      • Rossetto V.
      • et al.
      A risk assessment model for the identification of hospitalized medical patients at risk for venous thromboembolism: the Padua Prediction Score.
      • Kucher N.
      • Koo S.
      • Quiroz R.
      • et al.
      Electronic alerts to prevent venous thromboembolism among hospitalized patients.
      • Spyropoulos A.C.
      • Anderson Jr., F.A.
      • Fitzgerald G.
      • et al.
      Predictive and associative models to identify hospitalized medical patients at risk for VTE.
      • Woller S.C.
      • Stevens S.M.
      • Jones J.P.
      • et al.
      Derivation and validation of a simple model to identify venous thromboembolism risk in medical patients.
      published risk assessment models have limited generalizability and validation.
      • Huang W.
      • Anderson F.A.
      • Spencer F.A.
      • Gallus A.
      • Goldberg R.J.
      Risk-assessment models for predicting venous thromboembolism among hospitalized non-surgical patients: a systematic review.
      Although a few recent studies have validated existing risk assessment models,
      • Woller S.C.
      • Stevens S.M.
      • Jones J.P.
      • et al.
      Derivation and validation of a simple model to identify venous thromboembolism risk in medical patients.
      • Mahan C.E.
      • Liu Y.
      • Turpie A.G.
      • et al.
      External validation of a risk assessment model for venous thromboembolism in the hospitalised acutely-ill medical patient (VTE-VALOURR).
      • Rosenberg D.
      • Eichorn A.
      • Alarcon M.
      • McCullagh L.
      • McGinn T.
      • Spyropoulos A.C.
      External validation of the risk assessment model of the International Medical Prevention Registry on Venous Thromboembolism (IMPROVE) for medical patients in a tertiary health system.
      these validation efforts lacked external validation,
      • Woller S.C.
      • Stevens S.M.
      • Jones J.P.
      • et al.
      Derivation and validation of a simple model to identify venous thromboembolism risk in medical patients.
      have had limited external validation,
      • Nendaz M.
      • Spirk D.
      • Kucher N.
      • et al.
      Multicentre validation of the Geneva Risk Score for hospitalised medical patients at risk of venous thromboembolism. Explicit ASsessment of Thromboembolic RIsk and Prophylaxis for Medical PATients in SwitzErland (ESTIMATE).
      or have had large-scale external validation with a risk assessment model that has limitations in predicting at-admission venous thromboembolism risk.
      • Mahan C.E.
      • Liu Y.
      • Turpie A.G.
      • et al.
      External validation of a risk assessment model for venous thromboembolism in the hospitalised acutely-ill medical patient (VTE-VALOURR).
      • Rosenberg D.
      • Eichorn A.
      • Alarcon M.
      • McCullagh L.
      • McGinn T.
      • Spyropoulos A.C.
      External validation of the risk assessment model of the International Medical Prevention Registry on Venous Thromboembolism (IMPROVE) for medical patients in a tertiary health system.
      Additional large-scale, external validation studies are important to help confirm or refute the accuracy of available risk assessment models, especially during hospital admission and based on detailed clinical data, as the performance of existing risk assessment models has been moderate at best in this setting.
      The Michigan Hospital Medicine Safety Consortium (HMS) is a state-wide quality collaborative focused on preventing adverse events in hospitalized medical patients. The consortium collects detailed patient-level data on venous thromboembolism risk factors and outcomes. The aim of the present study was to externally validate several existing risk assessment models using the large HMS cohort to determine which risk assessment model optimally predicts venous thromboembolism in acutely ill, hospitalized, medical patients.

      Methods

      Study Setting and Participants

      The setting and design of HMS have been previously described.
      • Flanders S.A.
      • Greene M.T.
      • Grant P.
      • et al.
      Hospital performance for pharmacologic venous thromboembolism prophylaxis and rate of venous thromboembolism: a cohort study.
      Although participation is voluntary, each hospital receives payments for participating in the consortium and for data collection.
      Eligible patients included those admitted to a medicine service for 2 days or longer. Patients were excluded if they met any of the following criteria: 1) under the age of 18 years; 2) pregnant; 3) underwent any surgical procedure during the admission; 4) direct admission to an intensive care unit (ICU); 5) direct admission for end-of-life care; 6) diagnosis of venous thromboembolism in the 6 months prior to admission; 7) admitted for presumed venous thromboembolism; 8) admitted under observation status; 9) re-admitted within 90 days of discharge from an admission included in the registry; or 10) received systemic anticoagulation on day 1 or day 2 of the index hospitalization.
      Detailed patient demographic, medical history, predefined risk factors for venous thromboembolism, and laboratory and medication data were collected through a standardized process using a trained medical record abstractor at each hospital. Patients discharged from each participating hospital were sampled on an 8-day rolling cycle to avert bias in selecting cases for review.
      American College of Surgeons
      American College of Surgeons National Surgical Quality Improvement Program (ACS NSQIP®).
      Data on the first 18 eligible cases discharged during each cycle were collected. Follow-up data were collected by medical record review and direct telephone follow-up at 90 days post-hospital discharge. Each hospital is audited on an annual basis by data quality coordinators to ensure completeness and accuracy of data abstraction. The University of Michigan is the HMS coordinating center.

      Ascertainment of Outcomes

      The primary outcome of interest was clinically diagnosed, image-confirmed hospital-associated venous thromboembolism, including proximal upper- or proximal lower-extremity deep vein thrombosis and pulmonary embolism. In order to be considered hospital acquired, venous thromboembolism events must have occurred on the third day after admission or later during an index hospitalization. The diagnosis of deep vein thrombosis was based on positive findings via compression Doppler ultrasound or venography, whereas pulmonary embolism was confirmed via computed tomography scan, ventilation perfusion scan, or pulmonary angiography. Venous thromboembolism outcomes were assessed out to 90 days from the date of index hospital admission. Patients transferred to an ICU or palliative care, and those who died during follow-up were censored; however, venous thromboembolism events that contributed to death or were the reason for transfer to the ICU were included. Patients who were alive and free of venous thromboembolism occurrence at 90 days following admission were right-censored. Telephone follow-up at 90 days was completed for 58% of patients. Medical record review at 90 days was completed for 100% of eligible patients.

      Statistical Analysis

      External Validation of Existing Risk Assessment Models

      All eligible patients were included as an external validation sample for existing risk assessment models. Numerous risk factors thought to increase risk for venous thromboembolism as specified by each risk assessment model were assessed for all patients. Bivariable Cox regression was used to assess the independent associations between putative risk factors and 90-day venous thromboembolism. Cumulative risk scores based on the presence of individual risk factor associated weights for the Kucher,
      • Kucher N.
      • Koo S.
      • Quiroz R.
      • et al.
      Electronic alerts to prevent venous thromboembolism among hospitalized patients.
      Padua,
      • Barbar S.
      • Noventa F.
      • Rossetto V.
      • et al.
      A risk assessment model for the identification of hospitalized medical patients at risk for venous thromboembolism: the Padua Prediction Score.
      predictive IMPROVE,
      • Spyropoulos A.C.
      • Anderson Jr., F.A.
      • Fitzgerald G.
      • et al.
      Predictive and associative models to identify hospitalized medical patients at risk for VTE.
      and Intermountain
      • Woller S.C.
      • Stevens S.M.
      • Jones J.P.
      • et al.
      Derivation and validation of a simple model to identify venous thromboembolism risk in medical patients.
      risk assessment models were calculated for each patient. In this manner, patients were assigned an “at-risk” status for each risk assessment model based on established cut points.
      • Barbar S.
      • Noventa F.
      • Rossetto V.
      • et al.
      A risk assessment model for the identification of hospitalized medical patients at risk for venous thromboembolism: the Padua Prediction Score.
      • Kucher N.
      • Koo S.
      • Quiroz R.
      • et al.
      Electronic alerts to prevent venous thromboembolism among hospitalized patients.
      • Spyropoulos A.C.
      • Anderson Jr., F.A.
      • Fitzgerald G.
      • et al.
      Predictive and associative models to identify hospitalized medical patients at risk for VTE.
      • Woller S.C.
      • Stevens S.M.
      • Jones J.P.
      • et al.
      Derivation and validation of a simple model to identify venous thromboembolism risk in medical patients.
      Cox regression models with gamma shared frailty by hospital were used to determine the hazards of developing venous thromboembolism for patients determined to be “at-risk” for each risk assessment model. Model discrimination was assessed via Harrell's C-index with 95% confidence intervals using the somersd package in Stata v.13 (StataCorp LP, College Station, Texas).
      • Harrell Jr., F.E.
      • Lee K.L.
      • Mark D.B.
      Multivariable prognostic models: issues in developing models, evaluating assumptions and adequacy, and measuring and reducing errors.
      • Newson R.
      Confidence intervals for rank statistics: Somers' D and extensions.
      Calibration was assessed by comparing estimated incidence rates for both “low-risk” and “at-risk” groups for each of the existing risk assessment models. Event rates and hazard ratios for each of the represented scores across the scales of each risk assessment model were also investigated to assess calibration.

      Sensitivity Analysis

      Given that the pathophysiology of upper extremity deep vein thrombosis is unique from those of pulmonary embolism and lower-extremity deep vein thrombosis,
      • Spyropoulos A.C.
      Upper vs. lower extremity deep vein thrombosis: outcome definitions of venous thromboembolism for clinical predictor rules or risk factor analyses in hospitalized patients.
      we conducted sensitivity analyses for all of the validation analyses described above, only considering pulmonary embolism or lower-extremity deep vein thrombosis as the venous thromboembolism outcome definition.

      Ethical and Regulatory Oversight

      As the purpose of the consortium is to measure and improve the quality of existing patient care practices, this project received a “not regulated” status by the University of Michigan Medical School's Institutional Review Board.

      Results

      From January 2011 through March 2014, data on 63,548 eligible patients were collected from 48 Michigan hospitals. A total of 38,723 (60.9%) patients received pharmacologic venous thromboembolism prophylaxis on admission. The patients' average age was 65.8 years, and 35,264 (55.5%) were female. The average length of stay was 4.5 days (median 4 days). The distribution of select putative risk factors for venous thromboembolism is displayed in Table 1. Factors most strongly associated with venous thromboembolism risk included recent cancer diagnosis or treatment within the last year, central venous catheter present on admission, history of prior venous thromboembolism, family history of venous thromboembolism, and surgery within the 30 days prior to admission (P <.001 for all comparisons). The Figure illustrates the risk factors and respective weights used for scaling and the percentage of patients determined to be “at-risk” for each of the risk assessment models. The percentage of “at-risk” patients was <20% for existing risk assessment models assessed.
      Table 1Patient Characteristics Associated with 90-day Venous Thromboembolism Risk, Unadjusted (N = 63,548)
      Risk Factorn (%)HR (CI)P-Value
      Cancer within last year
      Cancer diagnosis or acute cancer treatment within the last year.
      4990 (7.85)3.70 (3.08-4.45)<.001
      Central venous catheter present on admission
      Central venous catheters were considered present on admission if patient was admitted with device or if one was inserted on day 1 or day 2 of the admission.
      5011 (7.89)3.61 (3.01-4.33)<.001
      Prior venous thromboembolism4069 (6.40)2.96 (2.41-3.63)<.001
      Family history of venous thromboembolism441 (0.69)2.60 (1.47-4.60).001
      Postpartum (<1 mo)41 (0.06)2.24 (0.31-15.90).42
      Surgery (<1 mo)1695 (2.67)1.93 (1.36-2.74)<.001
      Leg edema (current)15,707 (24.72)1.76 (1.50-2.06)<.001
      Immobile/not ambulating
      Immobile/not ambulating defined as having at least one of the following: immobilizing plaster cast, paralysis, or bed rest for ≥72 hours prior to hospitalization.
      3371 (5.30)1.74 (1.32-2.29)<.001
      History of thrombophilia
      History of thrombophilia includes positive Factor V Leiden, positive prothrombin, G20210A gene mutation, positive lupus anticoagulant, other congenial or acquired thrombophilia.
      228 (0.36)1.66 (0.62-4.42).32
      Severe lung disease5726 (9.01)1.63 (1.30-2.04)<.001
      Age, ≥75 (y)22,660 (35.66)1.34 (1.15-1.56)<.001
      Pneumonia (<1 mo)9036 (14.22)1.23 (1.00-1.51).05
      Other acute infection
      Other acute infections include primary reason for admission based on the following Healthcare Cost and Utilization Project (HCUP) infection categories: 1.1 - Tuberculosis; 1.2 - Septicemia; 1.3 - Bacterial infection, unspecified site; 1.4 - Mycoses; 1.5 - HIV infection; 1.6 - Hepatitis; 1.7 - Viral infections; 1.8 - Other infections, including parasitic; 1.9 - Sexually transmitted infection (not HIV or hepatitis); 6.78 - Other CNS infection and poliomyelitis; 6.9 - Inflammation, infection of the eye; 8.126 - Other upper respiratory infection; 9.135 - Intestinal infection; 10.159 - Urinary tract infections; 12.197 - Skin and subcutaneous tissue infection; and 13.201 - Infective arthritis and osteomyelitis.
      8777 (13.81)1.23 (1.00-1.51).05
      Congestive heart failure5928 (9.33)1.13 (0.87-1.45).36
      Sepsis (<1 mo)6555 (10.32)1.10 (0.86-1.40).45
      Obesity (BMI >30)22,369 (35.20)0.82 (0.69-0.96).02
      Myocardial infarction (<1 mo)1061 (1.67)0.86 (0.44-1.65).65
      Inflammatory bowel disease2015 (3.17)0.78 (0.48-1.26).31
      Stroke3037 (4.78)0.40 (0.23-0.69).001
      BMI = body mass index; CI = confidence interval; HR = hazard ratio.
      Cancer diagnosis or acute cancer treatment within the last year.
      Central venous catheters were considered present on admission if patient was admitted with device or if one was inserted on day 1 or day 2 of the admission.
      Immobile/not ambulating defined as having at least one of the following: immobilizing plaster cast, paralysis, or bed rest for ≥72 hours prior to hospitalization.
      § History of thrombophilia includes positive Factor V Leiden, positive prothrombin, G20210A gene mutation, positive lupus anticoagulant, other congenial or acquired thrombophilia.
      Other acute infections include primary reason for admission based on the following Healthcare Cost and Utilization Project (HCUP) infection categories: 1.1 - Tuberculosis; 1.2 - Septicemia; 1.3 - Bacterial infection, unspecified site; 1.4 - Mycoses; 1.5 - HIV infection; 1.6 - Hepatitis; 1.7 - Viral infections; 1.8 - Other infections, including parasitic; 1.9 - Sexually transmitted infection (not HIV or hepatitis); 6.78 - Other CNS infection and poliomyelitis; 6.9 - Inflammation, infection of the eye; 8.126 - Other upper respiratory infection; 9.135 - Intestinal infection; 10.159 - Urinary tract infections; 12.197 - Skin and subcutaneous tissue infection; and 13.201 - Infective arthritis and osteomyelitis.
      Figure thumbnail gr1
      FigureRisk assessment model description. BMI = body mass index; CHF = congestive heart failure; MI = myocardial infarction; PICC = peripherally inserted central catheter; RAM = risk assessment model; VTE = venous thromboembolism. aHistory of thrombophilia includes positive Factor V Leiden, positive prothrombin, G20210A gene mutation, positive lupus anticoagulant, other congenial or acquired thrombophilia. bCurrent cancer defined as cancer diagnosed or treated within the year prior to hospital admission. cImmobile defined as having at least one of the following: immobilizing plaster cast, paralysis, or bed rest for ≥72 hours prior to hospitalization. dOther acute infections include primary reason for admission based on the following Healthcare Cost and Utilization Project (HCUP) infection categories: 1.1 - Tuberculosis, 1.2 - Septicemia, 1.3 - Bacterial infection, unspecified site, 1.4 - Mycoses, 1.5 - HIV infection, 1.6 - Hepatitis, 1.7 - Viral infections, 1.8 - Other infections, including parasitic, 1.9 - Sexually transmitted infection (not HIV or hepatitis), 6.78 - Other CNS infection and poliomyelitis, 6.9 - Inflammation, infection of the eye, 8.126 - Other upper respiratory infection, 9.135 - Intestinal infection, 10.159 - Urinary tract infections, 12.197 - Skin and subcutaneous tissue infection, and 13.201 - Infective arthritis and osteomyelitis. ePICC presence defined as admitted with PICC or inserted on day 1 or 2 of hospitalization.
      Among the eligible patients, 1183 (1.86%) were transferred to intensive care for reasons other than venous thromboembolism or transitioned to palliative/comfort care. Additionally, 3968 (6.24%) died between the index hospitalization admission and 90 days. Venous thromboembolism was associated with 3 transfers to the ICU and 10 deaths. A total of 670 (1.05%) patients developed a venous thromboembolism during follow-up. Two-by-two tables displaying the number of venous thromboembolism events and receipt of pharmacologic prophylaxis by “at-risk” status for each risk assessment model are presented in Table 2.
      Table 290-day Venous Thromboembolism and Receipt of Pharmacologic Prophylaxis by Risk Strata per Risk Assessment Model
      KucherPaduaIMPROVEIntermountain
      At-riskLow-riskAt-riskLow-riskAt-riskLow-riskAt-riskLow-risk
      Venous thromboembolism168502247423199471300370
      No venous thromboembolism640556,47310,33952,539724355,63511,85451,024
      Pharmacologic prophylaxis379634,927648132,242430434,419728531,438
      No pharmacologic prophylaxis277722,048410520,720313821,687486919,956
      Overall, 5,385,199 patient-days of follow-up were observed. The hazards of venous thromboembolism and model performance statistics estimated for each of the existing risk assessment models, by risk level and receipt of pharmacologic prophylaxis, are shown in Table 3. The hazards of 90-day venous thromboembolism were approximately threefold greater in the at-risk group relative to the low-risk group for all risk assessment models. Venous thromboembolism incidence in the low-risk group as defined by the Intermountain risk assessment model was the lowest of all risk assessment models assessed. Model discrimination was poor to moderate for all risk assessment models. The Intermountain risk assessment model also demonstrated the highest (albeit, moderate) degree of discrimination (Harrell's C-index = 0.61 [95% confidence interval, 0.60-0.62]). Receipt of pharmacologic prophylaxis during the index hospitalization was not found to be an independent predictor of 90-day venous thromboembolism, and model performance was essentially unchanged by receipt of pharmacologic prophylaxis for all risk assessment models assessed (Table 3). Additional calibration results based on the distributions of the 90-day venous thromboembolism event rates and hazard ratios for each score level per risk assessment model are illustrated in Supplementary Table 1 (available online).
      Table 390-day Venous Thromboembolism Incidence and Binary Risk Score Analysis by Risk Assessment Model
      KucherPaduaIMPROVEIntermountain
      All patients (n = 63,458)
       At-risk incidence
      Numerators in each “at-risk incidence” and “low-risk incidence” cells reflect number of VTEs and denominators within these cells represent total number of patient days. Incidence rates presented in parentheses per 10,000 patient-days.
      168/518,297 (3.24)247/832,082 (2.97)199/587,387 (3.39)300/958,318 (3.13)
       Low-risk incidence502/4,866,902 (1.03)423/4,553,117 (0.93)471/4,797,892 (0.98)370/4,426,881 (0.84)
       HR (95% CI)
      Hazard ratios and discrimination metrics based on binary risk models with “at-risk” per the respective scales of each risk assessment model coded as 1. As such, hazard ratios above 1 reflect increased hazard of VTE for patients in the respective “at-risk” categories. Harrell's C-index generated in somersd package in Stata.
      3.01 (2.52-3.59)3.08 (2.63-3.61)3.32 (2.81-3.92)3.59 (3.08-4.19)
       Harrell's C (95% CI)0.563 (0.558-0.568)0.600 (0.594-0.606)0.570 (0.565-0.576)0.611 (0.605-0.618)
      Pharmacologic prophylaxis
      Receipt of pharmacologic VTE prophylaxis was defined as receipt of any of the following treatments on day 1 or 2 of the index hospitalization: heparin 5000 units twice a day (BID); heparin 5000 units three times a day (TID); heparin 7500 units TID (for morbid obesity); enoxaparin 40 mg daily; enoxaparin 30 mg daily (for creatinine clearance <30 mL/min); enoxaparin 30 mg BID; dalteparin 5000 units daily; or fondaparinux 2.5 mg daily.
      during index hospitalization (n = 38,723)
       At-risk incidence98/302,398 (3.24)148/513,909 (2.88)120/342,432 (3.50)170/577,667 (2.94)
       Low-risk incidence301/2,980,961 (1.01)251/2,769,450 (0.91)279/2,940,927 (0.95)229/2,705,692 (0.85)
       HR (95% CI)3.08 (2.38-3.99)2.99 (2.37-3.76)3.71 (2.91-4.72)3.19 (2.54-4.00)
       Harrell's C (95% CI)0.551 (0.545-0.558)0.590 (0.582-0.597)0.560 (0.553-0.567)0.601 (0.593-0.609)
      No pharmacologic prophylaxis during index hospitalization (n = 24,825)
       At-risk incidence70/215,899 (3.24)99/318,173 (3.11)79/244,875 (3.23)130/380,651 (3.42)
       Low-risk incidence201/1,885,941 (1.07)172/1,783,667 (0.96)192/1,856,965 (1.03)141/1,721,189 (0.82)
       HR (95% CI)3.66 (2.70-4.97)3.61 (2.71-4.80)3.85 (2.87-5.17)4.16 (3.14-5.50)
       Harrell's C (95% CI)0.581 (0.572-0.590)0.616 (0.606-0.626)0.588 (0.578-0.597)0.625 (0.614-0.635)
      CI = confidence interval; HR = hazard ratio.
      Numerators in each “at-risk incidence” and “low-risk incidence” cells reflect number of VTEs and denominators within these cells represent total number of patient days. Incidence rates presented in parentheses per 10,000 patient-days.
      Hazard ratios and discrimination metrics based on binary risk models with “at-risk” per the respective scales of each risk assessment model coded as 1. As such, hazard ratios above 1 reflect increased hazard of VTE for patients in the respective “at-risk” categories. Harrell's C-index generated in somersd package in Stata.
      Receipt of pharmacologic VTE prophylaxis was defined as receipt of any of the following treatments on day 1 or 2 of the index hospitalization: heparin 5000 units twice a day (BID); heparin 5000 units three times a day (TID); heparin 7500 units TID (for morbid obesity); enoxaparin 40 mg daily; enoxaparin 30 mg daily (for creatinine clearance <30 mL/min); enoxaparin 30 mg BID; dalteparin 5000 units daily; or fondaparinux 2.5 mg daily.

      Sensitivity Analyses

      After modifying our venous thromboembolism outcome definition to include only pulmonary embolism or lower-extremity deep vein thromboses (or both), the 63,548 patients contributed a total of 5,391,673 patient-days of follow-up. A total of 512 (0.81%) patients had a 90-day venous thromboembolism using this modified outcome definition. Overall, risk assessment model performances for the modified venous thromboembolism definition were similar to those identified in our primary analyses (see Supplementary Table 2, Supplementary Table 3; available online).

      Discussion

      Several important findings emerged from our risk assessment model derivation and validation study of 90-day venous thromboembolism outcomes in a multi-site, large cohort of acutely ill hospitalized medical patients. First, overall venous thromboembolism incidence in this population was low (∼1%). Second, existing risk assessment models to be used at hospital admission classified <20% of our patient population as being “at-risk” for venous thromboembolism, suggesting either that: 1) risk stratification algorithms of existing risk assessment models may not adequately identify at-risk subsets of nonsurgical, hospitalized medical patients, or 2) using these risk-stratification methods the majority of this patient population is at low risk for venous thromboembolism. Third, although existing binary risk assessment models demonstrated good calibration, discrimination characteristics were uniformly poor to moderate across models. Taken together, these findings suggest that models for at-admission risk assessment in hospitalized medical patients are suboptimal and are deserving of further enquiry.
      The preexisting risk assessment models assessed in our external validation study have all been evaluated in medical populations. Many criticisms of the existing risk assessment models–derivation populations comprised mostly of surgical or cancer patients, risk factor selection based on expert consensus, and lack of external validation–have been noted.
      • Spyropoulos A.C.
      • McGinn T.
      • Khorana A.A.
      The use of weighted and scored risk assessment models for venous thromboembolism.
      Kucher et al
      • Kucher N.
      • Koo S.
      • Quiroz R.
      • et al.
      Electronic alerts to prevent venous thromboembolism among hospitalized patients.
      devised a computer-alert program to remind physicians to prescribe venous thromboembolism prophylaxis for high-risk patients. They developed their risk assessment model in a cohort of medical and surgical patients, with cancer present in 80% of their population. In our external validation analyses, the Kucher risk assessment model yielded the poorest discrimination characteristics, raising questions as to its generalizability to all medical patients.
      Barbar et al
      • Barbar S.
      • Noventa F.
      • Rossetto V.
      • et al.
      A risk assessment model for the identification of hospitalized medical patients at risk for venous thromboembolism: the Padua Prediction Score.
      developed the Padua risk assessment model, using a cohort of patients admitted to an internal medicine service in Italy. Limitations of this study include its single-center setting, non-evidence-derived data assessments, small number of venous thromboembolism events, no presentation of model performance characteristics, and suboptimal external validation.
      • Nendaz M.
      • Spirk D.
      • Kucher N.
      • et al.
      Multicentre validation of the Geneva Risk Score for hospitalised medical patients at risk of venous thromboembolism. Explicit ASsessment of Thromboembolic RIsk and Prophylaxis for Medical PATients in SwitzErland (ESTIMATE).
      Approximately 40% of the original derivation cohort patients used to develop the Padua score were classified as at risk,
      • Barbar S.
      • Noventa F.
      • Rossetto V.
      • et al.
      A risk assessment model for the identification of hospitalized medical patients at risk for venous thromboembolism: the Padua Prediction Score.
      compared with 17% in our validation cohort. In addition to a smaller percentage of cancer patients in our validation cohort, our operationalization of immobility classified only 5% of patients as immobile, compared with 23% in the original derivation cohort.
      • Barbar S.
      • Noventa F.
      • Rossetto V.
      • et al.
      A risk assessment model for the identification of hospitalized medical patients at risk for venous thromboembolism: the Padua Prediction Score.
      The IMPROVE multinational registry was designed to examine venous thromboembolism prophylaxis practices and clinical outcomes in hospitalized medical patients. The initial derivation of the predictive IMPROVE risk assessment model included 4 risk factors available at hospital admission (history of previous venous thromboembolism, known thrombophilia, cancer, and age >60 years) that were independently associated with venous thromboembolism events and yielded a c-statistic of 0.65.
      • Spyropoulos A.C.
      • Anderson Jr., F.A.
      • Fitzgerald G.
      • et al.
      Predictive and associative models to identify hospitalized medical patients at risk for VTE.
      Applying the predictive IMPROVE risk assessment model in our validation cohort classified approximately 12% of patients as at risk and demonstrated poorer discrimination characteristics compared with the original derivation cohort. Due to differences in clinical characteristics and uncertainty in the need for pharmacologic prophylaxis in upper- vs lower-extremity deep vein thrombosis,
      • Spyropoulos A.C.
      Upper vs. lower extremity deep vein thrombosis: outcome definitions of venous thromboembolism for clinical predictor rules or risk factor analyses in hospitalized patients.
      the risk assessment model developed in the original derivation study excluded upper-extremity deep vein thrombosis events from the combined venous thromboembolism outcome. In our sensitivity analysis excluding upper-extremity deep vein thrombosis, we found that the predictive IMPROVE risk assessment model did yield marginally better calibration, but discrimination characteristics were comparable with those detected in our primary analyses. Given that a patient's risk status may change during the course of hospitalization, and especially to determine a patient's continued venous thromboembolism risk after hospital discharge, the original IMPROVE study also developed an associative model, including the additional risk factors of immobility, lower limb paralysis, and ICU stay, which yielded an improved c-statistic of 0.69.
      • Spyropoulos A.C.
      • Anderson Jr., F.A.
      • Fitzgerald G.
      • et al.
      Predictive and associative models to identify hospitalized medical patients at risk for VTE.
      Recent large-scale external validation studies of the associative IMPROVE risk assessment model have shown good calibration and discrimination, suggesting that the IMPROVE associative venous thromboembolism risk assessment model may reliably stratify venous thromboembolism risk,
      • Mahan C.E.
      • Liu Y.
      • Turpie A.G.
      • et al.
      External validation of a risk assessment model for venous thromboembolism in the hospitalised acutely-ill medical patient (VTE-VALOURR).
      • Rosenberg D.
      • Eichorn A.
      • Alarcon M.
      • McCullagh L.
      • McGinn T.
      • Spyropoulos A.C.
      External validation of the risk assessment model of the International Medical Prevention Registry on Venous Thromboembolism (IMPROVE) for medical patients in a tertiary health system.
      • Cohen A.T.
      • Spiro T.E.
      • Spyropoulos A.C.
      • et al.
      D-dimer as a predictor of venous thromboembolism in acutely ill, hospitalized patients: a subanalysis of the randomized controlled MAGELLAN trial.
      and this model has been incorporated into a recent large multicenter, multinational trial of post-hospital-discharge thromboprophylaxis, the MARINER trial (ClinicalTrials.gov Identifier: NCT02111564). However, the main limitation of the IMPROVE associative risk assessment model is incomplete validation during hospital admission. We were unable to externally validate this associative model in our study, primarily due to data censoring when a patient is transferred to the ICU in our study cohort.
      Woller et al
      • Woller S.C.
      • Stevens S.M.
      • Jones J.P.
      • et al.
      Derivation and validation of a simple model to identify venous thromboembolism risk in medical patients.
      derived the Intermountain risk assessment model using administrative and electronic medical record data from 143,000 medical patient admissions. Four risk factors (history of previous venous thromboembolism, an order for bed rest, peripherally inserted central venous catheter, and cancer) were predictive of hospital-associated venous thromboembolism in up to 90 days of follow-up, with good discrimination in their derivation and validation cohorts (c-statistics = 0.87 and 0.84, respectively).
      • Woller S.C.
      • Stevens S.M.
      • Jones J.P.
      • et al.
      Derivation and validation of a simple model to identify venous thromboembolism risk in medical patients.
      Although the discrimination of the Intermountain risk assessment model in our validation cohort was the best of the risk assessment models assessed, underlying patient differences between the Intermountain cohorts and our validation cohort likely account for poorer performance in our validation sample. The prevalence of each of the 4 risk factors included in this risk assessment model was greater among the original derivation cohort compared with our validation cohort. Most notably, approximately 45% of patients in the Intermountain derivation cohort had a history of cancer.
      Our study has limitations. First, obtaining individual patient-level venous thromboembolism risk factors was constrained by each abstractor's ability to retrospectively find information in patient charts. However, data collection was standardized and quality ensured through use of annual audits at each hospital. Although we employed a reliable data-collection strategy and made every effort to match how venous thromboembolism risk factors were operationalized, differences in how we implemented venous thromboembolism risk factors (eg, immobility) relative to the derivation studies may exist. Second, although complete 90-day follow-up chart review was completed for 100% of the patients, the successful telephone follow-up rate to identify venous thromboembolism events without medical record documentation was 58%. Therefore, it is possible that some venous thromboembolism events were missed after hospital discharge, especially if these were treated at another institution. However, our reported venous thromboembolism incidence in medical patients is comparable with other estimates in the literature.
      • Lederle F.A.
      • Zylla D.
      • MacDonald R.
      • Wilt T.J.
      Venous thromboembolism prophylaxis in hospitalized medical patients and those with stroke: a background review for an American College of Physicians Clinical Practice Guideline.
      • Rothberg M.B.
      • Lindenauer P.K.
      • Lahti M.
      • Pekow P.S.
      • Selker H.P.
      Risk factor model to predict venous thromboembolism in hospitalized medical patients.
      Third, given that risk assessment strategies in medical patients often follow a binary approach (at-risk vs low-risk), we assessed model discrimination based upon binary risk prediction models for each risk assessment model. These performance metrics may not align with those reported in derivation studies that report (explicitly or otherwise) c-statistics based on the full regression models containing all respective risk factors as independent predictors.
      Despite these limitations, our study has several strengths. First, our study provides external validation of several existing risk assessment models within a large, representative, multi-site, cohort of hospitalized medical patients–a population that comprises the greatest proportion of patients in hospitals. Second, our standardized data collection by trained medical record abstractors coupled with centralized audits ensures data reliability and avoids problems often associated with retrospective datasets (eg, missing data, misclassification of International Classification of Diseases, 9th Revision coding) to determine venous thromboembolism risk factors and outcomes. Third, to ensure transparency of our results, we followed the recently published TRIPOD checklist for presenting results on risk prediction model validation studies.
      • Collins G.S.
      • de Groot J.A.
      • Dutton S.
      • et al.
      External validation of multivariable prediction models: a systematic review of methodological conduct and reporting.
      • Collins G.S.
      • Reitsma J.B.
      • Altman D.G.
      • Moons K.G.
      Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): the TRIPOD statement.
      • Moons K.G.
      • Altman D.G.
      • Reitsma J.B.
      • et al.
      Transparent Reporting of a multivariable prediction model for Individual Prognosis or Diagnosis (TRIPOD): explanation and elaboration.
      In doing so, we have shown that several existing risk assessment models have limited ability in identifying the greatest at-risk population for whom pharmacologic prophylaxis may benefit. Given the low incidence of venous thromboembolism in hospitalized medical (non-ICU, nonsurgical) patients, the risk assessment models described in this work may have greater utility in identifying subsets of this particular patient population for whom pharmacologic prophylaxis is not warranted.

      Acknowledgment

      The authors thank Jason Mann, MSA for his assistance with proofing and submitting this manuscript.

      Supplementary Data

      Supplementary Table 190-day Venous Thromboembolism Events and Hazard Ratios at Each Risk Level per Risk Assessment Model
      ScoreVTENo VTEEvent RateHR (95% CI)
      Kucher
       08715,7350.55%ref
       124928,9900.85%1.59 (1.25-2.03)
       27282540.86%1.64 (1.20-2.25)
       39434942.62%5.21 (3.89-6.97)
       410543402.36%4.76 (3.58-6.33)
       53013042.25%4.45 (2.94-6.74)
       6203894.89%10.04 (6.17-16.32)
       782712.87%6.07 (2.94-12.52)
       83793.66%8.32 (2.63-26.30)
       921611.11%28.78 (7.08-117.02)
       10060.00%-
      Padua
       05310,7000.49%ref
       115720,3010.77%1.59 (1.16-2.16)
       213014,9120.86%1.82 (1.32-2.50)
       38366261.24%2.68 (1.90-3.78)
       412048462.42%5.39 (3.90-7.45)
       56332101.92%4.36 (3.02-6.28)
       63112332.45%5.71 (3.67-8.90)
       7205673.41%7.93 (4.74-13.26)
       873012.27%5.49 (2.50-12.09)
       951213.97%9.91 (3.96-24.79)
       100430.00%-
       111156.25%20.55 (2.84-148.72)
       12020.00%-
       13010.00%-
      IMPROVE
       011520,2500.56%ref
       135635,3851.00%1.85 (1.50-2.29)
       28631482.66%5.65 (4.27-7.47)
       33010242.85%5.16 (3.45-7.71)
       47027182.51%4.71 (3.50-6.34)
       5112663.97%8.86 (4.77-16.45)
       61442.22%3.96 (0.55-28.37)
       71392.50%4.72 (0.66-33.77)
       8040.00%-
      Intermountain
       037051,0240.72%ref
       125310,6702.32%3.46 (2.95-4.06)
       24311323.66%5.82 (4.24-7.98)
       34507.41%14.25 (5.32-38.18)
       4020.00%-
      CI = confidence interval; HR = hazard ratio; VTE = venous thromboembolism.
      Supplementary Table 2Incidence and Binary Risk Score Analysis for 90-day Pulmonary Embolism or Lower-Extremity Deep Vein Thrombosis per Risk Assessment Model
      KucherPaduaIMPROVEIntermountain
      VTE events/patient-days (incidence per 10,000 patient-days)
       At-risk137/519,532 (2.64)193/834,306 (2.31)165/588,622 (2.80)226/961,466 (2.35)
       Low-risk375/4,872,141 (0.77)319/4,557,367 (0.70)347/4,803,051 (0.72)286/4,430,207 (0.65)
      Binary RAM performance
       HR (95% CI)3.29 (2.70-4.00)3.19 (2.67-3.82)3.75 (3.11-4.52)3.52 (2.95-4.20)
       Harrell's C (95% CI)0.563 (0.558-0.569)0.599 (0.594-0.606)0.571 (0.565-0.576)0.611 (0.604-0.617)
      CI = confidence interval; DVT = deep vein thrombosis; HR = hazard ratio; RAM = risk assessment model; VTE = venous thromboembolism.
      Supplementary Table 390-day Pulmonary Embolism or Lower-Extremity Deep Vein Thrombosis Events and Hazard Ratios at Each Risk Level per Risk Assessment Model
      ScoreVTENo VTEEvent RateHR (95% CI)
      Kucher
       06615,7560.42%ref
       118729,0520.64%1.58 (1.19-2.09)
       25382730.64%1.60 (1.11-2.29)
       36935191.92%5.03 (3.59-7.05)
       48843571.98%5.27 (3.83-7.25)
       52613081.95%5.08 (3.23-8.00)
       6123972.93%7.89 (4.27-14.60)
       772722.51%7.02 (3.22-15.30)
       82802.44%7.31 (1.79-29.87)
       921611.11%38.25 (9.36-156.31)
       10060.00%-
      Padua
       04110,7120.38%ref
       111820,3400.58%1.54 (1.08-2.20)
       29614,9460.64%1.73 (1.20-2.50)
       36466450.95%2.66 (1.80-3.94)
       49048761.81%5.23 (3.61-7.57)
       55232211.59%4.64 (3.08-6.99)
       62512391.98%5.96 (3.62-9.80)
       7165712.73%8.14 (4.57-14.52)
       873012.27%7.09 (3.18-15.81)
       921241.59%5.06 (1.22-20.93)
       100430.00%-
       111156.25%26.92 (3.70-195.94)
       12020.00%-
       13010.00%-
      IMPROVE
       06720,2980.33%ref
       128035,4610.78%2.50 (1.91-3.26)
       27131632.20%8.07 (5.77-11.27)
       32710272.56%7.97 (5.10-12.46)
       45627322.01%6.47 (4.54-9.22)
       5102673.61%13.96 (7.18-27.15)
       61442.22%6.82 (0.95-49.11)
       70400.00%-
       8040.00%-
      Intermountain
       028651,1080.56%ref
       119210,7311.76%3.39 (2.83-4.08)
       23311422.81%5.74 (4.00-8.23)
       31531.85%4.56 (0.64-32.50)
       4020.00%-
      CI = confidence interval; DVT = deep vein thrombosis; HR = hazard ratio; VTE = venous thromboembolism.

      References

        • Kahn S.R.
        • Lim W.
        • Dunn A.S.
        • et al.
        Prevention of VTE in nonsurgical patients: Antithrombotic Therapy and Prevention of Thrombosis, 9th ed: American College of Chest Physicians Evidence-Based Clinical Practice Guidelines.
        Chest. 2012; 141: e195S-e226S
        • Barbar S.
        • Noventa F.
        • Rossetto V.
        • et al.
        A risk assessment model for the identification of hospitalized medical patients at risk for venous thromboembolism: the Padua Prediction Score.
        J Thromb Haemost. 2010; 8: 2450-2457
        • Kucher N.
        • Koo S.
        • Quiroz R.
        • et al.
        Electronic alerts to prevent venous thromboembolism among hospitalized patients.
        N Engl J Med. 2005; 352: 969-977
        • Spyropoulos A.C.
        • Anderson Jr., F.A.
        • Fitzgerald G.
        • et al.
        Predictive and associative models to identify hospitalized medical patients at risk for VTE.
        Chest. 2011; 140: 706-714
        • Woller S.C.
        • Stevens S.M.
        • Jones J.P.
        • et al.
        Derivation and validation of a simple model to identify venous thromboembolism risk in medical patients.
        Am J Med. 2011; 124: 947-954.e2
        • Huang W.
        • Anderson F.A.
        • Spencer F.A.
        • Gallus A.
        • Goldberg R.J.
        Risk-assessment models for predicting venous thromboembolism among hospitalized non-surgical patients: a systematic review.
        J Thromb Thrombolysis. 2013; 35: 67-80
        • Mahan C.E.
        • Liu Y.
        • Turpie A.G.
        • et al.
        External validation of a risk assessment model for venous thromboembolism in the hospitalised acutely-ill medical patient (VTE-VALOURR).
        Thromb Haemost. 2014; 112: 692-699
        • Rosenberg D.
        • Eichorn A.
        • Alarcon M.
        • McCullagh L.
        • McGinn T.
        • Spyropoulos A.C.
        External validation of the risk assessment model of the International Medical Prevention Registry on Venous Thromboembolism (IMPROVE) for medical patients in a tertiary health system.
        J Am Heart Assoc. 2014; 3: e001152
        • Nendaz M.
        • Spirk D.
        • Kucher N.
        • et al.
        Multicentre validation of the Geneva Risk Score for hospitalised medical patients at risk of venous thromboembolism. Explicit ASsessment of Thromboembolic RIsk and Prophylaxis for Medical PATients in SwitzErland (ESTIMATE).
        Thromb Haemost. 2014; 111: 531-538
        • Flanders S.A.
        • Greene M.T.
        • Grant P.
        • et al.
        Hospital performance for pharmacologic venous thromboembolism prophylaxis and rate of venous thromboembolism: a cohort study.
        JAMA Intern Med. 2014; 174: 1577-1584
        • American College of Surgeons
        American College of Surgeons National Surgical Quality Improvement Program (ACS NSQIP®).
        American College of Surgeons, Chicago2013
        • Harrell Jr., F.E.
        • Lee K.L.
        • Mark D.B.
        Multivariable prognostic models: issues in developing models, evaluating assumptions and adequacy, and measuring and reducing errors.
        Stat Med. 1996; 15: 361-387
        • Newson R.
        Confidence intervals for rank statistics: Somers' D and extensions.
        Stata J. 2006; 6: 309-334
        • Spyropoulos A.C.
        Upper vs. lower extremity deep vein thrombosis: outcome definitions of venous thromboembolism for clinical predictor rules or risk factor analyses in hospitalized patients.
        J Thromb Haemost. 2009; 7: 1041-1042
        • Spyropoulos A.C.
        • McGinn T.
        • Khorana A.A.
        The use of weighted and scored risk assessment models for venous thromboembolism.
        Thromb Haemost. 2012; 108: 1072-1076
        • Cohen A.T.
        • Spiro T.E.
        • Spyropoulos A.C.
        • et al.
        D-dimer as a predictor of venous thromboembolism in acutely ill, hospitalized patients: a subanalysis of the randomized controlled MAGELLAN trial.
        J Thromb Haemost. 2014; 12: 479-487
        • Lederle F.A.
        • Zylla D.
        • MacDonald R.
        • Wilt T.J.
        Venous thromboembolism prophylaxis in hospitalized medical patients and those with stroke: a background review for an American College of Physicians Clinical Practice Guideline.
        Ann Intern Med. 2011; 155: 602-615
        • Rothberg M.B.
        • Lindenauer P.K.
        • Lahti M.
        • Pekow P.S.
        • Selker H.P.
        Risk factor model to predict venous thromboembolism in hospitalized medical patients.
        J Hosp Med. 2011; 6: 202-209
        • Collins G.S.
        • de Groot J.A.
        • Dutton S.
        • et al.
        External validation of multivariable prediction models: a systematic review of methodological conduct and reporting.
        BMC Med Res Methodol. 2014; 14: 40
        • Collins G.S.
        • Reitsma J.B.
        • Altman D.G.
        • Moons K.G.
        Transparent reporting of a multivariable prediction model for individual prognosis or diagnosis (TRIPOD): the TRIPOD statement.
        J Clin Epidemiol. 2015; 68: 134-143
        • Moons K.G.
        • Altman D.G.
        • Reitsma J.B.
        • et al.
        Transparent Reporting of a multivariable prediction model for Individual Prognosis or Diagnosis (TRIPOD): explanation and elaboration.
        Ann Intern Med. 2015; 162: W1-W73