Cohort characterisitics and dynamic risk
The cohort comprised 56 warfarin recipients, 38 edoxaban 60-mg recipients, 41 edoxaban 30-mg recipients, 35 rivaroxaban 15-mg recipients, and 32 rivaroxaban 10-mg recipients; results are expressed in Tables 1–4 and Supplementary Tables 1–9. The groups differed in body weight, CrCl, frailty, CHA₂DS₂-VASc, antiplatelet use, baseline PT, and baseline APTT (Table S1). Overall, 85 participants (42.1%) had baseline HAS-BLED scores of at least 3; all 202 had CHA₂DS₂-VASc scores of at least 2. CHA₂DS₂-VASc did not change during follow-up. HAS-BLED increased by 1 point in 16 participants, decreased by 1 point in 12, and remained unchanged in 174. At month 6, 92 participants (45.5%) had HAS-BLED scores of at least 3. Complete baseline data and longitudinal risk-score results are reported in Supplementary Tables 1 and 2.
Official-label dose concordance
At baseline, 72 of 146 DOAC recipients (49.3%) received a label-concordant dose, 72 (49.3%) received a lower-than-label dose, and 2 (1.4%) received a higher-than-label dose (Table 2). Among edoxaban recipients, 68 of 79 (86.1%) were label-concordant, 9 received a lower-than-label dose, and 2 received a higher-than-label dose. Among rivaroxaban recipients, 4 of 67 (6.0%) were label-concordant, and 63 received a lower-than-label dose; this included 31 of 35 recipients of rivaroxaban 15 mg and all 32 recipients of rivaroxaban 10 mg. Across 730 DOAC assessments during follow-up, 355 (48.6%) were label-concordant, 359 (49.2%) were lower than label, and 16 (2.2%) were higher than label. Follow-up results are provided in Supplementary Table 3.
Clinical outcomes
Eight participants experienced an embolic event: 3 ischemic strokes and 5 TIAs. Twenty-one participants had bleeding, all classified as ISTH minor; no major or clinically relevant non-major bleeding occurred. Neither embolic events nor minor bleeding differed significantly among groups (Fisher–Freeman–Halton exact p = 0.846 and p = 0.199, respectively; Table 3). Thirty-five participants experienced a non-bleeding adverse reaction: 18 had gastrointestinal symptoms and 17 had dizziness/vertigo. The overall distribution differed among groups (p = 0.014; Cramer’s V = 0.248). Pairwise risk ratios were imprecise; after Holm adjustment, only rivaroxaban 15 mg versus warfarin remained significant for other adverse reactions (risk ratio 0.19, 95% CI 0.05–0.77; adjusted p = 0.027). Detailed phenotypes and all unadjusted pairwise estimates are provided in Supplementary Tables 4 and 5.
Exploratory adjusted analysis
Generalized overlap weighting yielded an effective sample size of 89.9. The maximum absolute pairwise standardized mean difference decreased from 2.04 before weighting to 1.36 after weighting, driven mainly by persistent CrCl imbalance. Weighted point estimates generally favored the DOAC regimens for embolic and bleeding outcomes, but bootstrap CIs were wide and frequently included 1.00 (Supplementary Table 6). These estimates were exploratory and did not support causal or equivalence claims.
Laboratory trajectories
Hemoglobin showed a small time effect (p < 0.001; partial η2 = 0.024), with no treatment-by-time interaction (p = 0.888) or difference in baseline-to-month-6 change; it was not attributed to treatment. ALT, serum creatinine, and PT showed no time effects or interactions. Complete serial values are provided in Supplementary Table 7. APTT showed a small time effect (p = 0.0015; partial η2 = 0.022) and a treatment-by-time interaction (p < 0.001; partial η2 = 0.054; Table 4 and Figure 1). Baseline-to-month-6 APTT change differed among groups (H = 11.52, p = 0.021; ε2 = 0.038), whereas changes in the other four indicators did not (Supplementary Table 8). Within-group sensitivity analyses are reported in Supplementary Table 9. PT and APTT were interpreted descriptively, not as calibrated measures of drug exposure or clinical benefit24.
DATA AVAILABILITY:
All data are provided in the Raw Data Supplementary Package, named “raw data”.

Figure 1: Mean activated partial thromboplastin time across baseline, day 3, month 1, month 3, and month 6 by anticoagulant group. Error bars show 95% CIs. APTT is shown descriptively and is not a calibrated measure of comparative factor Xa inhibitor exposure, effectiveness, or safety. Please click here to view a larger version of this figure.
| Characteristic | Warfarin (n = 56) | Edoxaban 60 mg (n = 38) | Edoxaban 30 mg (n = 41) | Rivaroxaban 15 mg (n = 35) | Rivaroxaban 10 mg (n = 32) | p value |
| Age, years | 85 [83–87] | 86 [83–88] | 86 [83–88] | 85 [83–88] | 84 [82–87] | 0.448 |
| Body weight, kg | 64.8 [56.0–71.7] | 68.4 [64.1–75.4] | 61.5 [55.9–65.8] | 65.7 [61.5–73.9] | 58.1 [53.5–60.5] | <0.001 |
| Cockcroft–Gault CrCl, mL/min | 41.0 [32.1–55.2] | 66.3 [57.2–72.9] | 47.7 [39.7–55.2] | 65.7 [53.4–73.8] | 41.0 [32.8–48.0] | <0.001 |
| Clinical Frailty Scale | 5 [5–6] | 4 [4–5] | 5 [4–5] | 5 [4–6] | 5 [5–6] | <0.001 |
| CHA2DS2-VASc score | 4 [4–5] | 5 [3–6] | 4 [4–6] | 5 [4–5] | 6 [4–6] | 0.040 |
| HAS-BLED score | 2 [2–3] | 2 [2–3] | 2 [2–3] | 2 [2–3] | 2 [2–3] | 0.425 |
| Female sex | 28 (50.0%) | 22 (57.9%) | 27 (65.9%) | 20 (57.1%) | 18 (56.2%) | 0.654 |
| CrCl ≤50 mL/min | 38 (67.9%) | 2 (5.3%) | 23 (56.1%) | 4 (11.4%) | 25 (78.1%) | <0.001 |
| Heart failure | 18 (32.1%) | 18 (47.4%) | 13 (31.7%) | 10 (28.6%) | 15 (46.9%) | 0.265 |
| Hypertension | 24 (42.9%) | 19 (50.0%) | 16 (39.0%) | 12 (34.3%) | 19 (59.4%) | 0.252 |
| Prior stroke/TIA/systemic embolism | 14 (25.0%) | 12 (31.6%) | 15 (36.6%) | 12 (34.3%) | 11 (34.4%) | 0.767 |
| Vascular disease | 14 (25.0%) | 7 (18.4%) | 14 (34.1%) | 8 (22.9%) | 12 (37.5%) | 0.324 |
| Previous clinically important bleeding | 25 (44.6%) | 24 (63.2%) | 23 (56.1%) | 17 (48.6%) | 16 (50.0%) | 0.463 |
| Concomitant antiplatelet therapy | 28 (50.0%) | 14 (36.8%) | 9 (22.0%) | 21 (60.0%) | 9 (28.1%) | 0.003 |
| HAS-BLED ≥3 | 26 (46.4%) | 16 (42.1%) | 14 (34.1%) | 15 (42.9%) | 14 (43.8%) | 0.820 |
Table 1: Baseline demographic and clinical characteristics by anticoagulant group. Values are median [IQR] or n (%). Continuous variables were compared with Kruskal–Wallis tests; categorical variables were compared with Pearson chi-square or Fisher–Freeman–Halton exact tests, as appropriate. Abbreviations: APTT = activated partial thromboplastin time; The CHA₂DS₂-VASc score comprises congestive heart failure, hypertension, age 65–74 (1 point) or age ≥75 (2 points), diabetes, prior stroke/TIA (2 points), vascular disease, and female sex (1 point) to estimate annual stroke risk in non-valvular atrial fibrillation; CrCl = creatinine clearance; HAS-BLED, hypertension, abnormal renal/liver function, stroke, bleeding history, labile INR, elderly age, and drugs/alcohol score; IQR, interquartile range; PT, prothrombin time; TIA, transient ischemic attack.
| Dose assessment characteristic | Edoxaban 60 mg (n = 38) | Edoxaban 30 mg (n = 41) | Rivaroxaban 15 mg (n = 35) | Rivaroxaban 10 mg (n = 32) | DOAC total (n = 146) |
| CrCl label reduction criterion met | 2 (5.3%) | 23 (56.1%) | 4 (11.4%) | 25 (78.1%) | 54 (37.0%) |
| Weight label reduction criterion met | 0 (0.0%) | 16 (39.0%) | Not applicable | Not applicable | 16 (11.0%) |
| Relevant P gp inhibitor criterion met | 0 (0.0%) | 8 (19.5%) | Not applicable | Not applicable | 8 (5.5%) |
| Label concordant dose | 36 (94.7%) | 32 (78.0%) | 4 (11.4%) | 0 (0.0%) | 72 (49.3%) |
| Lower than label dose | 0 (0.0%) | 9 (22.0%) | 31 (88.6%) | 32 (100.0%) | 72 (49.3%) |
| Higher than label dose | 2 (5.3%) | 0 (0.0%) | 0 (0.0%) | 0 (0.0%) | 2 (1.4%) |
Table 2: Official-label dose concordance at baseline. Edoxaban was assessed against 60 mg once daily, reduced to 30 mg once daily for CrCl 15–50 mL/min, body weight 60 kg or less, or a relevant P-glycoprotein inhibitor. Rivaroxaban was assessed against 20 mg once daily for CrCl greater than 50 mL/min and 15 mg once daily for CrCl 15–50 mL/min. Warfarin was individualized by INR and excluded from the DOAC concordance denominator. CrCl, creatinine clearance; DOAC, direct oral anticoagulant; INR, international normalized ratio; P-gp, P-glycoprotein.
| Outcome | Warfarin (n = 56) | Edoxaban 60 mg (n = 38) | Edoxaban 30 mg (n = 41) | Rivaroxaban 15 mg (n = 35) | Rivaroxaban 10 mg (n = 32) | Primary p value |
| Embolic event | 4 (7.1%) | 1 (2.6%) | 1 (2.4%) | 1 (2.9%) | 1 (3.1%) | 0.846 |
| ISTH minor bleeding | 11 (19.6%) | 3 (7.9%) | 3 (7.3%) | 2 (5.7%) | 2 (6.2%) | 0.199 |
| Other adverse reaction | 17 (30.4%) | 8 (21.1%) | 5 (12.2%) | 2 (5.7%) | 3 (9.4%) | 0.014 |
| Gastrointestinal symptoms | 10 (17.9%) | 3 (7.9%) | 2 (4.9%) | 1 (2.9%) | 2 (6.2%) | 0.133 |
| Dizziness/vertigo | 7 (12.5%) | 5 (13.2%) | 3 (7.3%) | 1 (2.9%) | 1 (3.1%) | 0.327 |
Table 3: Six-month clinical outcomes by anticoagulant group. Values are n (%). Primary p values are Fisher–Freeman–Halton exact tests for embolic and bleeding outcomes and Pearson chi-square tests for non-bleeding adverse reactions. All recorded bleeding was ISTH minor; no ISTH major or clinically relevant non-major bleeding occurred. ISTH, International Society on Thrombosis and Haemostasis.
| Laboratory indicator | Effect | Complete n | F | df1 | df2 | p value | Partial eta-squared |
| Hemoglobin, g/dL | Treatment group | 202 | 4.75 | 4.00 | 197.0 | 0.001 | 0.088 |
| Hemoglobin, g/dL | Time | 202 | 4.94 | 3.82 | 751.6 | <0.001 | 0.024 |
| Hemoglobin, g/dL | Treatment × time | 202 | 0.59 | 15.26 | 751.6 | 0.888 | 0.012 |
| ALT, U/L | Treatment group | 202 | 5.14 | 4.00 | 197.0 | <0.001 | 0.095 |
| ALT, U/L | Time | 202 | 0.04 | 3.95 | 778.2 | 0.996 | 0.000 |
| ALT, U/L | Treatment × time | 202 | 0.00 | 15.80 | 778.2 | 1.000 | 0.000 |
| Serum creatinine, mg/dL | Treatment group | 202 | 85.55 | 4.00 | 197.0 | <0.001 | 0.635 |
| Serum creatinine, mg/dL | Time | 202 | 0.03 | 3.91 | 769.5 | 0.998 | 0.000 |
| Serum creatinine, mg/dL | Treatment × time | 202 | 0.02 | 15.62 | 769.5 | 1.000 | 0.000 |
| PT, seconds | Treatment group | 202 | 4309.19 | 4.00 | 197.0 | <0.001 | 0.989 |
| PT, seconds | Time | 202 | 0.93 | 3.88 | 763.6 | 0.446 | 0.005 |
| PT, seconds | Treatment × time | 202 | 0.24 | 15.50 | 763.6 | 0.999 | 0.005 |
| APTT, seconds | Treatment group | 202 | 15.74 | 4.00 | 197.0 | <0.001 | 0.242 |
| APTT, seconds | Time | 202 | 4.48 | 3.91 | 769.6 | 0.002 | 0.022 |
| APTT, seconds | Treatment × time | 202 | 2.80 | 15.63 | 769.6 | <0.001 | 0.054 |
Table 4: Omnibus longitudinal analysis of laboratory indicators. Greenhouse–Geisser-corrected degrees of freedom are shown for time and treatment-by-time effects. Treatment-group main effects represent between-group differences and do not establish treatment-induced change. F, F-statistic; df1, numerator degrees of freedom; df2, denominator degrees of freedom; ALT, alanine aminotransferase; APTT, activated partial thromboplastin time; PT, prothrombin time.
Supplementary Table 1: Complete baseline characteristics. Continuous variables are median [interquartile range] and were compared with Kruskal–Wallis tests. Categorical variables are n (%) and were compared with Pearson chi-square or Fisher–Freeman–Halton exact tests, as appropriate. ALT, alanine aminotransferase; APTT, activated partial thromboplastin time; CrCl, creatinine clearance; PT, prothrombin time; TIA, transient ischemic attack. Please click here to download this file.
Supplementary Table 2: Dynamic CHA₂DS₂-VASc and HAS-BLED assessment through month 6. CHA₂DS₂-VASc was unchanged in all 202 participants. HAS-BLED changed by one point in 28 participants. A score of at least 3 denotes validated high bleeding risk; inclusion in the cohort required at least one documented bleeding-risk factor and did not require HAS-BLED ≥3. Please click here to download this file.
Supplementary Table 3: Official-label dose assessment at baseline and throughout follow-up. Edoxaban was assessed against 60 mg once daily, reduced to 30 mg once daily for CrCl 15–50 mL/min, body weight ≤60 kg, or concomitant use of a relevant P-glycoprotein inhibitor. Rivaroxaban was assessed against 20 mg once daily for CrCl >50 mL/min and 15 mg once daily for CrCl 15–50 mL/min. Warfarin dosing was individualized according to INR and was excluded from the DOAC concordance denominator. CrCl, creatinine clearance; DOAC, direct oral anticoagulant; INR, international normalized ratio; P-gp, P-glycoprotein. Please click here to download this file.
Supplementary Table 4: Detailed clinical event phenotypes. Bleeding events were classified with ISTH definitions. Gastrointestinal symptoms and dizziness/vertigo were non-bleeding adverse reactions and were analyzed separately. ISTH, International Society on Thrombosis and Haemostasis. Please click here to download this file.
Supplementary Table 5: Unadjusted pairwise risk ratios versus warfarin. Risk ratios use the observed six-month patient-level risks. Wald 95% CIs are shown; raw two-sided Fisher exact p values were adjusted with the Holm procedure within each outcome. Wide CIs reflect sparse events and do not establish equivalence. Please click here to download this file.
Supplementary Table 6: Exploratory generalized overlap-weighted risk ratios versus warfarin. Weights used multinomial propensity scores for age, sex, body weight, Cockcroft–Gault CrCl, Clinical Frailty Scale, CHA₂DS₂-VASc, HAS-BLED, and antiplatelet therapy. The effective sample size was 89.9. The maximum absolute pairwise standardized mean difference was 2.04 before and 1.36 after weighting. Bootstrap CIs are exploratory because substantial residual imbalance remained. Please click here to download this file.
Supplementary Table 7: Laboratory indicators at each assessment. Values are mean ± standard deviation. Measurements were complete for all 202 participants at baseline, day 3, month 1, month 3, and month 6 after verification of one source creatinine value. ALT, alanine aminotransferase; APTT, activated partial thromboplastin time; PT, prothrombin time. Please click here to download this file.
Supplementary Table 8: Baseline-to-month-6 laboratory changes by treatment group. Values are mean change ± standard deviation. Between-group comparisons used Kruskal–Wallis tests. ALT, alanine aminotransferase; APTT, activated partial thromboplastin time; PT, prothrombin time. Please click here to download this file.
Supplementary Table 9: Within-group nonparametric longitudinal sensitivity analysis. Friedman tests evaluated change across five assessments. Holm adjustment was applied across 25 tests. Kendall W is the within-group effect size. Please click here to download this file.