$$\rightleftharpoonup{xx}$$
$$\longleftharp{xx}$$,
$$\longrightharp{xx}$$,
This study was conducted using de-identified, routinely collected clinical data. In accordance with institutional policies and national regulations, formal ethical approval and informed consent were waived as no identifiable patient information was used and no intervention was performed. The study adhered to the principles of the Declaration of Helsinki.
Study design
This multicenter retrospective cohort study was conducted using data from institutional joint replacement registries and electronic health records across several high-volume orthopedic centers. The study utilized prospectively collected registry data complemented by retrospective chart review. The study adhered to the Strengthening the Reporting of Observational Studies in Epidemiology (STROBE) guidelines for observational cohort studies12. Data quality was ensured through internal validation checks and selective audits. A schematic overview of the study design, including patient selection, grouping, and outcome assessment, is provided in Figure 1.

Figure 1: Study flow diagram and schematic overview of study design. Please click here to view a larger version of this figure.
Study population
Adults (≥18 years) undergoing elective primary THA or TKA for osteoarthritis were included. Exclusion criteria included preoperative anticoagulation therapy, known coagulopathy, revision surgery, or trauma-related arthroplasty. In cases of bilateral arthroplasty, only the first procedure was analyzed to avoid duplication. Patients were identified from multiple tertiary orthopedic centers to enhance generalizability.
Data collection
Data were extracted from electronic medical records and institutional databases using a standardized data collection form. Variables included demographics (age, sex, body mass index [BMI]), lifestyle factors (smoking status and alcohol use), comorbidities (history of VTE, diabetes mellitus, hypertension, dyslipidemia, chronic kidney disease, active cancer, and thrombophilia), and surgical variables (procedure type, operative duration, and tranexamic acid administration). The American Society of Anesthesiologists (ASA) classification and Charlson Comorbidity Index were calculated. Perioperative care variables included pharmacologic and mechanical prophylaxis, early mobilization (postoperative day 0–1), and hospital length of stay.
Thromboprophylaxis groups
Patients were categorized according to postoperative anticoagulation strategy. The rivaroxaban group received rivaroxaban 10 mg orally once daily, starting 6–10 h postoperatively and continued for 14 days (TKA) or 35 days (THA), in line with clinical guidelines13. The LMWH group received low molecular weight heparin (e.g., enoxaparin 40 mg once daily or 30 mg twice daily per institutional protocol), initiated 12–24 h postoperatively and continued for the same duration. Adherence was assessed using prescription records and patient self-reports. Rivaroxaban’s oral route was noted to potentially improve adherence in comparison to parenteral LMWH14.
Imaging and DVT/PE diagnostic protocols
All patients underwent standardized bilateral lower-extremity duplex ultrasonography for DVT screening between postoperative days 7 and 10 or at the time of hospital discharge, irrespective of symptom status. DVT diagnosis included both symptomatic and asymptomatic cases, identified through routine bilateral lower-extremity duplex ultrasonography performed according to the study protocol. Duplex studies were performed by certified vascular technologists using compression and Doppler flow criteria. Symptomatic DVT was defined based on clinical presentation with confirmatory imaging, whereas asymptomatic DVT was detected through scheduled screening ultrasonography. Symptomatic pulmonary embolism was diagnosed based on clinical suspicion and confirmed using computed tomography pulmonary angiography. Imaging protocols were harmonized across participating centers to ensure consistency in diagnostic criteria and timing.
Adherence and exposure assessment
Exposure to thromboprophylaxis was verified using a triangulated approach including electronic prescription records, pharmacy refill documentation, and structured patient self-report in the course of the postoperative follow-up visits. Adherence was defined as the patient having 80% or more of the prescription doses verified during the prophylaxis period. A per-protocol sensitivity analysis excluded patients with early discontinuation, non-adherence, or crossover between anticoagulant regimens.
Outcome measures
The primary endpoint was the 30-day incidence of postoperative DVT, confirmed by duplex ultrasonography. All patients underwent standardized ultrasound screening around postoperative day 7–10 or at discharge, regardless of symptoms. Secondary outcomes included pulmonary embolism, major bleeding, minor bleeding, laboratory parameters, wound complications, readmissions, and mortality.
Major bleeding was defined according to the International Society on Thrombosis and Haemostasis (ISTH) criteria15, including bleeding leading to reoperation, transfusion ≥ 2 units, a hemoglobin decrease of ≥2 g/dL within the first 5 postoperative days, or critical organ involvement. Minor bleeding included wound oozing, hematomas, and prolonged drainage. Laboratory parameters, including D-dimer, hemoglobin (Hb), and glucose levels, were recorded preoperatively and on postoperative days 1, 3, and 5. Elevated D-dimer has been associated with postoperative DVT, although its specificity is limited. Wound complications included prolonged drainage, superficial infection, or dehiscence.
Risk factor assessment
A comprehensive set of known or suspected DVT risk factors was recorded, including age, obesity, smoking, diabetes, prior VTE, and surgery-specific variables. Advanced age, obesity, and history of VTE are established risk factors16. Diabetes mellitus was examined because it has been found to be associated with enhanced DVT risk in joint replacement arthroplasty17. The influence of bilateral procedures, extended operating time16, and TXA use18 on thrombotic and bleeding outcomes was also compared.
Center-level variation handling
Random-effects modeling was used to test the heterogeneity, which might occur among the participating institutions, by including the treatment center as a clustering variable. Multilevel mixed-effects logistic regression models were applied to remove center-level variations in surgical volume, perioperative practice, and imaging practice. Sensitivity analyses incorporating center-specific random intercepts were performed to ensure robustness of the primary and secondary outcome estimates.
Statistical analysis
Baseline characteristics were compared using Student’s t-test or Mann–Whitney U test for continuous variables and chi-square or Fisher’s exact test for categorical variables. Multivariable logistic regression identified independent predictors of DVT and bleeding. Variables significant in univariate analysis (p < 0.10) or strongly supported by literature were included, including pre-existing venous thromboembolism (VTE), defined as any documented history of deep vein thrombosis or pulmonary embolism prior to the index surgery, and baseline anemia, defined according to World Health Organization (WHO) criteria as a preoperative hemoglobin level < 13 g/dL in males and <12 g/dL in females. Preoperative hemoglobin values were obtained from routine laboratory testing performed within 48 h prior to surgery.
The effect of prophylactic agents on DVT incidence was assessed through multivariate adjustment, where adjusted odds ratios (ORs) with 95% confidence intervals (CIs) were reported; propensity score matching (nearest-neighbor matching with caliper restriction), where patients were matched on key baseline covariates including age, sex, body mass index (BMI), smoking status, comorbidities (e.g., diabetes, prior VTE), ASA class, and type of surgery (THA/TKA) to minimize confounding, followed by paired statistical tests. Balance between groups was assessed using standardized mean differences, with values <0.1 indicating adequate covariate balance.
Survival analysis was performed using Kaplan-Meier curves and Cox proportional hazards models to assess DVT-free survival within 30 days. Sensitivity analyses included exclusion of asymptomatic DVT cases, per-protocol analysis (excluding patients with early discontinuation or crossover), and random-effects modeling to account for center-level clustering. Statistical significance was set at two-sided p < 0.05, and analyses were conducted using an appropriate statistical software.